hdu2068 RPG的错排

最新推荐文章于 2020-02-12 21:57:21 发布

钳子夹夹夹

最新推荐文章于 2020-02-12 21:57:21 发布

阅读量1.4k

点赞数

分类专栏： ACM 文章标签： hdu2068 hdu2068 RPG的错排错排问题错排组合数学

本文链接：https://blog.csdn.net/qq_31186409/article/details/50571182

版权

ACM 专栏收录该内容

20 篇文章 0 订阅

订阅专栏

RPG的错排

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

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

Sample Output

1
1

错排列，组合问题

错排公式：D(n) = n! [(-1)^2/2! + … + (-1)^(n-1)/(n-1)! + (-1)^n/n!] ， D(1) = 0， D(2) = 1

D(n)为错排列个数，n为元素个数

本题需要求0～n/2之间错排列的总个数

所以， D(1)*c(1, n) + ……+D(k)*c(k, n) +……+ D(n)*c(n, n)

#include <iostream>
typedef long long LL;

using namespace std;

LL f(int i, int n)
{
    LL ret = 1;
    while (i++ < n)
    {
        ret *= i;
    }
    return ret;
}

LL C(int k, int n)            //求C(k, n)
{
    LL ret = 1;
    for (int i = n-k+1; i <= n; i++)
        ret *= i;
    for (int i = 2; i <= k; i++)
        ret /= i;
    return ret;
}

int main()
{
    int n;
    LL sta[15] = {0};
    sta[1] = 0;
    sta[2] = 1;
    for (int i = 3; i < 14; i++)
    {
        for (int j = 2; j <= i; j++)
        {
            if (j & 1)
                sta[i] -= f(j, i);
            else
                sta[i] += f(j, i);
        }
    }
    while (cin >> n, n)
    {
        LL ans = 1;
        for (int i = 1; i <= n/2; i++)
            ans += C(i, n)*sta[i];
        cout << ans << endl;
    }
    return 0;
}