错位排序公式及理解

今天集训问错排是什么被hszx教练嘲讽了,痛下决心学一下竟然还挺简单的


错位排序递推公式:

设f[i]为i个数错位排序 (任意1<=i<=n a[i]!=i)

f[0]=1,f[1]=0;

f[i]=(f[i-1]*f[i-2])*(i-1)  (i>=2)

公式理解:

情况1:插入第i个元素时,前i-1个已经错位排好,则选择其中任意一个与第i个互换一定满足要求,选择方法共i-1种,前i-1位错排f[i-1]种,记f[i-1]*(i-1)

情况2:插入第i个元素时,前i-1个中恰有一个元素a[j]使得a[j]=j,其他i-2个错位排好,则将i与j交换,j在i-2位中的插入共i-1种,前i-2位错排f[i-2]种,记f[i-2]*(i-1)

以上两种情况求和可得

f[i]=(f[i-1]*f[i-2])*(i-1)  (i>=2)

递推代码

long long cp[maxn];
inline void get_cp()
{
	cp[0]=1ll;
	for(int i=2;i<=n;i++)
		cp[i]=(i-1)*(cp[i-2]+cp[i-1])%mod;
}
模板题:bzoj4517

附代码(写慢了,不预处理阶乘逆元能降一个log QAQ)

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define maxn 1000005
#define mod 1000000007
int n=maxn,m,T;
long long stair[maxn];
long long invs[maxn];
long long cp[maxn];
long long ans[maxn];
inline long long qpow(long long b,long long t)
{
	long long ans=1ll;
	while(t)
	{
		if(t&1) ans=(ans*b)%mod;
		t>>=1;
		b=(b*b)%mod;
	}
	return ans;
}
inline long long inv(long long x)
{
	return qpow(x,mod-2);
}
inline void init()
{
	stair[0]=1ll;
	invs[0]=1ll;
	for(int i=1;i<=n;i++)
	{
		stair[i]=stair[i-1]*i%mod;
		invs[i]=inv(stair[i]);
	}
	cp[0]=1ll;
	for(int i=2;i<=n;i++)
		cp[i]=(i-1)*(cp[i-2]+cp[i-1])%mod;
}
inline long long comb(int a,int b)
{
	return stair[a]*invs[b]%mod*invs[a-b]%mod;
}
int main()
{
 	scanf("%d",&T);
 	init();
 	while(T--)
 	{
 		scanf("%d%d",&n,&m);
 		printf("%lld\n",comb(n,m)*cp[n-m]%mod);
 	}
}



阅读更多
版权声明:本文为博主原创文章,转载请注明源网址blog.csdn.net/leo_h1104 https://blog.csdn.net/Leo_h1104/article/details/51508384
文章标签: 组合数学 OI
个人分类: 算法 组合数学 BZOJ
上一篇后缀数组模板
下一篇Link-Cut-Tree 的基本姿势
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭
关闭