結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 39 ms
52,736 KB
testcase_01 AC 105 ms
76,660 KB
testcase_02 AC 207 ms
77,512 KB
testcase_03 AC 107 ms
76,160 KB
testcase_04 AC 1,032 ms
77,816 KB
testcase_05 AC 1,035 ms
77,768 KB
testcase_06 AC 1,042 ms
78,208 KB
testcase_07 AC 1,032 ms
77,640 KB
testcase_08 AC 1,030 ms
77,516 KB
testcase_09 AC 1,055 ms
77,364 KB
testcase_10 AC 1,066 ms
77,820 KB
testcase_11 AC 1,043 ms
77,304 KB
testcase_12 AC 1,047 ms
77,560 KB
testcase_13 AC 1,041 ms
77,432 KB
testcase_14 AC 985 ms
77,360 KB
testcase_15 AC 984 ms
77,752 KB
testcase_16 AC 977 ms
77,452 KB
testcase_17 AC 993 ms
77,364 KB
testcase_18 AC 984 ms
77,560 KB
testcase_19 AC 434 ms
78,964 KB
testcase_20 AC 428 ms
78,068 KB
testcase_21 AC 437 ms
77,952 KB
testcase_22 AC 430 ms
77,808 KB
testcase_23 AC 429 ms
78,648 KB
testcase_24 AC 626 ms
77,720 KB
testcase_25 AC 612 ms
77,176 KB
testcase_26 AC 619 ms
77,112 KB
testcase_27 AC 624 ms
77,828 KB
testcase_28 AC 621 ms
77,436 KB
testcase_29 AC 379 ms
77,436 KB
testcase_30 AC 109 ms
76,544 KB
testcase_31 AC 110 ms
76,368 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