int sieve [MAX_N] ;
void init_sieve(int n)  {
    for (int i = 1; i <= n; i++)
        sieve[i] = i;
    for (int i = 2; i <= n; i++)
        for (int j = i; j <= n; j += i )
            sieve[j]--;
}