F. Vasilije Loves Number Theory

题解

$a, n$ 互质，所以 $d (n \cdot a) = d (a) d (n)$ ，即 $n m o d d (n) == 0$ 是否成立。（总是能构造出一与 $n$ 互质，且 $d (a)$ 任意的 $a$ )
由于 $n$ 会很大，所以我们将 $n$ 质因子分解， $n = p_{1}^{m_{1}} p_{2}^{m_{2}} . . . p_{k}^{m_{k}}$
这样一来 $d (n) = (m_{1} + 1) (m_{2} + 1) . . . (m_{k} + 1)$
乘 $x$ 的时候，将 $x$ 也进行质因子分解
最后在判断 $d (n)$ 进行质因子分解，然后判断因子是否足够

code

#include<bits/stdc++.h>
#define ll long long
using namespace std;

bool mark[1000006]={0};
int prime[1000006];


void solve()
{
    int n,q;
    cin>>n>>q;
    int tem=n;
    map<int,int> yz;

    while(n!=1)
    {
        yz[prime[n]]++;
        n/=prime[n];
    }

    while(q--)
    {
        int op;
        cin>>op;
        if(op==2)
        {
            yz.clear();
            n=tem;
            while(n!=1)
            {
                yz[prime[n]]++;
                n/=prime[n];
            }
        }
        else
        {
            int x;
            cin>>x;
            while(x!=1)
            {
                yz[prime[x]]++;
                x/=prime[x];
            }

            ll ans=1;
            for(auto it:yz)
            {
                ans*=(it.second+1);
            }

            bool flag=1;
            for(ll i=2;i*i<=ans;i++)
            {
                if(ans%i==0)
                {
                    int cnt=0;
                    while(ans%i==0)
                    {
                        ans/=i;
                        cnt++;
                    }
                    if(yz[i]<cnt)
                    {
                        flag=0;
                        break;
                    }
                }
            }
            if(ans>1)
            {
                if(yz[ans]<1) flag=0;
            }
            if(flag) cout<<"yes\n";
            else cout<<"no\n";
        }
    }
    cout<<'\n';
}

int main()
{
    ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    for(ll i=2;i<=1e6;i++)
    {
        if(!mark[i])
        {
            prime[i]=i;
            for(ll j=i*i;j<=1e6;j+=i)
            {
                if(!mark[i])
                {
                    prime[j]=i;
                    mark[j]=1;
                }
            }
        }
    }
    int t;
    cin>>t;

    while(t--) solve();
    return 0;
}