This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
read(transforms) : # ANDnos def'd here A353709 := proc(n) option remember; local c, i, known ; if n <= 2 then n; else for c from 1 do known := false ; for i from 1 to n-1 do if procname(i) = c then known := true; break ; end if; end do: if not known and ANDnos(c, procname(n-2)) = 0 and ANDnos(c, procname(n-1)) = 0 then return c; end if; end do: end if; end proc: # Following R. J. Mathar's program for A109812. [seq(A353709(n), n=0..256)] ; # second Maple program: b:= proc() false end: t:= 2: a:= proc(n) option remember; global t; local k; if n<2 then n else for k from t while b(k) or Bits[And](k, a(n-2))>0 or Bits[And](k, a(n-1))>0 do od; b(k):=true; while b(t) do t:=t+1 od; k fi end: seq(a(n), n=0..100); # Alois P. Heinz, May 06 2022
nn = 83; c[] = -1; a[0] = c[0] = 0; a[1] = c[1] = 1; u = 2; Do[k = u; While[Nand[c[k] == -1, BitAnd[a[n - 1], k] == 0, BitAnd[a[n - 2], k] == 0], k++]; Set[{a[n], c[k]}, {k, n}]; If[k == u, While[c[u] > -1, u++]], {n, 2, nn}]; Array[a, nn+1, 0] (* _Michael De Vlieger, May 06 2022 *)
from itertools import count, islice def A353709_gen(): # generator of terms s, a, b, c, ab = {0,1}, 0, 1, 2, 1 yield from (0,1) while True: for n in count(c): if not (n & ab or n in s): yield n a, b = b, n ab = a|b s.add(n) while c in s: c += 1 break A353709_list = list(islice(A353709_gen(),20)) # Chai Wah Wu, May 07 2022
A353709 has offset 0 and begins 0, 1, 2, 4, 8, 3, 16, 12, 32, 17, 6, 40, 64, 5, 10, ..., so a(1) = 2 (from A353709(2) = 2), and a(7) = 9 (from A353709(9) = 17 = prime(7)).
Accumulate[Flatten[Map[ConstantArray[LucasL[#],LucasL[#-1]]&,Range[15]]]] (* Peter J. C. Moses, May 02 2022 *)
from itertools import islice def A351628_gen(): # generator of terms a, b, c = 1, 3, 0 while True: yield from (c+i*a for i in range(1,b-a+1)) a, b, c = b, a+b, c + a*(b-a) A351628_list = list(islice(A351628_gen(),40)) # Chai Wah Wu, Jun 09 2022
Comments