A119614 a(1)=1. a(2^m +k) = a(2^m + 1 - k)*a(k) + 1, where 1 <= k <= 2^m, m>=0.

1, 2, 3, 3, 4, 7, 7, 4, 5, 15, 22, 13, 13, 22, 15, 5, 6, 31, 67, 40, 53, 155, 106, 21, 21, 106, 155, 53, 40, 67, 31, 6, 7, 63, 202, 121, 213, 1086, 743, 85, 106, 1591, 3411, 690, 521, 1475, 466, 31, 31, 466, 1475, 521, 690, 3411, 1591, 106, 85, 743, 1086, 213, 121, 202

Offset

1,2

Links

John Tyler Rascoe, Table of n, a(n) for n = 1..8192

Maple

A119614 := proc (mmax) local a, m, k, an; a := [1] ; for m from 0 to mmax do for k from 1 to 2^m do an := a[2^m+1-k]*a[k]+1 ; a := [op(a), an] ; od ; od ; RETURN(a) ; end: an := A119614(6) : for i from 1 to nops(an) do printf("%d, ", an[i]) ; od ; # R. J. Mathar, Aug 06 2006

Prog

(Python)
from itertools import count, islice
def a_gen():
    A = [1]
    yield(1)
    for m in count(0):
        for k in range(2**m):
            x = A[(2**m)-k-1]*A[k] + 1
            A.append(x)
            yield(x)
A119614_list = list(islice(a_gen(), 2**6)) # John Tyler Rascoe, May 10 2024

Crossrefs

Sequence in context: A269098 A119795 A207618 * A260167 A035540 A114863

Adjacent sequences: A119611 A119612 A119613 * A119615 A119616 A119617

Keyword

easy,nonn,look

Author

Leroy Quet, Jun 05 2006

Extensions

More terms from R. J. Mathar, Aug 06 2006

Status

approved