結果

問題 No.368 LCM of K-products
ユーザー YazatenYazaten
提出日時 2016-04-30 09:46:16
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 146 ms / 2,000 ms
コード長 1,053 bytes
コンパイル時間 1,875 ms
コンパイル使用メモリ 184,924 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-22 14:44:42
合計ジャッジ時間 3,166 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include "bits/stdc++.h"
#include <cassert>
using namespace std;
typedef long long ll;
typedef pair<int,int> pii;
#define rep(i,n) for(ll i=0;i<(ll)(n);i++)
#define all(a)  (a).begin(),(a).end()
#define vi vector<int>
#define pb push_back
#define INF 999999999
//#define INF (1LL<<59)

ll power(ll a,ll b/*,ll mod*/){
	ll ret=1;
	if(b>0){
		ret = power(a,b/2/*,mod*/);
		if(b%2==0)ret = (ret*ret)/*%mod*/;
		else ret = (((ret*ret)/*%mod*/)*a)/*%mod*/;
	}
	return ret;
}

map<int,int> prime_factorization(int n){
	map<int,int> ret;
	
	for(int div=2;div<=1+sqrt(n);div++){
		while(n%div==0){
			ret[div]++;
			n/=div;
		}
	}
	if(n!=1)ret[n]++;
	return ret;
}

int main(){
	map<int,vector<int>> mp;
	
	int n,k;
	cin>>n>>k;
	rep(i,n){
		int a;
		cin>>a;
		map<int,int> res = prime_factorization(a);
		for(auto &p : res){
			mp[p.first].pb(p.second);
		}
	}
	
	ll ans = 1;
	for(auto &elm : mp){
		sort(all(elm.second),greater<int>());
		rep(i,min<ll>(k,elm.second.size())){
			ans = (ans*power(elm.first,elm.second[i]))%1000000007;
		}
	}
	cout<<ans<<endl;
}
0