HDU - 4815(概率dp)

/*
    HDU - 4815 
    题意为
    n道题，老虎不低于P的概率下不会输
    求使老虎不会输于猴子的概率大于等于p的分数是多少


    转化下就是
    猴子得分小于等于k的概率大于p
    则老虎赢得的概率大于p的保险分数为k


    用概率dp
    dp[i][j]表示猴子前i道题得分为j的概率


    最后求和sum
    当sum加dp[n][k]之后概率大于等于p
    此时k即为老虎不低于p的概率不输于猴子的分数
*/
#include <iostream>
#include <algorithm>
#include <stdio.h>
#include <math.h>
#include <string.h>
#define mod p
#define MAX 3
#define ll long long
using namespace std;


double dp[45][40005];
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n;
        double p;
        cin>>n>>p;
        memset(dp,0,sizeof(dp));
        dp[0][0]=1;//猴子前0题得分为0的概率是1
        int max=0;
        for(int i=0;i<n;i++)
        {
            int w;
            cin>>w;
            for(int k=0;k<=max;k++)
            {
                dp[i+1][k+w]+=dp[i][k]*0.5;//有一半的概率答对
                dp[i+1][k]+=dp[i][k]*0.5;//有一半的概率答错
            }
            max+=w;
        }
        double sum=0.0;
        for(int i=0;i<=max;i++)
        {
            sum+=dp[n][i];
            if(sum>=p)//题目说的是不低于p不输
            {
                cout<<i<<endl;
                break;
            }
        }
    }
    return 0;
}