結果

問題 No.811 約数の個数の最大化
ユーザー peroonperoon
提出日時 2019-04-19 11:40:10
言語 C++14
(gcc 10.1.0 + boost 1.73.0)
結果
AC  
実行時間 107 ms / 2,000 ms
コード長 1,585 Byte
コンパイル時間 1,591 ms
使用メモリ 1,596 KB
最終ジャッジ日時 2020-08-07 04:44:32
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
使用メモリ
testcase_00 AC 2 ms
1,596 KB
testcase_01 AC 3 ms
1,596 KB
testcase_02 AC 31 ms
1,596 KB
testcase_03 AC 2 ms
1,584 KB
testcase_04 AC 2 ms
1,596 KB
testcase_05 AC 3 ms
1,596 KB
testcase_06 AC 4 ms
1,596 KB
testcase_07 AC 6 ms
1,592 KB
testcase_08 AC 15 ms
1,592 KB
testcase_09 AC 24 ms
1,596 KB
testcase_10 AC 12 ms
1,592 KB
testcase_11 AC 20 ms
1,592 KB
testcase_12 AC 8 ms
1,592 KB
testcase_13 AC 107 ms
1,588 KB
testcase_14 AC 29 ms
1,596 KB
権限があれば一括ダウンロードができます

ソースコード

diff #
#include<bits/stdc++.h>
using namespace std;
using ll = long long;

#define FOR(i,a,b) for(ll i=(a);i<(b);++i)
#define ALL(v) (v).begin(), (v).end()
#define p(s) cout<<(s)<<endl
#define p2(s, t) cout << (s) << " " << (t) << endl
#define br() p("")
#define pn(s) cout << (#s) << " " << (s) << endl
#define p_yes() p("Yes")
#define p_no() p("No")

const ll mod = 1e9 + 7;
const ll inf = 1e18;

ll gcd(ll a,ll b){
    if(b == 0) return a;
    return gcd(b,a%b);
}

map<ll, ll> factorize(ll n){
    map<ll, ll> mp;
    ll sq = sqrt(n);
    FOR(i, 2, sq+1){
        while(n%i==0){
            mp[i]++;
            n/=i;
        }
    }
    // 残り
    if(n!=1){
        mp[n]++;
    }
    return mp;
}

// 素因数の個数 (重複あり) ex. {2, 2, 2, 3}
ll prime_factor_num(ll a){
    ll count = 0;
    auto mp = factorize(a);
    for(auto p : mp){
        count += p.second;
    }
    return count;
}

// a, bの共通の素因数の数
ll common(ll a, ll b){
    ll g = gcd(a, b);
    return prime_factor_num(g);
}

// 約数の個数
ll divisor_num(ll a){
    auto mp = factorize(a);
    ll ret = 1;
    for(auto p : mp){
        ret *= (p.second + 1);
    }
    return ret;
}

int main(){
    cin.tie(0);
    ios::sync_with_stdio(false);

    // input
    ll N, K;
    cin >> N >> K;

    ll max_divisor_num = 0;
    ll max_M;
    FOR(i, 1, N){
        if(common(i, N)>=K){
            ll div = divisor_num(i);
            if(div > max_divisor_num){
                max_divisor_num = div;
                max_M = i;
            }
        }
    }

    p(max_M);
    
    return 0;
}
0