wikioi 1012最大公约数和最小公倍数【根据最大公约数和最小公倍数求原来的两个数a、b】

/*======================================================================
题目描述
输入二个正整数x0,y0(2<=x0<100000,2<=y0<=1000000),求出满足下列条件的P,Q的个数
条件:  1.P,Q是正整数
2.要求P,Q以x0为最大公约数,以y0为最小公倍数.
试求:满足条件的所有可能的两个正整数的个数.
输入描述 Input Description
二个正整数x0,y0
输出描述 Output Description
满足条件的所有可能的两个正整数的个数
样例输入 Sample Input
3 60
样例输出 Sample Output
4

思路：
　　正常的思维是枚举每一个p，x<=p<=y。对每一个p再枚举每一个q，x<=q<=y。
假如p和q满足gcd(p，q)==x&&lcm(p，q)==y则认为得到了一组p、q，解得个数N加1。
如此往复，直到统计完所有的解的个数即可输出N。
但是这个方法，循环复杂度是10^5*10^5==10^10，这个已经远远超出了10^8这个勉强能接受的范围，所以肯定会超时。
优化的思路：
    枚举每一个p，假如y%p==0，则这个p可能是解，但要根据p、x、y算出q再验证现在这一组p、q是否是解。
其中q=y*x/p。然后：如果gcd（p，q）==x&&lcm(p,q)==y，那这一组p、q就是一组解了。
这样的解法，循环只有一重，复杂度是10^5，很容易接受。
========================================================================*/

 #include<stdio.h>
 int gcd(int a,int b);//最大公约数
 int lcm(int a,int b);//最小公倍数
 int main()
 {
     int x,y,p,q;
     int N=;
     scanf("%d%d",&x,&y);
     for(p=x;p<=y;p++)
     {
         if(y%p==)
         {
             q=y/p*x;
             if(gcd(p,q)==x&&lcm(p,q)==y) N++;
         }
     }
     printf("%d\n",N);
     return ;
 }
 int gcd(int a,int b)//最大公约数，输入要求a>=b>0
 {
     int c;
     if(b==) return -;
     c=a%b;
     while(c!=)
     {
         a=b;
         b=c;
         c=a%b;
     }
     return b;
 }
 int lcm(int a,int b)
 {
     if(a==||b==) return -;
     return a*b/gcd(a,b);
 }

x<=p<=y