結果

問題 No.3187 Mingle
コンテスト
ユーザー LyricalMaestro
提出日時 2026-05-05 18:51:00
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
AC  
実行時間 1,921 ms / 2,500 ms
コード長 2,040 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 242 ms
コンパイル使用メモリ 84,884 KB
実行使用メモリ 176,832 KB
最終ジャッジ日時 2026-05-05 18:51:35
合計ジャッジ時間 30,354 ms
ジャッジサーバーID
（参考情報） 		judge2_0 / judge1_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

## https://yukicoder.me/problems/no/866


class BinaryIndexTree:
    """
    フェニック木(BinaryIndexTree)の基本的な機能を実装したクラス
    """
    def __init__(self, size, mod):
        self.size = size
        self.array = [0] * (size + 1)
        self.mod = mod
    
    def add(self, x, a):
        index = x
        while index <= self.size:
            self.array[index] += a
            self.array[index] %= self.mod
            index += index & (-index)
    
    def sum(self, x):
        index = x
        ans = 0
        while index > 0:
            ans += self.array[index]
            ans %= self.mod
            index -= index & (-index)
        return ans

    def least_upper_bound(self, value):
        if self.sum(self.size) < value:
            return -1
        elif value <= 0:
            return 0

        m = 1
        while m < self.size:
            m *= 2

        k = 0
        k_sum = 0
        while m > 0:
            k0 = k + m
            if k0 < self.size:
                if k_sum + self.array[k0] < value:
                    k_sum += self.array[k0]
                    k += m
            m //= 2
        if k < self.size:
            return k + 1
        else:
            return -1


def main():
    N, P = map(int, input().split())

    # イベント派生
    events = [[] for _ in range(N + 1)]
    for p in range(1, N + 1):
        x = p
        while x <= N:
            events[x].append(p)
            x += p
            
    bit = BinaryIndexTree(N, P)
    dp = [0] * (N + 1)
    prevs = [-1] * (N + 1)
    for p in range(3, N + 1):
        cnt = len(events[p])

        for x in events[p]:
            if prevs[x] == -1:
                continue
            bit.add(prevs[x], -dp[prevs[x]])

        inv_nx = pow(p - cnt, P - 2, P)
        dp[p] = (inv_nx * bit.sum(p - 1)) % P
        dp[p] += (p * inv_nx) % P
        dp[p] %= P

        for x in events[p]:
            prevs[x] = p
            bit.add(p, dp[p])

    print(dp[N])


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