測試地址：排列計數
做法：本題需要用到組合數學中的錯排問題。
首先，如果兩個序列中穩定的位置不同，那麼兩個序列肯定不同，因此我們列舉穩定的位置，有Cmn$C_n^m$種方案，然後對於剩下的元素，它們不能處在編號和它相同的位置上，學過組合數學的同學應該知道這就是經典的錯排問題，令長為n$n$的排列中，對所有1≤i≤n$1\le i\le n$都有ai≠i$a_i\ne i$的排列數為D(n)$D(n)$個，我們有遞推式：
D(n)=(n−1)(D(n−1)+D(n−2))$D(n)=(n-1)(D(n-1)+D(n-2))$
邊界條件為D(0)=1,D(1)=0$D(0)=1,D(1)=0$。那麼本題的答案就是Cmn⋅D(n−m)$C_n^m\cdot D(n-$

#include <bits/stdc++.h>
 

using namespace std;
typedef long long ll;
const ll mod=1000000007;
int T,n[500010],m[500010],maxn=0;
ll fac[1000010],inv[1000010],fi[1000010],D[1000010];

ll C(int n,int m)
{
    return fac[n]*fi[m]%mod*fi[n-m]%mod;
}

int main()
{
    scanf("%d",&T);
    for(int i=1;i<=T;i++)
    {
        scanf("%d%d",&n[i],&m[i]);
        maxn=max(maxn,n[i]);
    }

    fac[0
 
]=fac[1]=inv[0]=inv[1]=fi[0]=fi[1]=1;
    D[0]=1,D[1]=0;
    for(ll i=2;i<=(ll)maxn;i++)
    {
        fac[i]=fac[i-1]*i%mod;
        inv[i]=(mod-mod/i)*inv[mod%i]%mod;
        fi[i]=fi[i-1]*inv[i]%mod;
        D[i]=(i-1)*(D[i-1]+D[i-2])%mod;
    }

    for(int i=1;i<=T;i++)
        printf("%lld\n",C(n[i],m[i])*D[n[i]-m[i]]%mod);

    return 0;
}