upc 5212 Coins I

概率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);
    }
}