2016多校联合训练7&&HDU5816

最新推荐文章于 2017-07-29 20:02:21 发布

just_sort

最新推荐文章于 2017-07-29 20:02:21 发布

阅读量809

点赞数

CC 4.0 BY-SA版权

分类专栏： ACM/ICPC_多校联合训练

本文链接：https://blog.youkuaiyun.com/just_sort/article/details/52167968

ACM/ICPC_多校联合训练专栏收录该内容

42 篇文章

订阅专栏

本文介绍了一种使用状态压缩动态规划（状压DP）结合动态内存分配的方法来解决一个具体的概率计算问题：给定特定条件下的攻击牌和抽牌机会，计算击败怪物的概率。该文提供了详细的AC代码实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

【题意】一个怪初始有P的血量，并且这个人初始有n次抽牌机会，并且这个抽牌是可以换两次抽牌机会的，还有M张攻击牌，问用这些攻击牌和N次抽牌机会杀死怪物的概率。

【解题方法】状压DP+动态内存分配，这道题如果静态开会无线MLE。

【AC 代码】

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define LL long long
int hp,n,m;
int a[22];
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%d",&hp,&n,&m);
        long long **dp;
        dp = new LL*[n + 1];
        for(int i = 0; i <= n; i++)
        {
            dp[i] = new LL[1 << m];
            memset(dp[i], 0, sizeof(LL) * (1 << m));
        }
        dp[n][0]=1;
        long long sum=0;
        for(int i=0; i<m; i++){
            scanf("%d",&a[i]);
            sum+=1LL*a[i];
        }
        if(sum<hp){
            printf("0/1\n");
            continue;
        }
        for(int i=n; i>=0; i--){
            for(int j=0; j<(1<<m); j++){
                if(dp[i][j]){
                    int cursum=0,curcnt=0;
                    for(int k=0; k<m; k++){
                        if(j&(1<<k)) cursum+=a[k],curcnt++;
                    }
                    int nexcnt=n-i+1-curcnt;
                    if(cursum>=hp){
                        for(int k=0; k<m; k++){
                            if(!(j&(1<<k))) dp[i][j|(1<<k)]+=dp[i][j];
                        }
                        if(i) dp[i-1][j]=dp[i][j]*i;
                    }
                    else if(nexcnt>=1){
                        for(int k=0; k<m; k++){
                            if(!(j&(1<<k))) dp[i][j|(1<<k)]+=dp[i][j];
                        }
                        if(i) dp[i-1][j]=dp[i][j]*i;
                    }
                }
            }
        }
        long long ans1=dp[0][(1<<m)-1];
        long long ans2=1;
        for(int i=1; i<=n+m; i++){
            ans2=ans2*(1LL)*i;
        }
        long long ans=__gcd(ans1,ans2);
        ans1/=ans;
        ans2/=ans;
        printf("%I64d/%I64d\n",ans1,ans2);
    }
}