結果
問題 |
No.2262 Fractions
|
ユーザー |
![]() |
提出日時 | 2023-03-10 19:57:05 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 806 ms / 2,000 ms |
コード長 | 778 bytes |
コンパイル時間 | 378 ms |
コンパイル使用メモリ | 82,432 KB |
実行使用メモリ | 167,168 KB |
最終ジャッジ日時 | 2024-11-27 02:02:07 |
合計ジャッジ時間 | 27,237 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 45 |
ソースコード
from sys import stdin input=lambda :stdin.readline()[:-1] def calc(n,p): dp=[0]*(n+1) for i in range(1,n+1): s=(p*i)>>38 dp[i]+=s if dp[i]==0: continue for j in range(2*i,n+1,i): dp[j]-=dp[i] return sum(dp) 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()