結果
問題 |
No.2262 Fractions
|
ユーザー |
![]() |
提出日時 | 2023-03-20 05:09:48 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 349 ms / 2,000 ms |
コード長 | 920 bytes |
コンパイル時間 | 353 ms |
コンパイル使用メモリ | 82,304 KB |
実行使用メモリ | 114,580 KB |
最終ジャッジ日時 | 2024-11-27 02:02:34 |
合計ジャッジ時間 | 14,196 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 45 |
ソースコード
from sys import stdin input=lambda :stdin.readline()[:-1] M=3*10**5+100 sieve=[0]*2+[1]*M primes=[] for p in range(2,M): if sieve[p]: primes.append(p) for i in range(2*p,M,p): sieve[i]=0 def calc(n,p): dp=[0]*(n+1) for i in range(1,n+1): dp[i]=(p*i)>>38 for p in primes: if p>n: return sum(dp) for i in range(n//p,0,-1): dp[i*p]-=dp[i] m=1<<38 def solve(): n,k=map(int,input().split()) c=calc(n,m-1) if k==c+1: print('1/1') return if k<=c: rev=False elif k<=2*c+1: rev=True k=2*c+2-k else: print(-1) return ng,ok=0,1<<38 for _ in range(38): mid=(ok+ng)>>1 if calc(n,mid)>=k: ok=mid else: ng=mid ans=-1 for p in range(1,n+1): q=(ok*p)>>38 if m*q>ng*p: if rev: p,q=q,p print(str(q)+'/'+str(p)) return raise Exception for _ in range(int(input())): solve()