結果
問題 |
No.1322 Totient Bound
|
ユーザー |
![]() |
提出日時 | 2025-03-20 20:53:13 |
言語 | PyPy3 (7.3.15) |
結果 |
TLE
|
実行時間 | - |
コード長 | 1,881 bytes |
コンパイル時間 | 166 ms |
コンパイル使用メモリ | 82,748 KB |
実行使用メモリ | 85,100 KB |
最終ジャッジ日時 | 2025-03-20 20:53:56 |
合計ジャッジ時間 | 7,495 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 10 TLE * 1 -- * 25 |
ソースコード
import sys sys.setrecursionlimit(1 << 25) def is_prime(n): if n < 2: return False for p in [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37]: if n % p == 0: return n == p d = n - 1 s = 0 while d % 2 == 0: d //= 2 s += 1 for a in [2, 325, 9375, 28178, 450775, 9780504, 1795265022]: if a >= n: continue x = pow(a, d, n) if x == 1 or x == n - 1: continue for _ in range(s - 1): x = pow(x, 2, n) if x == n - 1: break else: return False return True def next_prime(start): n = start + 1 while True: if is_prime(n): return n n += 1 def count_totient_le_N(N): if N == 0: return 1 # Only x=1 count = 1 # x=1 def dfs(current_phi, last_p): p = last_p while True: p = next_prime(p) new_phi_part = (p - 1) if current_phi == 0: new_phi = new_phi_part else: new_phi = current_phi * new_phi_part if new_phi > N: break # Add counts for exponents starting from 1 nonlocal count count += 1 temp_phi = new_phi # Check higher exponents for p (e >=2) while True: if temp_phi > N // p: break temp_phi *= p count += 1 dfs(temp_phi, p) # Recursively check new primes with this phi dfs(new_phi, p) try: dfs(0, 1) # Start with current_phi=0 for the first prime except: pass # Handles cases where recursion might end early return count # Read input N = int(sys.stdin.readline()) print(count_totient_le_N(N))