HDU 2049 不容易系列之(4)——考新郎

啊。。。

突然发现我只有讲课才会写博客

我以后会改进的

---

原题戳

别戳了，我粘过来了

---

Problem Description

国庆期间,省城HZ刚刚举行了一场盛大的集体婚礼,为了使婚礼进行的丰富一些,司仪临时想出了有一个有意思的节目,叫做”考新郎”,具体的操作是这样的:

首先,给每位新娘打扮得几乎一模一样,并盖上大大的红盖头随机坐成一排;
然后,让各位新郎寻找自己的新娘.每人只准找一个,并且不允许多人找一个.
最后,揭开盖头,如果找错了对象就要当众跪搓衣板…

看来做新郎也不是容易的事情…

假设一共有N对新婚夫妇,其中有M个新郎找错了新娘,求发生这种情况一共有多少种可能.

Input

输入数据的第一行是一个整数C,表示测试实例的个数，然后是C行数据，每行包含两个整数N和M($1<N<M<20$)。

Output

对于每个测试实例，请输出一共有多少种发生这种情况的可能，每个实例的输出占一行。

Sample Input

2
2 2
3 2

Sample Output

1
3

---

这是一道组合+错排的公式模板题。。。
总之就是很水很水很短很短就对了。。。
我也不知道为什么要讲它。。。

---

当n个编号元素放在n个编号位置，元素编号与位置编号各不对应的方法数用D(n)表示，那么D(n-1)就表示n-1个编号元素放在n-1个编号位置，各不对应的方法数，其它类推.
第一步，把第n个元素放在一个位置，比如位置k，一共有n-1种方法；
第二步，放编号为k的元素，这时有两种情况：
⑴把它放到位置n，那么，对于剩下的n-1个元素，由于第k个元素放到了位置n，剩下n-2个元素就有D(n-2)种方法；
⑵第k个元素不把它放到位置n，这时，对于这n-1个元素，有D(n-1)种方法；
综上得到
D(n) = (n-1) [D(n-2) + D(n-1)]
特殊地，D(1) = 0, D(2) = 1.

下面通过这个递推关系推导通项公式：
为方便起见，设D(k) = k! N(k), k = 1, 2, …, n,
则N(1) = 0, N(2) = 1/2.
n ≥ 3时，n! N(n) = (n-1) (n-1)! N(n-1) + (n-1)! N(n-2)
即 nN(n) = (n-1) N(n-1) + N(n-2)
于是有
N(n) - N(n-1) = - [N(n-1) - N(n-2)] / n = (-1/n) [-1/(n-1)] [-1/(n-2)]…(-1/3) [N(2) - N(1)] = (-1)^n / n!.
因此
N(n-1) - N(n-2) = (-1)^(n-1) / (n-1)!,
N(2) - N(1) = (-1)^2 / 2!.
相加，可得
N(n) = (-1)^2/2! + … + (-1)^(n-1) / (n-1)! + (-1)^n/n!
因此
D(n) = n! [(-1)^2/2! + … + (-1)^(n-1)/(n-1)! + (-1)^n/n!].
此即错排公式。

---

嗯组合公式我就不讲了吧

$C^m_n=\frac {n!}{m!*(n-m)!}$

---

代码

#include<bits/stdc++.h>
using namespace std;

int c,n,m;
long long a[30],b[30];
long long ans;

int main()
{
    a[0]=1;a[1]=1;a[2]=2;
    for(int i=3;i<=20;i++)
        a[i]=a[i-1]*i;
    b[0]=0;b[1]=0;b[2]=1;b[3]=2;
    for(int i=4;i<=20;i++)
        b[i]=(i-1)*(b[i-1]+b[i-2]);
    scanf("%d",&c);
    while(c--)
    {
        scanf("%d%d",&n,&m);
        ans=(a[n]/a[m]/a[n-m])*b[m];
        printf("%I64d\n",ans);
    }
    return 0;
}

---

最后注意俩坑
一个是初始化
还有一个是long long （所以别想着把俩东西直接乘一块儿去。。。劲那。。。）