这题主要难点是高维前缀和,当然还有关于子图的联通性问题,其实是个很简单的dp,如果不太了解这方面的知识可以做hdu5713
现在就说一下高维前缀和:
首先我们想一个简单的二位前缀和,这个在字符串hash种很常见,
比如(0,0)~(1,1)
可以是
(0,0) (1,0)
(0,1) (1,1)
所以dp[1][1]=value[0][0]+value[1][0]+value[0]][1]+value[1][1];
dp[1][0]=value[0][0]+value[1][0];dp[0][1]=value[0][0]+value[0][1];
dp[0][0]=value[0][0];
回忆一下我们是怎么用进行加的,首先我们是先dp第列
dp[0][0]=value[0][0],dp[1][0]=value[1][0],
dp[0][1]=value[0][0]+value[0][1],
dp[1][1]=value[1][0]+value[1][1];
接着dp行
dp[1][1]+=dp[0][1];
dp[1][0]+=dp[0][0];
同样我们可以把这个方法进行拓展
在二位的时候先x轴,再y轴
那么三维就可以是x轴,y轴,然后z轴。。
如此继续下去,就是高维前缀和了
前缀和的性质是
dp[1][1][x][1][1][1]
这个dp[1][1][x][1][1][1]值包含了dp[][][x][][]的所有dp值的和
#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int t, n, m;
bool isconnected[2000200];
struct edgee
{
int from, to;
};
edgee rem[500];
int edge[500];
int dp[2000200],value[1000];
int main()
{
scanf("%d", &t);
int k = 1;
while (t--)
{
scanf("%d%d", &n, &m);
for (int i = 0; i <400; i++)
value[i] = 0;
for (int i = 0; i < n; i++)
edge[i] = 0;
for (int i = 0; i < m; i++)
{
int a, b;
scanf("%d%d", &a, &b);
edge[a] |= (1 << b);
edge[b] |= (1 << a);
rem[i].from = a; rem[i].to = b;
}
int sum = 1 << n;
for (int i = 1; i < sum; i++)
isconnected[i] = 0,
dp[i]=0;
for (int i = 0; i < n; i++)
{
int ss = 1 << i;
isconnected[ss] = 1;
}
for (int i = 1; i < sum; i++)
{
for (int j = 0; j < n; j++)
{
if (isconnected[i] == 1)
break;
int kk = 1 << j;
if ((i&kk) == 0)
continue;
int tempp = kk^i;
if (isconnected[tempp] == 0)
continue;
if (edge[j] & tempp)
isconnected[i] = 1;
}
}
for (int i = 1; i < sum; i++)
{
int another = (sum-1) - i;
if (isconnected[i] && isconnected[another])
dp[i] = 1;
}
for (int i = 0; i < n; i++)
{
for (int j = 1; j < sum; j++)
{
int kk = 1 << i;
if (kk&j)
continue;
int tempp = j | kk;
dp[tempp] += dp[j];
}
}
for (int i = 0; i < m; i++)
{
int a = rem[i].from; int b = rem[i].to;
int temp1 = 1 << a; int temp2 = 1 << b;
int temp3 = 0;
for (int j = 0; j < n; j++)
{
if (j == a || j == b)
continue;
temp1 |= (1 << j); temp2 |= (1 << j);
temp3 |= (1 << j);
}
value[i] = dp[temp1] + dp[temp2] - 2 * dp[temp3];
}
printf("Case #%d: ", k++);
for (int i = 0; i <m; i++)
{
if (i != m - 1)
printf("%d ", value[i]/2);
else
printf("%d\n", value[i]/2);
}
}
return 0;
}