FZU2282--错排公式

http://acm.fzu.edu.cn/problem.php?pid=2282

F[i]=(i-1)(F[i-2]+F[i-1])

#include<cstdio>
#include<iostream>
#define LL long long 
#define Mod 1000000007
#define maxn 10010
using namespace std;
int n,m,k;
LL C[maxn][120],F[maxn],ans,tot[maxn];
int main(){
    F[1]=0,F[2]=1,tot[1]=1,tot[2]=2,tot[0]=1;
    for(int i=3;i<=10000;i++){
        F[i]=(i-1)*((F[i-2]+F[i-1])%Mod);
        F[i]%=Mod;
        tot[i]=tot[i-1]*i;
        tot[i]%=Mod;
    //    cout<<i<<' '<<tot[i]<<endl;
    }
    C[0][0]=1;
    for(int i=1;i<=10000;i++){
        C[i][0]=1;
        for(int j=1;j<=min(i,101);j++){
            C[i][j]=(C[i-1][j]+C[i-1][j-1])%Mod;
        //    cout<<C[i][j]<<' ';
        }
        //cout<<endl;
    }
    int T;
    cin>>T;
    while(T--){
        cin>>n>>k;
        ans=0;
        for(int i=0;i<k;i++){
            ans+=C[n][i]*F[n-i];
            ans%=Mod;    
        }
        LL res=tot[n]-ans;
        while(res<0) res+=Mod;
        cout<<(res%Mod)<<endl;
    }
    return 0;
}