【原创】【NOIP2009】Hankson的趣味题

【NOIP2009】Hankson

题目描述

题目描述

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的个数。请你帮助他编程求解这个问题。

输入

输入文件名为son.in。第一行为一个正整数n，表示有n组输入数据。接下来的n行每行一组输入数据，为四个正整数a0，a1，b0，b1，每两个整数之间用一个空格隔开。输入数据保证a0能被a1整除，b1能被b0整除。

输出

输出文件名为son.out。输出共n行。每组输入数据的输出结果占一行，为一个整数。 对于每组数据：若不存在这样的x，请输出0； 若存在这样的x，请输出满足条件的x的个数；

样例输入

2
41 1 96 288
95 1 37 1776

样例输出

6
2

题目分析

这道题目的大概意思就是求满足gcd(x,a0)=a1，lcm(x,b0)=b1，的数x有几个。

这道题的正确思路是分解质因数：

设a,b,c,d,v分别是a0,a1,b0,b1,x含有质因数p的数量。
如果a小于b，即a0的因数a1居然含有比a0还多的p！这是不可能的！不用算了，肯定无解了。
如果a=b，那么a0和a1含有同样多的p，那么x只需要含有>=b个p就可以，满足条件了。
如果a>b,即a0含有的p比a1多，怎么能让它小下来，小到b呢？那么v必须为b。

类似论述可以得到：
如果c小于d，v=d;
如果c=d，v<=d
如果c>d，无解。
//好像遇到了bug，如果上面两个‘小于’换成‘<’，那么小于之后的内容都会消失不见！

我们就可以规定出v的取值范围，也就是说，x可以取这么多种的p，所有用乘法原理乘起来就是了。

此题还有暴力解法，枚举1~sqrt(b1)，(因为b1总是最大的←此言论没有经过验证)，看这些数满不满足条件。但是注意，计算lcm的时候（如果你用的是乘积除以gcd），小心爆int！

详见代码

正解：

#include<cstdio>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;

void Read(int &p)
{
    p=0;
    char c=getchar();
    while(c<'0' || c>'9') c=getchar();
    while(c>='0' && c<='9')
        p=p*10+c-'0',c=getchar();
}

int pr[5517],alls;
bool vis[50017];
int a0,a1,b0,b1,_;

void search_p()
{
    int n=50000;
    for(int i=2;i<=n;i++)
    {
        if(!vis[i]) pr[++alls]=i;
        for(int j=1;j<=alls && i*pr[j]<=n;j++)
        {
            vis[i*pr[j]]=1;
            if(i%pr[j]==0) break;
        }
    }
}

void cd(int p,int &x,int &cnt)
{
    while(x%p==0) x/=p,cnt++;
}

int ans;

int dow(int p)
{
    int a=0,b=0,c=0,d=0;
    cd(p,a0,a); cd(p,a1,b); cd(p,b0,c); cd(p,b1,d);

    if(a<b || c>d) return 0;
    if(a==b)
    {
        if(c==d) return d-b+1>0?d-b+1:0;
        return d>=b?1:0;
    }
    //if(a>b) 
    if(c==d) return b<=d?1:0;
    return b==d?1:0;
}

bool check()
{
    return a0!=1 || a1!=1 || b0!=1 || b1!=1;
}

int main()
{
    search_p();
    Read(_);
    while(_--)
    {
        Read(a0); Read(a1); Read(b0); Read(b1);

        ans=1;
        for(int i=1;i<=alls && ans && check();i++)
            ans*=dow(pr[i]);
        if(a0>1) ans*=dow(a0);
        if(b1>1 && b1!=a0) ans*=dow(b1);
        cout<<ans<<endl;
    }   
    return 0;
}

暴解(小心超时！此代码耗时超高！但勉勉强强能过题)：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>

void Read(int &p)
{
    p=0;
    char c=getchar();
    while(c<'0' || c>'9') c=getchar();
    while(c>='0' && c<='9')
        p=p*10+c-'0',c=getchar();
} 

int n,a0,a1,b0,b1,ans;

int gcd(int a,int b)
{
    return b==0?a:gcd(b,a % b);
}

long long lcm(int a,int b)
{
    return (long long)a*b/gcd(a,b);
}

void solve(int x)
{
    if(gcd(x,a0)==a1 && lcm(x,b0)==b1)
        ans++;
}

int main()
{
   Read(n);
   while(n--)
   {
       ans=0;
       Read(a0); Read(a1); Read(b0); Read(b1);
       int sss=sqrt(b1);
       for(int i=1;i<=sss;i++)
       if(b1%i==0)
       {
           solve(i);
           if (i*i!=b1) solve(b1/i);
       }
       printf("%d\n",ans);
   }    
 }