結果
問題 | No.1250 汝は倍数なりや? |
ユーザー |
![]() |
提出日時 | 2020-10-15 23:04:11 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 122 ms / 1,000 ms |
コード長 | 1,776 bytes |
コンパイル時間 | 1,626 ms |
コンパイル使用メモリ | 170,808 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-20 20:05:12 |
合計ジャッジ時間 | 4,775 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 49 |
ソースコード
#include <bits/stdc++.h> //#include <chrono> //#pragma GCC optimize("O3") using namespace std; #define reps(i,s,n) for(int i = s; i < n; i++) #define rep(i,n) reps(i,0,n) #define Rreps(i,n,e) for(int i = n - 1; i >= e; --i) #define Rrep(i,n) Rreps(i,n,0) #define ALL(a) a.begin(), a.end() #define fi first #define se second using ll = long long; using vec = vector<ll>; using mat = vector<vec>; ll N,M,H,W,Q,K,A,B; string S; typedef pair<ll, ll> P; const ll INF = (1LL<<48); P extgcd(ll a, ll b, bool gcd = false){ // return p : a * p.fi + b * p.se = gcd(a, b); struct eq_elem{ // ax + by = u ll u, x, y; eq_elem(ll _a, ll _b, ll _c){ u = _a; x = _b; y = _c; } eq_elem& operator *=(const ll &t){ u *= t; x *= t; y *= t; return *this; } eq_elem operator * (const ll &t){ return eq_elem(*this) *= t; } eq_elem& operator -=(const eq_elem &t){ u -= t.u; x -= t.x; y -= t.y; return *this; } eq_elem operator - (const eq_elem &t){ return eq_elem(*this) -= t; } }; bool flag = a < b; if(flag) swap(a, b); eq_elem s(a, 1, 0), t(b, 0, 1); while(t.u != 0){ s -= t * (s.u / t.u); swap(s, t); } if(gcd){ return make_pair(s.u, 0LL); }else { if(flag) swap(s.x, s.y); return make_pair(s.x, s.y); } } int main() { cin>>N>>K; vec a(N); rep(i, N) {cin>>a[i]; a[i] = abs(a[i]);} bool flag(false); if(*min_element(ALL(a)) == 0 || K == 1) flag = true; rep(i, N){ K /= extgcd(K, a[i], true).fi; if(K == 1) flag = true; if(flag) break; } cout<<(flag ? "YES" : "NO")<<endl; }