HDU 4320 Arcane Numbers 1 (质因子分解)

#include <iostream>
#include <math.h>
#include <string.h>
#include <algorithm>
#include <stdio.h>
#include <stdlib.h>
const int maxn=1000100;
bool isprime[maxn];
int prime[maxn];
using namespace std;
int i,j,k=0;
void pr()
{
    for(i=2;i<=maxn;i++)
    {
        if(!isprime[i])
        {
            for(j=i+i;j<=maxn;j=j+i)
            isprime[j]=true;
            prime[k++]=i;
        }
    }
}
int main()
{
    pr();
    int t,c=1;
    bool flag=true;
    long long a,b;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lld%lld",&a,&b);
        flag=true;
        for(i=0;i<k&&prime[i]<=a;i++)
        {
            if(a%prime[i]==0)
            {
                if(b%prime[i]!=0)
                {
                    flag=false;
                    break;
                }
            }
            while(a%prime[i]==0)
            a/=prime[i];
        }
        if(b%a!=0)
        flag=false;
        if(flag)
        cout<<"Case #"<<c++<<": YES"<<endl;
        else
        cout<<"Case #"<<c++<<": NO"<<endl;
    }
    return 0;
}