結果

問題 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
権限があれば一括ダウンロードができます

ソースコード

diff #

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()
0