考虑 2s 容斥，那么限制就变成了不能经过某些点，求路径长度为 l 的方案 数，直接使用矩阵快速幂优化即可，时间复杂度 O(2sn3 log l)。