神、上帝以及老天爷


http://acm.hdu.edu.cn/showproblem.php?pid=2048



Time Limit: 2000/1000 MS (Java/Others)    

Memory Limit: 65536/32768 K (Java/Others)


三种思路:

1. N张字条的所有可能排列自然是N!(分母)。
现在的问题就是求N张字条的错排数f(N)(分子)。
首先我们考虑,如果前面N-1个人拿的都不是自己的票,即前N-1个人满足错排,现在又来了一个人,他手里拿的是自己的票。只要他把自己的票与其他N-1个人中的任意一个交换,就可以满足N个人的错排。这时有(N-1)*f(N-1)种方法。

Besides,我们考虑,如果前N-1个人不满足错排,而第N个人把自己的票与其中一个人交换后恰好满足错排。
这种情况发生在原先N-1人中,N-2个人满足错排,有且仅有一个人拿的是自己的票,而第N个人恰好与他做了交换,这时候就满足了错排。
又因为前N-1个人中,每个人都有机会拿着自己的票。所以有N-1种交换的可能。故这时有(N-1)*f(N-2)种方法。

综上所述:f(N)=(N-1)*[f(N-1)+f(N-2)] 

2. 另一种推导思路是先通过容斥定理直接计算出错排数f(n)=n!(1/2!-1/3!+…..+((-1)^n)/n!)

所以f(n)/n!-f(n-1)/(n-1)!=((-1)^n)/n!

两边同时乘上n!有f(n)-n*f(n-1)=(-1)^n

所以f(n-1)-(n-1)*f(n-2)=(-1)^(n-1)=(-1)*[f(n)-n*f(n-1)]

化简得f(n)=(n-1)*[f(n-1)+f(n-2)] 

3. 通过错位排列数的简化公式

来计算(参考维基百科


完整代码:

/*0ms,224KB*/

#include<cstdio>
#include<cmath>
const double e = exp(1.0);

double p[21];

int main()
{
	double temp = 1.0;
	int t, n, i;
	for (i = 2; i <= 20; ++i)
	{
		temp *= i;
		p[i] = round(temp / e) / temp;
	}
	scanf("%d", &t);
	while (t--)
	{
		scanf("%d", &n);
		printf("%.2f%%\n", p[n] * 100.0);
	}
	return 0;
}