hdu5812Distance(数学)

题目链接：点这里！！！

题意：

我们定义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

题解：

1、我们很容易知道d(x,y)=f(x/gcd(x,y))+f(y/gcd(x,y))。f(a)表示a质因子的个数。f数组我们可以预处理出来。

2、对于每个询问x，我们知道gcd(x,y)一定是x的约数，我们就去枚举x的约数i(gcd(x,y))，很容易得到f(x/gcd(x,y))的值，但是如何求f(y/gcd(x,y))??

3、其实对于每个在集合里的数，我们一样去枚举他的约数，然后对于每个约数i，我用“某个东西“存下f(x/i)。删除的话，逆着做就可以了。

4、有了”3“知道，我们就可以在2的情况下找到”某个东西“下gcd(x,y)存的最小值就ok啦！！

至于某个东西的话，我用了multset和(set+map)两种方法来实现，multset优了很多！！

代码：(multset实现)

#include<cstdio>
#include<cstring>
#include<iostream>
#include<sstream>
#include<algorithm>
#include<vector>
#include<bitset>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<cstdlib>
#include<cmath>
#define pb push_back
#define pa pair<int,int>
#define clr(a,b) memset(a,b,sizeof(a))
#define lson lr<<1,l,mid
#define rson lr<<1|1,mid+1,r
#define bug(x) printf("%d++++++++++++++++++++%d\n",x,x)
#define key_value ch[ch[root][1]][0]
#pragma comment(linker, "/STACK:102400000000,102400000000")
typedef  long long LL;
const LL  MOD = 1000000007;
const int N = 1e6+15;
const int maxn = 1e6+15;
const int letter = 130;
const int INF = 1e9;
const double pi=acos(-1.0);
const double eps=1e-10;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int vis[N],prime[N/10],num,cnt[N],mark[N];
multiset<int>ss[N];
void init(){
    int sum=0;
    cnt[1]=0;
    vis[1]=1;
    for(int i=2;i<N;i++){
        if(!vis[i]) prime[++num]=i,cnt[i]=1;
        for(int j=1;j<=num&&prime[j]*i<N;j++){
            vis[prime[j]*i]=1;
            cnt[i*prime[j]]=cnt[i]+1;
            if(i%prime[j]==0) break;
        }
    }
}
void add(int x,int y){
    ss[x].insert(y);
}
void del(int x,int y){
    auto it=ss[x].find(y);
    if(it!=ss[x].end()) ss[x].erase(it);
}
int get(int x){
    if(ss[x].empty()) return INF;
    return *ss[x].begin();
}
int main(){
    init();
    int q,cas=0;
    while(~scanf("%d",&q)&&q){
        clr(mark,0);
        for(int i=0;i<N;i++) ss[i].clear();
        char s[3];
        int x,sum=0;
        printf("Case #%d:\n",++cas);
        while(q--){
            scanf("%s%d",s,&x);
            if(s[0]=='I'){
                if(mark[x]) continue;
                mark[x]=1,sum++;
                int i;
                for(i=1;i*i<x;i++){
                    if(x%i==0) add(i,cnt[x/i]),add(x/i,cnt[i]);
                }
                if(i*i==x) add(i,cnt[i]);
            }
            else if(s[0]=='D'){
                if(!mark[x]) continue;
                mark[x]=0,sum--;
                int i;
                for(i=1;i*i<x;i++){
                    if(x%i==0) del(i,cnt[x/i]),del(x/i,cnt[i]);
                }
                if(i*i==x) del(i,cnt[i]);
            }
            else {
                if(sum==0) {puts("-1");continue;}
                int ans=INF;
                for(int i=1;i*i<=x;i++){
                    if(x%i==0) ans=min(ans,min(get(i)+cnt[x/i],get(x/i)+cnt[i]));
                }
                printf("%d\n",ans);
            }
        }
    }
    return 0;
}

代码：(set+map实现)

#include<cstdio>
#include<cstring>
#include<iostream>
#include<sstream>
#include<algorithm>
#include<vector>
#include<bitset>
#include<set>
#include<queue>
#include<stack>
#include<map>
#include<cstdlib>
#include<cmath>
#define pb push_back
#define pa pair<int,int>
#define clr(a,b) memset(a,b,sizeof(a))
#define lson lr<<1,l,mid
#define rson lr<<1|1,mid+1,r
#define bug(x) printf("%d++++++++++++++++++++%d\n",x,x)
#define key_value ch[ch[root][1]][0]
#pragma comment(linker, "/STACK:102400000000,102400000000")
typedef  long long LL;
const LL  MOD = 1000000007;
const int N = 1e6+15;
const int maxn = 1e6+15;
const int letter = 130;
const int INF = 1e9;
const double pi=acos(-1.0);
const double eps=1e-10;
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int vis[N],prime[N/10],num,cnt[N],mark[N];
set<int>ss[N];
map<int,int>mp[N];
void init(){
    int sum=0;
    cnt[1]=0;
    vis[1]=1;
    for(int i=2;i<N;i++){
        if(!vis[i]) prime[++num]=i,cnt[i]=1;
        for(int j=1;j<=num&&prime[j]*i<N;j++){
            vis[prime[j]*i]=1;
            cnt[i*prime[j]]=cnt[i]+1;
            if(i%prime[j]==0) break;
        }
    }
}
void add(int x,int y){
    if(ss[x].find(y)==ss[x].end()) ss[x].insert(y);
    mp[x][y]++;
}
void del(int x,int y){
    mp[x][y]--;
    if(mp[x][y]==0) ss[x].erase(y);
}
int get(int x){
    if(ss[x].empty()) return INF;
    return *ss[x].begin();
}
int main(){
    init();
    int q,cas=0;
    while(~scanf("%d",&q)&&q){
        clr(mark,0);
        for(int i=0;i<N;i++) ss[i].clear(),mp[i].clear();
        char s[3];
        int x,sum=0;
        printf("Case #%d:\n",++cas);
        while(q--){
            scanf("%s%d",s,&x);
            if(s[0]=='I'){
                if(mark[x]) continue;
                mark[x]=1,sum++;
                int i;
                for(i=1;i*i<x;i++){
                    if(x%i==0) add(i,cnt[x/i]),add(x/i,cnt[i]);
                }
                if(i*i==x) add(i,cnt[i]);
            }
            else if(s[0]=='D'){
                if(!mark[x]) continue;
                mark[x]=0,sum--;
                int i;
                for(i=1;i*i<x;i++){
                    if(x%i==0) del(i,cnt[x/i]),del(x/i,cnt[i]);
                }
                if(i*i==x) del(i,cnt[i]);
            }
            else {
                if(sum==0) {puts("-1");continue;}
                int ans=INF;
                for(int i=1;i*i<=x;i++){
                    if(x%i==0) {
                            ans=min(ans,min(get(i)+cnt[x/i],get(x/i)+cnt[i]));
                    }
                }
                printf("%d\n",ans);
            }
        }
    }
    return 0;
}