結果

問題 No.1428 PeRmutation Question
ユーザー lam6er
提出日時 2025-04-15 23:01:24
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,294 bytes
コンパイル時間 290 ms
コンパイル使用メモリ 81,604 KB
実行使用メモリ 90,272 KB
最終ジャッジ日時 2025-04-15 23:02:33
合計ジャッジ時間 3,372 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 18 WA * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

MOD = 10**9 + 7

def main():
    import sys
    input = sys.stdin.read
    data = input().split()
    n = int(data[0])
    P = list(map(int, data[1:n+1]))
    
    # Convert to 0-based index
    P = [x - 1 for x in P]
    
    visited = [False] * n
    from collections import defaultdict
    counter = defaultdict(int)
    
    for i in range(n):
        if not visited[i]:
            cycle_length = 0
            j = i
            while not visited[j]:
                visited[j] = True
                j = P[j]
                cycle_length += 1
            counter[cycle_length] += 1
    
    # Precompute factorial and inverse factorial
    max_n = n
    fact = [1] * (max_n + 1)
    for i in range(1, max_n + 1):
        fact[i] = fact[i-1] * i % MOD
    inv_fact = [1] * (max_n + 1)
    inv_fact[max_n] = pow(fact[max_n], MOD-2, MOD)
    for i in range(max_n-1, -1, -1):
        inv_fact[i] = inv_fact[i+1] * (i+1) % MOD
    
    denominator_inv = 1
    for c, m in counter.items():
        term_c = pow(c, m, MOD)
        term_m_fact = fact[m]
        term = term_c * term_m_fact % MOD
        term_inv = pow(term, MOD-2, MOD)
        denominator_inv = denominator_inv * term_inv % MOD
    
    ans = fact[n] * denominator_inv % MOD
    print(ans)

if __name__ == '__main__':
    main()
0