【动态规划02】bzoj3191卡牌游戏（概率）

最新推荐文章于 2021-09-17 14:29:02 发布

Flanoc

最新推荐文章于 2021-09-17 14:29:02 发布

阅读量826

点赞数

分类专栏：动态规划（弱省省选~NOI）文章标签：动态规划

本文链接：https://blog.csdn.net/Flanoc/article/details/72814487

版权

动态规划（弱省省选~NOI）专栏收录该内容

25 篇文章 0 订阅

订阅专栏

对于一个数学渣，且文化课并没有讲到概率与期望的一般普通高中生来讲，dp中最难搞的方面就是概率与期望，所以今天继续概率/期望。

题目描述

N个人坐成一圈玩游戏。一开始我们把所有玩家按顺时针从1到N编号。首先第一回合是玩家1作为庄家。每个回合庄家都会随机（即按相等的概率）从卡牌堆里选择一张卡片，假设卡片上的数字为X，则庄家首先把卡片上的数字向所有玩家展示，然后按顺时针从庄家位置数第X个人将被处决即退出游戏。然后卡片将会被放回卡牌堆里并重新洗牌。被处决的人按顺时针的下一个人将会作为下一轮的庄家。那么经过N-1轮后最后只会剩下一个人，即为本次游戏的胜者。现在你预先知道了总共有M张卡片，也知道每张卡片上的数字。现在你需要确定每个玩家胜出的概率。
这里有一个简单的例子：
例如一共有4个玩家，有四张卡片分别写着3,4,5,6.
第一回合，庄家是玩家1，假设他选择了一张写着数字5的卡片。那么按顺时针数1,2,3,4,1，最后玩家1被踢出游戏。
第二回合，庄家就是玩家1的下一个人，即玩家2.假设玩家2这次选择了一张数字6，那么2,3,4,2,3,4，玩家4被踢出游戏。
第三回合，玩家2再一次成为庄家。如果这一次玩家2再次选了6，则玩家3被踢出游戏，最后的胜者就是玩家2.

输入输出格式

第一行包括两个整数N,M分别表示玩家个数和卡牌总数。
接下来一行是包含M个整数，分别给出每张卡片上写的数字。
输出一行包含N个百分比形式给出的实数，四舍五入到两位小数。分别给出从玩家1到玩家N的胜出概率，每个概率之间用空格隔开，最后不要有空格。

对于100%的数据，有1<=N<=50 1<=M<=50 1<=每张卡片上的数字<=50
从数据范围上来看，这道题能达到O(n^4)的时间复杂度。
首先对于每个人获胜的概率，只与庄家与自己相对位置与剩余人数相关。
当满足这两个条件时的状态，显然可以看作是同一状态。

那么考虑每次转移，枚举抽到每张牌之后的贡献，设f[i][j]为剩余i人时，距离庄家第j个位置的人获胜的概率。
注意一下题内距离庄家位置的描述，距离庄家为1的人事庄家自己。
由逆推可知，显然如果该次转移令j被淘汰，则贡献为0.
以外的转移，由于m张牌等概率抽取，故贡献为f[i-1][lastj]/m，需要注意的是此时的lasyj和之前的j的数值不同，因为即使代表的是同一位置，但由于庄家改变，距离庄家的位置也有所改变。
由于抽到牌后的状态时固定的，所以我们可以算出lastj的。
由于是顺时针计算，所以当庄家在j的逆时针方向时与庄家在j的顺时针方向时的值是不同的。
具体的差异自己随便举个栗子用手推一推就推出来了，直接放代码。

#include<bits/stdc++.h>
#define fer(i,j,n) for(int i=j;i<=n;i++)
#define far(i,j,n) for(int i=j;i>=n;i--)
#define ll unsigned long long
const int maxn=60;;
const int INF=1e9+7;
const int mod=100003;
using namespace std;
/*----------------------------------------------------------------------------*/
inline int read()
{
    char ls;int x=0,sng=1;
    for(;ls<'0'||ls>'9';ls=getchar())if(ls=='-')sng=-1;
    for(;ls>='0'&&ls<='9';ls=getchar())x=x*10+ls-'0';
    return x*sng;
}
/*----------------------------------------------------------------------------*/
int n,m,a[maxn];
double f[maxn][maxn];
int main()
{
    n=read();m=read();
    fer(i,1,m)a[i]=read();
    f[1][1]=1;
    fer(i,2,n)
        fer(j,0,i)
            fer(k,1,m)
            {
                int t=a[k]%i;
                if(t==j)continue;
                if(!t)t=i;
                if(t>j)f[i][j]+=f[i-1][i-t+j]/m;
                if(t<j)f[i][j]+=f[i-1][j-t]/m;  
            }
    fer(i,1,n-1)
        printf("%.2lf%% ",f[n][i]*100);
    printf("%.2lf%%",f[n][n]*100);
}