cp's OEIS Frontend

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.

A337120 Factor complexity (number of subwords of length n) of the regular paperfolding sequence (A014577), and all generalized paperfolding sequences.

Original entry on oeis.org

1, 2, 4, 8, 12, 18, 23, 28, 32, 36, 40, 44, 48, 52, 56, 60, 64, 68, 72, 76, 80, 84, 88, 92, 96, 100, 104, 108, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 180, 184, 188, 192, 196, 200, 204, 208, 212, 216, 220, 224, 228
Offset: 0

Views

Author

Kevin Ryde, Aug 17 2020

Keywords

Examples

			For n=4, all length 4 subwords except 0000, 0101, 1010, 1111 occur, so a(4) = 16-4 = 12.  (These words do not occur because odd terms in a paperfolding sequence alternate, so a subword wxyz must have w!=y or x!=z.)
		

Crossrefs

Cf. A014577, A214613 (Abelian complexity), A333994 (arithmetical complexity).
Cf. A005943 (GRS).

Programs

  • Mathematica
    LinearRecurrence[{2, -1}, {1, 2, 4, 8, 12, 18, 23, 28, 32}, 100] (* Paolo Xausa, Feb 29 2024 *)
  • PARI
    Vec((1 + x^2)*(1 + 2*x^3 - x^6) / (1 - x)^2 + O(x^50)) \\ Colin Barker, Sep 08 2020

Formula

a(1..6) = 2,4,8,12,18,23, then a(n) = 4*n for n>=7. [Allouche]
From Colin Barker, Sep 05 2020: (Start)
G.f.: (1 + x^2)*(1 + 2*x^3 - x^6) / (1 - x)^2.
a(n) = 2*a(n-1) - a(n-2) for n>8.
(End)