const int X=1000010;
bool isPrime[X+1];
int total;//计数
int prime[79000];
void getPrime()
{
total=0;
memset(isPrime,true,sizeof(isPrime));
memset(prime,0,sizeof(prime));
for(int i=2;i<=X;i++)
{
if(isPrime[i]) prime[total++]=i;
for(int j=0; j<total && i*prime[j]<=X; j++)
{
isPrime[i*prime[j]]=false;
if(i%prime[j]==0)
break;
}
}
}
int num[X];
void getCot()//每个数质因子的个数
{
memset(num,0,sizeof(num));
for(int i=0;prime[i]<=1000000;i++)
for(int j=prime[i];j<=1000000;j+=prime[i])
num[j]++;
}
//有可能会超内存,建议用上面的
//这个求的是每个数的质因子情况
#define X 1000005
int visited[X];
vector<int>a[X];
void init()
{
int i,j;
for(i=0;i<X;i++)
a[i].clear();//vector的清空
memset(visited,0,sizeof(visited));
for(i=2;i<=100000;i++)
{
if(visited[i]==0)//i是素数这是可以保证的
{
a[i].push_back(i);
for(j=i+i;j<=100000;j+=i)//筛选素数,其实这种方法没以前那种素数法快,但是这里用来求一个数的质因子就比较好了
{
visited[j]=1;
a[j].push_back(i);
}
}
}
}