洛谷 P2059 [JLOI2013]卡牌游戏

题目描述

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的胜出概率，每个概率之间用空格隔开，最后不要有空格。

输入输出样例

输入样例#1：

输入样例1：

5 5
2 3 5 7 11

输入样例2：

4 4
3 4 5 6

输出样例#1：

输出样例1：

22.72% 17.12% 15.36% 25.44% 19.36%

输出样例2：

25.00% 25.00% 25.00% 25.00%

说明

对于30%的数据，有1<=N<=10
对于50%的数据，有1<=N<=30
对于100%的数据，有1<=N<=50 1<=M<=50 1<=每张卡片上的数字<=50

分析

这是一个约瑟夫问题的变种。
一、约瑟夫问题
1、什么是约瑟夫问题：
n个人（编号0~(n-1）），从0开始报数，报到（m-1）的退出，剩下的人继续从0开始报数。求胜利者的编号。（来自百度百科）
2、约瑟夫问题的解法
记f(i)表示i个人玩约瑟夫游戏的最终胜利者。
显然f(1) = 0
如果n个人玩约瑟夫游戏，显然，第一轮游戏后，(m-1) mod n号玩家将会退出。问题变成了剩下的n-1个人玩约瑟夫游戏。
我们将这n-1个人重新编号，令k = m mod n（也就是退出游戏的玩家的下一位）。
k –> 0

k+1 –> 1

k+2 –> 2
… … k-2 –> n-2
假如我们知道f(n-1)的答案是x，那么x号玩家在n个人玩的游戏中的编号应该是(x+m) mod n。
也就是说f(n) = (f(n-1)+m) mod n
这样我们就得到了递推式，根据该式就可以得到约瑟夫问题的答案。
二、本问题
由于m的值在本问题中无法确定，那么游戏的胜利者也无法确定。
所以，我们记f(i,j)为i个人玩约瑟夫游戏，j号玩家胜利的概率。
由约瑟夫问题可知，f(i+1,j) = 求和(f(i,k)*1/m) 。设第x张卡片上的数字为px，则k为满足(k+x) mod (i+1) = j的所有数。
这样我们就求出了所有玩家获胜的概率。

代码

#include <bits/stdc++.h>

#define N 55

double f[N][N];
int num[N];
int n,m;
int main()
{
    scanf("%d%d",&n,&m);
    for (int i = 1; i <= m; i++)
        scanf("%d",&num[i]);
    f[1][1] = 1.0;
    for (int i = 2; i <= n; i++)
        for (int j = 1; j <= n; j++)
            for (int  k = 1; k <= m; k++)
            {
                int now = num[k] % i;
                if (now == 0)
                    now = i;
                if (now < j)
                    f[i][j] = f[i][j] + f[i - 1][j - now] / m;
                if (now > j)
                    f[i][j] = f[i][j] + f[i - 1][i - now + j] / m;
            }
    for (int i = 1; i <= n; i++)
    {
        printf("%.2lf%s",f[n][i]*100,"%");
        if (i != n)
            printf(" ");
    }
}