bzoj2560 串珠子 子集dp

       虽然我这么弱但做做这种sb题还是没问题的。。。都水了一个早上sb题了还不狗带

       令f[i]表示集合为i时的方案数,然后全集减去不合法的就好了。不合法的话可以用一个定点来枚举。时间复杂度O(3^N+N^2·2^N)。

AC代码如下:

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
#define mod 1000000007
using namespace std;

int n,bin[25],f[1<<16],g[1<<16],a[25][25];
int main(){
	scanf("%d",&n); int i,j,k;
	bin[0]=1; for (i=1; i<=n; i++) bin[i]=bin[i-1]<<1;
	for (i=1; i<=n; i++)
		for (j=1; j<=n; j++) scanf("%d",&a[i][j]);
	for (k=1; k<bin[n]; k++){
		f[k]=1;
		for (i=1; i<=n; i++) if (k&bin[i-1])
			for (j=1; j<i; j++) if (k&bin[j-1])
				f[k]=(ll)f[k]*(a[i][j]+1)%mod;
		g[k]=f[k]; i=k^(k&-k);
		for (j=i; j; j=(j-1)&i) f[k]=(f[k]-(ll)g[j]*f[k^j]%mod+mod)%mod;
	}
	printf("%d\n",f[bin[n]-1]);
	return 0;
}


by lych

2016.4.16

你可能感兴趣的:(动态规划,子集Dp)