HDU 5812 Distance 暴力+素数打表

/*
   I insert x;
   D delect x;
   Q query min(d(x,y);
     d(x,y) =x 乘或除一个素数 变成y的最小次数;

   d(x,y) =ans[x/gcd(x,y)]+ans[y/gcd(x,y]; ans[i] i的素因子数；
   保留之前插入 y时,每一个y因子对应的ans[y/i]值 ->因子，ans值较小，
   作为查询表的key 统计满足条件的y个数 ->C[i][ans[y/i]
   查询暴力每个因子，如果有gcd时,显然最小。 
*/ 
#include<iostream>
#include<algorithm> 
#include<cstdio>
#include<cmath> 
#include<cstring>
#include<set>
using namespace std;
const int Max=1000000+1; 
int ans[Max]={0};
int C[Max][21];
void init()
{
	for(int i=2;i<=1000000;i++)
	  {
	  	if(!ans[i]) 
		  {
		  ans[i]=1;
		  for(int j=i+i;j<=1000000;j +=i) ans[j]=ans[j/i]+1;
	      }
	  } 
}
int main()
{
	int n,x,cas=0;
	init();
	char q[2];
	set<int> s;
	while(scanf("%d",&n)&&n)
	{
		printf("Case #%d:\n",++cas);
		s.clear();
		memset(C,0,sizeof(C));
		for(int i=1;i<=n;i++)
		 {
		 	scanf("%s%d",q,&x);
		 	if(q[0]=='I')
		 	{
		 		if(s.count(x)) continue;
		 		else
		 		{
		 			s.insert(x);
					for(int i=1;i<=sqrt(x)+1;i++)
					{
						if(x%i==0) 
						{
						  C[i][ans[x/i]]++;
						  C[x/i][ans[i]]++;
					    }
					} 
				}
			}
			if(q[0]=='D')
			{
				if(!s.count(x)) continue;
				else
				{
					s.erase(x);
					for(int i=1;i<=sqrt(x)+1;i++)
					{
						if(x%i==0) 
						{
						  C[i][ans[x/i]]--;
						  C[x/i][ans[i]]--;
					    }
					}
				}
			} 
		 if(q[0]=='Q')
		 {
		 	int Min=1234;
		 	if(s.size()==0) 
			 {
			   printf("-1\n");
			   continue;
		     }
		 	for(int i=1;i<=sqrt(x)+1;i++)
		 	{
		 		if(x%i==0)
		 		{
		 			int j;
		 			for(j=0;j<=20;j++)
					  if(C[i][j])
					  {
					   Min=min(Min,ans[x/i]+j);
					   break;
				      }
		 		    
					for(j=0;j<=20;j++) 
					if(C[x/i][j])
					 {
					 Min=min(Min,ans[i]+j);
					 break;
				     }
				}
			}
		  	printf("%d\n",Min);
		 } 
		} 
	}
	return 0;
}