結果
| 問題 |
No.1428 PeRmutation Question
|
| コンテスト | |
| ユーザー |
convexineq
|
| 提出日時 | 2021-03-12 23:06:38 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,199 bytes |
| コンパイル時間 | 432 ms |
| コンパイル使用メモリ | 81,996 KB |
| 実行使用メモリ | 93,484 KB |
| 最終ジャッジ日時 | 2024-10-14 13:41:27 |
| 合計ジャッジ時間 | 3,791 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 2 WA * 28 |
ソースコード
SIZE=10**5+1; MOD=10**9+7 #998244353 #ここを変更する
inv = [0]*SIZE # inv[j] = j^{-1} mod MOD
fac = [0]*SIZE # fac[j] = j! mod MOD
finv = [0]*SIZE # finv[j] = (j!)^{-1} mod MOD
fac[0] = fac[1] = 1
finv[0] = finv[1] = 1
for i in range(2,SIZE):
fac[i] = fac[i-1]*i%MOD
finv[-1] = pow(fac[-1],MOD-2,MOD)
for i in range(SIZE-1,0,-1):
finv[i-1] = finv[i]*i%MOD
inv[i] = finv[i]*fac[i-1]%MOD
def choose(n,r): # nCk mod MOD の計算
if 0 <= r <= n:
return (fac[n]*finv[r]%MOD)*finv[n-r]%MOD
else:
return 0
def perm(n,r): # nPr mod MOD
if 0 <= r <= n:
return fac[n]*finv[n-r]%MOD
else:
return 0
n = int(input())
*p, = map(int,input().split())
for i in range(n):
p[i] -= 1
res = []
used = [0]*n
for i in range(n):
if used[i]: continue
c = 0
while used[i]==0:
used[i] = 1
i = p[i]
c += 1
res.append(c)
#print(res)
ans = fac[n]
for x in res:
ans *= finv[x]
ans %= MOD
from collections import Counter
d = Counter(res)
for v in d.values():
ans *= finv[v]
ans %= MOD
for x in res:
if x%2:
ans *= fac[x]
else:
ans *= fac[x//2]
ans %= MOD
print(ans)
convexineq