【C++心路历程12】1133 因数分解，公因数公倍数综合题

#include<cstdio>
#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
vector<int>fac;
void getfac(int a,vector<int>&fac)     //取因数
{
 int len=(int)sqrt(a+0.5);
 for(int i=1;i<=len;i++)
 {
  if(a%i==0)
  {
   fac.push_back(i);
   if(a/i>i) fac.push_back(a/i);
  }
 } 
}
int gcd(int x,int y)          //最大公因数
{
 if(y==0) return x;
 else return gcd(y,x%y);
}
int lcm(int x,int y)       //最小公倍数
{
 return x/gcd(x,y)*y;
}
int main()
{
 int a0,a1,b0,b1,n,j,cnt;
 scanf("%d",&n);
 for(int i=1;i<=n;i++)
 {
  scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
  getfac(b1,fac);      //把b1因数存入数组中
  cnt=0;
  int sz=fac.size();
  for(j=0;j<sz;j++)
  {
   int g=gcd(a0,fac[j]);
   int d=lcm(b0,fac[j]);
   if(g==a1 && d==b1)
   cnt++;
  }
  printf("%d\n",cnt);
  fac.clear();     //非常重要
 }
 return 0;
}