HDU 2068 RPG的错排

RPG的错排

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

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

 

利用错排和组合公式

错排：A[i]=(i-1)X(A[i-1]+A[i-2]);A[0]=0;A[1]=1;A[2]=2;

组合：C(m,n)=m!/(m-n)!/n!

import java.util.*;
import java.io.*;
import java.math.BigInteger;

public class Main {
    public static void main(String[] args) {
        Scanner sc = new Scanner(new BufferedInputStream(System.in));
        while (sc.hasNextInt()) {
            int n = sc.nextInt();
            int m=0;
            BigInteger sum=BigInteger.ONE;
            if (n == 0)
                System.exit(0);
            
            
            for(int i=0;i<=n/2;i++){
                sum=sum.add(zhuHe(n).divide(zhuHe(n-i)).divide(zhuHe(i)).multiply(cuoPai(i)));
                
            }
            System.out.println(sum);
        }
    }
    public static BigInteger zhuHe(int m){
        BigInteger big[]=new BigInteger[25+1];
        big[0]=BigInteger.ONE;
        big[1]=BigInteger.ONE;
        for(int i=1;i<=m;i++){
            big[i]=big[i-1].multiply(BigInteger.valueOf(i));
        }
        return big[m];
    }
    public static BigInteger cuoPai(int m){
        BigInteger big[]=new BigInteger[26];
        big[0]=BigInteger.ZERO;
        big[1]=BigInteger.ZERO;
        big[2]=BigInteger.valueOf(1);
        for(int i=3;i<=m;i++){
            big[i]=BigInteger.valueOf((i-1)).multiply(big[i-1].add(big[i-2]));
        }
        return big[m];
    }
}