結果

問題 No.2066 Simple Math !
ユーザー GOTKAKOGOTKAKO
提出日時 2024-04-26 00:26:11
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
RE  
実行時間 -
コード長 1,536 bytes
コンパイル時間 2,383 ms
コンパイル使用メモリ 196,980 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-04-26 00:26:24
合計ジャッジ時間 13,067 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 83 ms
6,940 KB
testcase_02 AC 107 ms
6,944 KB
testcase_03 RE -
testcase_04 AC 355 ms
6,944 KB
testcase_05 AC 373 ms
6,940 KB
testcase_06 AC 372 ms
6,940 KB
testcase_07 AC 362 ms
6,940 KB
testcase_08 AC 358 ms
6,940 KB
testcase_09 AC 374 ms
6,940 KB
testcase_10 AC 365 ms
6,944 KB
testcase_11 AC 355 ms
6,944 KB
testcase_12 AC 378 ms
6,940 KB
testcase_13 AC 360 ms
6,940 KB
testcase_14 AC 310 ms
6,944 KB
testcase_15 AC 312 ms
6,940 KB
testcase_16 AC 345 ms
6,940 KB
testcase_17 AC 313 ms
6,944 KB
testcase_18 AC 309 ms
6,940 KB
testcase_19 RE -
testcase_20 RE -
testcase_21 RE -
testcase_22 RE -
testcase_23 RE -
testcase_24 AC 297 ms
6,940 KB
testcase_25 AC 314 ms
6,940 KB
testcase_26 AC 334 ms
6,944 KB
testcase_27 AC 310 ms
6,940 KB
testcase_28 AC 313 ms
6,940 KB
testcase_29 AC 85 ms
6,944 KB
testcase_30 RE -
testcase_31 RE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

__int128_t safemod(__int128_t a,__int128_t m){
    a %= m; a += m;
    return a>=m?a-m:a;
}
__int128_t floor(__int128_t a,__int128_t b){
    if(b < 0) a *= -1,b *= -1;
    return a<0?(a+1)/b-1:a/b; 
}
__int128_t ceil(__int128_t a,__int128_t b){
    if(b < 0) a *= -1,b *= -1;
    return a>0?(a-1)/b+1:a/b; 
}
__int128_t floorsum(__int128_t N,__int128_t M,__int128_t a,__int128_t b){
    //return Σ(i=0~N-1)floor((a*i+b)/M).
    __int128_t ret = 0;
    if(a >= M) ret += a/M*(N*(N-1)/2),a %= M;
    if(b >= M) ret += b/M*N,b %= M;
    __int128_t maxy = a*N+b;
    if(maxy >= M) ret += floorsum(maxy/M,a,M,maxy%M);
    return ret;
}
long long floorsum2(__int128_t N,__int128_t M,__int128_t a,__int128_t b){
    if(N <= 0) return 0;
    __int128_t ret = 0;
    if(a < 0 || a >= M) ret += floor(a,M)*(N*(N-1)/2),a = safemod(a,M);
    if(b < 0 || b >= M) ret += floor(b,M)*N,b = safemod(b,M);
    __int128_t maxy = a*N+b;
    ret += floorsum(maxy/M,a,M,maxy%M);
    return ret;
}

int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);

    int T; cin >> T;
    while(T--){
        long long P,Q,K; cin >> P >> Q >> K;
        long long g = gcd(P,Q); P /= g,Q /= g;
        long long low = 0,high = 1e18;
        while(high-low > 1){
            long long mid = (high+low)/2;
            long long now = floorsum2(min(Q,mid/P+1),Q,-P,mid)+min(Q,mid/P+1);
            if(now >= K+1) high = mid;
            else low = mid;
        }
        cout << high*g << "\n";
    }    
}
0