hdu 5985 期望 概率dp

题目链接：https://vjudge.net/problem/HDU-5985

题目大意：有n种硬币，每种硬币在一次抛掷之后正面朝上的概率为p[i]，每次都将所有硬币都抛掷一次，反面的硬币出局，直至最后只剩一种硬币存活或者没有硬币存活，存活下来的硬币为lucky coin 问每种硬币成为lucky coin的概率。

, 

一眼概率dp。。

解题思路：第i种硬币在第j轮成为lucky coin的概率为 ：

第i种硬币恰好活到第j轮存活(在第j+1轮死了)的情况下，其他种类的硬币在前j轮全部死光。

结果就是对j从1到无穷求和。

根据上述的结果，我们现在需要知道两个变量，一个是第i种硬币在前k轮死掉的概率，和第i种硬币恰好在第k轮存活的概率，恰好的情况可以用（在第k轮存活的概率-在第k+1轮存活的概率）

这个地方一定是恰好在第k轮存活的概率，因为如果只用第k轮存活的概率来算，会包含在k+1.....轮成为lucky coin的情况，一旦在第k轮成为lucky coin，那么游戏就在第k轮结束了，不会包含之后的情况。

设die[][]和alive[][]两个数组来表示上述的两个结果

为一个第i种的硬币活到第j轮的概率，那么   

最后当n==1时直接输出1

注意当数字特别接近1的时候pow速度会很慢，这个地方需要快速幂，由于限制了所以步数的上限不会很大，几百就够了。

#include <iostream>
#include <cstring>
#include <algorithm>
#include <string>
#include <vector>
#include <unordered_set>
#include <unordered_map>
#include <queue>
using namespace std;
double eps=1e-10;
const int MAXD=500;
double die[12][MAXD+5];
double alive[12][MAXD+5];

double mul(double a,int b)
{
    double res=1.0;
    while(b)
    {
        if(b&1)
        {
            res=res*a;
        }
        a=(a*a);
        b>>=1;
    }
    return res;
}

int main()
{
    int a[15];
    double p[15];
    int t;
    scanf("%d", &t);
    while(t--)
    {
        int n;
        scanf("%d", &n);
        for(int i=1;i<=n;i++)
        {
            scanf("%d%lf", &a[i], &p[i]);
        }
        if(n==1){ printf("1.000000\n"); continue;}
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=MAXD;j++)
            {
                die[i][j]=mul((1.0-mul(p[i],j)),a[i]); // TODO 1
                alive[i][j]=1.0-die[i][j];
            }
        }
        double ans;
        double now,tem;
        for(int i=1;i<=n;i++)//shu zu
        {
            ans=0;
            now=1.0;
            for(int j=1;j<MAXD;j++)//time
            {
                tem=1.0;
                for(int k=1;k<=n;k++)//other die
                {
                    if(k==i)continue;
                    tem*=die[k][j];
                }
                now=(alive[i][j]-alive[i][j+1])*tem;
                ans+=now;
                //cout<<j<<endl;
            }
            printf(i==n?"%.6lf\n":"%.6lf ", ans);
        }
    }

    return 0;
}