login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A007824 a(n) = f(a(n-1)), with f(m) = Sum i*b(i)*2^(i-1), m = Sum b(i)*2^i, and starting value 16. 1
16, 32, 80, 224, 720, 2976, 15120, 92448, 704592, 6455520, 69518032, 897158048, 12700875536, 203746961696, 3674124288080, 74060581880032, 1692860294726352, 42020514867170208, 1144626896435067664, 33222330682700081440 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,1

COMMENTS

Previous name was : For N=Sum a(i).2^i, a(i)=0,1, set D(N)/D(2)=Sum i.a(i).2^(i-1); sequence gives n-th derivative of 16.

Actually, it appears that c(n) for n=0, 1, 2, ... gives 0, 0, 1, 1, 4, 4, 5, 5, 12, 12, 13, 13, 16, 16, 17, 17, 32, 32, 33, 33, 36, ... that is A136013. - Michel Marcus, Jul 18 2013

LINKS

Table of n, a(n) for n=0..19.

PROG

(PARI) D(n) = {b = binary(n); lb = #b; rb = vector(lb, i, b[lb-i+1]); return(sum (i=0, lb-1, i*2^(i-1)*rb[i+1])); }

lista(nn) = {m = 16; for (n=0, nn, print1(m, ", "); m = D(m); ); } \\ Michel Marcus, Jul 18 2013

CROSSREFS

Sequence in context: A271634 A031463 A045071 * A235056 A115686 A088112

Adjacent sequences:  A007821 A007822 A007823 * A007825 A007826 A007827

KEYWORD

nonn

AUTHOR

Ralph Buchholz [ ralph(AT)defcen.gov.au ], Leisa Condie

EXTENSIONS

Edited by Michel Marcus, Jul 18 2013

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 24 14:11 EDT 2019. Contains 326282 sequences. (Running on oeis4.)