4517: [Sdoi2016]排列计数

Description

求有多少种长度为 n 的序列 A，满足以下条件：

1 ~ n 这 n 个数在序列中各出现了一次

若第 i 个数 A[i] 的值为 i，则称 i 是稳定的。序列恰好有 m 个数是稳定的

满足条件的序列可能很多，序列数对 10^9+7 取模。

题解：

序列里任意m个数稳定都是等价的，考虑不稳定的那几个就是裸错排

错排问题的递推式f[n]=(f[n-1]+f[n-2])*(n-1)。

我们考虑这n个人中编号最小的人i，设他站到了j号位置。

1.如果j号人站到了i号位置，问题转化成一个n-2的错排问题。

2.如果j号人没有站到i号位置，那么除去i号人所有人都站错了，问题转化成一个n-1的错排问题。

代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define mod 1000000007 
#define ll long long
const int N=1000010;
ll t;
ll n,m;
ll fac[N],ine[N],f[N];
int main()
{
	
	fac[0]=ine[0]=ine[1]=1;
	
	for(ll i=1;i<1000001;i++)fac[i]=(fac[i-1]*i)%mod;
	 for(int i=2;i<1000001;i++){  
        ine[i]=(mod-mod/i)*ine[mod%i]%mod;  
    }  
	for(ll i=2;i<1000001;i++)ine[i]=(ine[i]*ine[i-1])%mod;
	f[0]=1;
	f[1]=0;
	f[2]=1;
	for(ll i=3;i<=1000000;i++)
	f[i]=((i-1)*(f[i-1]+f[i-2]))%mod;
	scanf("%lld",&t);
	while(t--)
	{
		scanf("%lld%lld",&n,&m);
		printf("%lld\n",fac[n]*ine[m]%mod*ine[n-m]%mod*f[n-m]%mod);
	}
}