D1. Domino (easy version)（分解质因数+欧拉筛）

D1. Domino (easy version)

题目入口

思路
设a=p1 ^ a1 * p2 ^ a2 * p3 ^ a3…pn ^ an (p1,p2…为a的质因数），b=k1 ^ b1 * k2 ^ b2 * k3 ^ b3…kn ^ bn (k1,k2…为b的质因数），则a可以在cnt1=a1+a2+a3…an的次数内变成1，b能够在 cnt2=b1+b2+b3…+bn的次数内变成1，于是可以发现，若k>cnt1+cnt2，两个数在cnt1+cnt2操作后无法操作，输出NO，在k<cnt1+cnt2时，特判k=1的情况即可。本题分解质因数可用欧拉筛预处理。

代码

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e7+5;
int prime[N];
bool vis[N];
void euler()
{
    vis[1]=true;
    int cnt=0;
    for(int i=2;i<=N-5;i++)
    {
        if(!vis[i]) prime[++cnt]=i;
        for(int j=1;i*prime[j]<=N-5;j++)
        {
            vis[i*prime[j]]=true;
            if(i%prime[j]==0) break;
        }
    }
}
signed main()
{
    int t;
    cin>>t;
    euler();
    while(t--)
    {
        int a,b,k;
        scanf("%lld %lld %lld",&a,&b,&k);
        int tmpa=a,tmpb=b;
        int cnt1=0,cnt2=0;
        if(a>N-5||vis[a])
        for(register int j=1; prime[j]*prime[j]<=a; j++)
        {
            if(a%prime[j]==0)
            {
                while(a%prime[j]==0)
                {
                    a/=prime[j];
                    cnt1++;
                }
            }
        }
        if(a>1) cnt1++;
        if(b>N-5||vis[b])
        for(register int j=1; prime[j]*prime[j]<=b; j++)
        {
            if(b%prime[j]==0)
            {
                while(b%prime[j]==0)
                {
                    b/=prime[j];
                    cnt2++;
                }
            }
        }
        if(b>1) cnt2++;
        if(cnt1+cnt2<k) printf("No\n");
        else if(k==1)
        {
            if((tmpa%tmpb==0||tmpb%tmpa==0)&&tmpa!=tmpb) printf("Yes\n");
            else if(tmpa==tmpb) printf("No\n");
            else printf("No\n");
        }
        else printf("Yes\n");
    }
    return 0;
}