1. 程式人生 > >lightoj 1095 - Arrange the Numbers(dp+組合數)

lightoj 1095 - Arrange the Numbers(dp+組合數)

鏈接 簡單 turn amp using ace ring ++ ber

題目鏈接:http://www.lightoj.com/volume_showproblem.php?problem=1095

題解:其實是一道簡單的組合數只要推導一下錯排就行了。在這裏就推導一下錯排

dp[i]=(i-1)*dp[i-2](表示新加的那個數放到i-1中的某一個位置然後那個被放位置的數放在i這個位置就是i-2的錯排)+(i-1)*dp[i-1](表示新加的那個數放到i-1中的某一個位置然後用那個位置被占的數代替i這個位置的數就是i-1的錯排)

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
#define mod 1000000007
using namespace std;
typedef long long ll;
const int M = 1e3 + 10;
ll dp[M];
ll up[M] , down[M];
ll inv(ll a) {
    return a == 1 ? 1 : (ll)(mod - mod / a) * inv(mod % a) % mod;
}
void fk() {
    dp[0] = 1 , dp[1] = 0 , dp[2] = 1;
    for(int i = 3 ; i < M ; i++) dp[i] = (i - 1) * ((dp[i - 1] + dp[i - 2]) % mod) , dp[i] %= mod;
}
ll C(ll n , ll m)
{
    if(m < 0)return 0;
    if(n < m)return 0;
    if(m > n-m) m = n-m;
    ll up = 1, down = 1;
    for(ll i = 0 ; i < m ; i++){
        up = up * (n-i) % mod;
        down = down * (i+1) % mod;
    }
    return up * inv(down) % mod;
}
int main() {
    fk();
    int t , Case = 0;
    scanf("%d" , &t);
    while(t--) {
        int n , m , k;
        scanf("%d%d%d" , &n , &m , &k);
        ll ans = 0;
        ll gg = C(m , k);
        up[0] = 1 , down[0] = 1;
        for(int i = 1 ; i <= (n - m) / 2 ; i++) up[i] = up[i - 1] * ((n - m) - i + 1) % mod , down[i] = down[i - 1] * i % mod;
        for(int i = (n - m) / 2 + 1 ; i <= (n - m) ; i++) up[i] = up[(n - m) - i] , down[i] = down[(n - m) - i];
        for(int i = n - k ; i >= (m - k) ; i--) {
            ans += dp[i] * (up[n - k - i] * (inv(down[n - k - i]) % mod) % mod);
            ans %= mod;
        }
        ans *= gg;
        ans %= mod;
        printf("Case %d: %lld\n" , ++Case , (ans + mod) % mod);
    }
    return 0;
}

lightoj 1095 - Arrange the Numbers(dp+組合數)