简单错排——HDU-2049-考新郎

链接：http://acm.hdu.edu.cn/showproblem.php?pid=2049

题目：

国庆期间,省城HZ刚刚举行了一场盛大的集体婚礼,为了使婚礼进行的丰富一些,司仪临时想出了有一个有意思的节目,叫做"考新郎",具体的操作:
首先,给每位新娘打扮得几乎一模一样,并盖上大大的红盖头随机坐成一排;
然后,让各位新郎寻找自己的新娘.每人只准找一个,并且不允许多人找一个.
最后,揭开盖头,如果找错了对象就要当众跪搓衣板...
看来做新郎也不是容易的事情...
假设一共有N对新婚夫妇,其中有M个新郎找错了新娘,求发生这种情况一共有多少种可能. 

Input

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

 Output

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

Sample Input

22 23 2

 Sample Output

13

本来不知道错排是什么，看到题解公式真的懵，看来我的数学是真的差。QAQ

错排：考虑一个有n个元素的排列，若一个排列中所有的元素都不在自己原来的位置上，那么这样的排列就称为原排列的一个错排。

公式推导：

第一步，把第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)]

又因为，题目要求总共有m对，其中m对找错了，也就是m对进行错排，剩下的（n-m）对找对了，答案也就是从n中挑出来m对，然后将其错排，也就是：D（m）*C(n,m)=(m-1)*(D(m-2)+D(m-1))*n*(n-1)*(n-2)*...*(n-m+1)/(m!)

代码：

#include<bits/stdc++.h>

using namespace std;

int main(){
int t, n, m;
long long C;//排列组合
long long D[22]={0, 0, 1};//错排//比较方便的赋值
for(int i = 3; i <= 20; i++)
D[i]=( D[i-1] + D[i-2] )*( i-1 );//预处理
cin>>t;
while(t-- ) {
cin>>n>>m;
C=1;
for(int i=n;i>=n-m+1;i--)
C*=i;
for(int i=m;i>1;i--)
C/=i;
cout<<D[m]*C<<endl;
}
return 0;
}