[BZOJ] 2560:串珠子

 

子集枚举,DP,容斥

今天晚上突然不太舒服。。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define MOD 1000000007
#define MAXN 20
#define M 101000
using namespace std;

int n,up,a[MAXN][MAXN];
long long f[M],g[M];

int main(){
    cin>>n;
    for(int i=0;i<n;i++)for(int j=0;j<n;j++)cin>>a[i][j];
    up=1<<n;
    for(int i=1;i<up;i++){
        g[i]=1;
        for(int j=0;j<n;j++)if(i&(1<<j))
            for(int k=j+1;k<n;k++)if(i&(1<<k))
                g[i]=g[i]*(a[j][k]+1)%MOD;
        f[i]=g[i];
        int cur=i^(i&-i);
        for(int j=cur;j;j=(j-1)&cur)
            f[i]=(f[i]-g[j]*f[i^j]%MOD+MOD)%MOD;
    }
    cout<<f[up-1];
}

 

转载于:https://www.cnblogs.com/ghostcai/p/9383954.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值