結果

問題 No.1250 汝は倍数なりや?
ユーザー mugen_1337mugen_1337
提出日時 2020-10-24 20:22:12
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
MLE  
実行時間 -
コード長 3,900 bytes
コンパイル時間 2,613 ms
コンパイル使用メモリ 217,624 KB
実行使用メモリ 814,592 KB
最終ジャッジ日時 2024-07-21 15:36:03
合計ジャッジ時間 5,400 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 2 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 2 ms
5,376 KB
testcase_14 AC 2 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 2 ms
5,376 KB
testcase_17 AC 3 ms
5,376 KB
testcase_18 AC 2 ms
5,376 KB
testcase_19 AC 2 ms
5,376 KB
testcase_20 AC 3 ms
5,376 KB
testcase_21 MLE -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
testcase_46 -- -
testcase_47 -- -
testcase_48 -- -
testcase_49 -- -
testcase_50 -- -
testcase_51 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>
using namespace std;
#define ALL(x) begin(x),end(x)
#define rep(i,n) for(int i=0;i<(n);i++)
#define debug(v) cout<<#v<<":";for(auto x:v){cout<<x<<' ';}cout<<endl;
#define mod 1000000007
using ll=long long;
const int INF=1000000000;
const ll LINF=1001002003004005006ll;
int dx[]={1,0,-1,0},dy[]={0,1,0,-1};
// ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
template<class T>bool chmax(T &a,const T &b){if(a<b){a=b;return true;}return false;}
template<class T>bool chmin(T &a,const T &b){if(b<a){a=b;return true;}return false;}

struct IOSetup{
    IOSetup(){
        cin.tie(0);
        ios::sync_with_stdio(0);
        cout<<fixed<<setprecision(12);
    }
} iosetup;
 
template<typename T>
ostream &operator<<(ostream &os,const vector<T>&v){
    for(int i=0;i<(int)v.size();i++) os<<v[i]<<(i+1==(int)v.size()?"":" ");
    return os;
}
template<typename T>
istream &operator>>(istream &is,vector<T>&v){
    for(T &x:v)is>>x;
    return is;
}

// 2^127 = 170141183460469231731687303715884105728 ~ 10^38
using int128=__int128;
 
ostream &operator<<(ostream &os,const __int128 n){
    if(n==0){
        os<<0;
        return os;
    }
    __int128 num=n;
    bool neg=false;
    if(num<0)neg=true,num=-num;
    string res="";
    while(num>0){
        res.push_back('0'+num%10);
        num/=10;
    }
    if(neg) res.push_back('-');
    reverse(begin(res),end(res));
    os<<res;
    return os;
}
 
istream &operator>>(istream &is,__int128 &n){
    string s;
    is>>s;
    int idx=0;
    bool neg=false;
    if(s[0]=='-')neg=true,idx++;
    n=0;
    for(;idx<(int)s.size();idx++){
        n=10*n+s[idx]-'0';
    }
    if(neg) n=-n;
    return is;
}
 
__int128 abs128(__int128 x){
    return x<0?-x:x;
}
 
__int128 gcd128(__int128 a,__int128 b){
    if(a==0) return b;
    if(b==0) return a;
    return b>0?gcd128(b,a%b):a;
}
 
/*
pollard rho 素因数分解
乱択の一種，計算量はよくわかっていない(未解決問題)
*/
namespace fastprime{
    int128 pow_mod(int128 a,int128 n,int128 m){
        int128 ret=1;
        while(n){
            if(n&1)ret=ret*a%m;
            a=a*a%m;
            n>>=1;
        }
        return ret;
    }
 
    bool miller_rabin(int128 n){
        if(n<=1) return false;
        if(n==2) return true;
        if(n%2==0) return false;
        int128 d=n-1,s=0;
        while(!(d&1))d>>=1,s++;// n-1=2^s*d
        for(int128 a:{2,3,5,7,11,13,17,19,23,29,31,37}){
            if(n<=a) break;
            int128 y=pow_mod(a,d,n);// y=a^d (mod n)
            int128 t;
            for(t=d;t<n-1 and y!=1 and y!=n-1;t<<=1) y=y*y%n;
            if(y!=n-1 and t%2==0) return false;
        }
        return true;
    }
 
    int128 pollard_rho_single(int128 n){
        auto f=[&](int128 x){return (x*x+1)%n;};
        if(miller_rabin(n)) return n;
        if(!(n&1)) return 2;
        int128 a=0;
        while(true){
            a++;
            int128 x=a;
            int128 y=f(x);
            while(true){
                int128 g=gcd128(y-x+n,n);
                if(g==0 or g==n) break;
                if(g!=1) return g;
                x=f(x);
                y=f(f(y));
            }
        }
    }
 
    vector<int128> pollard_rho(int128 n){
        if(n==1) return {};
        int128 x=pollard_rho_single(n);
        if(x==n) return {x};
        vector<int128> l=pollard_rho(x);
        vector<int128> r=pollard_rho(n/x);
        l.insert(l.end(),r.begin(),r.end());
        return l;
    }
}


signed main(){
    int n;ll h;cin>>n>>h;
    vector<ll> a(n);
    cin>>a;
    for(auto &x:a) x=abs(x);

    auto vh=fastprime::pollard_rho(h);
    map<int128,int> cnt;
    for(auto &x:vh) cnt[x]++;

    rep(i,n){
        auto v=fastprime::pollard_rho(a[i]);
        for(auto &x:v)if(cnt.count(x)) cnt[x]--;
    }

    for(auto [u,v]:cnt)if(v>0){
        cout<<"NO"<<endl;
        return 0;
    }
    cout<<"YES"<<endl;
    return 0;
}
0