結果

問題 No.129 お年玉(2)
ユーザー rlangevin
提出日時 2023-02-11 01:55:49
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 404 ms / 5,000 ms
コード長 826 bytes
コンパイル時間 248 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 72,832 KB
最終ジャッジ日時 2024-07-07 18:58:50
合計ジャッジ時間 8,410 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 46
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

from math import *
from collections import *
def Sieve(n):
    lst = [True] * (n + 1)
    S = set()
    for i in range(2, ceil(sqrt(n)) + 1):
        if lst[i]:
            for j in range(2 * i, n + 1, i):
                lst[j] = False
    for i in range(2, n + 1):
        if lst[i]:
            S.add(i)
    return S
N = int(input())//1000
M = int(input())
r = N % M
ans = 1
mod = 10 ** 9
inv = 1
seen = [0] * (r + 1)
P = Sieve(M)
D = defaultdict(int)
for p in P:
    now = p
    while True:
        D[p] += r//now
        now *= p
        if M // now == 0:
            break
for i in range(r):
    now = M - i
    for k, v in D.items():
        while D[k]:
            if now % k == 0:
                now //= k
                D[k] -= 1
            else:
                break
    ans *= now
    ans %= mod
print(ans)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0