結果

問題 No.368 LCM of K-products
ユーザー face4
提出日時 2019-02-13 10:59:01
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 40 ms / 2,000 ms
コード長 1,707 bytes
コンパイル時間 740 ms
コンパイル使用メモリ 80,136 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-22 14:58:49
合計ジャッジ時間 2,271 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

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