login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A212500 a(n) is the difference between multiples of 5 with even and odd digit sum in base 4 in interval [0,4^n). 29
1, 4, 7, 36, 65, 340, 615, 3220, 5825, 30500, 55175, 288900, 522625, 2736500, 4950375, 25920500, 46890625, 245522500, 444154375, 2325622500, 4207090625, 22028612500, 39850134375, 208658012500, 377465890625, 1976437062500, 3575408234375 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Let the term "transform" mean the operation of summing the products of the numbers in the n-th row of an m-nomial triangle (m-nomial T(n,k)) and  the ascending numbers of a sequence. And let T(0,0) be the top entry (0th row, 0th column) in an m-nomial triangle.  Then starting with a(1)=1, the bisection of this sequence (1,7,65,615,5825...) is the quadrinomial (4-nomial) transform of A000045 (Fibonacci sequence, F(j)), starting at T(0,0)=1, F(1)=1. - Bob Selcoe, May 24 2014

LINKS

Vincenzo Librandi, Table of n, a(n) for n = 1..1000

V. Shevelev, On monotonic strengthening of Newman-like phenomenon on (2m+1)-multiples in base 2m, arXiv:0710.3177 [math.NT], 2007.

Index entries for linear recurrences with constant coefficients, signature (0,10,0,-5).

FORMULA

For n>=5, a(n) = 10*a(n-2)-5*a(n-4).

a(n) = 0.4*((5+2*sqrt(5))^(n/2)+ (5-2*sqrt(5))^(n/2)) , if n is even, and

a(n) = 0.1*((5+2*sqrt(5))^((n-1)/2)*sqrt(30+10*sqrt(5))+(5-2*sqrt(5))^((n-1)/2)*sqrt(30-10*sqrt(5))), if n is odd.

a(2n+1) = Sum_{j=0..3n+1} fibonacci(j+1)*A008287(n,j). - Bob Selcoe, May 28 2014

G.f.: -x*(-1-4*x+3*x^2+4*x^3) / ( 1-10*x^2+5*x^4 ). - R. J. Mathar, Jun 16 2014

EXAMPLE

Let n=3. In interval [0,4^3) we have 13 multiples of 5,from which in base 4 only three (namely, 35,50,55) have odd digit sums. Thus a(3)=(13-3)-3=7.

From Bob Selcoe, May 28 2014: (Start)

n=2: a(5)=65 because T(2,k) {k=0..6} is {1,2,3,4,3,2,1} and {j=1..7} is {1,1,2,3,5,8,13}: 1*1 + 2*1 + 3*2 + 4*3 + 3*5 + 2*8 + 1*13 = 65.

n=3: a(7)=615 because T(3,k) {k=0..9} is {1,3,6,10,12,12,10,6,3,1} and {j=1..10} is {1,1,2,3,5,8,13,21,34,55}: 1*1 + 3*1 + 6*2 + 10*3 + 12*5 + 12*8 + 10*13 + 6*21 + 3*34 + 1*55 = 615. (End)

MATHEMATICA

CoefficientList[Series[-(-1 - 4 x + 3 x^2 + 4 x^3)/(1 - 10 x^2 + 5 x^4), {x, 0, 30}], x] (* Vincenzo Librandi, Jun 17 2014 *)

LinearRecurrence[{0, 10, 0, -5}, {1, 4, 7, 36}, 30] (* Harvey P. Dale, Apr 07 2019 *)

PROG

(MAGMA) I:=[1, 4, 7, 36]; [n le 4 select I[n] else 10*Self(n-2)-5*Self(n-4): n in [1..40]]; // Vincenzo Librandi, Jun 17 2014

(PARI) Vec(-x*(-1-4*x+3*x^2+4*x^3)/(1-10*x^2+5*x^4) + O(x^30)) \\ Michel Marcus, Feb 06 2016

CROSSREFS

Cf. A038754, A084990, A189334 (bisection).

Sequence in context: A124624 A209480 A209339 * A302198 A279830 A245305

Adjacent sequences:  A212497 A212498 A212499 * A212501 A212502 A212503

KEYWORD

nonn,base,easy

AUTHOR

Vladimir Shevelev and Peter J. C. Moses, May 19 2012

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 14 01:33 EDT 2020. Contains 336473 sequences. (Running on oeis4.)