最小公倍数 (Benefit,UVa 11889)

#include <iostream>
#include <string.h>
#include <string>
#include <fstream>
#include <algorithm>
#include <stdio.h>
#include <vector>
#include <queue>
#include <set>
#include <cmath>
using namespace std;
const double eps = 1e-8;
const int INF=0x7fffffff;
#define MAXN 10000002
typedef long long LL;
int vis[MAXN];
int prime[700000];

void sieve(int n)
{
    int m=(int)sqrt(n+0.5);
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=m;i++)if(!vis[i])
    for(int j=i*i;j<=n;j+=i)vis[j]=1;
}

int gen_prime(int n)
{
    sieve(n);
    int c=0;
    for(int i=2;i<=n;i++)if(!vis[i])
        prime[c++]=i;
    return c;
}

LL gcd(LL a,LL b)
{
    return b==0?a:gcd(b,a%b);
}

int main()
{
    int T;scanf("%d",&T);
    LL a,b,c;
    while(T--)
    {
        cin>>a>>c;
        if(c%a!=0)printf("NO SOLUTION
");
        else
        {
        b=c/a;
        LL temp=b;
        while(temp!=1)
        {
            temp=gcd(b,a);
            b*=temp;
            a/=temp;
        }
        cout<<b<<endl;
        }
    }
    return 0;
}

c/a=b/gcd(a,b) 找出a与b/gcd(a,b)所用相同的质因子 乘上b/gcd(a,b)即是最小的b