結果

問題 No.2262 Fractions
ユーザー Carpenters-CatCarpenters-Cat
提出日時 2023-04-07 23:36:44
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 1,745 bytes
コンパイル時間 2,228 ms
コンパイル使用メモリ 205,440 KB
実行使用メモリ 44,096 KB
最終ジャッジ日時 2024-10-02 20:44:43
合計ジャッジ時間 25,038 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 TLE -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
int main () {
    int T;
    cin >> T;
    vector<vector<int>> X(300030);
    for (int i = 1; i <= 300000; i ++) {
        for (int j = 1; i * j <= 300000; j ++) {
            X[i * j].push_back(i);
        }
    }
    int meb[300030];
    fill(meb, meb + 300030, -2);
    meb[1] = 1;
    for (int i = 2; i <= 300000; i ++) {
        if (meb[i] == -2) {
            meb[i] = -1;
        }
        for (int j = 2; i * j <= 300000; j ++) {
            meb[i * j] *= -1;
            if (j % i == 0) {
                meb[i * j] = 0;
            }
        }
    }
    while (T --) {
        ll N, K;
        cin >> N >> K;
        if (N == 1) {
            puts("1/1");
            continue;
        }
        double l = 0, r = N + 1;
        bool ok = false;
        while(true) {
            double mu = (l + r) / 2;
            ll sum = 0;
            for (ll m = 1; m <= N; m ++) {
                ll n = min(N, (ll)(mu * m));
                for (auto a : X[m]) {
                    sum += (n / a) * meb[a];
                }
            }
            if (sum <= K) {
                l = mu;
            } else {
                r = mu;
            }
            if (sum == K) {
                ok = true;
                break;
            }
            if (l > N) {
                break;
            }
        }
        if (ok) {
            ll s = 0, b = 1;
            for (ll m = 1; m <= N; m ++) {
                ll a = (ll)(m * l);
                a = min(a, N);
                if (b * a >= m * s) {
                    s = a; b = m;
                }
            }
            printf("%lld/%lld\n", s, b);
        } else {
            puts("-1");
        }
    }
}
0