首页 > 代码库 > POJ 3071 Football
POJ 3071 Football
概率,$dp$。
设dp[i][j]表示第i轮结束后,编号为j的人还活着的概率。$dp[i][j]=sum(dp[i-1][j]*dp[i-1][k]*p[j][k])$,即在第$i$轮所有与$i$可能交锋的人$k$的概率之和。
$a$与$b$能在第$i$轮交锋的条件是:$a/2^{i-1}!=b/2^{i-1}且a/2^i==b/2^i$。
#pragma comment(linker, "/STACK:1024000000,1024000000")#include<cstdio>#include<cstring>#include<cmath>#include<algorithm>#include<vector>#include<map>#include<set>#include<queue>#include<stack>#include<iostream>using namespace std;typedef long long LL;const double pi=acos(-1.0),eps=1e-6;void File(){ freopen("D:\\in.txt","r",stdin); freopen("D:\\out.txt","w",stdout);}template <class T>inline void read(T &x){ char c = getchar(); x = 0; while(!isdigit(c)) c = getchar(); while(isdigit(c)) { x = x * 10 + c - ‘0‘; c = getchar(); }}int n;double p[200][200],dp[10][200];int POW[10];int main(){ POW[0]=1; for(int i=1;i<=8;i++) POW[i]=2*POW[i-1]; while(~scanf("%d",&n)) { if(n==-1) break; for(int i=0;i<(1<<n);i++) for(int j=0;j<(1<<n);j++) scanf("%lf",&p[i][j]); for(int i=0;i<(1<<n);i++) dp[0][i]=1; for(int i=1;i<=n;i++) { for(int j=0;j<(1<<n);j++) { double sum=0; for(int k=0;k<(1<<n);k++) { int a=j,b=k; if(a/POW[i-1]!=b/POW[i-1]&&a/POW[i]==b/POW[i]) sum=sum+dp[i-1][j]*dp[i-1][k]*p[j][k]; } dp[i][j]=sum; } } double ans=dp[n][0]; int pos=1; for(int i=1;i<(1<<n);i++) if(ans<dp[n][i]) ans=dp[n][i], pos=i+1; printf("%d\n",pos); } return 0;}
POJ 3071 Football
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。