結果

問題 No.1815 K色問題
ユーザー 259_Momone
提出日時 2022-09-18 03:44:03
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
WA  
実行時間 -
コード長 1,335 bytes
コンパイル時間 458 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 63,900 KB
最終ジャッジ日時 2024-09-29 22:46:42
合計ジャッジ時間 13,752 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 11 WA * 2
権限があれば一括ダウンロードができます

ソースコード

diff # 

N, M, K = map(int, input().split())

ans = [0] * (K + 1)

if N == 1:
    ans = [i * pow(i, M - 1, 1000000007) % 1000000007 for i in range(K + 1)]
elif N == 2:
    ans = [i * (i - 1) * pow((i - 3) * i + 3, M - 1, 1000000007) % 1000000007 for i in range(K + 1)]
else:
    import numpy as np
    ans = np.array([(i - 2) * ((i - 3) * i + 5) % 1000000007 for i in range(K + 1)])
    b = np.array([(((i - 6) * i + 13) * i - 11) % 1000000007 * (i - 1) % 1000000007 for i in range(K + 1)])
    c = np.array([i * (i - 1) * (i - 1) % 1000000007 for i in range(K + 1)])
    d = np.array([((i - 3) * i + 1) % 1000000007 * i * (i - 1) % 1000000007 for i in range(K + 1)])
    n = M - 1
    while n:
        if n & 1:
            c = ans * c - d
            c = c % 1000000007
            d = b * d
        else:
            d = ans * d - b * c
        d = d % 1000000007
        ans = ans * ans - 2 * b
        ans = ans % 1000000007
        b = b * b
        b = b % 1000000007
        n >>= 1
    ans = c.tolist()

coef = [1] * (K + 1)

for i in range(K):
    coef[i + 1] = coef[i] * (K - i) % 1000000007

c = 1
for i in range(K + 1, 0, -1):
    c = -c * i % 1000000007
    coef[i - 1] = coef[i - 1] * c % 1000000007

ret = 0
for i in range(K + 1):
    ret += ans[i] * coef[i]

ret *= pow(coef[-1], 1000000005, 1000000007)
print(ret % 1000000007)
0