#include <stdio.h>
#include <math.h>
#define N 50
#define MAX_PRIME (3*N*N*N)
int c3[N];
int c2[N];
int not_prime[MAX_PRIME];
#define L (3000) //L must be more than sqrt(MAX_PRIME)
#define IS_PRIME(x) (!not_prime[x])
void init_prime(){
int i;
not_prime[0]=not_prime[1]=1;
for(i=2;i<L;i++){
if(IS_PRIME(i)){
int j;
for(j=i+i;j<MAX_PRIME;j+=i)not_prime[j]=1;
}
}
}