>>655
その場合、1024bitでモジュラ演算をしなくちゃいけないでしょ？
ビットマップがわからないなら、フラグ配列を考えよう、

FlagArray[n] に対応する RandomMap + n を篩うようにするわけ n はM迄の小さな数ね
素数 kの倍数を篩う時は

(RandomMap + Nk )% k =0 となる最小のNkを求めておいて

for(i=0;i<M/3;i++) FlagArray[n3+i*3] = 素数でない　
for(i=0;i<M/5;i++) FlagArray[n5+i*5] = 素数でない　
　　・・・
とやってくわけ