結果
問題 |
No.2262 Fractions
|
ユーザー |
![]() |
提出日時 | 2023-04-07 01:11:41 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 444 ms / 2,000 ms |
コード長 | 907 bytes |
コンパイル時間 | 224 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 114,440 KB |
最終ジャッジ日時 | 2024-11-27 02:03:23 |
合計ジャッジ時間 | 17,041 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
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)//(n*n) for p in primes: if p>n: return sum(dp) for i in range(n//p,0,-1): dp[i*p]-=dp[i] def solve(): n,k=map(int,input().split()) c=calc(n,n*n-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=-1,n*n while ok-ng>1: mid=(ok+ng)>>1 if calc(n,mid)>=k: ok=mid else: ng=mid for p in range(1,n+1): q=(ok*p)//(n*n) if n*n*q>ng*p: if rev: p,q=q,p print(str(q)+'/'+str(p)) return raise Exception for _ in range(int(input())): solve()