hdu 5985 Lucky Coins (概率)

好有趣的一个题，真是让人进一步理解的概率问题，如果太小选择忽略，这个题目中很简单就能找到小的那个，0.5的30次方在double中已经为0了，为了保险起见可以选择100.

令两个函数 f[i][k] 表示第i种硬币第K步死光的概率，所以

f[i][k]=(1−pki)ni

revf[i][k] 表示第i种硬币第k步至少有一个活着所以

revf[i][k]=1−f[i][k]

ans[i]=∑k=1100∑j=1,j!=in(rev[i][k]−rev[i][k+1])∗f[j][k−1]

答案看起来很简单，但是想的是一个个算，然后没办法算，没想怎么去反着搞。。。

#include <bits\stdc++.h>
using namespace std;
double f[20][120],s[20][160],p[20],ans[20];
int a[20];
int main()
{
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){cout<<1.000000<<endl;continue;}
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=99;j++)
        {
            double temp=pow(p[i],j);
            temp=pow(1-temp,a[i]);
            f[i][j]=temp;
            s[i][j]=1-temp;
        }
    }
    for(int i=1;i<=n;i++)
    {
        ans[i]=0;
        for(int j=1;j<99;j++)
        {
            double temp=1;
            for(int k=1;k<=n;k++)
            {
                if(i==k)continue;
                temp*=f[k][j];
            }
            ans[i]+=(s[i][j]-s[i][j+1])*temp;
        }
    }
    for(int i=1;i<n;i++)
    {
        printf("%.6f ",ans[i]);
    }
    printf("%.6f\n",ans[n]);
}
    return 0;
}