結果

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

ソースコード

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::map<INT, int> prime_factorize(INT n, std::vector<INT> lp){
    std::map<INT, int> res;
    while (n != 1){
        res[lp[n]]++;
        n /= lp[n];
    }

    return res;
}

int main(){
    int x;
    cin >> x;
    auto lp = get_min_prime_fact_table(x);
    vector<int> A(x + 1);
    int ma = 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 &[k, v] : pa)
            ca *= (v + 1);
        for (const auto &[k, v] : pb)
            cb *= (v + 1);
        auto s = abs(a - ca - (b - cb));
        A[a] = s;
        ma = min(ma, s); 
    }

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