|
|
A033485
|
|
a(n) = a(n-1) + a(floor(n/2)), a(1) = 1.
(Formerly N0236)
|
|
35
|
|
|
1, 2, 3, 5, 7, 10, 13, 18, 23, 30, 37, 47, 57, 70, 83, 101, 119, 142, 165, 195, 225, 262, 299, 346, 393, 450, 507, 577, 647, 730, 813, 914, 1015, 1134, 1253, 1395, 1537, 1702, 1867, 2062, 2257, 2482, 2707, 2969, 3231, 3530, 3829, 4175, 4521, 4914, 5307, 5757
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Sequence gives the number of partitions of 2n into "strongly decreasing" parts (see the function s*(n) in the paper by Bessenrodt, Olsson, and Sellers); see the example in A040039.
Partial sums of the sequence a(1)=1, a(1), a(1), a(2), a(2), a(3), a(3), a(4), a(4), a(5), a(5), a(6), ...; example: a(1) = 1, a(2) = 1+1 = 2, a(3) = 1+1+1 = 3, a(4) = 1+1+1+2 = 5, a(5) = 1+1+1+2+2 = 7, ... - Philippe Deléham, Jan 02 2004
The number of odd numbers before the n-th even number in this sequence is A003156(n). - Philippe Deléham, Mar 27 2004
There are no terms divisible by 4 and there are infinitely many terms divisible by {2,3,5,7} (see Computational Complexity link). - Ivan N. Ianakiev, Aug 06 2022
|
|
REFERENCES
|
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
Conjecture: lim_{n->infinity} a(2n)/a(n)*log(n)/n = c = 1.64.... and a(n)/A(n) is bounded where A(n)=1 if n is a power of 2, otherwise A(n) = sqrt(n)*Product_{k<log_2(n)} n/2^k/log(n/2^k)). - Benoit Cloitre, Jan 26 2003
G.f.: (1/2)*(((1-x)*Product_{n>=0} (1-x^(2^n)))^(-1)-1). a(n) modulo 4 = A007413(n). - Philippe Deléham, Feb 28 2004
|
|
MAPLE
|
a:= proc(n) option remember;
`if`(n<2, n, a(n-1)+a(iquo(n, 2)))
end:
|
|
MATHEMATICA
|
b[1]=1; b[n_] := b[n]=Sum[b[k], {k, 1, n/2}]; Table[b[n], {n, 3, 105, 2}] (* Robert G. Wilson v, Apr 22 2001 *)
|
|
PROG
|
(PARI) a(n)=if(n<2, 1, a(floor(n/2))+a(n-1))
(Haskell)
import Data.List (transpose)
a033485 n = a033485_list !! (n-1)
a033485_list = 1 : zipWith (+)
a033485_list (concat $ transpose [a033485_list, a033485_list])
(Magma) [n le 1 select 1 else Self(n-1) + Self(Floor(n/2)) : n in [1..60]]; // Vincenzo Librandi, Nov 20 2015
(Python)
from itertools import islice
from collections import deque
def A033485_gen(): # generator of terms
aqueue, f, b, a = deque([2]), True, 1, 2
yield from (1, 2)
while True:
a += b
yield a
aqueue.append(a)
if f: b = aqueue.popleft()
f = not f
|
|
CROSSREFS
|
Also half of A000123 (with first term omitted).
|
|
KEYWORD
|
nonn,nice,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|