\(n\)个技能 , 打\(k\)轮 , 每轮按顺序可能打出某个技能 , 第\(i\)个技能有\(p_i\)概率发动并打出伤害 , 且\(k\)轮内每个技能只能发动一次 , 求期望伤害
\(19.4.1\)
设第\(i\)张牌被使用的概率为\(f[i]\) , 则\(ans=\sum{f[i]*d[i]}\)
有\(f[1]=1-(1-p[i])^m\)
设\(g[i][j]\)表示前\(i\)张牌中使用\(j\)张的概率 , 则
\[f[i]=\sum{g[i-1][j]*(1-(1-p[i])^{m-j})}\]
即有\(j\)轮不会考虑到\(i\) , 还有\(m-j\)轮会考虑到\(i\) , 这些轮中又要选\(i\)的概率
转化为求\(g[i][j]\)的转移方程
\(1.\)取了第\(i\)张牌 \(g[i][j]=\sum{g[i-1][j-1]*(1-(1-p[i])^{m-j+1})}\)
\(2.\)没取第\(i\)张牌 \(g[i][j]=\sum{g[i-1][j]*(1-p[i])^{m-j}}\)
\(19.4.4\)
枚举当前点的思路 :
使用第\(i\)张牌的概率为\(f[i]\)
设\(g[i][j]\)表示前i轮中选\(j\)张牌的概率 , 那么就可以推到到当前第\(i\)张有多少轮是考虑过了的
(利用好题目意思可以想到 , 因为如果轮数没打完是一定要继续打的)
一般复杂的纯期望题都要转化一步 , 求另一个式子的转移方程 , 但也一定要考虑好题目要求的信息的转移方程
求\(g[i][j]\)的时候还是枚举当前点状态(取没取第\(i\)张) : 那么前面的要怎么样(已经取了\(j\)张或\(j-1\)张) , 转移条件是怎么样(剩下的轮怎么样)
这里也体现了阶段的重要性 , 即依次考虑每个点
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
const int INF=1e9+7;
inline LL read(){
register LL x=0,f=1;register char c=getchar();
while(c<48||c>57){if(c=='-')f=-1;c=getchar();}
while(c>=48&&c<=57)x=(x<<3)+(x<<1)+(c&15),c=getchar();
return f*x;
}
const int N=227;
const int M=143;
double g[N][M],pw[N][M],f[N],p[N];
int d[N],n,m;
inline double solve(){
n=read(),m=read();
for(int i=1;i<=n;i++){
scanf("%lf%d",&p[i],&d[i]);
f[i]=0,pw[i][0]=1;
for(int j=0;j<=m;j++) g[i][j]=0;
}
if(m==0) return 0;
for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++) pw[i][j]=pw[i][j-1]*(1.0-p[i]);
}
g[1][0]=pw[1][m],g[1][1]=f[1]=1.0-pw[1][m];
for(int i=2;i<=n;i++){
for(int j=0;j<=min(i,m);j++){
if(j) g[i][j]+=g[i-1][j-1]*(1.0-pw[i][m-j+1]);
if(i!=j) g[i][j]+=g[i-1][j]*pw[i][m-j];
}
}
for(int i=2;i<=n;i++){
for(int j=0;j<=min(i,m);j++) f[i]+=g[i-1][j]*(1.0-pw[i][m-j]);
}
double ans=0;
for(int i=1;i<=n;i++) ans+=f[i]*d[i];
return ans;
}
int main(){
for(int i=read();i;i--) printf("%.10lf\n",solve());
}