• hdu 4649 Professor Tian 多校联合训练的题


    这题起初没读懂题意,悲剧啊,然后看了题解写完就AC了

    题意是给一个N,然后给N+1个整数 接着给N个操作符(只有三种操作  即  或 ,与 ,和异或 |   &  ^ )这样依次把操作符插入整数之间就可以得到一个表达式

    接着给出 N 给浮点数(在0~1之间表示概率 )表示的是 操作符和他右边的整数丢失的概率。 例如下面这组数据

    1

    1 2

    &

    0.5

    整数与操作符间可以组成一个表达式即 1&2 但是由于某些原因表达式的操作符和他操作的右边的

    那个数有一定的概率会丢失这组数据就是 &2有 0.5的概率会丢失 ,要你算这个表达式的期望值是

    多少 这组数据可以这样算 当&2不丢失的时候1&2=0

    当他丢失的时候表达式就变成了1 这样这个表达式的期望值就是 0*0.5+1*0.5=0.5 ;

    貌似没思路啊。。。。

    思路如下:

    先反状态压缩——把数据转换成20位的01来进行运算因为只有20位,而且&,|,^都不会进位,那么一位一位地看,每一位不是0就是1,这样求出每一位是1的概率,再乘以该位的十进制数,累加,就得到了总体的期望。对于每一位,状态转移方程如下:

     

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

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

    f[i][0]同理。

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

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

    最后的结果就是把每一位为一的概率乘上这位的位权相加就可以了

    下面是我的代码:

    #include<cstdio>
    #include<bitset>
    #include<iostream>
    using namespace std;
    bitset<20> team[205];//强烈推荐用bitset,这个比较快而且简单好用,关键还省空间 
    char ch[205];
    int pow[20],n,cas=0;
    double p[205], dp[20][205][2];
    int cal(int i,char c,int j,int f=0)//f=0表示 默认是计算通过运算为0的概率,f=1则是去计算为1的概率 
    {//这个函数计算的运算结果概率只能是1或0,所以可以返回真假值
    	int ans;
    	if(c=='&') ans=i&j;
    	if(c=='^') ans=i^j;
    	if(c=='|') ans=i|j;
    	if(!f)	return ans==0;//返回通过计算为0的概率 
    	else  return ans==1;//返回通过计算1的概率 
    }
    int read()
    {
    	if(!(cin>>n)) return 0;//没有读入了就返回假值 
    	for(int a,i=0;i<=n;i++)
    	{
    		cin>>a;
    		team[i]=a;//用bitset转化为二进制 
    	}
    	for(int i=0;i<n;i++)
    		cin>>ch[i];//读操作符 
    	for(int i=0;i<n;i++)
    		cin>>p[i];//读概率 
    	return 1;
    }
    void deal()
    {
    	double ans=0;
    	for(int i=0;i<20;i++)//初始化第一位的概率 
    	if(team[0][i]) dp[i][0][1]=1,dp[i][0][0]=0;
    	else dp[i][0][1]=0,dp[i][0][0]=1;
    	
    	for(int i=0;i<19;i++)
    	for(int j=1;j<=n;j++)
    	{
    		dp[i][j][1]=dp[i][j-1][1]*p[j-1]+(1-p[j-1])*(dp[i][j-1][0]*cal(0,ch[j-1],team[j][i],1)+dp[i][j-1][1]*cal(1,ch[j-1],team[j][i],1));//这个dp式如果还不懂就看一下,文章尾部我的说明。。
    		//计算当前为1的概率 
    		dp[i][j][0]=dp[i][j-1][0]*p[j-1]+(1-p[j-1])*(dp[i][j-1][0]*cal(0,ch[j-1],team[j][i])+dp[i][j-1][1]*cal(1,ch[j-1],team[j][i]));
    		//计算当前为0的概率 和为1的概率是一回事的
    	}
    	
    	for(int i=0;i<20;i++)
    		ans+=pow[i]*dp[i][n][1];//统计结果 
    	printf("Case %d:
    ",++cas);
    	printf("%.6lf
    ",ans);
    }
    int main()
    {
    	for(int i=pow[0]=1;i<20;i++)
    		pow[i]=pow[i-1]*2;//计算二进制的位权 
    	while(read()) deal();
    	return 0;
    }

    dp[i][j][1]=dp[i][j-1][1]*p[j-1]+(1-p[j-1])*(dp[i][j-1][0]*cal(0,ch[j-1],team[j][i],1)+dp[i][j-1][1]*cal(1,ch[j-1],team[j][i],1));

    其中dp[i][j][1]表示前j个操作后使第i位变成1的概率,dp[i][j-1][1]*p[j-1],表示当前的操作符和其右边的数消失后,使这位为1的概率

    (1-p[j-1])*(dp[i][j-1][0]*cal(0,ch[j-1],team[j][i],1)+dp[i][j-1][1]*cal(1,ch[j-1],team[j][i],1));这整串是在算当前的操作符和其右边的数不消失的概率,所以有公因子(1-p[j-1])即表示,表达式不消失的概率,dp[i][j-1][0]*cal(0,ch[j-1],team[j][i],1)这句是在算如果上一步当前这位计算的结果是0,那么通过这步运算得到1的概率 ,dp[i][j-1][0] 这是上一步算出这位结果为0的概率,cal(0,ch[j-1],team[j][i],1)这是算0和当前的操作符及其右边的数进行运算为1的概率,这个概率要么是一要么是0,dp[i][j-1][1]*cal(1,ch[j-1],team[j][i],1)这句就是算若上次运算的结果是1那么和当前这一步运算后得到的是1的概率,终于说完这个恶心的式子了。。。下面那个算0的式子就可以对应的去理解了,即dp[i][j][0]=dp[i][j-1][0]*p[j-1]+(1-p[j-1])*(dp[i][j-1][0]*cal(0,ch[j-1],team[j][i])+dp[i][j-1][1]*cal(1,ch[j-1],team[j][i]));


  • 相关阅读:
    SVG PATH d参数的 ace
    如果你想动态创建一个iframe
    canvas构建一个平面直角坐标系
    JavaScript版几种常见排序算法
    【分享】 Stip,让表单验证轻松愉快。
    【分享】javascript合并混淆压缩
    【分享】页面提示插件更新
    nodejs 的 session 简单实现
    页面模块之间的通信依赖解决方案
    在平面旋转一个点
  • 原文地址:https://www.cnblogs.com/riskyer/p/3244061.html
Copyright © 2020-2023  润新知