結果

問題 No.368 LCM of K-products
ユーザー face4face4
提出日時 2019-02-13 10:59:01
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 41 ms / 2,000 ms
コード長 1,707 bytes
コンパイル時間 765 ms
コンパイル使用メモリ 79,108 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-23 21:29:03
合計ジャッジ時間 2,303 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 28 ms
4,348 KB
testcase_01 AC 36 ms
4,348 KB
testcase_02 AC 35 ms
4,348 KB
testcase_03 AC 36 ms
4,348 KB
testcase_04 AC 36 ms
4,348 KB
testcase_05 AC 41 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 3 ms
4,348 KB
testcase_09 AC 3 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 2 ms
4,348 KB
testcase_12 AC 2 ms
4,348 KB
testcase_13 AC 20 ms
4,348 KB
testcase_14 AC 31 ms
4,348 KB
testcase_15 AC 34 ms
4,348 KB
testcase_16 AC 30 ms
4,348 KB
testcase_17 AC 25 ms
4,348 KB
testcase_18 AC 35 ms
4,348 KB
testcase_19 AC 11 ms
4,348 KB
testcase_20 AC 23 ms
4,348 KB
testcase_21 AC 7 ms
4,348 KB
testcase_22 AC 33 ms
4,348 KB
testcase_23 AC 3 ms
4,348 KB
testcase_24 AC 3 ms
4,348 KB
testcase_25 AC 2 ms
4,348 KB
testcase_26 AC 2 ms
4,348 KB
testcase_27 AC 3 ms
4,348 KB
testcase_28 AC 3 ms
4,348 KB
testcase_29 AC 3 ms
4,348 KB
testcase_30 AC 2 ms
4,348 KB
testcase_31 AC 3 ms
4,348 KB
testcase_32 AC 3 ms
4,348 KB
testcase_33 AC 15 ms
4,348 KB
testcase_34 AC 35 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
typedef long long ll;

const ll mod = 1000000007;

ll modpow(ll a, ll b, ll p = 1e9+7){
    if(b == 0)  return 1;

    if(b % 2 == 0){
        ll d = modpow(a, b/2, p);
        return (d*d) % p;
    }else{
        return (a%p * modpow(a, b-1, p)) % p;
    }
}

int main(){
    vector<int> p;
    bool nonp[100001] = {};
    nonp[0] = nonp[1] = true;
    for(int i = 2; i < 100001; i++){
        if(nonp[i]) continue;
        p.push_back(i);
        for(int j = i+i; j < 100001; j += i)    nonp[j] = true;
    }

    int n, k;
    cin >> n >> k;
    vector<int> a(n);
    for(int i = 0; i < n; i++)  cin >> a[i];

    ll ans = 1;
    // 素因数毎に独立に考える
    for(int prime : p){
        vector<int> num;
        for(int i = 0; i < n; i++){
            int cnt = 0;
            while(a[i]%prime == 0)  a[i] /= prime, cnt++;
            if(cnt) num.push_back(cnt);
        }
        sort(num.rbegin(), num.rend());
        ll sum = 0;
        for(int i = 0; i < min((int)num.size(), k); i++)    sum += num[i];
        ans *= modpow((ll)prime, sum, mod);
        ans %= mod;
    }
    
    for(int i = 0; i < n; i++){
        if(a[i] == 1)   continue;
        vector<int> num;
        ll prime = a[i];
        for(int j = i; j < n; j++){
            int cnt = 0;
            while(a[j]%prime == 0)  a[j] /= prime, cnt++;
            if(cnt) num.push_back(cnt);
        }
        sort(num.rbegin(), num.rend());
        ll sum = 0;
        for(int i = 0; i < min((int)num.size(), k); i++)    sum += num[i];
        ans *= modpow((ll)prime, sum, mod);
        ans %= mod;
    }

    cout << ans << endl;

    return 0;
}
0