# Python program for A284462 # Michael S. Branicky, Aug 18 2021 print(data) from numba import njit @njit def ok(k, bits): for i in range(bits-1, 0, -1): left1 = (1 << i) mask = left1 - 1 suffix = (mask & k) + left1 count = 1 for j in range(1, bits - i + 1): sj = (((mask << j) & k) >> j) + left1 if sj == suffix: count += 1 if count > 1: return count == 2 return 0 @njit def a(n): if n == 1: return 2 c = 0 for i in range(2**(n-1)): k = (1 << (n-1)) + i if ok(k, n): c += 1 return 2*c print([a(n) for n in range(1, 23)]) # about 1 second on Google Colab # print n a(n) elapsed time, AND # write b-file from time import time time0 = time() alst = [] for n in range(1, 10001): an = a(n) alst.append(an) print(n, an, time()-time0, flush=True) if n < len(data): print(" ", data[n-1]==an, flush=True) with open('b284462.txt', 'a') as f: f.write(f"{n} {an}\n") # NEW TERMS from a single Colab run # 33 6131253398 3336.0645339488983 # 34 12266876820 7025.415410518646 # 35 24542430776 14766.55267739296 # 36 49101773008 30945.675590991974 # 37 98235992952 64642.07109642029