HDU5812 Distance 构造，预处理

分析：怎么看都是超时，但是可以先筛一遍1e6以内的每个数的最小素数

算出每个数由多少个素数组成，然后应用，c[1e6][20]

就是题解的那一套，参照题解，比赛的时候没有想到好的办法筛一个数的因子，醉了

然后赛后发现，预处理因子肯定超时，虽然是O(nlogn)的，但是n是1e6啊，常数太大

而且单组操作只有5e4,所以暴力sqrt(x)即可

#include <iostream>
#include <vector>
#include <queue>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <algorithm>
using namespace std;
const int N = 1e6+;
int c[N][],tot,prime[],mn[N],cnt[N],q,kase;
bool check[N];
char op[];
void getprime()
{
    for(int i=; i<=N-; ++i)
    {
        if(!check[i])mn[i]=prime[++tot]=i;
        for(int j=; j<=tot; ++j)
        {
            if(i*prime[j]>N-)break;
            check[i*prime[j]]=true;
            mn[i*prime[j]]=prime[j];
            if(i%prime[j]==)break;
        }
    }
}
void getcnt()
{
    for(int i=; i<=N-; ++i)
    {
        int tmp=i;
        while(tmp!=)++cnt[i],tmp/=mn[tmp];
    }
}
int main()
{
    getprime();
    getcnt();
    memset(mn,,sizeof(mn));
    while(~scanf("%d",&q),q)
    {
        printf("Case #%d:\n",++kase);
        memset(c,,sizeof(c));
        int ttt=;
        for(int i=; i<q; ++i)
        {
            int x;
            scanf("%s%d",op,&x);
            if(op[]=='I')
            {
                if(mn[x]==kase)continue;
                mn[x]=kase;++ttt;
                for(int j=; j*j<=x; ++j)
                {
                    if(x%j)continue;
                    ++c[j][cnt[x/j]];
                    if(x/j!=j)++c[x/j][cnt[j]];
                }
            }
            else if(op[]=='D')
            {
                if(mn[x]!=kase)continue;
                mn[x]=;--ttt;
                for(int j=; j*j<=x; ++j)
                {
                    if(x%j)continue;
                    --c[j][cnt[x/j]];
                    if(x/j!=j)--c[x/j][cnt[j]];
                }
            }
            else
            {
                if(ttt==){printf("-1\n");continue;}
                int ans=;
                for(int j=; j*j<=x; ++j)
                {
                    if(x%j)continue;
                    for(int k=; k<=; ++k)
                    {
                        if(c[j][k])
                        {
                            ans=min(ans,k+cnt[x/j]);
                            break;
                        }
                    }
                    if(x/j!=j)
                    {
                        for(int k=; k<=; ++k)
                        {
                            if(c[x/j][k])
                            {
                                ans=min(ans,k+cnt[j]);
                                break;
                            }
                        }
                    }
                }
                if(ans==)ans=-;
                printf("%d\n",ans);
            }
        }
    }
    return ;
}