HDUOJ-----2068RPG的错排

RPG的错排

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

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专场练习赛

 

错排公式  fn=（n-1）*（fn-1+fn-2）；

 

此题是错排公式的一个扩展，就是第K错牌的和...

特别要注意算法的格式。。。

这里介绍一下比较规范的组合数的求解

/*@coder 龚细军*/

int  cm( int n ,int c)
{
  //先比较 c=min{n-c,c}
  if(*c>n) c=n-c;  //由于cm(n,c)==cm(n,n-c);
  int i,j,ret=;
  for(i=,j=;i<n;i++)
 {
   ret*=(n-i);
   while(j<=c&&(j%c==))
   {
     ret/=j;
      j++;
   }
 }
return ret;
}

//时间复杂度0(n*c)

此题的代码。。。。简答了啦。。

 #include<stdio.h>
 int work(int n,int c)    /*组合数求c(n,c)*/
 {
   if(n<*c) c=n-c;   //c(n,c)=c(n,n-c)
   int i,j,s=;
   for(i=,j=;i<c;i++)
   {
     s*=(n-i);  //精简一下避免出现超出int型
     for( ; j<=c&&s%j==;++j)    s/=j;
   }
   return s;
 }
 int main()
 {
     int i,n;;
     _int64 sav[]={,,},ans;
      for(i=;i<;i++)
       sav[i]+=(i-)*(sav[i-]+sav[i-]);
     while(scanf("%d",&n),n)
     {
         ans=;
         for(i=(n+)/;i<=n;i++)
         {
             ans+=work(n,i)*sav[n-i];
         }
         printf("%I64d\n",ans);
     }
     return ;
 }