lightoj 1095 - Arrange the Numbers(dp+組合數)
阿新 • • 發佈:2017-07-08
鏈接 簡單 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+組合數)