<SCOI2008>奖励关

weixin_30483013

于 2018-10-31 21:11:00 发布

阅读量60

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/pile8852/p/9885964.html

本文探讨了一种结合期望动态规划（Expectation DP）与状态压缩（State Compression）的算法解决策略，通过C++实现了一个具体示例。该算法主要用于处理涉及概率与状态空间大的问题，如背包问题的变种等。文章详细展示了算法的代码实现，包括状态转移方程的推导与优化，为读者提供了深入理解与实践的机会。

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

emmm第一道期望dp+个状压

真有趣..

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#define max(a,b) (a>b?a:b)
#define rint register int
template <class T>inline void read(T &X)
{
    X=0;int W=0;char ch=0;
    while(!isdigit(ch))W|=ch=='-',ch=getchar();
    while(isdigit(ch))X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    X=W?-X:X;return;
}

int K,n,maxn,need[55];
double ans=0,p[20],dp[55][1<<16];
int main()
{
    read(K),read(n);
    maxn=(1<<n)-1;
    for(rint i=1;i<=n;++i)
    {
        scanf("%lf",&p[i]);
        while(1)
        {
            int x;
            read(x);
            if(!x)break;
            need[i]=need[i]|(1<<x-1);
        }
    }
    
    for(rint i=K;i>0;--i)
    {
        for(rint s=0;s<=maxn;++s)//枚举状态 
        {
            for(rint j=1;j<=n;++j)//枚举可能的物品种类 
            {
                if((s&need[j])!=need[j]) {dp[i][s]+=dp[i+1][s];continue;}
                dp[i][s]+=max(dp[i+1][s],dp[i+1][s|(1<<j-1)]+p[j]);
            }
            dp[i][s]/=n;
        }
    }
    printf("%.6lf\n",dp[1][0]);
return 0;
}