HDU 2068：RPG的错排（错排问题）

HDU 2068：RPG的错排

题目链接

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 Input1
2
0
Sample Output1
1

解决这个问题首先得了解题名中所说的错排问题：假设现在有N个数字，按照一定的次序排列，现在要使得每个数字都不放在原来的位置，有多少种摆法。那么首先将第一个数字错排，那么它将有N-1种选择，那么现在假设放在了K上，现在再处理原本K位置的数字，假设K放在了第一个位置，那么剩下的问题就等同于将剩下的N-2个数字错排，假设K没有放在第一个位置，那么剩下的问题就等同于错排N-1个数字，因为K的位置虽然已经被占用了，但是多出来一个第一个位置的空位，现在第一个位置的空位的效果是等同于K的那个已经被占用的空位的，都是除了K以外都可以放的空位，效果等同。所以这里就可以推出递推式x[N]= (N-1) * ( x[N-1] + x[N-2] )
很明显，这样定义的递推式是有限制的，那么就是N>=3，但是前两个是错排是很好算的，x[1]=0,x[2]=1;

那么现在我们就可以把这样的方法运用到题目中去了，题目要求答对一半或以上算过关，那么相反我们错排一半或以下就自然可以可以过关了，当然此题数据量有些大，最后几个可能会超范围，如果不用高精度计算的话可以用计算器算好打表。当然此题还有一个很重要的地方，x[0]=1，一个都不错排也是一种正确情况！
代码：

#include<iostream>
#include<cstdio>
long long x[30];
using namespace std;
int main()
{
    int n;
    while(cin>>n)
    {
        if(n==0)
            break;
        if(n==25)
             cout<<986478419687006<<endl;
        else if(n==24)
            cout<<515971556258510<<endl;
        else if(n==23)
             cout<<21498814787768<<endl;
        else if(n==22)
            cout<<11293696692640<<endl;
        else if(n==21)
            cout<<513349865698<<endl;
        else{
        int m;
        m=n/2;
        long long sum=0,r,t;
        for(int i=0;i<=m;++i)
        {
            r=t=1;
            x[0]=1,x[1]=0,x[2]=1;
            for(int j=3;j<=i;++j)
            {
                x[j]=(j-1)*(x[j-1]+x[j-2]);
            }
            if(n>2)
            {
                for(int k=n;k>=i+1;--k)
                {
                    r*=k;
                }
                for(int k=1;k<=n-i;++k)
                {
                    t*=k;
                }
                x[i]*=r/t;
            }
            sum+=x[i];
        }
        cout<<sum<<endl;}
    }
}