結果

問題 No.2066 Simple Math !
ユーザー とりゐとりゐ
提出日時 2022-08-02 16:10:53
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,027 ms / 2,000 ms
コード長 1,009 bytes
コンパイル時間 196 ms
コンパイル使用メモリ 82,416 KB
実行使用メモリ 78,652 KB
最終ジャッジ日時 2024-04-27 15:40:54
合計ジャッジ時間 22,293 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 34 ms
52,880 KB
testcase_01 AC 97 ms
76,116 KB
testcase_02 AC 191 ms
77,236 KB
testcase_03 AC 105 ms
76,216 KB
testcase_04 AC 1,027 ms
77,648 KB
testcase_05 AC 986 ms
77,616 KB
testcase_06 AC 987 ms
77,988 KB
testcase_07 AC 981 ms
77,616 KB
testcase_08 AC 976 ms
77,512 KB
testcase_09 AC 1,005 ms
77,256 KB
testcase_10 AC 1,025 ms
78,020 KB
testcase_11 AC 994 ms
77,000 KB
testcase_12 AC 1,000 ms
77,252 KB
testcase_13 AC 1,011 ms
76,976 KB
testcase_14 AC 952 ms
77,380 KB
testcase_15 AC 971 ms
77,380 KB
testcase_16 AC 958 ms
77,120 KB
testcase_17 AC 971 ms
77,092 KB
testcase_18 AC 967 ms
77,252 KB
testcase_19 AC 421 ms
78,652 KB
testcase_20 AC 401 ms
77,756 KB
testcase_21 AC 415 ms
77,896 KB
testcase_22 AC 405 ms
77,740 KB
testcase_23 AC 413 ms
78,640 KB
testcase_24 AC 640 ms
77,380 KB
testcase_25 AC 591 ms
77,120 KB
testcase_26 AC 586 ms
76,976 KB
testcase_27 AC 587 ms
77,388 KB
testcase_28 AC 591 ms
77,252 KB
testcase_29 AC 336 ms
77,112 KB
testcase_30 AC 95 ms
76,164 KB
testcase_31 AC 104 ms
76,280 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

def floor_sum(n,m,a,b):
  ans=0
  while True:
    if a>=m:
      ans+=(n-1)*n*(a//m)//2
      a%=m
    if b>=m:
      ans+=n*(b//m)
      b%=m
    y_max=(a*n+b)//m
    x_max=(y_max*m-b)
    if y_max==0:
      return ans
    ans+=(n-(x_max+a-1)//a)*y_max
    n,m,a,b=y_max,a,m,(a-x_max%a)%a

def calc1(p,q,n):
  x=(q*n)//p
  return (x+1)*n-floor_sum(x+1,q,p,0)

def calc2(p,q,n,m):
  ng,ok=q,0
  while abs(ok-ng)>1:
    mid=(ok+ng)//2
    if floor_sum(n,q,p,q-mid)-floor_sum(n,q,p,0)>=n+1-m:
      ok=mid
    else:
      ng=mid
  return ok

import math

def solve(p,q,k):
  g=math.gcd(p,q)
  p,q=min(p,q)//g,max(p,q)//g
  m=q*(p-1)-(p-1)*(q-1)//2
  if k>m:
    res=k-m-1+q*(p-1)
    return res*g

  ng,ok=0,p+1
  while abs(ok-ng)>1:
    mid=(ok+ng)//2
    if k<=calc1(p,q,mid):
      ok=mid
    else:
      ng=mid
  
  if ok==0:
    x=k
  else:
    x=k-calc1(p,q,ok-1)
  
  return (q*(ok-1)+calc2(p,q,(ok*q)//p+1,x))*g

t=int(input())
for _ in range(t):
  p,q,k=map(int,input().split())
  print(solve(p,q,k+1))
0