hdu 4649——Professor Tian

来源：互联网发布：js怎么判断undefined 编辑：程序博客网时间：2024/06/05 02:51

因为只有20位，而且&,|,^都不会进位，那么一位一位地看，每一位不是0就是1，这样求出每一位是1的概率，再乘以该位的十进制数，累加，就得到了总体的期望。

f[i][j]表示该位取前i个数，运算得到j(0或1)的概率是多少。

f[i][1]=f[i-1][1]*p[i]+根据不同运算符和第i位的值运算得到1的概率。

f[i][0]同理。

初始状态：f[0][0~1]=0或1（根据第一个数的该位来设置）

每一位为1的期望 f[n][1]

#include<iostream>#include<cstdio>using namespace std;int num[300];double p[300];char c[300];double dp[300][2];int main(){int n;int cnt=0;while(~scanf("%d",&n)){for(int i=1;i<=n+1;i++)scanf("%d",&num[i]);for(int i=2;i<=n+1;i++)cin>>c[i];for(int i=2;i<=n+1;i++)scanf("%lf",&p[i]);double ans=0.0;for(int i=1;i<=20;i++){dp[1][1]=num[1]&1;dp[1][0]=!dp[1][1];num[1]>>=1;for(int j=2;j<=n+1;j++){int k=num[j]&1;num[j]>>=1;if(c[j]=='|'){if(k==1){dp[j][1]=1*(1-p[j])+p[j]*dp[j-1][1];dp[j][0]=1-dp[j][1];}else{dp[j][1]=dp[j-1][1]*(1-p[j])+p[j]*dp[j-1][1];dp[j][0]=1-dp[j][1];}continue;}if(c[j]=='&'){if(k==1){dp[j][1]=dp[j-1][1]*(1-p[j])+p[j]*dp[j-1][1];dp[j][0]=1-dp[j][1];}else{dp[j][1]=0+p[j]*dp[j-1][1];dp[j][0]=1-dp[j][1];}continue;}if(c[j]=='^'){if(k==1){dp[j][1]=dp[j-1][0]*(1-p[j])+p[j]*dp[j-1][1];dp[j][0]=1-dp[j][1];}else{dp[j][1]=dp[j-1][1]*(1-p[j])+p[j]*dp[j-1][1];dp[j][0]=1-dp[j][1];}continue;}}ans+=(1<<(i-1))*dp[n+1][1];}printf("Case %d:\n%.6f\n",++cnt,ans);}return 0;}