結果

問題 No.1978 Permutation Repetition
ユーザー shino16shino16
提出日時 2022-06-03 20:08:13
言語 PyPy3
(7.3.13)
結果
AC  
実行時間 249 ms / 2,000 ms
コード長 975 bytes
コンパイル時間 184 ms
コンパイル使用メモリ 81,812 KB
実行使用メモリ 68,504 KB
最終ジャッジ日時 2023-10-21 04:38:42
合計ジャッジ時間 6,510 ms
ジャッジサーバーID
(参考情報)
judge14 / judge15
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 38 ms
53,672 KB
testcase_01 AC 39 ms
53,672 KB
testcase_02 AC 43 ms
59,648 KB
testcase_03 AC 44 ms
59,644 KB
testcase_04 AC 58 ms
59,656 KB
testcase_05 AC 65 ms
59,652 KB
testcase_06 AC 198 ms
63,756 KB
testcase_07 AC 59 ms
59,652 KB
testcase_08 AC 233 ms
63,756 KB
testcase_09 AC 171 ms
63,768 KB
testcase_10 AC 175 ms
63,772 KB
testcase_11 AC 135 ms
63,756 KB
testcase_12 AC 44 ms
59,644 KB
testcase_13 AC 87 ms
61,704 KB
testcase_14 AC 219 ms
63,756 KB
testcase_15 AC 115 ms
63,768 KB
testcase_16 AC 50 ms
61,708 KB
testcase_17 AC 58 ms
59,652 KB
testcase_18 AC 221 ms
63,756 KB
testcase_19 AC 98 ms
61,708 KB
testcase_20 AC 38 ms
53,672 KB
testcase_21 AC 70 ms
59,652 KB
testcase_22 AC 150 ms
61,704 KB
testcase_23 AC 110 ms
63,772 KB
testcase_24 AC 95 ms
61,724 KB
testcase_25 AC 120 ms
61,704 KB
testcase_26 AC 90 ms
61,704 KB
testcase_27 AC 155 ms
61,704 KB
testcase_28 AC 174 ms
63,756 KB
testcase_29 AC 138 ms
61,704 KB
testcase_30 AC 183 ms
63,768 KB
testcase_31 AC 101 ms
61,704 KB
testcase_32 AC 129 ms
63,768 KB
testcase_33 AC 127 ms
63,772 KB
testcase_34 AC 63 ms
59,656 KB
testcase_35 AC 148 ms
61,704 KB
testcase_36 AC 80 ms
61,700 KB
testcase_37 AC 58 ms
61,704 KB
testcase_38 AC 80 ms
61,720 KB
testcase_39 AC 39 ms
53,672 KB
testcase_40 AC 58 ms
59,656 KB
testcase_41 AC 182 ms
63,756 KB
testcase_42 AC 249 ms
68,504 KB
testcase_43 AC 210 ms
68,504 KB
testcase_44 AC 40 ms
53,672 KB
testcase_45 AC 38 ms
53,672 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