結果

問題 No.368 LCM of K-products
コンテスト
ユーザー latte0119
提出日時 2016-04-29 22:34:55
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,614 bytes
コンパイル時間 1,937 ms
コンパイル使用メモリ 171,056 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-10-04 18:15:15
合計ジャッジ時間 7,010 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 32 WA * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;

#define int long long
typedef 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;
    set<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.insert(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*(*it)%mod;
    cout<<ans<<endl;
    return 0;
}
0