最大公约数和最小公倍数问题

#include <iostream>
#include <cstdio>
#include <cmath>
#define LL long long
using namespace std;
int read()
{
	int x=0,f=1;
	char c=getchar();
	while(c<'0'||c>'9')
	{
		if(c=='-')  f=-1;
		c=getchar();
	}
	while(c>='0'&&c<='9')
	{
		x=x*10+c-'0';
		c=getchar();
	}
	return f*x;
}
int gcd(LL a,LL b)
{
    if(a<b)  a^=b,b^=a,a^=b;
	if(!b)  return a;
	return gcd(b,a%b);
}
LL a,b,n,ans;
int main()
{
	a=read();  b=read();
	n=sqrt(a*b);
	for(int i=1;i<=n;i++)
	  if(a*b%i==0&&gcd(i,a*b/i)==a)  ans++;
	cout<<ans*2;
	return 0;
}