素數打表，質因數分解

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);
            }
        }
    }
}