結果

問題 No.368 LCM of K-products
ユーザー 🐬hec🐬hec
提出日時 2016-04-04 13:16:45
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 78 ms / 2,000 ms
コード長 2,430 bytes
コンパイル時間 2,302 ms
コンパイル使用メモリ 177,800 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-23 21:07:39
合計ジャッジ時間 3,809 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 58 ms
4,348 KB
testcase_01 AC 77 ms
4,348 KB
testcase_02 AC 78 ms
4,348 KB
testcase_03 AC 78 ms
4,348 KB
testcase_04 AC 78 ms
4,348 KB
testcase_05 AC 78 ms
4,348 KB
testcase_06 AC 2 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 2 ms
4,348 KB
testcase_09 AC 2 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 2 ms
4,348 KB
testcase_12 AC 2 ms
4,348 KB
testcase_13 AC 39 ms
4,348 KB
testcase_14 AC 63 ms
4,348 KB
testcase_15 AC 69 ms
4,348 KB
testcase_16 AC 61 ms
4,348 KB
testcase_17 AC 51 ms
4,348 KB
testcase_18 AC 72 ms
4,348 KB
testcase_19 AC 18 ms
4,348 KB
testcase_20 AC 47 ms
4,348 KB
testcase_21 AC 12 ms
4,348 KB
testcase_22 AC 68 ms
4,348 KB
testcase_23 AC 2 ms
4,348 KB
testcase_24 AC 2 ms
4,348 KB
testcase_25 AC 2 ms
4,348 KB
testcase_26 AC 2 ms
4,348 KB
testcase_27 AC 2 ms
4,348 KB
testcase_28 AC 2 ms
4,348 KB
testcase_29 AC 2 ms
4,348 KB
testcase_30 AC 2 ms
4,348 KB
testcase_31 AC 2 ms
4,348 KB
testcase_32 AC 3 ms
4,348 KB
testcase_33 AC 28 ms
4,348 KB
testcase_34 AC 77 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

#define _overload(_1,_2,_3,name,...) name
#define _rep(i,n) _range(i,0,n)
#define _range(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload(__VA_ARGS__,_range,_rep,)(__VA_ARGS__)

#define _rrep(i,n) _rrange(i,n,0)
#define _rrange(i,a,b) for(int i=int(a)-1;i>=int(b);--i)
#define rrep(...) _overload(__VA_ARGS__,_rrange,_rrep,)(__VA_ARGS__)

#define _all(arg) begin(arg),end(arg)
#define uniq(arg) sort(_all(arg)),(arg).erase(unique(_all(arg)),end(arg))
#define getidx(ary,key) lower_bound(_all(ary),key)-begin(ary)
#define clr(a,b) memset((a),(b),sizeof(a))
#define bit(n) (1LL<<(n))
#define popcount(n) (__builtin_popcountll(n))

template<class T>bool chmax(T &a, const T &b) { return (a<b)?(a=b,1):0;}
template<class T>bool chmin(T &a, const T &b) { return (b<a)?(a=b,1):0;}

using namespace std;
using ll=long long;
const ll mod=1000000007LL;

inline ll extgcd(ll a,ll b,ll& x,ll& y){x=1,y=0;ll g=a;if(b!=0) g=extgcd(b,a%b,y,x),y-=a/b*x;return g;}
inline ll ADD(const ll &a, const ll &b,const ll &mod) { return a+b<mod?a+b:a+b-mod;}
inline ll SUB(const ll &a, const ll &b,const ll &mod) { return a-b>=0?a-b:a-b+mod;}
inline ll MUL(const ll &a, const ll &b,const ll &mod) { return (1LL*a*b)%mod;}
inline ll INV(ll a,ll mod){ll x,y;extgcd(a,mod,x,y);return (x%mod+mod)%mod;}
inline ll DIV(const ll &a, const ll &b,const ll &mod) {return MUL(a,INV(b,mod),mod);}
inline ll POW(ll a,ll n,ll mod){ll b=1LL;for(a%=mod;n;a=MUL(a,a,mod),n>>=1)if(n&1) b=MUL(b,a,mod); return b;}

const int limit=40000;
bool prime[limit+1];

void init(){
	rep(i,limit) prime[i]=true;
	prime[0]=prime[1]=false;
	for(ll i=2;i*i<=limit;++i){
		if(prime[i])
			for(ll j=i*i;j<=limit;j+=i)
				prime[j]=false;
	}
	return;
}

ll a[1010];


int main(void){
	init();

	int n,k;
	cin >> n >> k;
	rep(i,n) cin >> a[i];
	ll ans=1LL;

	ll factor[1010];
	for(int f=2;f<=limit;++f){
		if(prime[f]==false) continue;
		rep(i,n){
			factor[i]=0;
			while(a[i]%f==0) a[i]/=f,factor[i]++;
		}
		sort(factor,factor+n);
		reverse(factor,factor+n);
		ll sum=0LL;
		rep(i,k) sum+=factor[i];
		ans=MUL(ans,POW(f,sum,mod),mod);
	}

	map<ll,int> remain;
	rep(i,n){
		if(a[i]==1) continue;
		remain[a[i]]++;
	} 

	for(auto &it:remain){
		ll index=min<ll>(it.second,k);
		ans=MUL(ans,POW(it.first,index,mod),mod);
	}

	cout << ans << endl;
	return 0;
}

// 考えられる危険なケース
// big prime numbers
// ソートの順番のミス
0