错排问题系列之 RPG的错排 hdoj 2068

/*
RPG的错排

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 7296    Accepted Submission(s): 2963

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
 

Source
RPG专场练习赛
 

Recommend
lcy   |   We have carefully selected several similar problems for you:  2065 2067 2066 2063 2062
 

Statistic | Submit | Discuss | Note
*/

/*

错排：

考虑一个有n个元素的排列，若一个排列中所有的元素都不在自己原来的位置上，那么这样的排列就称为原排列的一个错排

  a[n];

  a[0]=a[1]=0;

  a[2]=1;

  a[n]=(n-1)*(a[n-1]+a[n-2]);

组合：

 C(m,n)=(n!)/(m!*(n-m)!)

*/
#include<cstdio>
#include<string.h>
#include<stdlib.h>
#define M 1000
__int64 a[M];
__int64 C(int n,int m)
{
  __int64 sum1,sum2;
  int i;
  for(sum1=1,i=1;i<=m;i++)
     sum1*=i;
   for(sum2=1,i=n;i>=n-(m-1);i--)
     sum2*=i;
     return sum2/sum1;
}
int main()
{
 int i,j;
 memset(a,0,sizeof(a));
 a[0]=a[1]=0;
 a[2]=1;
 for(i=3;i<=25;i++)
   a[i]=(i-1)*(a[i-1]+a[i-2]);
  
   int n;
   while(scanf("%d",&n),n)
   {
       int m=n/2;
       __int64 sum=1;
       for(i=1;i<=m;i++)
          sum+=C(n,i)*a[i];
          printf("%I64d\n",sum);
   }
   return 0;
}