码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ4008】【HNOI2015】亚瑟王 概率DP

时间:2015-06-11 21:17:07      阅读:130      评论:0      收藏:0      [点我收藏+]

标签:bzoj4008   hnoi2015   亚瑟王   概率dp   

链接:

#include <stdio.h>
int main()
{
    puts("转载请注明出处[辗转山河弋流歌 by 空灰冰魂]谢谢");
    puts("网址:blog.csdn.net/vmurder/article/details/46461649");
}

题解:

f(i,j) 表示分配给第 [i,n] 张牌 j 次机会的期望。
然后 f(i,j)=f(i?1,j)?(1?pi?1)j)+f(i?1,j+1)?(1?(1?pi?1)j+1)

总结:

有的时候反过来想,比如求选的概率,记录不选的概率可能会更好一些。
再比如这道题,不从前面递推,而考虑 f(i) 表示 [i,n]。

代码:

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

int n,m;
long double ans;
long double p[233],d[233],f[233][333];

long double power(long double x,int p)
{
    long double ret=1.0;
    while(p)
    {
        if(p&1)ret*=x;
        x*=x,p>>=1;
    }
    return ret;
}

int main()
{
    freopen("test.in","r",stdin);
    int i,j,T;
    double a,b;
    for(scanf("%d",&T);T--;)
    {
        memset(f,0,sizeof f);
        scanf("%d%d",&n,&m);
        for(i=1;i<=n;i++)
        {
            scanf("%lf%lf",&a,&b);
            p[i]=a,d[i]=b;
        }
        f[0][m]=1.0,ans=0.0;
        for(i=1;i<=n;i++)for(j=1;j<=m;j++)
        {
            f[i][j]=f[i-1][j]*power(1-p[i-1],j)+f[i-1][j+1]*(1-power(1-p[i-1],j+1));
            ans+=f[i][j]*(1-power(1-p[i],j))*d[i];
        }
        printf("%.10lf\n",(double)ans);
    }
    return 0;
}

【BZOJ4008】【HNOI2015】亚瑟王 概率DP

标签:bzoj4008   hnoi2015   亚瑟王   概率dp   

原文地址:http://blog.csdn.net/vmurder/article/details/46461649

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!