結果

問題 No.3441 Sort Permutation 2
コンテスト
ユーザー 👑 loop0919
提出日時 2026-02-06 21:54:11
言語 PyPy3
(7.3.17)
結果
AC  
実行時間 850 ms / 2,000 ms
コード長 2,349 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 274 ms
コンパイル使用メモリ 82,292 KB
実行使用メモリ 128,844 KB
最終ジャッジ日時 2026-02-06 21:54:30
合計ジャッジ時間 17,715 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 41
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

from itertools import pairwise
from math import gcd, isqrt


class dsu:
    n = 1
    parent_or_size = [-1 for i in range(n)]

    def __init__(self, N):
        self.n = N
        self.parent_or_size = [-1 for i in range(N)]

    def merge(self, a, b):
        assert 0 <= a < self.n, "0<=a<n,a={0},n={1}".format(a, self.n)
        assert 0 <= b < self.n, "0<=b<n,b={0},n={1}".format(b, self.n)
        x = self.leader(a)
        y = self.leader(b)
        if x == y:
            return x
        if -self.parent_or_size[x] < -self.parent_or_size[y]:
            x, y = y, x
        self.parent_or_size[x] += self.parent_or_size[y]
        self.parent_or_size[y] = x
        return x

    def same(self, a, b):
        assert 0 <= a < self.n, "0<=a<n,a={0},n={1}".format(a, self.n)
        assert 0 <= b < self.n, "0<=b<n,b={0},n={1}".format(b, self.n)
        return self.leader(a) == self.leader(b)

    def leader(self, a):
        assert 0 <= a < self.n, "0<=a<n,a={0},n={1}".format(a, self.n)
        if self.parent_or_size[a] < 0:
            return a
        self.parent_or_size[a] = self.leader(self.parent_or_size[a])
        return self.parent_or_size[a]

    def size(self, a):
        assert 0 <= a < self.n, "0<=a<n,a={0},n={1}".format(a, self.n)
        return -self.parent_or_size[self.leader(a)]

    def groups(self):
        leader_buf = [0 for i in range(self.n)]
        group_size = [0 for i in range(self.n)]
        for i in range(self.n):
            leader_buf[i] = self.leader(i)
            group_size[leader_buf[i]] += 1
        result = [[] for i in range(self.n)]
        for i in range(self.n):
            result[leader_buf[i]].append(i)
        result2 = []
        for i in range(self.n):
            if len(result[i]) > 0:
                result2.append(result[i])
        return result2


def divisors(n: int) -> set[int]:
    div = set()
    for i in range(1, isqrt(n) + 1):
        if n % i == 0:
            div.update([i, n // i])

    return div


N = int(input())
P = [int(s) - 1 for s in input().split()]

uf = dsu(N)

for i, p in enumerate(P):
    uf.merge(i, p)

answers = [0] * N

for group in uf.groups():
    if len(group) <= 1:
        continue

    g = gcd(*[abs(x - y) for x, y in pairwise(group)])
    for div in divisors(g):
        answers[div] += len(group) - 1

print(*answers[1:], sep="\n")
0