結果

問題 No.584 赤、緑、青の色塗り
コンテスト
ユーザー kidodesu
提出日時 2026-04-16 19:13:50
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 1,090 ms / 2,000 ms
コード長 1,110 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 171 ms
コンパイル使用メモリ 84,864 KB
実行使用メモリ 96,256 KB
最終ジャッジ日時 2026-04-16 19:14:00
合計ジャッジ時間 3,525 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

mod = 10**9+7
comb_num = 10**6 + 2
F = [1] * comb_num
F_ = [1] * comb_num
for i in range(2, comb_num): F[i] = F[i-1]*i%mod
F_[-1] = pow(F[-1], -1, mod)
for i in range(comb_num-1, 0, -1): F_[i-1] = F_[i]*i%mod
def comb(a, b):
    if 0 <= b <= a: return F[a] * F_[b] * F_[a-b] % mod
    else: return 0

n, r, g, b = list(map(int, input().split()))
m = r + g + b
r, b = min(r, g, b), max(r, g, b)
g = m - r - b
ans = 0
if n * 2 + 3 < m * 3:
    print(0)
else:
    for a2 in range(m//2+1):
        a1 = m - a2 * 2
        ans0 = comb(n-m+1, a1+a2) * comb(a1+a2, a1) * pow(2, a2, mod) % mod
        if not ans0: continue
        for b01 in range(min(a2, r)+1):
            for b20 in range(min(r-b01, a2-b01)+1):
                b12 = a2-b01-b20
                b0 = r-b01-b20
                b1 = g-b01-b12
                b2 = b-b12-b20
                if min(b12, b1, b2) < 0: continue
                ans1 = F[a2] * F_[b01] % mod * F_[b12] % mod * F_[b20] % mod
                ans1 = ans1 * F[a1] % mod * F_[b0] % mod * F_[b1] % mod * F_[b2] % mod
                ans = (ans + ans0 * ans1) % mod
    print(ans)
0