def polymul(f,g): lf = len(f) lg = len(g) res = [0]*(lf+lg-1) for i in range(lf): for j in range(lg): res[i+j] += f[i]*g[j] res[i+j] %= MOD return res def fps_nth_term(f,g,N): assert g[0] != 0 while N: h = g[:] for i in range(1,len(g),2): h[i] = -h[i] f = polymul(f,h)[N%2:N+1:2] g = polymul(g,h)[:N+1:2] N //= 2 return f[0]*pow(g[0],MOD-2,MOD)%MOD MOD = 10**9+7 n,p,c = map(int,input().split()) MAX = 13*p+1 f = [[0]*MAX for _ in range(p+1)] f[0][0] = 1 for x in [2,3,5,7,11,13]: for i in range(p): for j in range(MAX-x): f[i+1][j+x] += f[i][j] MAX = 12*c+1 g = [[0]*MAX for _ in range(c+1)] g[0][0] = 1 for x in [4,6,8,9,10,12]: for i in range(c): for j in range(MAX-x): g[i+1][j+x] += g[i][j] h = polymul(f[-1],g[-1]) h1 = [-hi for hi in h] h2 = h1[:] h1[0] += sum(h) h2[0] += 1 h2 = polymul(h2,[1,-1]) ans = fps_nth_term(h1,h2,n-1) print(ans)