歐拉函數 素數篩選法模板

long long p[maxn],e[maxn];
bool hash[maxn];
void init()
{
		memset(hash,0,sizeof(hash));
		memset(hash,0,sizeof(e));
    long long i, j;
    
    p[0] = 1; //記錄素數個數
    p[1] = 2;
    for (i=3; i<N; i+=2)
    {
        if (hash[i])
            continue;
        p[++p[0]] = i;
        for (j=i*i; j<N; j+=i)
            hash[j] = true;
    } //篩素數
    
    e[1] = 1;

    for (i=1; i<=p[0]; i++)
        e[p[i]] = p[i] - 1; //初始化素數的phi

    for (i=2; i<N; i++)
    {
        if(!e[i])
        {
            for (j=1; j<=p[0]; j++)
                if (i % p[j]==0)
                {
                    if (i / p[j] % p[j])
                        e[i] = e[i / p[j]] * e[p[j]];
                    else
                        e[i] = e[i / p[j] ]* p[j];
                    break;
                } // 利用上述性質求解
        }        
    }
    return ;
}