题目链接: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;
}