HDU 2068 错排公式

HDU 2068
RPG的错排
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 17473 Accepted Submission(s): 7016

Problem Description
今年暑假杭电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输入结束。

Sample Input

1
2
0

Sample Output

1
1

Author
Rabbit

问题分析：
套用错排公式。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;

long long  a[30];

double cmn(int x,int y)
{
	double sum1=1,sum2=1,sum3=1,sum,z=x-y;
	for(int i=2;i<=x;i++)
	sum1*=i;
	for(int i=2;i<=y;i++)
	sum2*=i;
	for(int i=2;i<=z;i++)
	sum3*=i;
	sum=sum1/sum2/sum3;
	return sum;
}

int main()
{
	a[1]=0;
	a[2]=1;
	for(int i=3;i<=25;i++)
	{
		a[i]=(i-1)*(a[i-1]+a[i-2]);
	}
	int n;
	while(cin>>n&&n)
	{
         double sum=0;int l;
         if(n%2==0)
         l=n/2;
         else l=n/2+1;
        
		 for(int i=l;i<n;i++)
         {
         	sum+=cmn(n,i)*a[n-i];
		 }
		long long s= sum+1;
		cout<<s<<endl;
	}
	
}