HDU 2068 RPG的错排 (错排+组合数)

RPG的错排

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

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

 

POINT：

用错排递归求出错排数。

答对一半或以上就算过关，我们就求答错(0-一半)的数量。

就是错排数加上组合数。都不大，但是组合数别直接用阶乘求，会爆LL。

如果组合数很大，肯定要取模，取模一般用质数，我们就可以用费马小定理和Lucas来求大的组合数，显然这题不需要。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define  LL long long
LL k[30];
LL cp[30];
void init()//C(m,n)=C(m-1,n)+C(m-1,n-1)；
{
    k[0]=1;
    for(int i=1;i<=25;i++)
    {
        k[i]=k[i-1]*i;
    }
    cp[1]=0;
    cp[2]=1;
    for(int i=3;i<=25;i++)
    {
        cp[i]=(i-1)*(cp[i-1]+cp[i-2]);
    }
}
LL c(int n,int m){
    int a=1;
    if(m==0)
        return 1;
    for(int i=1;i<=m;i++){
        a=a*(n-i+1);
        a=a/i;//求组合数的一个技巧
    }
    return a;
}
int main()
{
    init();
    int n;
    while(~scanf("%d",&n)&&n)
    {
        int x=n/2;
        LL ans=1;
        for(int i=0;i<=x;i++)
        {
            ans+=cp[i]*c(n,i);
        }
        printf("%lld\n",ans);
    }
    return 0;
}