Hdu 5812 Distance(三种操作. 1.插入x 2.删除x 3.在集合中找到一个数y,使得y->x,每次能除以一个素数或者乘上一个素数,问最少的操作次数)

传送门:Hdu 5812 Distance

---

维护一个集合，有三种操作.
1.插入x 2.删除x 3.在集合中找到一个数y,使得y->x,每次能除以一个素数或者乘上一个素数,问最少的操作次数

---

思路:记f(x)为从1变到x所需要使用的最小的次数(可在线性筛的时候预处理)
要使y->x,即相当于对于每个x的因子z,从集合中找出一个y使得f[y/z]+f[x/z]最小,对于f[x/z]很好求解,
接下来问题的关键便转化为了f[y/z]的求解

C[x][s]表示集合Y中的数变到x需要处理s次的数有几个,D[x]表示现在能变到x所需要的次数的二进制表示
对于每次插入y,我们可以处理出它的因子,更新现在C[x][s],从而更新D[x]

---

#include<bits/stdc++.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int N=1e6;
void fre() { freopen("1004.in", "r", stdin); freopen("output.out", "w", stdout); }
int cnt[N+10],prime[N+10],tot,Hash[N+10],T=0;
int C[N+10][21],Vis[N+10][21],D[N+10],Min[N+10];
vector<int>G[N+10];
bool check[N+10],vis[N+10];

void init(){
    cnt[1]=0,tot=0;
    for(int i=2;i<=N;i++){
        if(!check[i]){
            prime[tot++]=i;
            cnt[i]=1;
        }
        for(int j=0;j<tot;j++){
            if(i*prime[j]>N)    break;
            check[i*prime[j]]=true;
            if(i%prime[j]==0){
                cnt[i*prime[j]]=cnt[i]+1;
                break;
            }
            else
                cnt[i*prime[j]]=cnt[i]+cnt[prime[j]];
        }
    }
}

char s[2];
void add(int num,int x){
    if(Vis[num][ cnt[x/num] ]!=T)
        C[num][ cnt[x/num] ]=1,Vis[num][ cnt[x/num] ]=T;
    else
        C[num][ cnt[x/num] ]++;
    if(C[num][ cnt[x/num] ]==1)
        D[num]|=(1<<cnt[x/num]);
}

void remove(int num,int x){
    C[num][ cnt[x/num] ]--;
    if(C[num][ cnt[x/num] ]==0)
        D[num]^=(1<<cnt[x/num]);
}

void insert(int x){
    if(Hash[x]==T)
        return ;
    Hash[x]=T;
    if(vis[x])
        for(int i=0;i<G[x].size();i++){
            int num=G[x][i];
            add(num,x);
        }
    else{
        int num=sqrt(x+0.5);
        for(int i=1;i<=num;i++)
            if(x%i==0){
                G[x].push_back(i),add(i,x);
                if(i*i!=x)
                    G[x].push_back(x/i),add(x/i,x);
            }
    }
    vis[x]=1;
}

void del(int x){
    if(Hash[x]!=T)
        return ;
    Hash[x]=T-1;
    for(int i=0;i<G[x].size();i++){
        int num=G[x][i];
        remove(num,x);
    }
}
int ans;

void getans(int num,int x){
    if(D[num]!=0)
        ans=min(ans,cnt[x/num]+Min[ D[num]&-D[num] ]);
}

void solve(int x){
    ans=INF;
    if(vis[x])
        for(int i=0;i<G[x].size();i++){
            int num=G[x][i];
            getans(num,x);
        }
    else{
        int num=sqrt(x+0.5);
        for(int i=1;i<=num;i++)
            if(x%i==0){
                G[x].push_back(i);
                getans(i,x);
                if(i*i!=x)
                    G[x].push_back(x/i),getans(x/i,x);
            }
    }
    vis[x]=1;
    if(ans==INF)
        printf("-1\n");
    else
        printf("%d\n",ans);
}

int main(){
    int n,x,case1=1;
    init();
    for(int Cnt=0,i=1;Cnt<=21&&i<=N;i<<=1,Cnt++)
        Min[i]=Cnt;
    while(scanf("%d",&n)!=EOF){
        if(n==0)
            break;
        printf("Case #%d:\n",case1++);
        ++T;
        memset(D,0,sizeof(D));
        for(int i=1;i<=n;i++){
            scanf("%s%d",s,&x);
            if(s[0]=='I')
                insert(x);
            else if(s[0]=='D')
                del(x);
            else
                solve(x);
        }
    }
    return 0;
}