[NOIP2009] $Hankson$ 的趣味题 (数论,gcd)

题目链接

Solution

此题,用到的结论都是比较浅显的,但是,我竟然没想到反过来枚举...
只有50分... 被自己蠢哭...

结论比较浅显:

1.对于两个正整数(a),(b),设 (gcd(a,b)=k),则存在(gcd(a/k,b/k)=1).

也就是说 (x=k_1*a_1),(a_0=k_2*a_1),它们最大公约数为(a_1),那么要求 (k_1) 与 (k_2) 必须互质,否则它们的最大公约数会是 (gcd(k_1,k_2)*a_1).

2.对于两个正整数(a),(b),设(lcm(a,b)=k),则存在(gcd(k/a,k/b)=1).

比较浅显,可以由 (a*b=gcd(a,b)*lcm(a,b)) 推出来.


然后通过分析题意结论,便可以分析出 (x) 满足 (x) 是 (b_1) 的因子,并且满足是 (a_1) 的倍数.
所以我们直接 (sqrt{b_1}) 枚举其因子,并且判断是否满足上述条件即可.

### Code ### 100 分做法 ```cpp #include #define ll long long using namespace std; ll n,a1,a0,b0,b1;

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

int main()
{
scanf("%lld",&n);
while(n--)
{
scanf("%lld%lld%lld%lld",&a0,&a1,&b0,&b1);
if(b1%a1!=0){printf("0 ");continue;}
ll ans=0,maxx=sqrt(b1);
for(int x=1;x<=maxx;x++)
{
if(b1%x!=0)continue;
if(x%a10)
if(gcd(b1/b0,b1/x)1)
if(gcd(x/a1,a0/a1)1)
ans++;
if(b1/xx)continue;
ll y=b1/x;
if(y%a10)
if(gcd(b1/b0,b1/y)1)
if(gcd(y/a1,a0/a1)==1)
ans++;
}
printf("%lld ",ans);
}
}

### 50 分做法(暴力枚举 $a_1$ 的倍数,然后判断)
```cpp
#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll n,a1,a0,b0,b1;

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

int main()
{
    scanf("%lld",&n);
    while(n--)
    {
        scanf("%lld%lld%lld%lld",&a0,&a1,&b0,&b1);
        if(b1%a1!=0){printf("0
");continue;}
        ll tt=0,ans=0;
        while(1)
        {
            tt++;
            if(tt*a1>b1)break;
            ll x=tt*a1;
            if(b1%x!=0)continue;
            if(gcd(x,a0)!=a1)continue;
            if(x*b0!=gcd(b0,x)*b1)continue;
            ans++;
        }
        printf("%lld
",ans);
    }
}