ZOJ-1577 求素因子个数 + gcd和lcm的性质?(或者你也可以暴力?)

题意：

给出两个数x,y，求有多少组p,q，满足gcd(p,q) = x 且 lcm(p,q) = y。

注意：pair(p, q) 和 pair(q, p)是两种方案。

思路：

两个很易得的结论：
lcm（x，y）/ gcd(x,y)  =  (x / gcd(x,y)) * (y / gcd(x,y))
lcm（x，y）* gcd(x,y)  =  x*y

假设g=gcd(p,q)，l=lcm(p,q)，key=l/g，如果l%g!=0则无解。

l是p和q公共的素因子的并集，而g则是交集(图中红色部分)，于是key的意思就是l刚好比g多的一部分素因子(图中蓝色部分)

于是我们可以选择这些素因子和g组合得到了p，另一部分和g组合得到了q。

假设key的素因子个数（重复的不算）有n个，那么选择方案就有：C(n,0)+C(n,1)+C(n,2)+...+C(n,n)，这个式子的和为2^n。C(n, 0)和C(n, n)就对应着上面要求注意的一对pair(p, q)和pair(q, p)。

重复的素因子不算个数，因为同1个素因子必须全部在p中或者全部在q中，如果同时属于p和q，GCD将增大。

以上摘自数论学习小计其之三 Gcd与Lcm - 时雨晴天

下面是代码：

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const int maxn = 1500;
int isprime[maxn], prime[maxn], np;
void init(){
	np = 0;
	for(int i = 2; i <= maxn; ++i)
	if(!isprime[i]){
		prime[np++] = i;
		for(int j = i*i; j <= maxn; j+=i) isprime[j] = 1;
	}
}
int calcu(int n){
	int ans = 0;
	for(int i = 0; prime[i]*prime[i] <= n && i < np; ++i){
		if(n%prime[i] == 0){
			++ans;
			while(n%prime[i] == 0) n /= prime[i];
		}
	}
	if(n > 1) ++ans;
	return ans;
}
int main(){
	int x, y, ans;
	init();
	while(cin >> x >> y){
		ans = 0;
		if(y%x == 0){
			int n = calcu(y/x);
			ans += pow(2, n); 
		}
		cout << ans << endl;
	}
	return 0;
}

不过还有一个暴力的做法：

就是求出key=l/g之后，然后枚举得两个互素且相乘为key的个数，再加一个特判就是当x == y时，ans为1即可。

#include <iostream>
#include <cstdio>
using namespace std;
int gcd(int a, int b){
	if(a%b == 0) return b;
	return gcd(b, a%b);
}
int main(){
	int x, y, ans;
	while(cin >> x >> y){
		ans = 0;
		if(y%x == 0){
			int n = y/x;
			for(int i = 1; i <= n; ++i){
				if(n/i <= i) break;	//之后就是重复枚举了 
				if(n%i == 0 && gcd(i, n/i) == 1) ans += 2;
				//每次枚举满足的都会增加两个ans，即p,q和q,p 
			}
		}
		if(x == y) ans = 1;
		cout << ans << endl;
	}
	return 0;
}

继续加油~