結果
問題 | No.1396 Giri |
ユーザー | buey_t |
提出日時 | 2023-04-16 15:54:10 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 536 ms / 2,000 ms |
コード長 | 3,260 bytes |
コンパイル時間 | 191 ms |
コンパイル使用メモリ | 82,228 KB |
実行使用メモリ | 125,952 KB |
最終ジャッジ日時 | 2024-04-20 05:12:21 |
合計ジャッジ時間 | 7,563 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 111 ms
84,972 KB |
testcase_01 | AC | 110 ms
85,088 KB |
testcase_02 | AC | 526 ms
125,920 KB |
testcase_03 | AC | 107 ms
84,756 KB |
testcase_04 | AC | 109 ms
84,912 KB |
testcase_05 | AC | 536 ms
125,932 KB |
testcase_06 | AC | 111 ms
84,940 KB |
testcase_07 | AC | 110 ms
85,244 KB |
testcase_08 | AC | 113 ms
84,932 KB |
testcase_09 | AC | 115 ms
84,896 KB |
testcase_10 | AC | 113 ms
85,124 KB |
testcase_11 | AC | 110 ms
85,064 KB |
testcase_12 | AC | 109 ms
84,988 KB |
testcase_13 | AC | 108 ms
84,928 KB |
testcase_14 | AC | 121 ms
89,316 KB |
testcase_15 | AC | 121 ms
89,364 KB |
testcase_16 | AC | 140 ms
90,012 KB |
testcase_17 | AC | 146 ms
90,744 KB |
testcase_18 | AC | 176 ms
92,760 KB |
testcase_19 | AC | 337 ms
108,320 KB |
testcase_20 | AC | 404 ms
116,048 KB |
testcase_21 | AC | 474 ms
122,244 KB |
testcase_22 | AC | 519 ms
125,464 KB |
testcase_23 | AC | 513 ms
125,920 KB |
testcase_24 | AC | 505 ms
125,952 KB |
testcase_25 | AC | 531 ms
125,732 KB |
ソースコード
from math import sqrt,sin,cos,tan,ceil,radians,floor,gcd,exp,log,log10,log2,factorial,fsum from heapq import heapify, heappop, heappush from bisect import bisect_left, bisect_right from copy import deepcopy import copy import random from collections import deque,Counter,defaultdict from itertools import permutations,combinations from decimal import Decimal,ROUND_HALF_UP #tmp = Decimal(mid).quantize(Decimal('0'), rounding=ROUND_HALF_UP) from functools import lru_cache, reduce #@lru_cache(maxsize=None) from operator import add,sub,mul,xor,and_,or_,itemgetter INF = 10**18 mod1 = 10**9+7 mod2 = 998244353 #DecimalならPython #再帰ならPython!!!!!!!!!!!!!!!!!!!!!!!!!! def sieve(n): is_prime = [True for _ in range(n+1)] is_prime[0] = False is_prime[1] = False for i in range(2, n+1): if is_prime[i]: j = 2 * i while j <= n: is_prime[j] = False j += i table = [ i for i in range(1, n+1) if is_prime[i]] return is_prime, table class Eratosthenes: def __init__(self, n): self.isprime = [True]*(n+1) self.minfactor = [-1]*(n+1) self.isprime[1] = False self.minfactor[1] = 1 for p in range(2, n + 1): if self.isprime[p]: self.minfactor[p] = p q = p while q <= n: self.isprime[q] = False q += p if q > n: break if self.minfactor[q] == -1: self.minfactor[q] = p def factorize(self, n): factors = [] while n > 1: p = self.minfactor[n] exp = 0 while self.minfactor[n] == p: n //= p exp += 1 factors.append((p, exp)) return factors def divisors(self, n): res = [1] pf = self.factorize(n) for p in pf: s = len(res) for i in range(s): v = 1 for j in range(p[1]): v *= p[0] res.append(res[i] * v) return res ''' エラトステネスっぽくやるのが思いつくけど ダメだな こんなん愚直にやればよくね gcdだけ取り除くべき 最大の素数はかけないほうが良い 逆元がない可能性があるじゃん 困った 逆元持つまで素因数分解するみたいなのもありか? いや、modの話だわ 素数じゃないといけないのは ansがmodとってるから、gcdが機能してない エラトステネスで、先に割ってあげる 後ろからかけていったらどう それの約数分解をmemoにぶちこむ N = 8のとき、6をかけるのではなく3をかけたほうがよい ''' N = int(input()) prime,t = sieve(N+1) er = Eratosthenes(N+1) pr = -1 memo = [0]*(N+1) ans = 1 for i in range(1,N+1): for p,e in er.factorize(i): if e-memo[p] <= 0: continue ans *= pow(p,e-memo[p],mod2) ans %= mod2 memo[p] += e-memo[p] if prime[i]: pr = i ans *= pow(pr,mod2-2,mod2) ans %= mod2 print(ans)