codeforces833 A题解

话说看各种提交的时候发现还有求导的……简直可怕

---

题的大意是这样的：假设有两个数，使得：

a=X21×X22×…×X2n×Y1×Y2×…×Yn

b=Y21×Y22×…×Y2n×X1×X2×…×Xn

则称a,b具有关系α。
现在给出q组询问，每组询问涉及两个数，要求判断两个数之间是否具有关系α。
(1 ≤ n ≤ 350000，1 ≤ a, b ≤ 109)
看见这个数据范围，大概是一道数论题。
我首先想到的是gcd，如果我忽略了一个问题：求gcd(a,b)不一定能求出
X1×X2×…×Xn×Y1×Y2×…×Yn
因为可能其中有些因子是重合的。看了某位神犇的题解之后才恍然大悟：其实完全不需要gcd。
如果该条件满足的话 a×b所求出的值一定是上面式子的三次方。两数之积的立方根为整数就是a,b具有关系α的充分条件。
再来看看如何满足必要性。如果一个数是一个立方，那么这个数必然能拆解为n个立方数的积。我们考虑其中一个，如k3=q是一个因子，那么有四种情况：
1.q3是a的因子
2.q2是a的因子，q是b的因子
3.q是a的因子，q2是b的因子
4.q3是b的因子
显然中间的两者是我们所需要的。那么我们就应该设置一个限定条件筛去1和4.怎么筛呢？非常简单：我们求出的ab−−√3中一定含有因子q，换言之，如果：

a,b≡0(modab−−√3)

就说明q在a,b之中。所以，我们得到的结论是：

a,b∃关系α⇔ab−−√3为整数且a,b≡0(modab−−√3)

实现：（注意读入优化，不然会T）

#include <bits/stdc++.h>using namespace std;#define eps 0.000001typedef long long ll;void read(ll &x){    x=0;    char c=getchar();    while(!isdigit(c))  c=getchar();    while(isdigit(c)){        x=x*10+c-'0';        c=getchar();    }}bool havealpha(ll a,ll b){    ll t=round(pow((long double)a*b,1.0/3));    cout<<(t*t*t-a*b);    if(fabs(t*t*t-a*b)<=eps&&a%t==0&&b%t==0)        return 1;    else return 0;    }int main(){    ll q;    read(q);    ll m,n;    for(int i=0;i<q;++i){        read(m),read(n);        if(havealpha(m,n))              cout<<"Yes"<<"\n";        else                cout<<"No"<<"\n";    }}