結果
| 問題 |
No.2755 行列の共役類
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2023-08-07 21:12:43 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 3,099 ms / 3,500 ms |
| コード長 | 1,043 bytes |
| コンパイル時間 | 266 ms |
| コンパイル使用メモリ | 82,140 KB |
| 実行使用メモリ | 79,556 KB |
| 最終ジャッジ日時 | 2024-11-08 02:37:45 |
| 合計ジャッジ時間 | 23,394 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| other | AC * 66 |
ソースコード
import sys, time, random
from collections import deque, Counter, defaultdict
input = lambda: sys.stdin.readline().rstrip()
ii = lambda: int(input())
mi = lambda: map(int, input().split())
li = lambda: list(mi())
inf = 2 ** 63 - 1
mod = 998244353
B, C = map(int,input().split())
P = []
b = B
e = 1
for i in range(2,B):
if i*i>b:break
while b % i == 0:
b //= i
P.append(i)
e *= (i - 1)
while b % i == 0:
b //= i
e *= i
if b>1:
P.append(b)
e *= (b - 1)
n = len(P)
u = [pow(i,e-1,B)for i in range(B)]
d = B // C
b = d * B
found = [0] * b
answer = 0
for i in range(b):
if found[i]:continue
q = i // d
r = i % d
c = 1
for p in P:
c *= q % p
if c == 0:
continue
answer +=1
if answer >100:
print("100+")
exit()
for j in range(b):
s = j // d
t = j % d
c = 1
for p in P:
c *=s % p
found[(u[s] * q * s % B) * d + (u[s] * (q * t + r - t)) % d] |= c
print(answer)