lightoj 1226 - One Unit Machine（dp+大组合数去摸）

题目链接：http://www.lightoj.com/volume_showproblem.php?problem=1226

题解：由于这些任务完成是有先后的所以最后一个完成的肯定是最后一个任务的子任务，不妨设dp[i]表示第几个任务完成后总共有几种方案，这里要逆着来至于为什么想想也是挺好理解的。于是有这么一个方程式dp[i]=dp[i + 1]*C(sum-1,k[i]-1)，这样列出来就更好理解了。就是最后一个位置肯定是确定的之后就靠组合来凑。

#include <iostream>
#include <cstring>
#include <cstdio>
#define mod 1000000007
using namespace std;
typedef long long ll;
const int M = ;
ll k[M];
ll dp[M] , Po[];
ll mod_pow(ll a , ll b) {
    ll res = ;
    while(b) {
        if(b & ) res = (res * a) % mod;
        a = (a * a) % mod;
        b >>= ;
    }
    return res;
}
ll C(ll n , ll m) {
    if(n < m) return ;
    if(!m || !n) return ;
    ll res = (Po[n] * mod_pow(Po[m] * Po[n - m] % mod , mod - )) % mod;
    return res;
}
int main() {
    int t , n;
    scanf("%d" , &t);
    int Case = ;
    Po[] = ;
    for(int i =  ; i <  ; i++) {
        Po[i] = Po[i - ] * i % mod;
        Po[i] %= mod;
    }
    while(t--) {
        memset(dp ,  , sizeof(dp));
        scanf("%d" , &n);
        ll sum = ;
        for(int i =  ; i <= n ; i++) {
            scanf("%lld" , &k[i]);
            sum += k[i];
        }
        dp[n] = C(sum -  , k[n] - );
        dp[n] %= mod;
        sum -= k[n];
        for(int i = n -  ; i >=  ; i--) {
            dp[i] = dp[i + ] * C(sum -  , k[i] - ) % mod;
            dp[i] %= mod;
            sum -= k[i];
        }
        printf("Case %d: %lld\n" , ++Case , dp[]);
    }
    return ;
}