8.20T3 无损加密（线性代数转LGV+状压dp+高维前缀和）

http://cplusoj.com/d/senior/p/NODSX2301C

对于式子：

这个神秘的线性代数形式比较难处理，但我们可以考虑其组合意义。行列式现存的可用组合意义之一就是LGV（矩阵式不太可用）

先把原先的矩阵转化为一个有向图。现在我们要构造一个图，满足 $B_{i,j}$ 代表从 $a_i$ 走到 $b_j$ 所有路径权值积的和为 $B_{i,j}$，而上面行列式求的就是从 $a_{1\to n}$ 走到 $b$ 的任意 $\left(\genfrac{}{}{0ex}{}{m}{n}\right)$ 个点的不想交路劲数。

现在我们已经理清条件，考虑构造图了。我们先弄一行 $a$ 和 $b$，

现在是满足条件的。

我们考虑在 $[l,r]$ 内加入 $d,c$，不妨令其为 $[2,3]$，在外面加没问题：

而对于 [l,r] 内的情况，我们相当于在原先基础上多了一个不断往右增 $c$ 的边：

这样子建边能完美解决我们的问题。

现在只需要对最终的图统计不相交路径数即可。

我们考虑第 $i$ 列的链：

根据题目，这条链的总入度和总出度不超过 $s$，其实我们状压一下。

我们可以先钦定哪些入度有流，从上一步的 $f[S]$ 转移过来，当然必须满足红色点也有流。

然后我们直接让这些流先走到最近的出度。

当走到最近的出度后，他们还可以往下走，此时我们可以拿一个类似高维前缀和的东西来实现。而这里，我们要按顺序，先走最上面，再往下一个一个考虑，那样子才能保证路径不交。

我们可以设 $g(s,i)$ 表示当前状态是 $s$，考虑到第 $i$ 个1的路径权值积之和、

最后我们把黄色点的东西去掉，就可以转移到新 $f[T]$ 了

答案为 $f[0]$

复杂度我实现有点劣，为 $O(n{2}^s{s}^2)$，但很多情况会直接停止，所以跑不满。

4k

#include<bits/stdc++.h>
using namespace std;
#ifdef LOCAL
 #define debug(...) fprintf(stdout, ##__VA_ARGS__)
 #define debag(...) fprintf(stderr, ##__VA_ARGS__)
#else
 #define debug(...) void(0)
 #define debag(...) void(0)
#endif
#define int long long
inline int read(){int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;
ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<1)+
(x<<3)+(ch^48);ch=getchar();}return x*f;}
#define Z(x) (x)*(x)
#define pb push_back
#define fi first
#define se second
#define M 9
#define mo (int)(1e9 + 7)
#define N 200010
int pw(int a, int b) {
	int ans = 1; 
	while(b) {
		if(b & 1) ans *= a; 
		a *= a; b >>= 1; 
		ans %= mo; a %= mo; 
	}
	return ans; 
}
int pw(int a) { return pw(a, mo - 2); }
inline void Mod(int &a) { if(a >= mo || a <= -mo) a %= mo; if(a < 0) a += mo; }
inline void Add(int &a, int b) { a += b; Mod(a); }
inline void Mul(int &a, int b) { Mod(b); a *= b; Mod(a); } 
int n, m, i, j, k, T;
int C[N * M], D[N * M], q, l, r, Ds[N * M]; 
vector<int>L[N], R[N]; 
int nl, nr, s, t, f[1 << 10], g[1 << 10][11], nxt[15]; 
int ans, c[15], d[15], del[15]; 

signed main()
{
	#ifdef LOCAL
	  freopen("in.txt", "r", stdin);
	  freopen("out.txt", "w", stdout);
	#endif
//	srand(time(NULL));
//	T=read();
//	while(T--) {
//
//	}
	n = read(); m = read(); q = read(); Ds[0] = C[q + 1] = 1; 
	for(T = 1; T <= q; ++T) {
		l = read(); r = read(); C[T] = read(); D[T] = read(); 
		Ds[T] = D[T] * Ds[T - 1] % mo; 
		for(j = l + 1; j <= r; ++j) L[j].pb(T); 
		for(j = l; j < r; ++j) R[j].pb(T); 
	}
	Ds[q + 1] = Ds[q]; 
	debug("Ds : "); for(i = 0; i <= q + 1; ++i) debug("%lld ", Ds[i]); debug("\n"); 
	f[0] = 1; 
	for(T = 1; T <= m; ++T) {
		nl = L[T].size(); j = 0; 
		R[T].pb(q + 1); nr = R[T].size(); 
		for(auto t : L[T]) debug("%lld ", t); debug("\n"); 
		for(auto t : R[T]) debug("%lld ", t); debug("\n"); 
		for(auto t : L[T]) {
			for(k = 0; k < nr; ++k) 
				if(R[T][k] >= t) break; 
			nxt[j] = k; del[j] = Ds[R[T][k]] * pw(Ds[t]) % mo; ++j; 
			debug("[%lld %lld]%lld ", t, R[T][k], nxt[j - 1]); 
		}
		debug("\n"); 
		memset(g, 0, sizeof(g)); 
		for(s = 0; s < (1 << nl); ++s) {
			debug("f[%lld] = %lld\n", s, f[s]); 
			t = (T <= n ? 1 : 0); 
			int G = (T <= n ? Ds[R[T][0]] : 1); 
			for(j = 0; j < nl; ++j) {
				if(!(s & (1 << j))) continue; 
				if(t & (1 << nxt[j])) break; 
				t |= (1 << nxt[j]); G = G * del[j] % mo; 
			}
			if(j < nl) continue; 
			debug("--> %lld[%lld %lld] += %lld\n", t, t, 1ll, f[s]); 
//			if(!t) Add(g[t][0], f[s]); 
//			else 
			Add(g[t][1], f[s] * G % mo); 
		}
		for(i = 0; i < nr; ++i) c[i] = C[R[T][i]]; 
		for(i = 0; i < nr - 1; ++i) d[i] = Ds[R[T][i + 1]] * pw(Ds[R[T][i]]) % mo; 
		debug("# C : "); for(i = 0; i < nr; ++i) debug("%lld ", c[i]); debug("\n"); 
		debug("# D : "); for(i = 0; i < nr - 1; ++i) debug("%lld ", d[i]); debug("\n"); 
		for(i = 1; i <= nl + 1; ++i) //第 
			for(s = 0; s < (1 << nr); ++s) {
//				debug("(%lld %lld)\n", s, i); 
				if(!g[s][i]) continue; 
				if(__builtin_popcount(s) < i) continue; 
				for(k = j = 0; k < nr; ++k) {
					if((s >> k) & 1) ++j; 
					if(j == i) break; 
				}
				debug("(%lld %lld)%lld * %lld [%lld] => (%lld %lld)\n", s, i, g[s][i], c[k], k, s, i + 1); 
				Add(g[s][i + 1], g[s][i] * c[k] % mo); 
				if(k + 1 >= nr || ((s >> k + 1) & 1)) continue; 
				t = s - (1 << k) + (1 << k + 1); 
				debug("(%lld %lld)%lld * %lld [%lld] => (%lld %lld)\n", s, i, g[s][i], d[k], k, t, i); 
				Add(g[t][i], g[s][i] * d[k] % mo); 
			}
		memset(f, 0, sizeof(f)); 
		for(s = 0; s < (1 << nr); ++s) {
			t = s & (1 << nr - 1) - 1; 
			int cnt = __builtin_popcount(s); 
			debug(">>> %lld ---> %lld (%lld %lld)[%lld]\n", s, t, s, cnt + 1, g[s][cnt + 1]); 
			Add(f[t], g[s][cnt + 1]); 
		}
		debug("------------\n"); 
	}
	ans = f[0]; 
	printf("%lld", ans); 
//	for(i = 1; i <= m; ++i) assert(L[i].size() <= 8 && R[i].size() <= 8); 
	return 0;
}