結果

問題 No.368 LCM of K-products
ユーザー odan3240
提出日時 2016-05-01 11:23:05
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 343 ms / 2,000 ms
コード長 1,522 bytes
コンパイル時間 895 ms
コンパイル使用メモリ 89,496 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-22 14:44:53
合計ジャッジ時間 2,641 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string>
#include <vector>
#include <sstream>
#include <map>
#include <set>
#include <queue>
#include <algorithm>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>

using namespace std;
using ll = long long;

#define all(c) (c).begin(), (c).end()
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define pb(e) push_back(e)
#define mp(a, b) make_pair(a, b)
#define fr first
#define sc second

const ll INF = 1e9;
const ll MOD = 1e9 + 7;
int dx[4] = {1, 0, -1, 0};
int dy[4] = {0, 1, 0, -1};

int N, K;
int a[1003];
map<int, vector<int>> mp;

void factor(int n) {
    ll i = 2;
    while (n != 1) {
        int cnt = 0;
        while (n % i == 0) {
            cnt++;
            n /= i;
        }
        if (cnt) {
            mp[i].push_back(cnt);
        }
        if (i * i > n) {
            mp[n].push_back(1);
            return;
        }
        i++;
    }
}

ll mod_pow(ll a, int n, int m) {
    if (n == 0)
        return 1;
    else if (n % 2 == 1)
        return (a * mod_pow((a * a) % m, n >> 1, m)) % m;
    else
        return mod_pow((a * a) % m, n >> 1, m);
}

int main() {
    cin >> N >> K;
    rep(i, N) cin >> a[i];
    rep(i, N) factor(a[i]);

    ll ans = 1;
    for (auto &e : mp) {
        sort(all(e.second));
        reverse(all(e.second));
        int cnt = 0;
        rep(i, min(K, (int)e.second.size())) cnt += e.second[i];
        ans *= mod_pow(e.first, cnt, MOD);
        ans %= MOD;
    }

    cout << ans << endl;

    return 0;
}
0