HDU2068(错排公式)

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

问提分析：随机抽出一半人，剩下的错排，得出多少种，然后猜对的人加一，剩下的错排数加上去，直到求出一到n/2的组合数与错排数的积的和即为答案。错排数推导公式：d[n]=(n-1)*(d[n-1]+d[n-2])

#include <iostream>
using namespace std;
__int64 c(int m, int n)//排列组合
{
	__int64 a, b, i;
	if (n > m / 2)n = m - n;//
	for (a = b = i = 1; i <= n; i++)
	{
		a *= (m - i + 1);
		b *= i;
	}
	return a / b;
}
int main()
{
	int n;
	
		__int64 a[26] = {1,0,1};//至少有1种
		for (int i = 3; i <= 13; i++)//求错排数
		{
			a[i] = (i - 1)*(a[i - 1] + a[i - 2]);
		}
		while (cin >> n && n)
		{
			__int64 ans=0;
			for (int i = 0; i <= n / 2; i++)
			{
				ans += c(n,n-i)*a[i];//求和
			}
			cout << ans << endl;
		}
	
}