結果

問題 No.2045 Two Reflections
ユーザー tktk_snsntktk_snsn
提出日時 2022-08-20 00:15:31
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 213 ms / 2,000 ms
コード長 1,913 bytes
コンパイル時間 267 ms
コンパイル使用メモリ 12,800 KB
実行使用メモリ 15,616 KB
最終ジャッジ日時 2024-10-08 12:03:13
合計ジャッジ時間 4,448 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 27
権限があれば一括ダウンロードができます

ソースコード

diff # 

from collections import defaultdict
mod = 998244353


def prime_factorization(n):  # 素因数分解
    res = []
    for i in range(2, int(n ** 0.5) + 1):
        if n % i == 0:
            cnt = 0
            while n % i == 0:
                cnt += 1
                n //= i
            res.append((i, cnt))
    if n > 1:
        res.append((n, 1))
    return res


class UnionFind:
    def __init__(self, n):
        self.root = [-1] * (n + 1)

    def find(self, x):
        stack = []
        while self.root[x] >= 0:
            stack.append(x)
            x = self.root[x]
        for i in stack:
            self.root[i] = x
        return x

    def same(self, x, y):
        return self.find(x) == self.find(y)

    def unite(self, x, y):
        x = self.find(x)
        y = self.find(y)
        if x == y:
            return False
        if -self.root[x] > -self.root[y]:
            x, y = y, x
        self.root[y] += self.root[x]
        self.root[x] = y
        return True

    def size(self, x):
        return -self.root[self.find(x)]


def main(N, p, q):

    if p == 1 and q == 1:
        return 1
    if p == 1 or q == 1:
        return 2
    if p + q <= N:
        return 4

    A = list(range(N))
    for i in range(p):
        j = p - 1 - i
        if i < j:
            A[i], A[j] = A[j], A[i]
    for i in range(q):
        j = q - 1 - i
        i = N - 1 - i
        j = N - 1 - j
        if i < j:
            A[i], A[j] = A[j], A[i]

    uf = UnionFind(N)
    for i, a in enumerate(A):
        uf.unite(i, a)

    prime = defaultdict(int)
    for i in range(N):
        if i == uf.find(i):
            for k, v in prime_factorization(uf.size(i)):
                if prime[k] < v:
                    prime[k] = v

    res = 1
    for k, v in prime.items():
        res = res * pow(k, v, mod) % mod
    return res * 2 % mod


N, p, q = map(int, input().split())
print(main(N, p, q))
0