BZOJ 4517: [Sdoi2016]排列计数 组合数 错排公式

4517: [Sdoi2016]排列计数

Time Limit: 60 Sec  Memory Limit: 128 MB
Submit: 1190  Solved: 727
[Submit][Status][Discuss]

Description

求有多少种长度为 n 的序列 A,满足以下条件:
1 ~ n 这 n 个数在序列中各出现了一次
若第 i 个数 A[i] 的值为 i,则称 i 是稳定的。序列恰好有 m 个数是稳定的
满足条件的序列可能很多,序列数对 10^9+7 取模。

Input

第一行一个数 T,表示有 T 组数据。
接下来 T 行,每行两个整数 n、m。
T=500000,n≤1000000,m≤1000000

Output

输出 T 行,每行一个数,表示求出的序列数

Sample Input

5
1 0
1 1
5 2
100 50
10000 5000

Sample Output

0
1
20
578028887
60695423

思路很简单

预处理阶乘、逆元和错排公式:D(n)=(n-1)(D(n-1)+D(n-2))D(1)=0,D(2)=1

ans=c(x,y)*er[x-y]      括弧:er是错排公式

省的大家百度 附上错排公式证明

n 个不同元素的一个错排可由下述两个步骤完成: 
第一步,“错排” 1 号元素(将 1 号元素排在第 2 至第 n 个位置之一),有 n - 1 种方法。 
第二步,“错排”其余 n - 1 个元素,按如下顺序进行。视第一步的结果,若1号元素落在第 k 个位置,第二步就先把 k 号元素“错排”好, k 号元素的不同排法将导致两类不同的情况发生:
1、 k 号元素排在第1个位置,留下的 n - 2 个元素在与它们的编号集相等的位置集上“错排”,有 f(n -2) 种方法;
2、 k 号元素不排第 1 个位置,这时可将第 1 个位置“看成”第 k 个位置(也就是说本来准备放到k位置为元素,可以放到1位置中),于是形成(包括 k 号元素在内的) n - 1 个元素的“错排”,有 f(n - 1) 种方法。据加法原理,完成第二步共有 f(n - 2)+f(n - 1) 种方法。 
根据乘法原理, n 个不同元素的错排种数 
f(n) = (n-1)[f(n-2)+f(n-1)] (n>2) 。 


#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<queue>
#include<set>
#include<map>
using namespace std;
typedef long long ll;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
	return f*x;
}
const int N=1001000;
const int mod=1000000007;
int fac[N],inv[N],er[N];
void initial()
{
	fac[0]=fac[1]=1;for(int i=2;i<=N;i++)fac[i]=int(ll(fac[i-1])*ll(i)%ll(mod));
	inv[1]=inv[0]=1;for(int i=2;i<=N;i++)inv[i]=int(ll(mod-mod/i)*inv[mod%i]%ll(mod));
	for(int i=2;i<=N;i++)inv[i]=int(ll(inv[i])*ll(inv[i-1])%ll(mod));
	er[0]=1;er[1]=0;er[2]=1;for(int i=3;i<=N;i++)er[i]=int(ll(i-1)*ll(er[i-1]+er[i-2])%ll(mod));
}
int main()
{
	int T=read(),x,y;
	initial();
	while(T--)
	{
		x=read();y=read();
		printf("%d\n",int(ll(er[x-y])%ll(mod)*ll(fac[x])%ll(mod)*ll(inv[x-y])%ll(mod)*ll(inv[y])%ll(mod)));
	}
	
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值