洛谷 1072

      这道题我们需要自己来推算，首先，我们有gcd(x,b0)*b1=x*b0，将b0除过去，有gcd(x,b0)*(b1/b0)=x，令t=b1/b0，则有gcd(x,b0)*t=x，那么我们可以发现x一定是能表示成kt的形式，那么再整理一下有gcd(kt,b0)*t=kt，消去t，有gcd(kt,b0)=k，同时除以k，有gcd(t,b0/k)=1。再看令一个式子，gcd(x,a0)=a1，我们同时除以a1，有gcd(x/a1,a0/a1)=1，令w=a0/a1，则有gcd(x/a1,w)=1。又有x=kt，那么最后的式子是gcd(kt/a1,w)=1，这个时候我们发现条件的变量只剩下了k，那么我们可以考虑枚举这个k，这样的复杂度是O(b*logb)的，会TLE，再观察发现k一定是b0的因子，于是我们可以根号b的枚举k，O(logb)的判断，这样就可以AC了。

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int a0,a1,b0,b1,t,x,T,ans,w    ;

int gcd(int x,int y)
{
    if (y==0) return x;
    return gcd(y,x%y);
}

int main()
{
    scanf("%d",&T);
    while (T--) 
    {
        ans=0;
        scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
        t=b1/b0;w=a0/a1;
        for (int k=1;k*k<=b0;k++) 
        {
            if (b0%k) continue;
            //printf("%d\n",k*t);
            if (k*t%a1==0&&gcd(t,b0/k)==1&&gcd(k*t/a1,w)==1) ans++;
            if (b0/k==k) break;
            int k1=b0/k;
            if (k1*t%a1) continue;
            //printf("%d\n",k1*t);
            if (gcd(t,b0/k1)!=1) continue;
            if (gcd(k1*t/a1,w)!=1) continue;
            ans++;
        }
        printf("%d\n",ans);
    }
    return 0;    
}