HDU-2048 神、上帝以及老天爷

 题目链接：HDU-2048

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

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-1个签中随便抽取，有n-1种抽发......最后一个人只有最后一种抽法

分母就是n！

就下来看分母怎么求

就这道题来说，我们可以看成有n个人进行抽签，为了方便叙述，我们就看最后一个人的抽签状态，首先最后一个人不会抽到自己的签，那他就有n-1种抽法

还是觉着网上普遍的叙述方式比较好说，就是第n个人如果抽到第k位置的人的名字，那就看第k个人抽什么签，如果说第k个人抽的是第n个人的签，那这两个人是不是可以除外了，因为他们抽的都是对方的签，不会再影响剩下n-2个人的抽签结果，即f(n-2)

再看另一种情况，如果说第k个人抽的不是第n个人的签，那看第n个人的签，它不会被第k个人抽到，因为上面已经讨论过相互抽到对方签的情况，也就是说第k个人不能抽第n个人的签，这是不是就像第n个人不能抽第n个人的签、自己不能抽自己的签
第n个人的签在这种情况下就相当于位置k的人签，也就是说这变成了n-1个人的游戏，即f(n-1)

于是就能得到递归关系(n-1)*(f(n-1)+f(n-2))

递归边界是只有一个人和只有两个人的情况

只有一个人，只会抽到自己的签，结果为0

两个人的话，只有一种情况，结果为1 

代码很好写，写个递归函数，调用就行了

我下面写一下我遇到的问题

用C++写，输出小数点后固定的位数，需要调用setprecision()函数，引用的头文件是iomanip

另一个问题是变量的类型问题，我一开始给求阶乘的变量设的是int类型的，结果超范围
后来改成double类型的就可以通过

 

AC代码：

#include<iostream>
#include<iomanip>
using namespace std;
double f(int n){
//	if(n==0){
//		return 0; 
//	}
	if(n==1||n==2){
		return n-1;
	}
	else if(n==3){
		return 2;
	}
	else 
	return (n-1)*(f(n-1)+f(n-2));
} 
int main(){
	int C,n;
	cin>>C;
	for(int i=0;i<C;i++){
		cin>>n;
//		if(n>=7){
//			cout<<"36.79%"<<endl;
//			continue; 
//		}
		double d=1;
		for(int j=1;j<=n;j++){
			d*=j;
			cout<<d<<endl;
		}
		//cout<<f(n)*100<<endl;
		//cout<<f(n)<<endl;
		cout<<fixed<<setprecision(2)<<(f(n)*100)/d<<'%'<<endl;
	}
}