bzoj2560串珠子 状压dp+容斥(?)

Posted wsy

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj2560串珠子 状压dp+容斥(?)相关的知识,希望对你有一定的参考价值。

2560: 串珠子

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 515  Solved: 348
[Submit][Status][Discuss]

Description

  铭铭有n个十分漂亮的珠子和若干根颜色不同的绳子。现在铭铭想用绳子把所有的珠子连接成一个整体。
  现在已知所有珠子互不相同,用整数1到n编号。对于第i个珠子和第j个珠子,可以选择不用绳子连接,或者在ci,j根不同颜色的绳子中选择一根将它们连接。如果把珠子看作点,把绳子看作边,将所有珠子连成一个整体即为所有点构成一个连通图。特别地,珠子不能和自己连接。
  铭铭希望知道总共有多少种不同的方案将所有珠子连成一个整体。由于答案可能很大,因此只需输出答案对1000000007取模的结果。

Input

 标准输入。输入第一行包含一个正整数n,表示珠子的个数。接下来n行,每行包含n个非负整数,用空格隔开。这n行中,第i行第j个数为ci,j。

Output

 标准输出。输出一行一个整数,为连接方案数对1000000007取模的结果。

Sample Input

3
0 2 3
2 0 4
3 4 0

Sample Output

50

HINT

  对于100%的数据,n为正整数,所有的ci,j为非负整数且不超过1000000007。保证ci,j=cj,i。每组数据的n值如下表所示。
编号 1 2 3 4 5 6 7 8 9 10  
n      8 9 9 10 11 12 13 14 15 16 

Source

2012国家集训队Round 1 day1

 

感觉还是有点不懂,这个dp来的莫名其妙
g[s]表示状态为s的点之间任意连边方案数
f[s]表示s点集构成连通块方案数
g[]很容易求

现在考虑f[s]应该怎么求
枚举s子集 把s划分成2个部分s1,s2,f[s]+=f[s1]*g[s2]
为什么这么转移可以得到答案呢?
因为这样可以保证s1和s2之间相互没有连边且可以枚举完所有子集情况又不会重复
感性理解。。

 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define mod 1000000007
#define N 20
#define M 101000
using namespace std;
int n,S,a[N][N];
long long f[M],g[M];

int main()
{
    int i,j,k,now;
    scanf("%d",&n);
    for(i=0;i<n;i++)
    for(j=0;j<n;j++)
    scanf("%d",&a[i][j]);
    S=1<<n;
    for(i=1;i<S;i++)
    {
        g[i]=1;
        for(j=0;j<n;j++)if(i&(1<<j))
            for(k=j+1;k<n;k++)if(i&(1<<k))
                g[i]=g[i]*(a[j][k]+1)%mod;
        f[i]=g[i];
        for(j=n-1;j+1;j--)
            if(i&(1<<j)){now=j;break;}
        now=i^(1<<now);
        for(j=now;j;j=(j-1)&now)
            f[i]=(f[i]-g[j]*f[i^j]%mod+mod)%mod;
    }
    cout<<f[S-1];
}

以上是关于bzoj2560串珠子 状压dp+容斥(?)的主要内容,如果未能解决你的问题,请参考以下文章

$bzoj2560$ 串珠子 容斥+$dp$

题解bzoj2560串珠子

题解-bzoj2560 串珠子

bzoj2560串珠子(子集dp)

bzoj 2560 串珠子

bzoj 2560 串珠子