BZOJ4517——组合数学+错排+乘法逆元

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

---

这题是一道错排题，对于错排，我们有公式：

f[i]=(i-1)*(f[i-1]+f[i-2]);

可以证明，当第i个数排在k时，若k在i，则剩下i-2个数有f[i-2]种错排方法；若k不在i，则这i-1个数有f[i-1]种可能性。由于k可能是1~i-1,所以可以推出上式。

对于本题，要求一个n序列，m个数固定的错排方案。则m个数固定，总共可能性有C（m,n）种，用C（m,n）乘上剩下的数（n-m）的错排方案f[n-m]数，即为答案。

然而数据范围十分大，必须使用mod,然而mod对于+，-，*（加减乘）都试用，但对于/（除法）就不适用了。

这里就需要用到乘法逆元，如果原本的结果中有除法，比如除以a,那就可以乘以a的逆元替代。那么该怎么求乘法逆元呢?乘法逆元的求法有很多，这里介绍快速幂求法。

在介绍快速幂求法之前，必须先了解费马小定理。

费马小定理(Fermat’s little theorem)是数论中的一个重要定理，在1636年提出，其内容为： 假如p是质数，且gcd(a,p)=1，那么 a(p-1)≡1（mod p），例如：假如a是整数，p是质数，则a,p显然互质(即两者只有一个公约数1)，那么我们可以得到费马小定理的一个特例，即当p为质数时候， a^(p-1)≡1(mod p)。

简单的说，若p为质数，a为整数，则a^(p-1)≡1(mod p)。若将其变形，则为:

a*a^(p-2)≡1(mod p),则可以看出x=a^(p-2)。直接快速幂得出。

Code：

#include<bits/stdc++.h>
#define ll long long
#define mod 1000000007
using namespace std;
int read(){
    char c;int x;while(c=getchar(),c<'0'||c>'9');x=c-'0';
    while(c=getchar(),c>='0'&&c<='9') x=x*10+c-'0';return x;
} 
ll f[1000005],n,T,x,y,ans,p[1000005];
ll k(int x){
    return p[x];
}
ll qpow(ll x,int y){
    ll ans=1;
    while(y)
    {
        if(y&1) ans=(ans*x)%mod;
        x=(x*x)%mod,y>>=1;
    }
    return ans;
}
ll C(int x,int y){
    return k(y)*qpow(k(x),mod-2)%mod*qpow(k(y-x),mod-2)%mod;
}
int main()
{
    T=read();
    f[0]=1;f[1]=0;f[2]=1;p[1]=1;p[0]=1;
    for(int i=2;i<=1000000;i++) p[i]=p[i-1]*i%mod;
    for(int i=3;i<=1000000;i++) 
      f[i]=(i-1)*(f[i-2]+f[i-1])%mod;
    while(T--){
        x=read(),y=read();
        ans=C(y,x)*f[x-y]%mod;
        printf("%lld\n",ans);
    }
    return 0;
}