#include #include namespace my{ using namespace std; #define eb emplace_back #define done(...) return pp(__VA_ARGS__) #define LL(...) ll __VA_ARGS__;lin(__VA_ARGS__) #define ST(...) string __VA_ARGS__;lin(__VA_ARGS__) #define FO(n) for(ll ij=n;ij--;) #define FOR(i,...) for(auto[i,i##stop,i##step]=range(0,__VA_ARGS__);isync_with_stdio(0);cout<auto range(bool s,A...a){arrayr{0,0,1};ll I=0;((r[I++]=a),...);if(!s&&I==1)swap(r[0],r[1]);r[0]-=s;return r;} constexpr char nl=10; constexpr char sp=32; constexpr auto Yes(bool y=1){return y?"Yes":"No";} constexpr auto No(){return Yes(0);} templatestruct pair{ A a;B b; pair()=default; pair(A a,B b):a(a),b(b){} pair(const std::pair&p):a(p.first),b(p.second){} auto operator<=>(const pair&)const=default; friend ostream&operator<<(ostream&o,const pair&p){return o<>auto&sort(auto&a,const F&f={}){ranges::sort(a,f);return a;} templateostream&operator<<(ostream&o,const std::pair&p){return o<ostream&operator<<(ostream&o,const array&a){fo(i,n)o<concept vectorial=is_base_of_v,V>; templateistream&operator>>(istream&i,vector&v){fe(v,e)i>>e;return i;} templateostream&operator<<(ostream&o,const vector&v){fe(v,e)o<?nl:sp);return o;} templatestruct vec:vector{ using vector::vector; vec(const vector&v){vector::operator=(v);} vec&operator^=(const vec&u){this->insert(this->end(),u.begin(),u.end());return*this;} vec operator^(const vec&u)const{return vec{*this}^=u;} vec&operator++(){fe(*this,e)++e;return*this;} vec&operator--(){fe(*this,e)--e;return*this;} }; void lin(auto&...a){(cin>>...>>a);} templatevoid pp(const auto&...a){ll n=sizeof...(a);((cout<0,c)),...);cout<auto rle(const vec&a){vec>r;fe(a,e)r.size()&&e==r.back().a?++r.back().b:r.eb(e,1).b;return r;} templateauto rce(veca){return rle(sort(a));} ll rand(auto...a){arrayv{0,0};ll I=0;((v[I++]=a),...);auto[l,r]=v;if(I==1)swap(l,r);static ll t=495;t^=t<<7,t^=t>>9;return l>64; return r>=N?r-N:r; } auto&operator+=(const modular&b){if((a+=b.a)>=N)a-=N;return*this;} auto&operator-=(const modular&b){if(i64(a-=b.a)<0)a+=N;return*this;} auto&operator*=(const modular&b){a=reduce(u128(a)*b.a);return*this;} auto&operator/=(const modular&b){*this*=b.inv();return*this;} friend auto operator+(const modular&a,const modular&b){return modular{a}+=b;} friend auto operator-(const modular&a,const modular&b){return modular{a}-=b;} friend auto operator*(const modular&a,const modular&b){return modular{a}*=b;} friend auto operator/(const modular&a,const modular&b){return modular{a}/=b;} friend bool operator==(const modular&a,const modular&b){return a.a==b.a;} auto operator-()const{return modular{}-modular{*this};} modular pow(u128 n)const{ modular r{1},x{*this}; while(n){ if(n&1)r*=x; x*=x; n>>=1; } return r; } modular inv()const{u64 a=val(),b=N,u=1,v=0;assert(gcd(a,b)==1);while(b)swap(u-=a/b*v,v),swap(a-=a/b*b,b);return u;} u64 val()const{return reduce(a);} friend istream&operator>>(istream&i,montgomery64&b){ll t;i>>t;b=t;return i;} friend ostream&operator<<(ostream&o,const montgomery64&b){return o<as){ ll d=n-1; while(~d&1)d>>=1; using modular=montgomery64; auto pre_mod=modular::mod(); modular::set_mod(n); modular one=1,minus_one=n-1; fe(as,a){ if(a%n==0)continue; ll t=d; modular y=modular(a).pow(t); while(t!=n-1&&y!=one&&y!=minus_one)y*=y,t<<=1; if(y!=minus_one&&~t&1)return modular::set_mod(pre_mod),0; } return modular::set_mod(pre_mod),1; } bool is_prime(ll n){ if(~n&1)return n==2; if(n<=1)return 0; if(n<4759123141LL)return miller_rabin(n,{2,7,61}); return miller_rabin(n,{2,325,9375,28178,450775,9780504,1795265022}); } ll pollard_rho(ll n){ if(~n&1)return 2; if(is_prime(n))return n; using modular=montgomery64; auto pre_mod=modular::mod(); modular::set_mod(n); modular R,one=1; auto f=[&](const modular&x){return x*x+R;}; while(1){ modular x,y,ys,q=one; R=rand(2,n),y=rand(2,n); ll g=1; constexpr ll m=128; for(ll r=1;g==1;r<<=1){ x=y; fo(r)y=f(y); for(ll k=0;g==1&&k0); auto f=[](auto&f,ll m){ if(m==1)return vec{}; ll d=pollard_rho(m); return d==m?vec{d}:f(f,d)^f(f,m/d); }; return rce(f(f,n)); } single_testcase void solve(){ ST(sa,sb); ll A; if(sb[0]!='-'){ if(sa.substr(sa.size()-4,sa.size())!="0000")done(No()); A=stoll(sa.substr(0,sa.size()-4)); }else{ fo(i,sa.size()-5,sa.size()-1)sa[i]=sa[i+1]; A=stoll(sa.substr(0,sa.size()-1)); if(10000%A)done(No()); A=10000/A; } fo(i,sb.size()-5,sb.size()-1)sb[i]=sb[i+1]; ll B=stoll(sb.substr(0,sb.size()-1)); fe(factorize(A),p,q)if(B*q%10000)done(No()); pp(Yes()); }}