hdu-5812- Distance-（数学+枚举约数+素筛法(素因子数量表)）

http://acm.hdu.edu.cn/showproblem.php?pid=5812

题目的询问是给一个Q,X ，在给出的集合里找到一个y满足d(x,y)最小。

1、首先素筛预处理出ans[i]，表示i能分解成的质因数个数。 

2、要从x变成y，就是从x出发去掉x独有的素因子（gcd），再加上所有y的独有的素因子。

3、所以题意的d(x,y)可以简要地表示为 d(x,y)= ans[x/gcd] + ans[y/gcd]， 其中的gcd是指的 gcd(x,y)即她俩的公约数

4、那么对于每一次询问 Q X， 我们可以大胆放心地枚举X的所有约数i，也即当作i即为d(x,y)的gcd，  对于答案 d(x,y)= ans[x/gcd] + ans[y/gcd]，既然我们枚举了gcd，那么接下来就是要找 对于这个gcd下，枚举y咯（是不是很暴力啊啊啊）

5、我们需要（暴力）维护一个表int c[MAX][20];//c[i][j]表示i的倍数Y,满足ans[Y/i]==j的数量

由于数最大也就是1000000，不超过2^20,所以第二维开20够了。

每次操作都要分解质因素，复杂度sqrt(x)，总复杂度q*sqrt(x)

#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <algorithm>
#include <iostream>
#include <queue>
#include <deque>
#include <set>
#include <vector>
using namespace std;
typedef long long ll;
const int MAX=1000000;
bool f[MAX+50];
int ans[MAX+50];
int c[MAX][20];//c[i][j]表示i的倍数Z,满足[Z/i]==j的数量
set<int >sb;
set<int >::iterator it;
void pre1()
{
    f[1]=true;
    for (ll i=2; i<=MAX; i++)
    {
        if (f[i]==false)
        {
            ans[i]=1;
            for (ll j= i+i; j<=MAX; j=j+i)
            {
                f[j]=true;
                if (j%i==0)
                    ans[j]=ans[j/i]+1;
            }
        }
    }
}
void add(int x,int times)
{
    for (int i=1; i*i<=x; i++)
    {
        if (x%i==0)
        {
            c[i][ans[x/i]]+=times;
            if (i!=x/i)
                c[x/i][ans[i]]+=times;
        }
    }
}
int main()
{
    pre1();
    int n;
    int cnt=1;
    while(scanf("%d",&n)!=EOF&&n)
    {
        for (it=sb.begin();it!=sb.end();it++)  add(*it,-1);
        sb.clear();
        printf("Case #%d:\n",cnt++);
        char op[5];
        int x;
        for (int i=1; i<=n; i++)
        {
            scanf("%s%d",op,&x);
            if (op[0]=='I')
            {
                if (!sb.count(x))
                    add(x,1),sb.insert(x);
            }
            else if (op[0]=='D')
            {
                if (sb.count(x))
                    add(x,-1),sb.erase(x);
            }
            else
            {
                if(sb.empty())
                {
                    printf("-1\n");continue;
                }
                int num=1e9;
                for (int j=1; j*j<=x; j++) //枚举gcd-约数
                {
                    if (x%j==0)
                    {
                        for (int k=0; k<20; k++)
                        {
                            if (c[j][k])//表示j(gcd)的倍数Z且 ans[z/j]==k
                                num=min(num,ans[x/j]+k);
                            if (c[x/j][k])
                                num=min(num,ans[j]+k);
                        }
                    }
                }
                printf("%d\n",num);
            }
        }


    }



    return 0;

}