HDU打怪升级之RPG的错排(2068)

问题链接：http://acm.hdu.edu.cn/showproblem.php?pid=2068

问题描述：

今年暑假杭电ACM集训队第一次组成女生队,其中有一队叫RPG,但做为集训队成员之一的野骆驼竟然不知道RPG三个人具体是谁谁。RPG给他机会让他猜猜，第一次猜：R是公主，P是草儿，G是月野兔；第二次猜：R是草儿，P是月野兔，G是公主；第三次猜：R是草儿，P是公主，G是月野兔；......可怜的野骆驼第六次终于把RPG分清楚了。由于RPG的带动，做ACM的女生越来越多，我们的野骆驼想都知道她们，可现在有N多人，他要猜的次数可就多了，为了不为难野骆驼，女生们只要求他答对一半或以上就算过关，请问有多少组答案能使他顺利过关。

输入：

 输入的数据里有多个case，每个case包括一个n，代表有几个女生，(n<=25)，n=0输入结束。

栗子输入：

1
2
0

栗子输出：

1
1

 

问题分析：

根据题目描述可知需要找到n个女生当中能够猜对一半或以上人数的名字，才能够过关。所以可以等同于看成在n个女生当中，有小于或等于一半的人数的真实名字与猜测名字不符合的错排问题。

对于错排问题，假设有n个人，每个人按进行编号，然后进行排队。则n个人错排表示这n个人编号与所站的位置不同。设f(n)为n个人错排的情况总数。考虑下面两种情况：
1. n-1个人错排，1个人编号和所占的位置相同，要使得能够达到n个人错排，只需要剩下的1个人与n-1个人中的某个人交换位置即可。则该种情况的总数为f(n-1)*(n-1)。
2. n-1个人中有1个人编号和位置相同，其余的n-2个人错排，而剩下的1个人编号和所占位置不同，要使得能够达到n个人错排，只需要n-1中那个编号与位置相同的人与剩下的1个人交换位置即可。则该种情况的总数为f(n-2)*(n-1)。
所以 f(n) = (n-1)*(n-1) + f(n-1)*(n-2)。

对于n中人有m人错排的问题，则等同于从n个人中选取m个人，然后这m个人错排，其余n-m个人未错排，则情况的总数为*f(m)。

对于该问题，n个人中有小于或等于一半的人数错排等同于0、1、2、3、…、n/2个人错排情况的总和。

 

参考代码：

#include<iostream>
using namespace std;

long long C(int n,int m)
{
	long long a,b;
	a=1;b=1;
	for(int i=0;i<m;i++)
	{
		a*=(n-i);
		b*=(m-i);		
	}
	return a/b;
}

int main()
{
	int n,m;
	long long a[26];
	long long total;
	a[0]=1; a[1]=0;	a[2]=1;
	for(int i=3;i<=25;i++)
		a[i] = (i-1)*(a[i-1]+a[i-2]);
	
	while(cin>>n)
	{
		if(!n)
			break;
		total = 0;
		m = n/2;
		for(int i=0;i<=m;i++)
			total += a[i]*C(n,i);
		cout<<total<<endl;
	}
	return 0;
} 

学习记录：

对于求的函数，如果先分别算出分子、分母，然后再相除，需要考虑在相乘时，分子和分母是否会溢出。如果是先相除后再相乘，需要考虑在相除的时候会不会出现算出整数结果、精度丢失而导致结果错误的问题。

1.
long long C(int n,int m)
{
	long long a,b;
	a=1;b=1;
	for(int i=0;i<m;i++)
	{
		a*=(n-i);
		b*=(m-i);		
	}
	return a/b;
}

2.
long long C(int n,int m)
{
	long long res=1;
	for(int i=0;i<m;i++)
	{
		res*=(n-i);
		res/=(i+1);
	}
	return res;
}