P1072 Hankson 的趣味题

Posted 2022-01-06 xcg123

原题链接  https://www.luogu.org/problemnew/show/P1072

这个题是数论+暴力吧，需要一些优化。

先来看看怎么求最大公约数和最小公倍数：

欧几里德算法 ------求最大公约数

概述

欧几里德算法又称辗转相除法，用于计算两个整数a,b的最大公约数。其计算原理依赖于下面的定理：

gcd函数就是用来求(a,b)的最大公约数的。

gcd函数的基本性质：

gcd(a,b)=gcd(b,a)=gcd(-a,b)=gcd(|a|,|b|)

公式表述

gcd(a,b)=gcd(b,a mod b)    

gcd(a,b)=gcd(b,a%b)    ----c++语言

证明：a可以表示成a = kb + r，则r = a mod b

假设d是a,b的一个公约数，则有d|a, d|b，而r = a - kb，因此d|r；

因此d是(b,a mod b)的公约数

假设d 是(b,a mod b)的公约数，则

d | b , d |r ，但是a = kb +r

因此d也是(a,b)的公约数

因此(a,b)和(b,a mod b)的公约数是一样的，其最大公约数也必然相等，得证。

 

然后我们还有一个定理来求最小公倍数：

两个数的乘积等于其最大公约数与最小公倍数的乘积

证明：

这两数记为A，B，设gcd(A,B)=d, 那么存在整数a=A/d，b=B/d， gcd(a,b)=1

则A=ad，B=bd；

lcm(A,B)=lcm(ad,bd) ===将d提出来===>d*lcm(a,b)

因为gcd(a,b)=1（a和b互质），所以lcm(a*b)=ab，所以lcm(A,B)=dab；
gcd(A,B)*lcm(A,B)=d*dab=da*db=AB，得证。

 

这样我们就得到一个求最小公倍数的方法：lcm(a,b)=a*b/gcd(a,b)；

我们就觉得暴力枚举，但也要有个上界啊，不能没完没了的枚举QwQ，因为x和b0的最小公倍数是b1，也就是说x是b1的因子，那么x一定小于等于b1，那么我们就找到了最大上界：

 

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int read()

    char ch=getchar();
    int a=0,x=1;
    while(ch<‘0‘||ch>‘9‘)
    
        if(ch==‘-‘) x=-x;
        ch=getchar();
    
    while(ch>=‘0‘&&ch<=‘9‘)
    
        a=(a<<3)+(a<<1)+(ch-‘0‘);
        ch=getchar();
    
    return x*a;

int t;
int a0,a1,b0,b1,ans;
int gcd(int a,int b)

    if(b==0) return a;
    else return gcd(b,a%b);

int main()

    t=read();
    for(int i=1;i<=t;i++)
    
        ans=0;
        a0=read();     //x和a0的最大公约数是a1 
        a1=read();     //gcd(x,a0)=a1
        b0=read();     //x和b0的最大公约数是b1 
        b1=read();     //x*b0=b1*gcd(x,b0)  =>  x=b*gcd(x,b0)
        int b=b1/b0;
        for(int i=1;i<=sqrt(b1);i++)        //枚举到sqrt(b1)就行了 
            
                if(b1%i==0)
                
                    if(gcd(i,a0)==a1&&b*gcd(i,b0)==i)  ans++;
                    int j=b1/i;
                       if(i!=j)
                     
                           if(gcd(j,a0)==a1&&b*gcd(j,b0)==j)  
                             ans++; 
                       
                    
                
            
        printf("%d\\n",ans);
    
    return 0;