【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);
}
}