概率dp,不要用pow会超时
代码如下
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
int T,n,m,t;
double dp[105][105];
double c[105][105];
double p[105];
int main()
{
p[0]=1;
for(int i=1;i<=100;i++) {
p[i]=p[i-1]*0.5;
}
for(int i=0;i<=100;i++) {
c[i][0]=1; c[i][i]=1;
}
for(int i=2;i<=100;i++) {
for(int j=1;j<i;j++) {
c[i][j]=c[i-1][j]+c[i-1][j-1];
}
}
scanf("%d",&T);
while(T--) {
scanf("%d%d%d",&n,&m,&t);
memset(dp,0,sizeof(dp));
dp[0][0]=1;
for(int i=1;i<=m;i++) {
for(int j=0;j<=n;j++) {
for(int k=0;k<=t;k++) {
if(n-j>=t) {
dp[i][j+k]+=dp[i-1][j]*c[t][k]*p[t];
}
else{
dp[i][n-t+k]+=dp[i-1][j]*c[t][k]*p[t];
}
}
}
}
double ans=0;
for(int i=0;i<=n;i++) {
ans+=dp[m][i]*i;
}
printf("%.3f\n",ans);
}
}