結果

問題 No.2066 Simple Math !
ユーザー とりゐとりゐ
提出日時 2022-08-02 16:10:53
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 1,126 ms / 2,000 ms
コード長 1,009 bytes
コンパイル時間 1,225 ms
コンパイル使用メモリ 87,192 KB
実行使用メモリ 80,612 KB
最終ジャッジ日時 2023-08-09 21:12:46
合計ジャッジ時間 26,925 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 72 ms
71,404 KB
testcase_01 AC 136 ms
78,304 KB
testcase_02 AC 241 ms
80,000 KB
testcase_03 AC 135 ms
78,212 KB
testcase_04 AC 1,094 ms
79,100 KB
testcase_05 AC 1,088 ms
78,564 KB
testcase_06 AC 1,099 ms
79,164 KB
testcase_07 AC 1,090 ms
79,664 KB
testcase_08 AC 1,088 ms
78,596 KB
testcase_09 AC 1,126 ms
79,940 KB
testcase_10 AC 1,114 ms
79,044 KB
testcase_11 AC 1,106 ms
78,524 KB
testcase_12 AC 1,111 ms
79,128 KB
testcase_13 AC 1,117 ms
78,880 KB
testcase_14 AC 1,058 ms
79,044 KB
testcase_15 AC 1,067 ms
78,884 KB
testcase_16 AC 1,050 ms
78,736 KB
testcase_17 AC 1,053 ms
78,904 KB
testcase_18 AC 1,060 ms
79,188 KB
testcase_19 AC 475 ms
80,032 KB
testcase_20 AC 470 ms
79,736 KB
testcase_21 AC 477 ms
80,612 KB
testcase_22 AC 474 ms
79,784 KB
testcase_23 AC 488 ms
79,896 KB
testcase_24 AC 693 ms
78,876 KB
testcase_25 AC 668 ms
78,768 KB
testcase_26 AC 674 ms
78,320 KB
testcase_27 AC 670 ms
79,048 KB
testcase_28 AC 669 ms
78,668 KB
testcase_29 AC 406 ms
79,336 KB
testcase_30 AC 136 ms
77,928 KB
testcase_31 AC 132 ms
78,148 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