結果
問題 | No.368 LCM of K-products |
ユーザー |
![]() |
提出日時 | 2016-04-29 22:37:28 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
WA
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,627 bytes |
コンパイル時間 | 1,422 ms |
コンパイル使用メモリ | 171,256 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-10-04 18:18:35 |
合計ジャッジ時間 | 6,630 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 34 WA * 1 |
ソースコード
#include<bits/stdc++.h>using namespace std;#define int long longtypedef pair<int,int>pint;typedef vector<int>vint;typedef vector<pint>vpint;#define pb push_back#define mp make_pair#define fi first#define se second#define all(v) (v).begin(),(v).end()#define rep(i,n) for(int i=0;i<(n);i++)#define reps(i,f,n) for(int i=(f);i<(n);i++)#define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)template<class T,class U>inline void chmin(T &t,U f){if(t>f)t=f;}template<class T,class U>inline void chmax(T &t,U f){if(t<f)t=f;}int N,K;int A[1000];vint prime;int f[111111];const int mod=1000000007;int mpow(int n,int m){int ret=1;while(m){if(m&1)ret=ret*n%mod;n=n*n%mod;m>>=1;}return ret;}signed main(){fill_n(f,111111,true);f[0]=f[1]=false;for(int i=2;i<111111;i++){if(!f[i])continue;prime.pb(i);for(int j=i+i;j<111111;j+=i)f[j]=false;}cin>>N>>K;map<int,int>hoge;rep(i,N){cin>>A[i];int a=A[i];rep(j,prime.size()){while(a%prime[j]==0)a/=prime[j];}if(a!=1)hoge[a]++;}int ans=1;rep(i,prime.size()){vint v;rep(j,N){int a=A[j];v.pb(0);while(a%prime[i]==0){v[v.size()-1]++;a/=prime[i];}}sort(all(v));int sum=0;rep(j,K)sum+=v[N-1-j];ans=ans*mpow(prime[i],sum)%mod;}each(it,hoge)ans=ans*mpow(it->fi,it->se)%mod;cout<<ans<<endl;return 0;}