P1072 Hankson 的趣味题

emmm这道题看到题目首先要打一下草稿，首先暴力的算法可能是会T掉了，最开始我就是用的暴力的算法，T掉了5个点，还是得了50分（如果想不到好办法暴力就暴力吧）
附上暴力代码~

#include <iostream>
using namespace std;
int gcd(int a,int b);//a>=b
int lcm(int a,int b); //a<=b
    int main()
    {
        int a0,a1,b0,b1;
        int n;
        cin>>n;
        int ans=0;
        int res1,res2;
        for(int i=1;i<=n;i++)
        {
            cin>>a0>>a1>>b0>>b1;
            ans=0;
            for(int i=1;i<=b1;i++)
            {
                if(i%a1!=0||b1%i!=0)
                {
                    continue;
                }

                if(a0<i)
                {
                    res1=gcd(i,a0);
                }
                else
                {
                    res1=gcd(a0,i);
                }
                if(b0<i)
                {
                    res2=lcm(b0,i);
                }
                else
                {
                    res2=lcm(i,b0);
                }
                if(res1==a1&&res2==b1)
                {
                    ans++;
                }
            }
            cout<<ans<<endl;
        }

        return 0;
    }

    int gcd(int a,int b)//a>=b;
    {
        if(a%b==0)
        {
            return b;
        }
        else
        {
            return gcd(b,a%b);
        }
    }

    int lcm(int a,int b) //a<=b
    {
        int i=1;
        while(a*i%b!=0)
        {
            i++;
        }
        return a*i;
    }

这道题的正解是要推理的：根据题意我们可以知道gcd(x,a0)=a1 lcm(b0,x)=b1

变一下形:
gcd(x/a1,a0/a1)=1
gcd(b0,x)=b0*x/b1 ---->gcd(b1/x,b1/b0)=1

上面的推理部分用到一个知识点：某两个数的最小公倍数乘以最大公约数等于那两个数的乘积

根据上面的两个式子
gcd(x/a1,a0/a1)=1
gcd(b1/x,b1/b0)=1
我们可以得出结论 a1是x的因子，x是 b1的因子，所以我们只需要枚举1—sqrt(b1)这些数字，康康哪些数既是a1的倍数，同时也是b1的因子.

咳咳一定要注意对称性：假如枚举的那个数字是x 那么还要注意对应的对称的一个因子b/x康康是否满足条件，还要注意，如果x==b/x那么就不用判断b/x，因为这样就重复了~

下面是ACcode~

#include <iostream>
#include <stdio.h>
#include <cmath>
using namespace std;
int gcd(int a,int b);//a>=b

    int main()
    {

        int a0,a1,b0,b1;
        int n;
        scanf("%d",&n);
        int ans;
        for(int i=1;i<=n;i++)
        {
            ans=0;
            scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
            for(int x=1;x<=sqrt(b1);x++)
            {
                if(b1%x==0)
                {
                    if(x%a1==0&&gcd(x/a1,a0/a1)==1&&gcd(b1/b0,b1/x)==1)
                    {
                        ans++;
                    }
                    int x0=b1/x;
                    if(x0%a1==0&&x0!=x)
                    {
                        if(gcd(x0/a1,a0/a1)==1&&gcd(b1/b0,b1/x0)==1)
                        {
                            ans++;
                        }
                    }
                }

            }
            cout<<ans<<endl;
        }
        return 0;
    }

    int gcd(int a,int b)
    {
       if(a<b)
       {
           int t=a;
           a=b;
           b=t;
       }
       while(a%b!=0)
       {
           int r=a%b;
           a=b;
           b=r;
       }
       return b;
    }