結果

問題 No.1815 K色問題
ユーザー 259_Momone259_Momone
提出日時 2022-09-18 03:44:51
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 1,695 ms / 2,000 ms
コード長 1,339 bytes
コンパイル時間 361 ms
コンパイル使用メモリ 12,928 KB
実行使用メモリ 61,892 KB
最終ジャッジ日時 2024-09-29 22:47:02
合計ジャッジ時間 13,062 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 449 ms
44,276 KB
testcase_01 AC 26 ms
11,008 KB
testcase_02 AC 439 ms
44,420 KB
testcase_03 AC 26 ms
11,008 KB
testcase_04 AC 26 ms
11,136 KB
testcase_05 AC 448 ms
44,168 KB
testcase_06 AC 1,673 ms
26,068 KB
testcase_07 AC 396 ms
14,360 KB
testcase_08 AC 1,527 ms
58,072 KB
testcase_09 AC 348 ms
18,712 KB
testcase_10 AC 507 ms
22,584 KB
testcase_11 AC 988 ms
58,012 KB
testcase_12 AC 26 ms
11,136 KB
testcase_13 AC 26 ms
11,136 KB
testcase_14 AC 1,695 ms
61,892 KB
testcase_15 AC 1,687 ms
61,776 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

ans = [0] * (K + 1)

if N == 1:
    ans = [i * pow(i - 1, 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