BZOJ4517 & 洛谷4071：[SDOI2016]排列计数——题解

https://www.lydsy.com/JudgeOnline/problem.php?id=4517

https://www.luogu.org/problemnew/show/P4071

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

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

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

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

sb题，有C(n,m)种可能稳定，剩下的就是使n-m错排即可。

错排公式d[i]=(i-1)*(d[i-1]+d[i-2])

（虽然这个公式是我现查的……）

#include<cmath>
#include<queue>
#include<vector>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int p=1e9+;
const int N=1e6+;
inline int read(){
    int X=,w=;char ch=;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
int qpow(int k,int n){
    int res=;
    while(n){
    if(n&)res=(ll)res*k%p;
    k=(ll)k*k%p;n>>=;
    }
    return res;
}
int jc[N],inv[N],d[N];
inline int C(int n,int m){
    return (ll)jc[n]*inv[m]%p*inv[n-m]%p;
}
void init(int n){
    jc[]=;
    for(int i=;i<=n;i++)jc[i]=(ll)jc[i-]*i%p;
    inv[n]=qpow(jc[n],p-);
    for(int i=n-;i;i--)inv[i]=(ll)inv[i+]*(i+)%p;
    inv[]=;
    d[]=;d[]=;d[]=;
    for(int i=;i<=n;i++)d[i]=(ll)(i-)*(d[i-]+d[i-])%p;
}
int main(){
    init(N-);
    int T=read();
    while(T--){
    int n=read(),m=read();
    printf("%lld\n",(ll)C(n,m)*d[n-m]%p);
    }
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++