結果

問題 No.1978 Permutation Repetition
ユーザー shino16shino16
提出日時 2022-06-03 20:08:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 234 ms / 2,000 ms
コード長 975 bytes
コンパイル時間 212 ms
コンパイル使用メモリ 82,184 KB
実行使用メモリ 67,456 KB
最終ジャッジ日時 2024-09-21 05:44:42
合計ジャッジ時間 6,165 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
51,968 KB
testcase_01 AC 39 ms
52,352 KB
testcase_02 AC 43 ms
58,240 KB
testcase_03 AC 42 ms
58,112 KB
testcase_04 AC 56 ms
59,008 KB
testcase_05 AC 63 ms
59,392 KB
testcase_06 AC 191 ms
63,104 KB
testcase_07 AC 57 ms
59,008 KB
testcase_08 AC 224 ms
62,336 KB
testcase_09 AC 162 ms
62,848 KB
testcase_10 AC 167 ms
62,848 KB
testcase_11 AC 130 ms
61,824 KB
testcase_12 AC 44 ms
58,368 KB
testcase_13 AC 86 ms
60,032 KB
testcase_14 AC 205 ms
62,592 KB
testcase_15 AC 111 ms
62,336 KB
testcase_16 AC 48 ms
59,904 KB
testcase_17 AC 56 ms
59,648 KB
testcase_18 AC 218 ms
62,720 KB
testcase_19 AC 96 ms
61,696 KB
testcase_20 AC 38 ms
51,840 KB
testcase_21 AC 69 ms
60,160 KB
testcase_22 AC 141 ms
61,440 KB
testcase_23 AC 104 ms
62,464 KB
testcase_24 AC 90 ms
61,824 KB
testcase_25 AC 120 ms
60,800 KB
testcase_26 AC 89 ms
60,288 KB
testcase_27 AC 150 ms
61,568 KB
testcase_28 AC 169 ms
62,336 KB
testcase_29 AC 134 ms
61,696 KB
testcase_30 AC 179 ms
62,592 KB
testcase_31 AC 101 ms
61,312 KB
testcase_32 AC 128 ms
62,208 KB
testcase_33 AC 122 ms
62,464 KB
testcase_34 AC 62 ms
60,288 KB
testcase_35 AC 139 ms
60,928 KB
testcase_36 AC 79 ms
59,392 KB
testcase_37 AC 56 ms
59,392 KB
testcase_38 AC 78 ms
60,800 KB
testcase_39 AC 40 ms
52,608 KB
testcase_40 AC 59 ms
59,648 KB
testcase_41 AC 175 ms
62,336 KB
testcase_42 AC 234 ms
67,456 KB
testcase_43 AC 205 ms
67,072 KB
testcase_44 AC 39 ms
51,968 KB
testcase_45 AC 40 ms
52,480 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

from math import gcd

MOD = 1000000007

N, M = map(int, input().split())
A = [int(s) - 1 for s in input().split()]

used = [False] * N
cycles = [0] * (N+1)
for i in range(N):
    if used[i]: continue
    cnt = 0
    while not used[i]:
        used[i] = True
        i = A[i]
        cnt += 1
    cycles[cnt] += 1

fact = [0] * (N+1)
ifact = [0] * (N+1)
fact[0] = 1
ifact[0] = 1
for i in range(N):
    fact[i+1] = fact[i] * (i+1) % MOD
    ifact[i+1] = pow(fact[i+1], MOD - 2, MOD)

def binom(n, k):
    return fact[n] * ifact[k] % MOD * ifact[n-k] % MOD

ans = 1
for i in range(1, N+1):
    dp = [0] * (cycles[i] + 1)
    ms = [m for m in range(1, N+1) if m // gcd(m, M) == i]
    dp[0] = 1
    for j in range(1, cycles[i] + 1):
        for m in ms:
            if j - m//i >= 0:
                dp[j] += dp[j - m//i] * binom(j-1, m//i - 1) % MOD * fact[m//i - 1] % MOD * pow(i, m//i - 1, MOD) % MOD
                dp[j] %= MOD
    ans = ans * dp[cycles[i]] % MOD
print(ans)
0