【HDOJ5447】Good Numbers（数论）

题意：

思路：From https://blog.csdn.net/qq_36553623/article/details/76683438

大概就是把1e6里面的质因子能除的都除光之后借助两者gcd中有最大质因子的条件分类讨论

其实也就5类，重点是如何设分类的标准，其实唯一的标准就是对答案的贡献

int128本机环境跑不了，只能交上去的时候改，总觉得有朝一日要掉坑里

 #include<bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 typedef unsigned int uint;
 typedef unsigned long long ull;
 typedef long double ld;
 typedef pair<int,int> PII;
 typedef pair<ll,ll> Pll;
 typedef vector<int> VI;
 typedef vector<PII> VII;
 typedef pair<ll,ll>P;
 #define N  1000000+10
 #define M  200000+10
 #define INF 1e9
 #define fi first
 #define se second
 #define MP make_pair
 #define pb push_back
 #define pi acos(-1)
 #define mem(a,b) memset(a,b,sizeof(a))
 #define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
 #define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
 #define lowbit(x) x&(-x)
 #define Rand (rand()*(1<<16)+rand())
 #define id(x) ((x)<=B?(x):m-n/(x)+1)
 #define ls p<<1
 #define rs p<<1|1
 #define fors(i) for(auto i:e[x]) if(i!=p)

 const int MOD=,inv2=(MOD+)/;
       //int p=1e4+7;
       double eps=1e-;
       int dx[]={-,,,};
       int dy[]={,,-,};

 int p[N],b[N];

 int read()
 {
    int v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 ll readll()
 {
    ll v=,f=;
    char c=getchar();
    while(c<||<c) {if(c=='-') f=-; c=getchar();}
    while(<=c&&c<=) v=(v<<)+v+v+c-,c=getchar();
    return v*f;
 }

 template<class T> void Read(T &x)
 {
     static int CH;
     while((CH=getchar())<''||CH>'');
     for(x=CH-'';(CH=getchar())>=''&&CH<='';x=(x<<)+(x<<)+(CH-''));
 }

 int is2(__int128 x)
 {
     __int128 t=pow(x,1.0/)+0.5;
     if(t*t==x) return ;
     return ;
 }

 int is3(__int128 x)
 {
     __int128 t=pow(x,1.0/)+0.5;
     if(t*t*t==x) return ;
     return ;
 }

 void calc(ll &ans,__int128 x,__int128 y)
 {
     if(y==)
     {
         if(x==) ans*=;
          else if(is3(x)) ans*=;
           else if(is2(x)) ans*=;
            else ans*=;
     }
      else if(y>&&is2(y))
      {
          if(y%x==) ans*=;
           else ans*=;
      }
       else if(x>&&is2(x))
       {
           if(x%y==) ans*=;
            else ans*=;
       }
        else if(y%x==) ans*=;
         else ans*=;
 }

 int main()
 {
     //freopen("1.in","r",stdin);
     //freopen("1.out","w",stdout);
     int n=;
     rep(i,,1e6)
     {
         if(!b[i]) p[++n]=i;
         rep(j,,n)
         {
             int t=p[j]*i;
             if(t>1e6) break;
             b[t]=;
             if(i%p[j]==) break;
         }
     }
     int cas=read();
     while(cas--)
     {
         __int128 s1,s2;
         Read(s1);
         Read(s2);
         ll ans1,ans2;
         ans1=ans2=;
         rep(i,,n)
         {
             if(s1%p[i]==)
             {
                 int s=;
                 while(s1%p[i]==)
                 {
                     s++;
                     s1/=p[i];
                 }
                 ans1*=s;
             }
             if(s2%p[i]==)
             {
                 int s=;
                 while(s2%p[i]==)
                 {
                     s++;
                     s2/=p[i];
                 }
                 ans2*=s;
             }
         }
         __int128 k=std::__gcd(s1,s2);
         calc(ans1,k,s1/k);
         calc(ans2,k,s2/k);
         printf("%I64d %I64d\n",ans1,ans2);
     }
     return ;
 }