hdu 1465 全错排 or 容斥 or 二项式反演

问题：n封信对应n个信封 求恰好全部装错了信封的方案数

解法一：全错排公式

要装第i封信的时候，先把前i-1个信全装错信封，然后随便选其中一个与第i封信交换，有i-1种选法
那么dp[i] = (i-1) * dp[i-1]
(即把新加入的一封和之前的任一一封交换，所得到的必然是错排。)
但是还有一种情况
要装第i封信的时候，先从i-1封信中任选i-2个信把他们全装错信封，然后把剩下的那个信与第i个交换，从i-1封信中任选i-2个信有i-1种选法
(即之前的 n - 1 封中有一封没有错排，把这封与新加入的一封交换进行错排。)
那么dp[i] = (i-1) * dp[i-2]
两个式子联合起来
就是那么dp[i] = (i-1) * (dp[i-1] + dp[i-2])
这就是全错排公式，递推，递归都可以做

#include<cstdio>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
int main()
{
	int n;
	ll ans=0;
	ll dp[50];
	dp[2]=1,dp[1]=0;
	for(int i=3;i<=23;i++)
	{
		dp[i]=(i-1)*(dp[i-1]+dp[i-2]);
	}
	while(cin>>n)
	{
		cout<<dp[n]<<endl; 
	}
	return 0;
}

解法二：容斥原理

首先，所有装信的总数是n!
（在n中任选一个信封放进一封信，然后在剩下的n-1中任选一个信封放进一封信，以此类推，所以是n*(n-1)*(n-2)... = n!）
假设
A1表示1封信装对信封，数量是(n-1)! （只有n-1个位置可以乱放）
A2表示2封信装对信封，数量是(n-2)! （只有n-2个位置可以乱放）
...
An表示n封信装对信封，数量是1
那么这题的答案就是
n! - C(n, 1)*|A1| + C(n, 2)*|A2| - C(n, 3)*|A3| + ... + (-1)^n * C(n, n)*|A4|

把C(n,m)用

表示 化简

n! - n! / 1! + n! / 2! - n! / 3! + ... + (-1)^n * n! / n!

#include<cstdio>
typedef long long LL;
int n, flag;
LL fac[25];
LL ans;
void init(){
    fac[0] = 1;
    for(int i = 1; i <= 20; i ++) fac[i] = fac[i-1] * i;    
}
int main(){
    init();
    while(~scanf("%d", &n)){
        ans = fac[n];
        flag = -1;//容斥的符号变化
        for(int i = 1; i <= n; i ++){
            ans += flag * fac[n] / fac[i];
            flag = -flag;  
        }
        printf("%I64d\n", ans);
    }
}

解法三：二项式反演

设g(i)表示正好有i封信装错信封

那么全部的C(n, i)*g(i)加起来正好就是所有装信的情况，总共n!种情况

n! = Σ C(n, i)*g(i) (i从0到n)

那么f(n) = n!，所以f(x) = x!

那么我们要求g(n)

根据公式

 

g(n) = Σ (-1)^(n-i) * C(n, i) * f(i)  (i从0到n)

#include<cstdio>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
ll fac[30];
int main()
{
	int n;
	fac[1]=1;
	fac[0]=1;
	for(int i=2;i<=25;i++)
	{
		fac[i]=fac[i-1]*i;
	}
	while(cin>>n)
	{
		int flag=(n&1)?-1:1;
		ll ans=0;
		for(int i=0;i<=n;i++)
		{
			ans+=flag*fac[n]/fac[n-i];
			flag=-flag;
		}
		cout<<ans<<endl;
	}
	return 0;
}