hdu2068-RPG的错排

原文链接：hdu2068-RPG的错排

原文：

RPG的错排

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

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

 

试题分析：全错排和组合的结合

AC代码：

#include "stdio.h"
#define LL long long
LL fac(long long a);
LL dp[25];
LL fac1(LL a,LL b);
int main(){
    LL i,n,s;
    LL k,sum;
    dp[0]=0;
    dp[1]=1;
    for(i=2;i<20;i++)
        dp[i]=i*(dp[i-2]+dp[i-1]);

    while(~scanf("%I64d",&n) && n){
        sum=0;
        k=n/2;
        if(n%2==1) k++;
        while(n-k>=2){
            sum+=fac1(n,n-k)*dp[n-k-1];
            k++;
        }
        printf("%I64d\n",sum+1);
    }
    return 0;
}
LL fac(long long a){
    int i;
    LL sum=1;
    for(i=1;i<=a;i++){
        sum*=i;
    }
    return sum;
}
LL fac1(LL a,LL b){
    long long i,k,t,sum,sum1;
    k=a-b;
    if(k>b){
        t=k;
        k=b;
        b=t;
    }
    sum=1;
    for(i=a;i>b;i--){
        sum*=i;
    }
    sum1=fac(k);
    return sum/sum1;
}