hdu4815 概率问题

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;

double dp[45][40005];///表示前i道题得分为j的概率
int a[45];

int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        int n;
        double p;
        cin>>n>>p;
        int sum=0;
        for(int i=1; i<=n; i++)
        {
            cin>>a[i];
            sum+=a[i];
        }
        memset(dp,0,sizeof(dp));
        dp[0][0]=1;  ///第0道题得分为0的概率为1
        for(int i=1; i<=n; i++)
        {
            for(int j=0; j<=sum; j++)
            {
                if(j>=a[i]) dp[i][j]+=dp[i-1][j-a[i]]/2;
                dp[i][j]+=dp[i-1][j]/2;
            }
        }
        double pp=0;
        int ans;
        for(int i=0; i<=sum; i++)
        {
            pp+=dp[n][i];
            if(pp>=p)
            {
                ans=i;
                break;
            }
        }
        cout<<ans<<endl;
    }
    return 0;
}