2018.11.07 hdu1465不容易系列之一（二项式反演）

传送门
其实标签只是搞笑的。
没那么难。
二项式反演只是杀鸡用牛刀而已。
这道题也只是让你$n\le 20$的错排数而已。
还记得那个$O(n)$的递推式吗？
没错那个方法比我今天用的要快一些。

---

言归正传。
回忆一下二项式反演的式子：
$f_n={\sum }_{i=0}^n\left(\genfrac{}{}{0ex}{}{n}{i}\right)g_i$
=>$g_n={\sum }_{i=0}^n((-1{)}^i\left(\genfrac{}{}{0ex}{}{n}{n-i}\right)f_i)$
证明很简单。
只用把第一个式子成立的条件带到第二个等式的右边就可以了。
然后这道题怎么用呢？
我们令$f_i$表示$i$张牌任意排列的总方案数。
$g_i$表示$i$张牌全部错排的方案数。
那么根据分类计数的原理显然有：
$f_n={\sum }_{i=0}^n{g}_i=n!$
于是$g_n={\sum }_{i=0}^n((-1{)}^i\left(\genfrac{}{}{0ex}{}{n}{i}\right)f_i)={\sum }_{i=0}^n((-1{)}^i\frac{n!}{(n-i)!})$
做完了。
代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=21;
ll fac[N];
int n;
int main(){
	fac[0]=1;
	for(int i=1;i<=20;++i)fac[i]=fac[i-1]*i;
	while(~scanf("%d",&n)){
		ll ans=0,tmp=1;
		for(int i=0;i<=n;++i,tmp*=-1)ans+=tmp*fac[n]/fac[i];
		cout<<ans<<'\n';
	}
	return 0;
}