CF599E Sandy and Nuts

这题挺难的，因为 $n\le 13$，很自然地想到状压（雾）。

题目

dp 式子

令 $dp[i][j]$ 表示在状态 $i$ 下，当前根为 $j$ 的方案数，显然，加第二维方便转移。

可以写出一下转移方程：$dp[i][j]=\sum dp[u][t]\times dp[uxori][j]$。

最后的答案就是 $\sum dp[2^n-1][k]$。

这显然对吧，意为将 $t$ 为根的子树加入 $j$ 的子树中。（即将 $j$ 与 $t$ 连边）

但是在这里会重复计算（显然），所以我们可以选一个 $i$ 中的特殊点，强制设定 $u$ 中必须有这个元素，这样就可以避免重复（不妨想想为什么），我选的是不等于 $j$ 的第 $1$ 个点。

p.s.在这里要满足以下条件 $t\in u$（此为dp1），$j\in uxori$（此为dp2），$j\in i$（此为dp3）。

LCA

现在我们考虑 $\mathcal{L}\mathcal{C}\mathcal{A}$。因为其本质为限制，不妨考虑它限制了什么。（废话）

考虑若 $j==c_k$，则 $a_k$ 和 $b_k$ 不能同时出现在 $u$。（这很显然，因为这样 $\mathcal{L}\mathcal{C}\mathcal{A}(a_k,b_k)$ 就不等于 $c_k$ 了）（此为LCA1）

另外，若 $c_k==t$，$a_{k}or{b}_k$ 有一个不在 $u$ 中，舍。（此为LCA2）

原有边

若存在一边使一个端点在 $i$ 里（不等于 $j$），另一个端点不在 $i$ 里，舍。（原因显然，要保证它是棵树）（此为e1）

若存在多个$t$ 与 $j$ 的连边大于等于 $2$，舍。（此为e2）（p.s.若连边数为 $0$，则随便选 $t$，若连边数为 $1$，$t$ 被固定）

到了这里，所有的元素考虑完了，时间复杂度 $\mathcal{O}(3^n\times n(n+m+q))$（枚举子集是 $3^n$）。由于这个时间复杂度跑不满，所以能 A。

感悟：剪枝能在外层循环剪就在外层循环剪。

另外，安利一个枚举子集的方法 $for(intk=i;k;k=((k-1)\&i))$。

Code

#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <climits>
#include <vector>
#define LL long long
using namespace std;
struct Node {
	int X, Y;
	Node() {}
	Node(int x, int y) { X = x; Y = y; }
};
const int MAXN = 15, MAXM = 105;
int n, m, q, X_[MAXN], Y_[MAXN], num, tmp;
Node u;
LL dp[1 << MAXN][MAXN], res;
bool vis[1 << MAXN][MAXN];
vector <Node> v[MAXN]; 
bool Check_In(int x, int y) { return ((x >> (y - 1)) & 1); }
int main() {
	int x, y, z, f;
	scanf("%d%d%d", &n, &m, &q);
	for(int i = 1; i <= m; i ++) scanf("%d%d", &X_[i], &Y_[i]);
	for(int i = 1; i <= q; i ++) scanf("%d%d%d", &x, &y, &z), v[z].push_back(Node(x, y));
	for(int i = 1; i <= n; i ++) dp[1 << (i - 1)][i] = 1;
	for(int i = 1; i < (1 << n); i ++) { // LCA2
		for(int j = 1; j <= n; j ++) {
			f = 0;
			for(unsigned int k = 0; k < v[j].size(); k ++) {
				if(!Check_In(i, v[j][k].X) || !Check_In(i, v[j][k].Y)) { f = 1; break; }
			}
			vis[i][j] = f;
		}
	}
	for(int i = 1; i < (1 << n); i ++) {
		for(int j = 1; j <= n; j ++) {
			if(!Check_In(i, j)) continue; // dp3
			f = 0;
			for(int t = 1; t <= m; t ++) if(j != X_[t] && Check_In(i, X_[t]) && !Check_In(i, Y_[t])) f = 1; // e1
			for(int t = 1; t <= m; t ++) if(j != Y_[t] && Check_In(i, Y_[t]) && !Check_In(i, X_[t])) f = 1; // e1
			int qwq = (i ^ (1 << (j - 1))), mst = 1;
			if(!qwq || f) continue;
			while(!(qwq & 1)) qwq >>= 1, mst ++;
			for(int k = i; k; k = ((k - 1) & i)) {
				if(!Check_In(k, mst)) continue; // 防止重复 
				f = 0;
				for(unsigned t = 0; t < v[j].size(); t ++) { // LCA1
					u = v[j][t];
					if(Check_In(k, u.X) && Check_In(k, u.Y)) { f = 1; break; }
				}
				if(f) continue;
				num = 0; tmp = 0;
				for(int t = 1; t <= m; t ++) { // e2
					if(Check_In(k, X_[t]) && Y_[t] == j) tmp ++, num = X_[t];
					else if(Check_In(k, Y_[t]) && X_[t] == j) tmp ++, num = Y_[t];
				}
				if(!tmp) {
					for(int t = 1; t <= n; t ++) {
						if(!Check_In(k, t) || vis[k][t] || !Check_In(i ^ k, j)) continue; // dp1/2 & LCA2
						dp[i][j] += dp[k][t] * dp[i ^ k][j];
					}
				}
				else if(tmp == 1) {
					if(!Check_In(k, num) || vis[k][num] || !Check_In(i ^ k, j)) continue; // dp1/2 & LCA2
					dp[i][j] += dp[k][num] * dp[i ^ k][j];
				}
			}
		}
	}
	printf("%lld", dp[(1 << n) - 1][1]); // 注意 1号点为根	 
	return 0;
}