結果

問題 No.1006 Share an Integer
コンテスト
ユーザー qLethon
提出日時 2020-03-06 23:27:14
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 717 ms / 2,000 ms
コード長 1,353 bytes
コンパイル時間 2,193 ms
コンパイル使用メモリ 196,492 KB
最終ジャッジ日時 2025-01-09 05:36:27
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

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

template<class INT>
std::vector<INT> get_min_prime_fact_table(INT n){
    std::vector<INT> lp(n + 1), P;
    for (INT i = 2; i <= n; i++){
        if (lp[i] == 0){
            lp[i] = i;
            P.push_back(i);
        }
        for (INT j = 0; j <= P.size() and P[j] <= lp[i] and i * P[j] <= n; j++)
            lp[i * P[j]] = P[j];
    }

    return lp;
}

template<class INT>
std::vector<int> prime_factorize(INT n, std::vector<INT> &lp){
    vector<int> res(1);
    int last = 2;
    while (n != 1){
        if (last != lp[n]){
            res.push_back(0);
            last = lp[n];
        }
        n /= lp[n];
        res[res.size() - 1]++;
    }

    return res;
}

int main(){
    int x;
    cin >> x;
    auto lp = get_min_prime_fact_table(x);
    vector<int> A(x + 1);
    int mi = 1e9;
    for (int a = 1; a < x; a++){
        auto b = x - a;
        auto pa = prime_factorize(a, lp);
        auto pb = prime_factorize(b, lp);
        int ca = 1, cb = 1;
        for (const auto &v : pa)
            ca *= (v + 1);
        for (const auto &v : pb)
            cb *= (v + 1);
        auto s = abs(a - ca - (b - cb));
        A[a] = s;
        mi = min(mi, s); 
    }

    for (int i = 1; i < x; i++){
        if (A[i] == mi){
            cout << i << " " << x - i << endl;
        }
    }
}
0