【数论】noip2009Hankson 的趣味题

【问题描述】
Hanks博士是BT（Bio-Tech，生物技术）领域的知名专家，他的儿子名叫Hankson。现在，刚刚放学回家的Hankson正在思考一个有趣的问题。 今天在课堂上，老师讲解了如何求两个正整数c1和c2的最大公约数和最小公倍数。现在Hankson认为自己已经熟练地掌握了这些知识，他开始思考一个“求公约数”和“求公倍数”之类问题的“逆问题”，这个问题是这样的：已知正整数a0,a1,b0,b1，设某未知正整数x满足：
1、 x和a0的最大公约数是a1；
2、 x和b0的最小公倍数是b1。
Hankson的“逆问题”就是求出满足条件的正整数x。但稍加思索之后，他发现这样的x并不唯一，甚至可能不存在。因此他转而开始考虑如何求解满足条件的x的个数。请你帮助他编程求解这个问题。

【问题分析】
根据题意可得

gcd(x,a0)=a1 <=> gcd(a0/a1,x/a1)=1
lcm(x,b0)=b1 <=> gcd(b1/x,b1/b0)=1

这是很好证明的，根据gcd是最大公约数，即找不出其他约数的情况下，所以a/gcd(a,b)与b/gcd(a,b)互质。那么lcm是俩个数的最小公倍数，怎么保证最小呢？就是对于俩个数重复的质因数只乘一遍就好，所以lcm(a,b)=a*b/gcd(a,b)。

所以根据列出的式子，枚举x为b1的约数，验证其他结论即可。需要注意的是gcd这里用的是辗转相除法，为什么对呢？因为gcd(a,b)是a的一个约数，同时也是b的约数，那么a%b一定也含有一个约数即gcd(a,b)。所以前后是等价的，这个代价是log级别的，为什么别问我，问欧拉。(他证明出来了，还有一个欧拉常数)。

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
int a0,a1,b0,b1,ans;
int gcd(int a,int b)
{   return b==0?a:gcd(b,a%b);   }
int work()
{
    if (a0%a1!=0||b1%b0!=0)
        return 0;
    int a=a0/a1,b=b1/b0;
    int i,x;  ans=0;
    int top=floor(sqrt(b1));
    for (i=1;i<=top;i++)
        if (b1%i==0)
        {
            if (i%a1==0)
            if (gcd(b,b1/i)==1)
                if (gcd(a,i/a1)==1)
                ans++;
            x=b1/i;
            if (x%a1==0)
            if (x!=i)
                if (gcd(b,b1/x)==1)
                    if (gcd(a,x/a1)==1)
                    ans++;
        }
    return ans;
}
int main()
{
    int t,i;
    scanf("%d",&t);
    for (i=1;i<=t;i++)
    {
        scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
        printf("%d\n",work());
    }
    return 0;
}