bool NotPrime[40005];
long long Prime[40005];
void init()//得到素數
{
    long long i,j,num=0;
    for(i=2;i<=40000;i++) //注意是從2開始的，直到你想要的范圍。
    {
        if(!NotPrime[i])  //如果不是不是素數（也就是說是素數了）
        {
            Prime[num++]=i; //將素數i存儲在prime數組中。
            for(j=i*i;j<=40000;j=j+i)
                NotPrime[j]=true;
        }
    }
}