E - 不容易系列之(4)――考新郎 错排数公式

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

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

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

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

Input输入数据的第一行是一个整数C,表示测试实例的个数，然后是C行数据，每行包含两个整数N和M(1<M<=N<=20)。 
Output对于每个测试实例，请输出一共有多少种发生这种情况的可能，每个实例的输出占一行。 
Sample Input

2

2 2

3 2

Sample Output

1

3

D(n) = (n-1) [D(n-2) + D(n-1)]
排公式的原形为D(n) = n! (1/0! - 1/1! + 1/2! - 1/3! - ..... + (-1)^n/n!)，当n很大时计算就很不方便。一个供参考的简化后的公式是D(n) = [n!/e+0.5] ，其中e是自然对数的底，[x]为x的整数部分。

#include <iostream>

#include <vector>

#include <algorithm>

#include<string>

#include<cmath>

#include<cstdio>

#include<cstring>

using namespace std;

#define MAXN 21

typedef long long LL;

/*

在N对中有M对找错了

错排数公式 Dn = (n-1)*(Dn-1 + Dn-2)

*/

LL D[MAXN],A[MAXN];

void init()

{

    D[] = ;D[] =;D[] = ;

    A[] = A[] = ;A[] =;A[] = ;

    for(int i=;i<MAXN;i++)

    {

        A[i] = i*A[i-];

        D[i] = (i-)*(D[i-]+D[i-]);

    }

}

void solve(LL n,LL m)

{

    LL ans;

    ans = A[n]/A[m]/A[n-m]*D[m];

    //ans = floor(ans*10000+0.5)/100;

    printf("%lld\n",ans);

}

int main()

{

    init();

    int c,n,m;

    scanf("%d",&c);

    while(c--)

    {

        scanf("%d%d",&n,&m);

        solve(n,m);

    }

}