一个集合中如果包含三个点相互连通或相互都不连通,那么这个集合算是合法,问有多少个这样的集合。
我们反过来想,如果正着求有多少个这样的集合,那么重复状况很难排除,反过来算有多少种不合法的集合,发现其实种类只有很少,直接用全局减去这些不合法集合的方案数就是答案了。
( 可以手动画一下,各种不合法的情况,而且由定理得到不合法情况点数不会达到6或以上)
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
const int maxn = 60;
const int MOD = 1000000007;
bool f[maxn][maxn];
int n, m, x, y, sum, ans;
int temp1, temp2, temp3;
int main() {
int tt;
scanf("%d", &tt);
for (int cases = 1; cases <= tt; cases++) {
scanf("%d%d", &n, &m);
memset(f, 0, sizeof(f));
for (int i = 1; i <= m; i++) {
scanf("%d%d", &x, &y);
f[x][y] = true;
f[y][x] = true;
}
sum = 1+n+n*(n-1)/2;
temp1 = 0;
temp2 = 0;
for (int i = 1; i <= n; i++)
for (int j = 1; j <= n; j++) if (i != j)
for (int k = 1; k <= n; k++) if (i != k && j != k) {
if (f[i][j] && !f[i][k] && !f[j][k]) temp1 = (temp1+1)%MOD;
else if (f[i][j] && f[j][k] && !f[i][k]) temp2 = (temp2+1)%MOD;
}
temp1 = temp1/2;
temp2 = temp2/2;
sum = sum+temp1+temp2;
temp1 = 0;
temp2 = 0;
temp3 = 0;
for (int i1 = 1; i1 <= n; i1++)
for (int i2 = 1; i2 <= n; i2++) if (i1 != i2)
for (int i3 = 1; i3 <= n; i3++) if (i1 != i3 && i2 != i3 && f[i1][i3] && !f[i2][i3])
for (int i4 = 1; i4 <= n; i4++) if (i1 != i4 && i2 != i4 && i3 != i4) {
if (f[i1][i2] && f[i1][i3] && !f[i1][i4] && !f[i2][i3] && f[i2][i4] && f[i3][i4]) temp1 = (temp1+1)%MOD;
else if (!f[i1][i2] && f[i1][i3] && !f[i1][i4] && !f[i2][i3] && f[i2][i4] && f[i3][i4]) temp2 = (temp2+1)%MOD;
else if (!f[i1][i2] && f[i1][i3] && !f[i1][i4] && !f[i2][i3] && f[i2][i4] && !f[i3][i4]) temp3 = (temp3+1)%MOD;
}
temp1 = temp1/8;
temp2 = temp2/2;
temp3 = temp3/8;
sum = (sum+temp1+temp2+temp3)%MOD;
temp1 = 0;
for (int i1 = 1; i1 <= n; i1++)
for (int i2 = 1; i2 <= n; i2++) if (i1 != i2 && f[i1][i2])
for (int i3 = 1; i3 <= n; i3++) if (i1 != i3 && i2 != i3 && !f[i1][i3] && f[i2][i3])
for (int i4 = 1; i4 <= n; i4++) if (i1 != i4 && i2 != i4 && i3 != i4 && !f[i1][i4] && !f[i2][i4] && f[i3][i4])
for (int i5 = 1; i5 <= n; i5++) if (i1 != i5 && i2 != i5 && i3 != i5 && i4 != i5 && f[i1][i5] && !f[i2][i5] && !f[i3][i5] && f[i4][i5]) {
temp1 = (temp1+1)%MOD;
}
temp1 = temp1/10;
sum = (sum+temp1)%MOD;
ans = 1;
for (int i = 1; i <= n; i++) ans = (ans*2)%MOD;
ans = ans-sum;
if (ans < 0) ans = ans+MOD;
printf("Case #%d: %d\n", cases, ans);
}
}