結果
| 問題 | No.368 LCM of K-products | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2016-04-29 22:52:54 | 
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 328 ms / 2,000 ms | 
| コード長 | 932 bytes | 
| コンパイル時間 | 1,651 ms | 
| コンパイル使用メモリ | 178,236 KB | 
| 実行使用メモリ | 6,944 KB | 
| 最終ジャッジ日時 | 2024-09-22 14:41:50 | 
| 合計ジャッジ時間 | 3,376 ms | 
| ジャッジサーバーID (参考情報) | judge4 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| other | AC * 35 | 
ソースコード
#include <bits/stdc++.h>
using namespace std;
const long long mod = 1e9 + 7;
map<long long, long long> primeFactors(long long n) {
	map<long long, long long> ret;
	for (long long i = 2; i * i <= n; i++)
		while (n % i == 0) ret[i]++, n /= i;
	if (n != 1) ret[n]++;
	return ret;
}
long long modpow(long long a, long long b, long long mod) {
	long long result = 1;
	for (; b > 0; a = a * a % mod, b /= 2) if (b & 1) result = result * a % mod;
	return result;
}
int main() {
	map<long long, vector<long long>> mp;
	int n, K;
	cin >> n >> K;
	for (int i = 0; i < n; i++) {
		int a;
		cin >> a;
		for (auto p : primeFactors(a)) {
			mp[p.first].push_back(p.second);
		}
	}
	long long ans = 1;
	for (auto p : mp) {
		sort(p.second.rbegin(), p.second.rend());
		int sum = 0;
		for (int i = 0; i < min<int>(p.second.size(), K); i++) {
			sum += p.second[i];
		}
		(ans *= modpow(p.first, sum, mod)) %= mod;
	}
	cout << ans << endl;
}
            
            
            
        