hdu—4497

GCD and LCM

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65535/65535 K (Java/Others)
Total Submission(s): 2624    Accepted Submission(s): 1150

Problem Description

Given two positive integers G and L, could you tell me how many solutions of (x, y, z) there are, satisfying that gcd(x, y, z) = G and lcm(x, y, z) = L? 
Note, gcd(x, y, z) means the greatest common divisor of x, y and z, while lcm(x, y, z) means the least common multiple of x, y and z. 
Note 2, (1, 2, 3) and (1, 3, 2) are two different solutions.

 

Input

First line comes an integer T (T <= 12), telling the number of test cases. 
The next T lines, each contains two positive 32-bit signed integers, G and L. 
It’s guaranteed that each answer will fit in a 32-bit signed integer.

 

Output

For each test case, print one line with the number of solutions satisfying the conditions above.

 

Sample Input

  
  

   
   2 
6 72 
7 33 
  
  

 

Sample Output

  
  

   
   72 
0
  
  

 

Source

2013 ACM-ICPC吉林通化全国邀请赛——题目重现

 

Recommend

liuyiding   |   We have carefully selected several similar problems for you:   6010  6009  6008  6007  6006

【分析】

题意很简单，让你求任意x,y,z满足他们的gcd是l，lcm是g的方案数。

首先肯定可以直接发现一个问题那就是，如果g%i！=0那么方案数就是0.

因为三个数的最大公约数一定是最小公倍数的因子，这个很容易证明。

那么麻烦的就是怎么求方案数了，首先可以知道，模拟xyz的组成一定是超时的因为数据范围到2^32-1.

那么就需要数学方法去算了，怎么算呢

首先我们可以考虑

三个数的gcd是l，需要满足什么？

显然，我们先把x,y,z分别质因数分解

可以得到

x=q1^p1*q2^p2*q3^p3....

y=q1^pp1*q2^pp2*q3^pp3....

z=q1^ppp1*q2^ppp2*q3^ppp3....

那么，x,y,z的gcd可以表示为

q1^min(p1,pp1,ppp1)*q2^min(p2,pp2,ppp2).....

lcm就是q1^max(p1,pp1,ppp1)*q2^max(p2,pp2,ppp2).....

那么显然我们可以发现的是，我们先假设x=y=z=l,那么进行怎样的变化可以让x,y,z的lcm变成g呢？显然就是对于每个质因数，只要三个数中有一个数的质因数的个数满足max(p,pp,ppp)就可以了，当然不能有数的质因数大于这个max

那么可以发现，我们要让x=y=z=l变化到lcm是g，只要找到g/l这个数每个质因数的个数，然后做一个组合就可以了..

比较考思路的一道题....不太适合我这种数学菜鸡0.0

【代码】

#include <stdio.h>
int prime[1000000];
int f[10000000]={0};
int len=0;
void init()
{
	for (int i=2;i<10000000;i++)
		if (!f[i])
		{
			prime[len++]=i;
			for (int j=i+i;j<10000000;j+=i) 
				f[j]=1;
		}
}
int main()
{
	init();
//	for (int i=0;i<10;i++) printf("%d ",prime[i]);
	int pp;scanf("%d",&pp);
	int l,g;
	while (pp--)
	{
		scanf("%d%d",&l,&g);
		if (g%l!=0)
		{
			printf("0\n");
			continue;
		}
		g/=l;
		int ans=1;
		for (int i=0;prime[i]<=g;i++)
			if (g%prime[i]==0)
			{
				int t=0;
				while (g%prime[i]==0)
				{
					g/=prime[i];
					t++;
				}
				ans=ans*t*6;
			}
		printf("%d\n",ans);
	}
	return 0;
}