https://code.google.com/codejam/contest/4214486/dashboard#s=p0

GCJ2014B_1.cpp

/*
f[i][n] = f[i-1][n-1]*C(n,1) + f[i-1][n02]*C(n,2) + ... + f[i-1][i-1] * C(n, n-(i-1));
*/

#include <iostream>
#include <vector>
using namespace std;

int base = 1000000007;
typedef long long llong;

llong combination[101][101];

void buildCombination() {
    for (int i = 0; i <= 100; i++) {
        for (int j = 0; j <= i; j++) {
            if (j == 0) {
                combination[i][j] = 1;
            } else {
                combination[i][j] = (combination[i-1][j] + combination[i-1][j-1]) % base;
            }
        }
    }
}

llong solve(int m, int n) {
    vector<vector<llong> > dp;
    dp.resize(m+1);
    for (int i = 0; i < m+1; i++) {
        dp[i].resize(n+1);
    }
    // i chars, len of j
    for (int i = 1; i <= m; i++) {
        for (int j = i; j <= n; j++) {
            if (i == 1) {
                dp[i][j] = 1;
                continue;
            }
            dp[i][j] = 0;
            for (int k = 1; j-k >= i-1; k++) {
                dp[i][j] = (dp[i][j] + dp[i-1][j-k] * combination[j][k]) % base;

            }
        }
    }
    return dp[m][n];
}

int main() {    
    int T;
    buildCombination();
    cin >> T;
    for (int i = 1; i <= T; i++) {
        int m, n;
        cin >> m >> n;
        llong r = solve(m, n);
        cout << "Case #" << i << ": " << r << endl;       
    }
    return 0;
}