#include #define GET_MACRO(a, b, c, NAME, ...) NAME #define rep(...) GET_MACRO(__VA_ARGS__, rep3, rep2)(__VA_ARGS__) #define rep2(i, a) rep3 (i, 0, a) #define rep3(i, a, b) for (int i = (a); i < (b); i++) #define repr(...) GET_MACRO(__VA_ARGS__, repr3, repr2)(__VA_ARGS__) #define repr2(i, a) repr3 (i, 0, a) #define repr3(i, a, b) for (int i = (b) - 1; i >= (a); i--) #define chmin(a, b) ((b) < a && (a = (b), true)) #define chmax(a, b) (a < (b) && (a = (b), true)) using namespace std; typedef long long ll; const ll mod = 1e9 + 7; struct Kitamasa { vector mul(vector &a, vector &b, vector &c) { int k = a.size(); vector res(k * 2); rep (i, k) rep (j, k) { (res[i + j] += a[i] * b[j]) %= mod; } repr (i, k, k * 2) rep (j, k) { (res[i - j - 1] += c[k - j - 1] * res[i]) %= mod; } res.resize(k); return res; } vector pow(vector &a, ll b) { vector p(a.size()), res(a.size()); p[1] = res[0] = 1; while (b) { if (b & 1) res = mul(res, p, a); p = mul(p, p, a); b /= 2; } return res; } // b[n+k]=a[0]b[n]+...+a[k-1]b[n+k-1] ll calc(vector a, vector b, ll n) { a.insert(a.begin(), 0); b.insert(b.begin(), 0); a = pow(a, n + 1); ll res = 0; rep (i, a.size()) { (res += a[i] * b[i]) %= mod; } return res; } }; ll dp[2][55][2020]; ll way[2020]; void calc(ll dp[55][2020], ll n, vector dice) { dp[0][0] = 1; rep (k, dice.size()) { rep (i, n) { rep (j, 1010) { (dp[i + 1][j + dice[k]] += dp[i][j]) %= mod; } } } } template ostream &operator <<(ostream &os, const vector &v) { rep (i, v.size()) { if (i) os << " "; os << v[i]; } return os; } int main() { ll N, P, C; cin >> N >> P >> C; calc(dp[0], P, {2, 3, 5, 7, 11, 13}); calc(dp[1], C, {4, 6, 8, 9, 10, 12}); rep (i, 1010) { rep (j, 1010) { (way[i + j] += dp[0][P][i] * dp[1][C][j]) %= mod; } } ll K = P * 13 + C * 12 + 1; vector a(K), b(K); rep (i, K) { a[i] = way[K - i]; b[i] = 1; } Kitamasa ktms; cout << ktms.calc(a, b, N + K - 1) << endl; return 0; }