HDU 5812 Distance

从a变到b，也就是将a一直除素因子，除到1为止，然后乘b的素因子，一直乘到b。

但是gcd(a,b)部分是不用除下去的。所以d(a,b)=a/gcd(a,b)的素因子个数+b/gcd(a,b)的素因子个数。

然后....脑洞开始......

枚举这个因子P，然后去计算a/P的素因子个数+b/P的素因子个数 最小值，也就是要在集合中寻找到一个b，使得式子取得最小值。

然后....可以搞一个数组 c[因子][素因子个数]=多少种情况 来存储集合中的信息。 然后删除，插入的时候维护这个数组即可。具体看看代码。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<iostream>
using namespace std;
typedef long long LL;
const double pi=acos(-1.0),eps=1e-;
void File()
{
    freopen("D:\\in.txt","r",stdin);
    freopen("D:\\out.txt","w",stdout);
}
inline int read()
{
    char c = getchar();  while(!isdigit(c)) c = getchar();
    int x = ;
    while(isdigit(c)) { x = x *  + c - ''; c = getchar(); }
    return x;
}

const int maxn=;
int f[maxn],c[maxn][],fac[maxn],sz;
bool m[maxn];
queue<int>Q;

bool prime(int x)
{
    if(x==) return ;
    for (int i=;i*i<=x;i++) if(x%i==) return ;
    return ;
}

void pre()
{
    memset(f,,sizeof f);
    for(int i=;i<=;i++)
    {
        if(!prime(i)) continue; f[i]=;
        int u=i+i; while(u<=) { int t=u; while(t%i==) f[u]++,t=t/i; u=u+i; }
    }
}

void get(int x)
{
    sz=;
    for(int i=;i*i<=x;i++)
    {
        if(x%i!=) continue; int x1=i,x2=x/i;
        fac[sz++]=x1; if(x2!=x1) fac[sz++]=x2;
    }
}

void Insert(int x) { m[x]=; get(x); for(int i=;i<sz;i++) c[fac[i]][f[x/fac[i]]]++; }
void Delete(int x) { m[x]=; get(x); for(int i=;i<sz;i++) c[fac[i]][f[x/fac[i]]]--; }

void Find(int x)
{
    int ans=; get(x);
    for(int i=;i<sz;i++) for(int j=;j<=;j++)
        if(c[fac[i]][j]) ans=min(ans,j+f[x/fac[i]]);
    if(ans==) ans=-; printf("%d\n",ans);
}

int main()
{
    pre(); int cas=,n;
    while(~scanf("%d",&n))
    {
        if(n==) break; printf("Case #%d:\n",cas++);
        memset(c,,sizeof c); memset(m,,sizeof m);
        for(int i=;i<=n;i++)
        {
            char op[]; int x; scanf("%s%d",op,&x);
            if(op[]=='I') { if(m[x]) continue; Insert(x); }
            else if(op[]=='D') { if(!m[x]) continue; Delete(x); }
            else Find(x);
        }
    }
    return ;
}