标签：cnt int 题解 Betting 状压 连向 集合 sum dp

题目链接

题目思路

数据范围很小，那么可以想到状压dp

设\(dp[i]\)表示胜利者的集合为\(i\)的概率，\(cal(i)\)为\(i\)的二进制数量

那么最后的答案为\(\large\sum_{i=0}^{i=(1<<n)-1}dp[i]*cal(i)\)

重点在于如何计算\(dp[i]\)

我们把\(n\)个元素分为两个集合，一个为\(A\)集合，一个为\(B\)集合

那么我们假设\(B\)集合中的所有边全部连向\(A\)集合

那么现在胜利者的集合是不是只能是\(B\)集合的所有子集的一种情况

假设\(B\)集合的所有边全部连向\(A\)集合的概率为\(prob\)

那么我们如何计算\(dp[B]\)

一种比较显然的想法就是减去B中所有真子集得的答案 即对于任意\(C\)集合为\(B\)的真子集

\(dp[B]=prob-\sum dp[C]\)

但这样不是正确的，为什么？

我们观察图，我们必须还要B的所有边全部指向\(C\)的

设\(rest=B \bigoplus C\)

设\(g\)为\(rest\)中所有的边全部连向\(A\)的概率

最后的答案为\(dp[B]=prob-\sum dp[C]\times g\)

时间复杂度为\(O(3^nn^2)\)

代码

#include<bits/stdc++.h>
#define fi first
#define se second
#define debug cout<<"I AM HERE"<<endl;
using namespace std;
typedef long long ll;
const int maxn=14+5,inf=0x3f3f3f3f,mod=1e9+7;
const double eps=1e-6;
int n;
int a[maxn];
int w[maxn][maxn];
bool vis[maxn];
int num[1<<14];
ll dp[1<<14];
ll qpow(ll a,ll b){
    ll ans=1,base=a;
    while(b){
        if(b&1) ans=ans*base%mod;
        base=base*base%mod;
        b=b>>1;
    }
    return ans;
}
int cal(int x){
    int cnt=0;
    while(x){
        if(x&1) cnt++;
        x=x/2;
    }
    return cnt;
}
signed main(){
    scanf("%d",&n);
    for(int i=0;i<n;i++){
        scanf("%d",&a[i]);
    }
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++){
            w[i][j]=a[i]*qpow(a[i]+a[j],mod-2)%mod;
        }
    }
    ll ans=0;
    for(int i=0;i<(1<<n);i++){
        for(int j=0;j<n;j++){
            if(i&(1<<j)) vis[j]=1;
            else vis[j]=0;
        }
        ll prob=1;
        for(int k=0;k<n;k++){
            if(!vis[k]) continue;
            for(int u=0;u<n;u++){
                if(vis[u]) continue;
                prob=prob*w[k][u]%mod;
            }
        }
        for(int j=i;j;j=(j-1)&i){
            ll g=1;
            int rest=(i^j);
            for(int k=0;k<n;k++){
                if(!(rest&(1<<k))) continue;
                for(int u=0;u<n;u++){
                    if(vis[u]) continue;
                    g=g*w[k][u]%mod;
                }
            }
            prob=((prob-g*dp[j])%mod+mod)%mod;
        }
        dp[i]=prob;
        ans=(ans+dp[i]*cal(i))%mod;
    }
    printf("%lld\n",ans);
    return 0;
}

标签：cnt,int,题解,Betting,状压,连向,集合,sum,dp
来源： https://www.cnblogs.com/hunxuewangzi/p/15233613.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。