ACM学习历程—HDU2068 RPG的错排(组合数学)

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输入结束。

 

Output

对应每组数据输出最小移动距离。

 

Sample Input

1

2

0

 

Sample Output

1

1

假设i个人在他本来位置，其余人错排的种数是f[i]，那么题目要求的就是所有大于等于(n+1)/2的f[i]的和，n+1是为了对奇数偶数情况统一。

假设k个人错排是p[k]，

那么就是n个人先取出i个人在自己位置C(n, i)，其余人再错排p[n-i]，然后控制i的范围就OK了。

这里需要注意的是由于题目没有模的情况，所以所有数可能会很大，所以在求p和c的时候最好是相邻项间递推。

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <set>
#include <map>
#include <queue>
#include <string>
#include <algorithm>
#define LL long long

using namespace std;

int n;
LL p[], c[], ans;

void init()
{
    p[] = ;
    p[] = ;
    p[] = ;
    for (int i = ; i < ; ++i)
        p[i] = (i-)*(p[i-]+p[i-]);
}

void cal()
{
    c[] = ;
    for (int i = ; i <= n; ++i)
        c[i] = c[i-]*(n-i+)/i;
}

void work()
{
    ans = ;
    int half = (n+)/;
    for (int i = ; i+half <= n; ++i)
        ans += c[i+half]*p[n-i-half];
    printf("%I64d\n", ans);
}

int main()
{
    //freopen("test.in", "r", stdin);
    init();
    while (scanf("%d", &n) != EOF && n)
    {
        cal();
        work();
    }
    return ;
}