zoukankan      html  css  js  c++  java
  • hdu 4649 Professor Tian 反状态压缩+概率DP

    思路:反状态压缩——把数据转换成20位的01来进行运算

    因为只有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[i-1][0]*p[i]+根据不同运算符和第i位的值运算得到0的概率。

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

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

    这题只要知道怎样表示状态就很简单了!!!以上是标程……

    代码如下:

     1 #include<iostream>
     2 #include<stdio.h>
     3 #include<algorithm>
     4 #include<iomanip>
     5 #include<cmath>
     6 #include<cstring>
     7 #include<vector>
     8 #define ll __int64
     9 #define pi acos(-1.0)
    10 #define MAX 50000
    11 using namespace std;
    12 double dp[202][2],p[202];
    13 char str[202];
    14 int f[202];
    15 double fun1(int i,int j){//位运算结果为1的期望值
    16     double ans=0.0,pp=1-p[i];
    17     int t=(f[i]>>j)&1;
    18     if(str[i]=='&'&&t){
    19         ans=dp[i-1][1]*pp;
    20     }
    21     else if(str[i]=='|'){
    22         ans=dp[i-1][1]*pp;
    23         if(t==1)
    24             ans+=dp[i-1][0]*pp;
    25     }
    26     else if(str[i]=='^'){
    27         if(t==1)
    28             ans=dp[i-1][0]*pp;
    29         else ans=dp[i-1][1]*pp;
    30     }
    31     return ans;
    32 }
    33 double fun0(int i,int j){//位运算结果为0的期望值
    34     double ans=0.0,pp=1-p[i];
    35     int t=(f[i]>>j)&1;
    36     if(str[i]=='&'){
    37         if(t==1)
    38             ans=dp[i-1][0]*pp;
    39         else ans=(dp[i-1][0]+dp[i-1][1])*pp;
    40     }
    41     else if(str[i]=='|'&&t==0)
    42         ans=dp[i-1][0]*pp;
    43     else if(str[i]=='^'){
    44         if(t==1) ans=dp[i-1][1]*pp;
    45         else ans=dp[i-1][0]*pp;
    46     }
    47     return ans;
    48 }
    49 int main(){
    50     int c,n,i,j,k=0;
    51     double ans;
    52     while(scanf("%d",&n)!=EOF){
    53         for(i=0;i<=n;i++) scanf("%d",&f[i]);
    54         for(i=1;i<=n;i++) scanf("%s",&str[i]);
    55         for(i=1;i<=n;i++) scanf("%lf",&p[i]);
    56         ans=0.0;
    57         for(i=0;i<20;i++){
    58             c=(f[0]>>i)&1;
    59             dp[0][1]=c; dp[0][0]=c^1;
    60             for(j=1;j<=n;j++){
    61                 dp[j][1]=dp[j-1][1]*(p[j])+fun1(j,i);
    62                 dp[j][0]=dp[j-1][0]*(p[j])+fun0(j,i);
    63             }
    64             ans+=(1<<i)*dp[n][1];
    65         }
    66         printf("Case %d:
    %.6lf
    ",++k,ans);
    67     }
    68     return 0;
    69 }
    View Code

  • 相关阅读:
    MFC基于对话框的程序添加菜单
    iOS远程推送原理及实现过程
    基于SolrCloud的内容搜索和热点推送
    Android 常用抓包工具介绍之Charles
    帮你快速实现全景应用性能可视化
    Android常用抓包工具之TcpDump
    JAVA IO 序列化与设计模式
    《2015中国移动应用性能管理白皮书》欢迎来看
    利用听云Server和听云Network实测Kubernetes和Mesos在高并发下的网络性能
    魅族电商运维之路
  • 原文地址:https://www.cnblogs.com/xin-hua/p/3248975.html
Copyright © 2011-2022 走看看