int  sprime()
 {
     int  i, j;
     int  b_size;
     int  MAXN  =   40000 ;
     for  (i = 2 ; i < MAXN; i ++ ) prime[i]  =   1 ;
     for  (i = 2 ; i < MAXN; i ++ )
     {
         if  (prime[i]  ==   1 )
             for  (j = 2 ; i * j < MAXN; j ++ )
                prime[i * j]  =   0 ;
    }
     for (i  =   0 , j  =   0 ; j  <  MAXN; j ++ )
     {
         if (prime[j]  ==   1 ) {
            b[i]  =  j;
            i ++ ;
        }
    }
    b_size  =  i;
     return  b_size;
}