結果
問題 | No.368 LCM of K-products |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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 secondconst 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;elsereturn 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;}