FZU -- 2282 和 DHU -- 2049 【简单错排 + 组合数学】

传送门
FZU — 2282
题意 ： 给你n个人, k 个数, 任意交换n个人的顺序，1号还在一号位就表示它的位置是对的. 则问你至少k个人还在对的位置的方案有几种.
//思路：思路这个应该就可以想到用错排公式加组合数学解决. 因为涉及到除法进行去mod运算, 所以要用到逆元. 首先预处理所有的数据,然后直接算就是了.

/** @Cain*/
#define ll long long int
const ll mod = 1e9+7;
const int maxn = 1e4+5;
int n,k;
ll qpow(ll x,ll y)
{
    ll res = 1;
    while(y){
        if(y&1) res = x * res % mod;
        x = x * x % mod;
        y >>= 1;
    }
    return res;
}

ll D[maxn],fac[maxn],inv[maxn];
void init()
{
    //因为这里我们是反着弄的,即0个放错时, 全部放对有一种方法的. 因为题目是求对的方法吗.
    //所以就应该初始化为1. 错的方法才是初始化为0.
    D[0] = 1; D[1] = 0; D[2] = 1;
    for(ll i=3;i<=maxn;i++){
        D[i] = (i-1)*(D[i-1]+D[i-2]) % mod;
    } //初始化D(n)错排数.
    fac[0] = 1;
    for(ll i=1;i<maxn;i++) fac[i] = (fac[i-1]*i) % mod;  // 初始化n!数组.
    inv[maxn-1] = qpow(fac[maxn-1],mod-2);//先求一最大数的逆元.
    for(ll i=maxn-2;i>=0;i--){
        inv[i] = (inv[i+1]*(i+1) % mod);//初始化每一个阶乘对应
   //的逆元. 这个看不太懂就把它写出来就行了. a*a的逆元==1,这样做可
   //以把每一个阶乘的逆元给预处理出来.

    }
}

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

void solve()
{
    scanf("%d%d",&n,&k);
    ll res = 0;
    for(int i=k;i<=n;i++){
        res += D[n-i]*C(n,n-i);    //C(n,i) 恒等 C(n,n-i),组合数学基本公式吗.
        res %= mod;
    }
    printf("%lld\n",res%mod);
}

int main()
{
    int t;
    scanf("%d",&t);
    init();
    while(t--){
        solve();
    }
}

HDU 2049
这个是中文题面就不说了.
思路 ：比较裸的错排题，直接上错排公式和组合数学公式就可以了．
AC Code

#define ll long long int
ll f[25],c[25];
void init()
{
    f[1] = 0; f[2] = 1;
    c[0] = 1;c[1] = 1; c[2] = 2;
    for(int i=3;i<=21;i++){
        f[i] = (i-1)*(f[i-1]+f[i-2]);
        c[i] = i*c[i-1];
    }
}
void solve()
{
    int n,m;
    scanf("%d%d",&n,&m);
    printf("%lld\n",f[m]*(c[n]/c[m]/c[n-m]));
}
int main()
{
    int t;
    scanf("%d",&t);
    init();
    while(t--){
        solve();
    }
}