hdu 5812 Distance 数论+思维枚举

题目连接

题意:

我们定义d(a,b)为a通过乘除最少的素数到达b。比如d(15,50)=3，15*2*5/3=50。
我们定义3个操作和一个空的集合s：
I x表示在集合s里插入一个x，如果存在x，忽略这个操作。
D x表示在集合x里删除x，如果不存在x，忽略这个操作。
Q x找到最小的z集合里面存在一个y，使得d(x,y)=z。
数据范围：
总操作数q<=50000,x<=100000

思路:

首先这个题目我们要发现怎么样才能快速的从x变到y，由于操作都是素数,我们可以发现我们每次*的都是x没有的只属于y的素因子,/的都是只属于x的素因子.
定义f(w)为w中质因子的个数,那么可以得到 d(x,y) = f(x/gcd(x,y)) + f(y/gcd(x,y)) .

有了上述的这个式子,我们对于每个x,可以$O(\sqrt(x))$ 枚举他的每一个约数i, 作为gcd(x,y) 。现在的问题变成了 已知f（x/gcd(x,y)）如何求最小的f(y/gcd(x,y)) .这里我们可以发现由于gcd已知,我们可以预处理出所有的 f（num/gcd） 表示的对应的质因子的个数 , 这里有一个技巧,我们将所有的f(num/gcd) 都放到对应的gcd组里,这样每次枚举一个gcd,我们就可以从对应组里拿出一个最小的来,这样就是f(y/gcd)了,维护这个过程从而得到最小的z。这个过程可以采用multiset来实现.

复杂度 $O(q * \sqrt(x) * logn)$ 常数大点,但是足够过这题了

#include<bits/stdc++.h>
#define pb push_back
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int maxn = 1e6+5;
bool book[maxn];
bool vis[maxn];
char s[10];
int  cnt,q,num[maxn];
vector<int>prime;
multiset<int>st[maxn];
multiset<int> ::iterator it;
void shai()
{
    vis[1] = 1;
    num[1] = 0;
    for(int i = 2;i < maxn;++i)
    {
        if(!vis[i])
        {
            prime.pb(i);
            vis[i] = 1;
            num[i] = 1;
        }
        for(int j = 0;j < prime.size()&&prime[j]*i < maxn;++j)
        {
            vis[i*prime[j]] = 1;
            num[i*prime[j]] = num[i] + 1;
            if(i % prime[j] == 0) break;
        }
    }
}
void init()
{
    cnt = 0;
    for(int i = 0;i < maxn;++i)
    {
        st[i].clear();
        book[i] = 0;
    }
    return ;
}
void add(int x,int y)
{
    st[x].insert(y);
}
void del(int x,int y)
{
    it = st[x].find(y);
    if(it != st[x].end())
    st[x].erase(it);
}
int get(int x)
{
    if(st[x].empty()) return inf;
    return *st[x].begin();
}
void Insert(int x)
{
    for(int i = 1;i*i <= x;++i)
    {
        if(x % i ==0)
        {
            add(i,num[x/i]);
            if(i != x/i)
            add(x/i,num[i]);
        }
    }
    return ;
}
void Del(int x)
{
    for(int i = 1;i*i <= x;++i)
    {
        if(x % i == 0)
        {
            del(i,num[x/i]);
            if(i != x / i)
            del(x/i,num[i]);
        }
    }
    return ;
}
int find(int x)
{
    int res = inf;
    for(int i = 1;i*i <= x;++i)
    {
        if(x % i == 0)
        {
            res = min(res,min(num[x/i]+get(i),num[i]+get(x/i)));
        }
    }
    return res;
}
int main()
{
    shai();
    int ca = 1;
    while(~scanf("%d",&q))
    {
        if(q == 0) break;
        init();
        int x;
        printf("Case #%d:\n",ca++);
        for(int i = 1;i <= q;++i)
        {
            scanf("%s %d",s,&x);
            if(s[0] =='I')
            {
                if(book[x]) continue;
                book[x] = 1;
                cnt++;
                Insert(x);
            }   
            else if(s[0] =='D')
            {
                if(!book[x])
                continue;
                book[x] = 0;
                cnt--;
                Del(x);
            }
            else
            {
                if(cnt ==  0)
                {
                    puts("-1");
                    continue;
                }
                ll ans = find(x);
                printf("%lld\n",ans);
            }
        }   
    }
    return 0;
}