【 HDU 2048】 神、上帝以及老天爷 递推

神、上帝以及老天爷

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 55657    Accepted Submission(s): 22235

Problem Description

HDU 2006'10 ACM contest的颁奖晚会隆重开始了！
为了活跃气氛，组织者举行了一个别开生面、奖品丰厚的抽奖活动，这个活动的具体要求是这样的：

首先，所有参加晚会的人员都将一张写有自己名字的字条放入抽奖箱中；
然后，待所有字条加入完毕，每人从箱中取一个字条；
最后，如果取得的字条上写的就是自己的名字，那么“恭喜你，中奖了！”

大家可以想象一下当时的气氛之热烈，毕竟中奖者的奖品是大家梦寐以求的Twins签名照呀！不过，正如所有试图设计的喜剧往往以悲剧结尾，这次抽奖活动最后竟然没有一个人中奖！

我的神、上帝以及老天爷呀，怎么会这样呢？

不过，先不要激动，现在问题来了，你能计算一下发生这种情况的概率吗？

不会算？难道你也想以悲剧结尾？！

Input

输入数据的第一行是一个整数C,表示测试实例的个数，然后是C 行数据，每行包含一个整数n(1<n<=20),表示参加抽奖的人数。

Output

对于每个测试实例，请输出发生这种情况的百分比，每个实例的输出占一行, 结果保留两位小数(四舍五入)，具体格式请参照sample output。

Sample Input

1 
2

Sample Output

50.00%

题意：有n个人将自己的名字写在纸条上，然后将纸条放进抽奖箱中，求每个人都抽不到自己名字的概率。

题解：其实这是一个错排问题，n个有序的元素应有n!个不同的排列，如若一个排列使得所有的元素不在原来的位置上，则称这个排列为错排。我们用dp[n]表示n个人都没有抽中自己的不同抽奖总数，还是来个递推三部曲，1：我们先考虑能不能容易的得到一组简单的解，一个人的时候肯定能抽中，两个人的时候只有一种情况，那就是互相交换。2：我们假设规模为n-1的问题已经被解决。3：当规模扩大到n时，如何枚举出所有情况，并且对于每一种子情况都能用已有的数据来解决。因为要每个人都不能抽到自己的名字，那么第n个人肯定要和前面n-1个人中的一个交换。一共需要考虑两种情况，前面n-1人抽到的都不是自己的名字和前面n-1个人中只有一个人抽到自己的名字，那么我们为什么只用考虑这两种情况呢，因为我们假设规模为n-1的问题合法，去考虑规模为n的问题。由于我们假设前n-1个人抽到的都不是自己的名字合法，当考虑第n个人的时候，为了让规模为n的问题合法，前n-1个人中最多有一个人抽到自己的名字，当然也可以每个人都没抽到自己的名字。

总结一下：递推问题中的情况考虑，规模为n-1的问题合法，由于n的加入，规模为n-1的部分不合法问题可以合法。

1：前面n-1个人中抽到的都不是自己的名字，那么第n个人和前面n-1个人中的任何一个人交换都合法。dp[n]=(n-1)*dp[n-1]

2：前面n-1个人中有一个人抽到自己的名字，即有n-2个人抽到的不是自己的名字，那么第n个人和抽到自己名字的那个人交换为合法,抽到自己名字的那个人可以是n-1个人中的任意一个人。dp[n]=(n-1)*dp[n-2]

递推式：dp[n]=(n-1)*dp[n-1]+(n-1)*dp[n-2]

概率：p=dp[n]/n!

#include<cstdio>
using namespace std;
const int maxn=20+2;
typedef long long int ll;
ll dp[maxn];
ll f[maxn];
int n;
int prize()                 //错排
{
    dp[1]=0,dp[2]=1;
    for(int i=3;i<maxn;i++)
    dp[i]=(i-1)*dp[i-1]+(i-1)*dp[i-2];
    return 0;
}
int factorial()                // 20! 没有爆long long 
{
    f[1]=1;
    for(int i=2;i<maxn;i++)
    f[i]=f[i-1]*i;
    return 0;
}
int main()
{
    prize();
    factorial();
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        printf("%.2f%%\n",(double)dp[n]/f[n]*100);    // 输出% 格式控制为 %%
    }
    return 0;
}