这道题我一开始的思路是枚举从a到b中的每一个数的组合,然后判断这个组合中的gcd和lcm是否分别等于a和b,然而这么做会超时
#include<iostream>
#include<algorithm>
using namespace std;
int a, b, ans;
int gcd(int x, int y) {
if (x < y)
swap(x, y);//让x一直是最大的那个数
return (y == 0) ? x : gcd(y, x % y);//算一下其最大公约数
}
int lcm(int x, int y) {
return x * y / a;
}
int main() {
cin >> a >> b;
for (int i = a; i <= b; i++)
for (int j = a; j <= b; j++)
if (gcd(i, j) == a && lcm(i, j) == b)
ans++;
cout << ans;
return 0;
}
于是我只好优化一下代码,我们其实不需要枚举每一个数的组合,对于不同的x,y交换它们的值一定可以等于另一个答案,因此我们只需要循环一次即可。
我们都知道,如果当x,y满足条件时,那么x*y=a*b;所以我们可以只进行一次循环,通过枚举x的值和这个式子求出另一个y,y=a*b/x;然后再通过相应判断来增加答案
所以我们可以让x枚举到sqrt(a*b),一开始时x<y,而当x=y时,x=sqrt(a*b),如果x的值再增大,那么x>y,我们的计算就会重复。
判断的时候我们进行a*b%i的判断其实就是在判断x*y=a*b。
值得注意的是,如果x,y的值相同,那么交换它们并不会得到另一个答案,所以当a和b相同时,ans要-1;
#include<iostream>
#include<algorithm>
#include<math.h>
using namespace std;
typedef long long ll;
ll a, b, ans;
int gcd(ll x, ll y) {
if (x < y)
swap(x, y);//让x一直是最大的那个数
return (y == 0) ? x : gcd(y, x % y);//算一下其最大公约数
}
int main() {
cin >> a >> b;
if (a == b) ans--;
for (ll i = a; i <= sqrt(a * b); i++)
if (a*b%i==0&&gcd(i, a * b / i) == a)
ans+=2;
cout << ans;
return 0;
}