CF16E Fish(概率dp+状压dp)

最新推荐文章于 2022-04-17 15:45:42 发布

Icefox_zhx

最新推荐文章于 2022-04-17 15:45:42 发布

阅读量594

点赞数

CC 4.0 BY-SA版权

分类专栏： codeforces ---------dp----------- 状压dp 概率与期望

本文链接：https://blog.youkuaiyun.com/Icefox_zhx/article/details/78275284

---------dp----------- 同时被 3 个专栏收录

304 篇文章

订阅专栏

codeforces

122 篇文章

订阅专栏

状压dp

44 篇文章

订阅专栏

本文介绍了一种使用概率转移矩阵解决特定问题的方法。通过构建概率矩阵并利用动态规划思想进行状态转移，最终实现了对所有可能状态的概率求解。代码中详细展示了如何初始化状态、读取输入及完成状态转移的过程。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这题好像思路还是很简单的。。。概率推了一会。。。太弱了。。。

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 20
inline int read(){
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x*f;
}
int n,bin[N];
double a[N][N],dp[300000];
int main(){
//  freopen("a.in","r",stdin);
    n=read();bin[0]=1;
    for(int i=1;i<=n;++i) bin[i]=bin[i-1]<<1;dp[bin[n]-1]=1;
    for(int i=0;i<n;++i)
        for(int j=0;j<n;++j) scanf("%lf",&a[i][j]);
    for(int s=bin[n]-1;s>=1;--s){
        int tot=0;for(int i=0;i<n;++i) if(s&bin[i]) tot++;
        if(tot==1) continue;double t=2.0/(tot*(tot-1));
        for(int i=0;i<n;++i){
            if((s&bin[i])==0) continue;
            for(int j=i+1;j<n;++j){
                if((s&bin[j])==0) continue;
                dp[s-bin[i]]+=a[j][i]*dp[s]*t;
                dp[s-bin[j]]+=a[i][j]*dp[s]*t;
            }
        }
    }
    for(int i=0;i<n;++i) printf("%.6f ",dp[bin[i]]);
    return 0;
}