哈理工OJ 1191 相亲（容斥原理—错排公式）

http://acm.hrbust.edu.cn/index.php?m=ProblemSet&a=showProblem&problem_id=1191
相亲
Time Limit: 1000 MS Memory Limit: 65536 K
Total Submit: 293(85 users) Total Accepted: 105(75 users) Rating: Special Judge: No
Description
为了追求邻家小妹，hrbust 隆重的举行了 我爱美女大型相亲活动，为了避免尴尬，找了N对幸福的小情侣，并且编号，首先，我们把清纯可爱的MM都躲在后台，在舞台上我们放置了一个大箱子，箱子里面有女生的编号，让他们的hero来抽取她们。如果抽到自己的princess的编号，我们就豪华大礼包—免挂金牌。
可是主持人想知道M个hero找错了princess一共有多少可能。
Input
输入数据的第一行是一个整数C,表示测试实例的个数，然后是C行数据，每行包含两个整数N和M(1

3

思路解析：
一共有m个人找错了自己的对象，那么有n-m个人找对了自己的对象，所以找对自己对象的可能数为C(n-m,n)，下面就要求m个人找错了总计的对象的种数，这明显是一个典型的伯努利装错信封问题。
伯努利装错信封问题可以由容斥原理来解，就是求至少有一个信封装对的对立事件。

下面就很简单了。

下面是AC代码：

#include<cstdio>
#define ll long long
#include<algorithm>
using namespace std;

ll calculatec(ll n,ll m)
{
    ll ans=1;
    for(ll i=n;i>=n-m+1;i--)
    {
        ans*=i;
    }
    for(ll i=1;i<=m;i++)
    {
        ans/=i;
    }
    return ans;
}

ll calculatea(ll n)
{
    ll ans=1;
    for(ll i=1;i<=n;i++)
    {
        ans*=i;
    }
    return ans;
}

int main()
{
    ll n,m,t;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&m);
        ll sum1=1,sum2=0;
        sum1=calculatec(n,n-m);
        ll aa=calculatea(m);
        for(ll i=1;i<=m;i++)
        {
            if(i%2==0)
            {
                sum2-=(aa/calculatea(i));
            }
            else
            {
                sum2+=(aa/calculatea(i));
            }
        }
        printf("%lld\n",sum1*(aa-sum2));
    }
    return 0;
}