BZOJ2560（dp难）

最新推荐文章于 2019-10-27 08:15:38 发布

sherlock_zhuang

最新推荐文章于 2019-10-27 08:15:38 发布

阅读量1.1k

点赞数

分类专栏： my-code 文章标签： dp

本文链接：https://blog.csdn.net/sherlock_zhuang/article/details/50811370

版权

my-code 专栏收录该内容

28 篇文章 0 订阅

订阅专栏

思路：
很明显的状压dp
一开始写的dp可能会出现重复统计的情况而且难以去重
假设一个状态s的随意连边集合是A;
那么 A应该是全部合法的方案(Ans)+sigma(某一部分合法(即某一部分是连通图)的方案*其他任意连边的方案);
那么可以把最终答案设置为f[i], 随意连边(也可以完全连边)设置成g[i];
先定一个基准点 x 和基准点相连的都是合法的, 其余集合 t=s^(1<<(x-1))可以随便连;
f[i]=g[i]-sigma((t的所有子集i)f[i]*g[s^i]);
为什么一个是f 一个是g 这样其实是要保证不重不漏

/**************************************************************
    Problem: 2560
    User: DtenSherlock
    Language: C++
    Result: Accepted
    Time:1784 ms
    Memory:2300 kb
****************************************************************/


#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
#include<vector>
using namespace std;
typedef long long LL;
const int imax=20;
const int kmax=1<<16;
const int mod=1e9+7;
int n,T;
LL a[imax][imax];
LL f[kmax],g[kmax];

void iread()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++) scanf("%lld",&a[i][j]);    
}

void iwork()
{
    memset(f,0,sizeof(f));
    int Max=(1<<n)-1;
    for(int s=1;s<=Max;s++)
    {
        g[s]=1;
        for(int i=1;i<=n;i++)
            if(s&(1<<(i-1)))
            {
                for(int j=i+1;j<=n;j++)//注意枚举起点 不要重复 
                    if(s&(1<<(j-1))) 
                        g[s]=(g[s]*(a[i][j]+1))%mod;//不连边也可以 +1 
            }
        int nows=0; 
        for(int i=1;i<=n;i++) 
            if(s&(1<<(i-1))) { nows=(s^(1<<(i-1))); break;}
        f[s]=g[s];
        for(int i=nows;i;i=nows&(i-1))//枚举子集的tips
        {
            f[s]=(f[s]-((g[i]*f[s^i])%mod)+mod)%mod;
        } 
    }   
    printf("%lld\n",f[Max]);
}

int main()
{
    iread(); 
    iwork();
    return 0;
}