hdu 2048 神、上帝以及老天爷 （错排公式，递推）

神、上帝以及老天爷

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

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%
  
  

错排：若一个排列式的所有的元素都不在原来的位置上，则称这个排列为错排；

1，递推的方法求错排公式：

当n个编号元素放在n个编号位置，元素编号与位置编号各不对应的方法数用M(n）表示，那么M(n-1）就表示n-1个编号元素放在n-1个编号位置，各不对应的方法数，其它类推.

第一步，把第n个元素放在一个位置，比如位置k，一共有n-1种方法；

第二步，放编号为k的元素，这时有两种情况⑴把它放到位置n，那么，对于剩下的n-1个元素，由于第k个元素放到了位置n，剩下n-2个元素就有M(n-2）种方法；⑵第k个元素不把它放到位置n，这时，对于这n-1个元素，有M(n-1）种方法；

综上得到

M(n)=(n-1）[M(n-2）+M(n-1）]

特殊地，M⑴=0,M⑵=1

下面通过这个递推关系推导通项公式：

为方便起见，设M(k)=k!N(k),(k=1,2，…，n)

则N⑴=0,N⑵=1/2

n>=3时，n!N(n)=(n-1）（n-1）！N(n-1）+(n-1）！N(n-2）

即 nN(n)=(n-1）N(n-1）+N(n-2）

于是有N(n)-N(n-1）=-[N(n-1）-N(n-2）]/n=(-1/n)[-1/(n-1）][-1/(n-2）]…（-1/3）[N⑵-N⑴]=(-1）^n/n!

因此

N(n-1）-N(n-2）=(-1）^(n-1）/(n-1）！

N⑵-N⑴=(-1）^2/2!

相加，可得

N(n)=(-1）^2/2!+…+(-1）^(n-1）/(n-1）！+(-1）^n/n!

因此

M(n)=n![(-1）^2/2!+…+(-1）^(n-1）/(n-1）！+(-1）^n/n!]

可以得到

错排公式为M(n)=n！（1/2!-1/3!+…..+(-1）^n/n!)

2，用容斥原理：

正整数1、2、3、……、n的全排列有n！种，其中第k位是k的排列有（n-1）！，当k取1、2、3、……、n时，共有n*（n-1）！种排列，由于是错排，这些排列应排除，但是此时把同时有两个数不错排的排列多排除了一次，应补上；在补上时，把同时有三个数不错排的排列多补上了一次，应排除；……；继续这一过程，得到错排的排列种数为

M(n)=n!-n!/1!+n!/2!-n!/3!+…+(-1）^n*n!/n!=sigma(k=2~n) (-1）^k*n!/k!

即M(n)=n![1/0!-1/1!+1/2!-1/3!+1/4!+..+(-1）^n/n!]

注：sigma表示连加符号，（k=2~n）是连加的范围

#include<iostream>
#include<cstdio>
using namespace std;

__int64 m[25];    //记录错排个数
double ans[25];

double f(double n)
{
    if(n==0||n==1)
        return 1;
    return n*f(n-1);
}

int main()
{
    m[1]=0;
    ans[1]=0.0;
    m[2]=1;
    ans[2]=1/2.0;
    for(int i=3;i<=23;i++)
    {
        m[i]=(i-1)*(m[i-1]+m[i-2]);
        ans[i]=m[i]/f(i);
    }
    int c;
    int n;
    cin>>c;
    while(c--)
    {
        scanf("%d",&n);
        printf("%.2lf",100*ans[n]);
        cout<<"%"<<endl;
    }
    return 0;
}

java版

import java.util.*;
//import java.text.DecimalFormat;
//import java.text.NumberFormat;
import java.text.*;

public class Main
{
	public static void main(String args[])
	{
		Scanner cin = new Scanner(System.in);
		long m[]=new long[25];
		double ans[]=new double[25];
		m[1]=0;
		ans[1]=0.0;
		m[2]=1;
		ans[2]=1/2.0;
		for(int i=3;i<=23;i++)
		{
			m[i]=(i-1)*(m[i-1]+m[i-2]);
			ans[i]=m[i]/f(i);
		}
		int c;
		int n;
		c=cin.nextInt();
		while(c-->0)
		{
			n=cin.nextInt();
			System.out.print(String.format("%.2f",100*ans[n]));  
			System.out.println("%");
		}
	}
	static double f(double n)
	{
		if(n==0||n==1)
			return 1;
		return n*f(n-1);
	}
}