[刷题笔记]HDU2049 错排公式

题目：HDU 2049 

地址：http://acm.hdu.edu.cn/showproblem.php?pid=2049

思路：组合数+错排

心得：

1. 所谓的错排公式 Dn=[n!/e+0.5]只是近似值！4次WA的代价啊！
2. 应该写成这样
3. 当n个编号元素放在n个编号位置,元素编号与位置编号各不对应的方法数用M(n)表示,那么M(n-1)就表示n-1个编号元素放在n-1个编号位置,各不对应的方法数,其它类推.
   第一步,把第n个元素放在一个位置,比如位置k,一共有n-1种方法;
   第二步,放编号为k的元素,这时有两种情况.1,把它放到位置n,那么,对于剩下的n-2个元素,就有M(n-2)种方法;2,不把它放到位置n,这时,对于这n-2个元素,有M(n-1)种方法;
   综上得到
   M(n)=(n-1)[M(n-2)+M(n-1)]
   特殊地，M(1)=0,M(2)=1
   
4. s[n] = {0, 0, 1}
   for (i=3;i<=20;i++) 
   s[i] = (i-1)*(s[i-1] + s[i-2]);
   n的错排数即为s[n]
   

代码：

#include<string>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<cmath>
using namespace std;
long long getNg(int m)
{
	int i = 2; 
	long long ans = m;
	for(;i<m;i++)
		ans *= i;
	return ans;
}
long long getC(int &n, int &m)
{
	int i = 1;
	long long ans = 1;
	for(;i<=m;i++)
		ans *= n--;
	return ans/getNg(m);
}
int main()
{
	int t, n, m, i; 
	long long s[21]={0,0,1};
	for (i=3;i<=20;i++) 
		s[i] = (i-1)*(s[i-1] + s[i-2]);//得到错排数
	cin>>t;
	while(t--)
	{
		cin>>n>>m;
		cout<<(getC(n, m)*s[m])<<endl;
	}
}