A010879 -id:A010879 - cp's OEIS frontend

A008959 Final digit of squares: a(n) = n^2 mod 10.

0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0, 1, 4, 9, 6, 5, 6, 9, 4, 1, 0

Offset: 0

N. J. A. Sloane, Mar 15 1996

a(m*n) = a(m)*a(n) mod 10; a(5*n+k) = a(5*n-k) for k <= 5*n. - Reinhard Zumkeller, Apr 24 2009
a(n) = n^6 mod 10. - Zerinvary Lajos, Nov 06 2009
a(n) = A002015(n) mod 10 = A174452(n) mod 10. - Reinhard Zumkeller, Mar 21 2010
Decimal expansion of 166285490/1111111111. - Alexander R. Povolotsky, Mar 09 2013

Harvey P. Dale, Table of n, a(n) for n = 0..1000
Index entries for sequences related to final digits of numbers
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,1).

Cf. A000290, A070431, A070435, A070438, A070442, A070452, A159852, A010879, A008960, A070514.

[0] cat [Intseq(n^2)[1]: n in [1..80]]; // Bruno Berselli, Feb 14 2013

[n^2 - 10*Floor(n^2/10): n in [0..80]]; // Vincenzo Librandi, Jun 16 2015

A008959:=n->(n^2 mod 10); seq(A008959(n), n=0..50); # Wesley Ivan Hurt, Jun 06 2014

Table[Mod[n^2,10],{n,0,200}] (* Vladimir Joseph Stephan Orlovsky, Apr 21 2011 *)
PowerMod[Range[0,80],2,10] (* or *) LinearRecurrence[{0,0,0,0,0,0,0,0,0,1},{0,1,4,9,6,5,6,9,4,1},120] (* Harvey P. Dale, Oct 16 2012 *)

a(n)=n^2%10 \\ Charles R Greathouse IV, Sep 24 2015

[power_mod(n,2,10) for n in range(0, 81)] # Zerinvary Lajos, Nov 06 2009

Periodic with period 10. - Franklin T. Adams-Watters, Mar 13 2006
a(n) = 4.5 - (1 + 5^(1/2))*cos(Pi*n/5) + (-1 - 3/5*5^(1/2))*cos(2*Pi*n/5) + (5^(1/2) - 1)*cos(3*Pi*n/5) + (-1 + 3/5*5^(1/2))*cos(4*Pi*n/5) - 0.5*(-1)^n. - Richard Choulet, Dec 12 2008
a(n) = A010879(A000290(n)). - Reinhard Zumkeller, Jan 04 2009
G.f.: (x^9+4*x^8+9*x^7+6*x^6+5*x^5+6*x^4+9*x^3+4*x^2+x)/(-x^10+1). - Colin Barker, Aug 14 2012
a(n) = n^2 - 10*floor(n^2/10). - Wesley Ivan Hurt, Jun 12 2013
a(n) = (n - 5*A002266(n + 2))^2 + 5*(5*A002266(n + 2) mod 2). - Wesley Ivan Hurt, Jun 06 2014
a(n) = A033569(n+3) mod 10. - Wesley Ivan Hurt, Dec 06 2014
a(n) = n^k mod 10; for k > 0 where k mod 4 = 2. - Doug Bell, Jun 15 2015

A002418 4-dimensional figurate numbers: a(n) = (5n-1)binomial(n+2,3)/4.

Original entry on oeis.org

0, 1, 9, 35, 95, 210, 406, 714, 1170, 1815, 2695, 3861, 5369, 7280, 9660, 12580, 16116, 20349, 25365, 31255, 38115, 46046, 55154, 65550, 77350, 90675, 105651, 122409, 141085, 161820, 184760, 210056, 237864, 268345, 301665, 337995

Offset: 0

N. J. A. Sloane

Partial sums of A002413.
Principal diagonal of the convolution array A213550, for n>0. - Clark Kimberling, Jun 17 2012
Convolution of A000027 with A000566. - Bruno Berselli, Dec 06 2012
Coefficients in the hypergeometric series identity 1 - 9*(x - 1)/(4*x + 1) + 35*(x - 1)*(x - 2)/((4*x + 1)*(4*x + 2)) - 95*(x - 1)*(x - 2)*(x - 3)/((4*x + 1)*(4*x + 2)*(4*x + 3)) + ... = 0, valid for Re(x) > 1. Cf. A000326 and A002412. Column 4 of A103450. - Peter Bala, Mar 14 2019

Albert H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 195.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992, arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Luis Verde-Star, A Matrix Approach to Generalized Delannoy and Schröder Arrays, J. Int. Seq., Vol. 24 (2021), Article 21.4.1.
Index entries for linear recurrences with constant coefficients, signature (5,-10,10,-5, 1).
Index to sequences related to pyramidal numbers.

Cf. A093562 ((5, 1) Pascal, column m=4).
Cf. A000332, A000566, A001622, A080852, A128064.
Cf. A220212 for a list of sequences produced by the convolution of the natural numbers with the k-gonal numbers.

List([0..40],n->(5*n-1)*Binomial(n+2,3)/4); # Muniru A Asiru, Mar 18 2019

[(5*n - 1)*Binomial(n + 2, 3)/4: n in [0..40]]; // Vincenzo Librandi, Oct 17 2012

/* A000027 convolved with A000566: */ A000566:=func; [&+[(n-i+1)*A000566(i): i in [0..n]]: n in [0..35]]; // Bruno Berselli, Dec 06 2012

Table[(5n-1) Binomial[n+2,3]/4,{n,0,40}] (* or *) LinearRecurrence[ {5,-10,10,-5,1},{0,1,9,35,95},40] (* Harvey P. Dale, Oct 16 2012 *)
CoefficientList[Series[x*(1 + 4*x)/(1 - x)^5, {x, 0, 40}], x] (* Vincenzo Librandi, Oct 17 2012 *)

a(n)=(5*n-1)*binomial(n+2,3)/4 \\ Charles R Greathouse IV, Sep 24 2015

[(5*n-1)*binomial(n+2,3)/4 for n in (0..40)] # G. C. Greubel, Jul 03 2019

G.f.: x*(1+4*x)/(1-x)^5. - Simon Plouffe in his 1992 dissertation.
Starting (1, 9, 35, 95, ...), = A128064 * A000332, (A000332 starting 1, 5, 15, 35, 70, ...), such that a(n) = n*C(n+3,4) - (n-1)*C(n+2,4). E.g., a(5) = 210 = 5*C(8,4) - 4*C(7,4) = 5*70 - 4*35. - Gary W. Adamson, Dec 28 2007
Unit digit, A010879(a(n)), is one of {0,1,9,5,6,4} [Eric Desbiaux] because a(n) mod 5 = 0,1,4,0,0, periodic with period 5. [Proof: A002413(n) mod 5 = 1,3,1,0,0 with period 5 and a(n) are the partial sums of A002413.] - R. J. Mathar, Mar 19 2008
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). - Harvey P. Dale, Oct 16 2012
a(n) = A080852(5,n-1). - R. J. Mathar, Jul 28 2016
a(n) = Sum_{i=0..n} (n-i) * Sum_{j=i..n} j. - J. M. Bergot, May 30 2017
E.g.f.: x*(24 + 84*x + 44*x^2 + 5*x^3)*exp(x)/4!. - G. C. Greubel, Jul 03 2019
Sum_{n>=1} 1/a(n) = (50*sqrt(5)*log(phi) + 125*log(5) - 50*sqrt(1+2/sqrt(5))*Pi - 26)/11, where phi is the golden ratio (A001622). - Amiram Eldar, Feb 11 2022

A106621 a(n) = numerator of n/(n+20).

Original entry on oeis.org

0, 1, 1, 3, 1, 1, 3, 7, 2, 9, 1, 11, 3, 13, 7, 3, 4, 17, 9, 19, 1, 21, 11, 23, 6, 5, 13, 27, 7, 29, 3, 31, 8, 33, 17, 7, 9, 37, 19, 39, 2, 41, 21, 43, 11, 9, 23, 47, 12, 49, 5, 51, 13, 53, 27, 11, 14, 57, 29, 59, 3, 61, 31, 63, 16, 13, 33, 67, 17, 69, 7, 71, 18, 73, 37, 15, 19, 77, 39, 79

Offset: 0

N. J. A. Sloane, May 15 2005

Contains as subsequences A026741, A017281, A017305, A005408, A017353, and A017377. - Luce ETIENNE, Nov 04 2018

Multiplicative and also a strong divisibility sequence: gcd(a(n),a(m)) = a(gcd(n,m)) for n, m >= 1. - Peter Bala, Feb 24 2019

Vincenzo Librandi, Table of n, a(n) for n = 0..2000
Peter Bala, A note on the sequence of numerators of a rational function, 2019.
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1).

Cf. A000010, A010879.

Cf. Sequences given by the formula numerator(n/(n + k)): A026741 (k = 2), A051176 (k = 3), A060819 (k = 4), A060791 (k = 5), A060789 (k = 6), A106608 thru A106612 (k = 7 thru 11), A051724 (k = 12), A106614 thru A106620 (k = 13 thru 19).

List([0..80],n->NumeratorRat(n/(n+20))); # Muniru A Asiru, Feb 19 2019

[Numerator(n/(n+20)): n in [0..100]]; // Vincenzo Librandi, Mar 06 2018

seq(numer(n/(n+20)),n=0..80); # Muniru A Asiru, Feb 19 2019

f[n_]:=Numerator[n/(n+20)];Array[f,100,0] (* Vladimir Joseph Stephan Orlovsky, Feb 17 2011 *)

a(n) = numerator(n/(n+20)); \\ Michel Marcus, Mar 07 2018

[lcm(20,n)/20 for n in range(0, 80)] # Zerinvary Lajos, Jun 12 2009

a(n) = lcm(20, n)/20. - Zerinvary Lajos, Jun 12 2009
a(n) = n/gcd(n, 20). - Andrew Howroyd, Jul 25 2018
From Luce ETIENNE, Nov 04 2018: (Start)
a(n) = 9*a(n-20) - 36*a(n-40) + 84*a(n-60) - 126*a(n-80) + 126*a(n-100) - 84*a(n-120) + 36*a(n-140) - 9*a(n-160) + a(n-180).
a(n) = (5*(119*m^9 - 4923*m^8 + 86250*m^7 - 832230*m^6 + 4807887*m^5 - 16882299*m^4 + 34770400*m^3 - 37855620m^2 + 16581744*m + 54432)*floor(n/10) + 72*m*(3*m^8 - 120*m^7 + 2030*m^6 - 18900*m^5 + 105329*m^4 - 356580*m^3 + 706220*m^2 - 733200*m + 300258) + ((19*m^9 - 855*m^8 + 15810*m^7 - 154350*m^6 + 849387*m^5 - 2597175*m^4 + 4037840*m^3 - 2600100*m^2 + 540144*m - 90720)*floor(n/10) - 72*m*(m^7 - 35*m^6 + 490*m^5 - 3500*m^4 + 13489*m^3 - 27335*m^2 + 26340*m - 9450))*(-1)^floor(n/10))/362880 where m = (n mod 10). (End)
From Peter Bala, Feb 24 2019: (Start)
a(n) = n/gcd(n,20) is a quasi-polynomial in n since gcd(n,20) is a purely periodic sequence of period 20.
O.g.f.: F(x) - F(x^2) - F(x^4) - 4*F(x^5) + 4*F(x^10) + 4*F(x^20), where F(x) = x/(1 - x)^2.
O.g.f. for reciprocals: Sum_{n >= 1} x^n/a(n) = Sum_{d divides 20} (phi(d)/d) * log(1/(1 - x^d)) = log(1/(1 - x)) + (1/2)*log(1/(1 - x^2)) + (2/4)*log(1/(1 - x^4)) + (4/5)*log(1/(1 - x^5)) + (4/10)*log(1/(1 - x^10)) + (8/20)*log(1/(1 - x^20)), where phi(n) denotes the Euler totient function A000010. (End)
From Amiram Eldar, Nov 25 2022: (Start)
Multiplicative with a(2^e) = 2^max(0, e-2), a(5^e) = 5^max(0,e-1), and a(p^e) = p^e otherwise.
Dirichlet g.f.: zeta(s-1)*(1 - 1/2^s - 1/4^s - 4/5^s + 4/10^s + 4/20^s).
Sum_{k=1..n} a(k) ~ (231/800) * n^2. (End)

Keyword:mult added by Andrew Howroyd, Jul 25 2018

A023961 First digit after decimal point of square root of n.

Original entry on oeis.org

0, 4, 7, 0, 2, 4, 6, 8, 0, 1, 3, 4, 6, 7, 8, 0, 1, 2, 3, 4, 5, 6, 7, 8, 0, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 0, 1, 2, 3, 4, 4, 5, 6, 7, 7, 8, 9, 0, 0, 1, 2, 2, 3, 4, 4, 5, 6, 6, 7, 8, 8, 9, 0, 0, 1, 1, 2, 3, 3, 4, 4, 5, 6, 6, 7, 7, 8, 8, 9, 0, 0, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 0, 0, 0, 1, 1, 2

Offset: 1

N. J. A. Sloane, Olivier Gérard

When n is a square, a(n) is equal to 0, but the converse is not true, see A034096. - Michel Marcus, Sep 21 2015

			sqrt(1) = 1.00000000... hence a(1) = 0.
sqrt(2) = 1.41421356... hence a(2) = 4.
sqrt(3) = 1.73205080... hence a(3) = 7.
sqrt(4) = 2.00000000... hence a(4) = 0.
Note that 26 = 2 * 13 and sqrt(26) = 5.09901951... so a(26) = 0 even though 26 is not a perfect square.

Nathaniel Johnston, Table of n, a(n) for n = 1..10000

Cf. A111862, A111850, A111851, A111852, A111853, A111854, A111855, A111856, A111857, A111858, A111859.

A023961 := proc(n) return floor(10*sqrt(n)) mod 10: end: seq(A023961(n),n=1..100); # Nathaniel Johnston, May 04 2011

Array[Function[n, RealDigits[N[Power[n, 1/2], 10], 10] // (#[[1, #[[2]] + 1]]) &], 110]
fd[n_] := Module[{rd = RealDigits[Sqrt[n], 10, 10]}, First[rd][[Last[rd] + 1]]]; Array[fd, 90] (* Harvey P. Dale, Jan 16 2014 *)

a(n) = floor(10*sqrt(n)) % 10; \\ Michel Marcus, Sep 21 2015

a(n) = A010879(A000196(100*n)). - Robert Israel, Jul 30 2015

A131242 Partial sums of A059995: a(n) = sum_{k=0..n} floor(k/10).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16, 18, 20, 22, 24, 26, 28, 30, 33, 36, 39, 42, 45, 48, 51, 54, 57, 60, 64, 68, 72, 76, 80, 84, 88, 92, 96, 100, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 156, 162, 168, 174, 180, 186, 192, 198

Offset: 0

Hieronymus Fischer, Jun 21 2007

Complementary with A130488 regarding triangular numbers, in that A130488(n)+10*a(n)=n(n+1)/2=A000217(n).

			As square array :
    0,   0,   0,   0,   0,   0,   0,   0,   0,    0
    1,   2,   3,   4,   5,   6,   7,   8,   9,   10
   12,  14,  16,  18,  20,  22,  24,  26,  28,   30
   33,  36,  39,  42,  45,  48,  51,  54,  57,   60
   64,  68,  72,  76,  80,  84,  88,  92,  96,  100
  105, 110, 115, 120, 125, 130, 135, 140, 145,  150
  156, 162, 168, 174, 180, 186, 192, 198, 204,  210
... - _Philippe Deléham_, Mar 27 2013

G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (2,-1,0,0,0,0,0,0,0,1,-2,1).

Cf. A008728, A059995, A010879, A002266, A130488, A000217, A002620, A130518, A130519, A130520, A174709, A174738, A118729, A218470.

Table[(1/2)*Floor[n/10]*(2*n - 8 - 10*Floor[n/10]), {n,0,50}] (* G. C. Greubel, Dec 13 2016 *)
Accumulate[Table[FromDigits[Most[IntegerDigits[n]]],{n,0,110}]] (* or *) LinearRecurrence[{2,-1,0,0,0,0,0,0,0,1,-2,1},{0,0,0,0,0,0,0,0,0,0,1,2},120] (* Harvey P. Dale, Apr 06 2017 *)

for(n=0,50, print1((1/2)*floor(n/10)*(2n-8-10*floor(n/10)), ", ")) \\ G. C. Greubel, Dec 13 2016

a(n)=my(k=n\10); k*(n-5*k-4) \\ Charles R Greathouse IV, Dec 13 2016

a(n) = (1/2)*floor(n/10)*(2n-8-10*floor(n/10)).
a(n) = A059995(n)*(2n-8-10*A059995(n))/2.
a(n) = (1/2)*A059995(n)*(n-8+A010879(n)).
a(n) = (n-A010879(n))*(n+A010879(n)-8)/20.
G.f.: x^10/((1-x^10)(1-x)^2).
From Philippe Deléham, Mar 27 2013: (Start)
a(10n)   = A051624(n).
a(10n+1) = A135706(n).
a(10n+2) = A147874(n+1).
a(10n+3) = 2*A005476(n).
a(10n+4) = A033429(n).
a(10n+5) = A202803(n).
a(10n+6) = A168668(n).
a(10n+7) = 2*A147875(n).
a(10n+8) = A135705(n).
a(10n+9) = A124080(n). (End)
a(n) = A008728(n-10) for n>= 10. - Georg Fischer, Nov 03 2018

A168184 Characteristic function of numbers that are not multiples of 10.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1

Offset: 0

Reinhard Zumkeller, Nov 30 2009

Reinhard Zumkeller, Table of n, a(n) for n = 0..1000
Index entries for characteristic functions
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,1).

Cf. A168185, A145568, A168182, A168181, A109720, A097325, A011558, A166486, A011655, A000035, A010690, A033442.

a168184 = (1 -) . (0 ^) . (`mod` 10)
a168184_list = cycle [0,1,1,1,1,1,1,1,1,1]
-- Reinhard Zumkeller, Oct 10 2012

Table[If[Mod[n,10]==0,0,1],{n,0,110}] (* or *) PadRight[{},110,{0,1,1,1,1,1,1,1,1,1}] (* Harvey P. Dale, Jun 03 2023 *)

a(n)=n%10>0 \\ Charles R Greathouse IV, Sep 24 2015

a(n+10) = a(n);
a(n) = A000007(A010879(n));
a(A067251(n)) = 1; a(A008592(n)) = 0;
not the same as A168046: a(n)=A168046 for n<=100;
A033442(n) = Sum_{k=0..n} a(k)*(n-k).
Dirichlet g.f.: (1-1/10^s)*zeta(s). - R. J. Mathar, Feb 19 2011
For the general case: the characteristic function of numbers that are not multiples of m is a(n)=floor((n-1)/m)-floor(n/m)+1, m,n > 0. - Boris Putievskiy, May 08 2013

A007584 9-gonal (or enneagonal) pyramidal numbers: a(n) = n(n+1)(7*n-4)/6.

Original entry on oeis.org

0, 1, 10, 34, 80, 155, 266, 420, 624, 885, 1210, 1606, 2080, 2639, 3290, 4040, 4896, 5865, 6954, 8170, 9520, 11011, 12650, 14444, 16400, 18525, 20826, 23310, 25984, 28855, 31930, 35216, 38720, 42449, 46410, 50610, 55056, 59755, 64714, 69940, 75440, 81221

Offset: 0

N. J. A. Sloane, R. K. Guy

For n > 1, the digital roots of this sequence A010888(A007584(n)) form the purely periodic 27-cycle 1, 1, 7, 8, 2, 5, 6, 3, 3, 4, 4, 1, 2, 5, 8, 9, 6, 6, 7, 7, 4, 5, 8, 2, 3, 9, 9. For n > 1, the units digits of this sequence A010879(A007584(n)) form the purely periodic 20-cycle 1, 0, 4, 0, 5, 6, 0, 4, 5, 0, 6, 0, 9, 0, 0, 6, 5, 4, 0, 0. - Ant King, Oct 30 2012

Partial sums of A001106. - Joerg Arndt, Jun 10 2013

A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index to sequences related to pyramidal numbers
Index entries for linear recurrences with constant coefficients, signature (4, -6, 4, -1).

Cf. A093564 ((7, 1) Pascal, column m=3).

Cf. similar sequences listed in A237616.

I:=[0, 1, 10, 34, 80]; [n le 5 select I[n] else 4*Self(n-1)-6*Self(n-2)+4*Self(n-3)-Self(n-4): n in [1..50]]; // Vincenzo Librandi, Jun 10 2013

a:=n->sum((n+j)^2-(n+j), j=0..n): seq(a(n)/2, n=0..30); # Zerinvary Lajos, May 26 2008

Table[n*(n+1)(7n-4)/6, {n, 0,100}] (* Vladimir Joseph Stephan Orlovsky, Jun 25 2009 *)
LinearRecurrence[{4,-6,4,-1},{1,10,34,80},30] (* Ant King, Oct 27 2012 *)
CoefficientList[Series[x (1 + 6 x) / (1 - x)^4, {x, 0, 50}], x] (* Vincenzo Librandi, Jun 10 2013 *)

A007584[n]:=n*(n+1)*(7*n-4)/6$
makelist(A007584[n],n,0,30); /* Martin Ettl, Oct 29 2012 */

a(n) = n*(n+1)*(7*n-4)/6; \\ Michel Marcus, Mar 04 2014

a(n) = (7*n-4)*binomial(n+1, 2)/3.
G.f.: x*(1+6*x)/(1-x)^4.
From Ant King, Oct 27 2012: (Start)
a(n) = a(n-1) + n*(7*n-5)/2.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 7.
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4).
a(n) = (n+1)*(2*A001106(n)+n)/6.
a(n) = A000292(n) + 6*A000292(n-1).
a(n) = A002414(n) + A000292(n-1).
a(n) = A000217(n) + 7*A000292(n-1).
a(n) = binomial(n+2,3) + 6*binomial(n+1,3). (End)
a(n) = Sum_{i = 0..n-1} (n-i)*(7*i+1) for n>0. - Bruno Berselli, Feb 10 2014
a(n) = A080851(7,n-1). - R. J. Mathar, Jul 28 2016
E.g.f.: (x/6)*(6 + 24*x + 7*x^2)*exp(x). - G. C. Greubel, Oct 29 2017

A007585 10-gonal (or decagonal) pyramidal numbers: a(n) = n(n + 1)(8*n - 5)/6.

Original entry on oeis.org

0, 1, 11, 38, 90, 175, 301, 476, 708, 1005, 1375, 1826, 2366, 3003, 3745, 4600, 5576, 6681, 7923, 9310, 10850, 12551, 14421, 16468, 18700, 21125, 23751, 26586, 29638, 32915, 36425, 40176, 44176, 48433, 52955, 57750

Offset: 0

N. J. A. Sloane, R. K. Guy

Binomial transform of [1, 10, 17, 8, 0, 0, 0, ...] = (1, 11, 38, 90, ...). - Gary W. Adamson, Mar 18 2009
This sequence is related to A000384 by a(n) = n*A000384(n) - Sum_{i=0..n-1} A000384(i) and this is the case d=4 in the identity n*(n*(d*n-d+2)/2) - Sum_{k=0..n-1} k*(d*k-d+2)/2 = n*(n+1)*(2*d*n - 2*d + 3)/6. - Bruno Berselli, Apr 21 2010
For n>0, (a(n)) is the principal diagonal of the convolution array A213750. - Clark Kimberling, Jun 20 2012
From Ant King, Oct 30 2012: (Start)
The partial sums of the figurate decagonal numbers A001107.
For n>1, the digital roots of this sequence A010888(A007585(n)) form the purely periodic 27-cycle {1,2,2,9,4,4,8,6,6,7,8,8,6,1,1,5,3,3,4,5,5,3,7,7,2,9,9}.
For n>1, the units’ digits of this sequence A010879(A007585(n)) form the purely periodic 20-cycle {1,1,8,0,5,1,6,8,5,5,6,6,3,5,0,6,1,3,0,0}.
(End)

A. H. Beiler, Recreations in the Theory of Numbers, Dover, NY, 1964, p. 194.
E. Deza and M. M. Deza, Figurate numbers, World Scientific Publishing (2012), page 93.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Harvey P. Dale, Table of n, a(n) for n = 0..1000
B. Berselli, A description of the recursive method in Comments lines: website Matem@ticamente (in Italian).
Index to sequences related to pyramidal numbers
Index entries for linear recurrences with constant coefficients, signature (4,-6,4,-1)

Cf. A000384.
Cf. A093565 ((8, 1) Pascal, column m=3). Partial sums of A001107.
Cf. similar sequences listed in A237616.

List([0..40], n-> n*(n+1)*(8*n-5)/6); # G. C. Greubel, Aug 30 2019

[n*(n+1)*(8*n-5)/6: n in [0..40]]; // G. C. Greubel, Aug 30 2019

seq(n*(n+1)*(8*n-5)/6, n=0..40); # G. C. Greubel, Aug 30 2019

Table[n(n+1)(8n-5)/6, {n,0,40}] (* Vladimir Joseph Stephan Orlovsky, Apr 18 2011 *)
LinearRecurrence[{4,-6,4,-1},{0,1,11,38},40] (* Harvey P. Dale, Dec 20 2011 *)

a(n)=(8*n^3+3*n^2-5*n)/6 \\ Charles R Greathouse IV, Sep 24 2015

[n*(n+1)*(8*n-5)/6 for n in (0..40)] # G. C. Greubel, Aug 30 2019

a(n) = (8*n-5)*binomial(n+1, 2)/3.
G.f.: x*(1+7*x)/(1-x)^4.
a(n) = (8*n^3 + 3*n^2 - 5*n)/6. - Vincenzo Librandi, Aug 01 2010
a(0)=0, a(1)=1, a(2)=11, a(3)=38, a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4). - Harvey P. Dale, Dec 20 2011
From Ant King, Oct 30 2012: (Start)
a(n) = a(n-1) + n*(4*n-3).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) + 8.
a(n) = (n+1)*(2*A001107(n) + n)/6.
a(n) = A000292(n) + 7*A000292(n-1).
a(n) = A007584(n) + A000292(n-1).
a(n) = A000217(n) + 8*A000292(n-1).
a(n) = binomial(n+2,3) + 7*binomial(n+1,3).
Sum_{n>=1} 1/a(n) = 6*(4*pi*(sqrt(2)-1) + 4*(8-sqrt(2))*log(2) + 8*sqrt(2)*log(2-sqrt(2))-5)/65 =  1.145932345...
(End)
a(n) = Sum_{i=0..n-1} (n-i)*(8*i+1), with a(0)=0. - Bruno Berselli, Feb 10 2014
E.g.f.: x*(6 + 27*x + 8*x^2)*exp(x)/6. - Ilya Gutkovskiy, May 12 2017

A008960 Final digit of cubes: n^3 mod 10.

Original entry on oeis.org

0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5, 6, 3, 2, 9, 0, 1, 8, 7, 4, 5

Offset: 0

N. J. A. Sloane

Decimal expansion of 208284810/1111111111. - Alexander R. Povolotsky, Mar 08 2013

Vincenzo Librandi, Table of n, a(n) for n = 0..5000
Index entries for sequences related to final digits of numbers
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,0,0,0,0,0,1).

Cf. A167176.

Cf. A010879, A008959, A070514. - Doug Bell, Jun 15 2015

[n^3 mod 10: n in [0..80]]; // Vincenzo Librandi, Mar 26 2013

Table[Mod[n^3,10],{n,0,200}] (* Vladimir Joseph Stephan Orlovsky, Apr 23 2011 *)
LinearRecurrence[{0, 0, 0, 0, 0, 0, 0, 0, 0, 1},{0, 1, 8, 7, 4, 5, 6, 3, 2, 9},81] (* Ray Chandler, Aug 26 2015 *)

a(n)=n^3%10 \\ Charles R Greathouse IV, Mar 08 2013

concat(0, Vec(x*(1+8*x+7*x^2+4*x^3+5*x^4+6*x^5+3*x^6+2*x^7+9*x^8) / ((1-x)*(1+x)*(1-x+x^2-x^3+x^4)*(1+x+x^2+x^3+x^4)) + O(x^100))) \\ Colin Barker, Nov 30 2015

[power_mod(n,3,10 ) for n in range(0, 81)] # Zerinvary Lajos, Oct 29 2009

Periodic with period 10. - Franklin T. Adams-Watters, Mar 13 2006
a(n) = 4.5 -cos(Pi*n/5) +(1/2*(-(5-5^(1/2))^(1/2) +(5+5^(1/2))^(1/2))*2^(1/2))*sin(Pi*n/5) -cos(2*Pi*n/5) +(-1/10*(-(5-5^(1/2))^(1/2)+3*(5+5^(1/2))^(1/2))*2^(1/2))*sin(2*Pi*n/5) -cos(3*Pi*n/5) +(-1/2*((5-5^(1/2))^(1/2) +(5+5^(1/2))^(1/2))*2^(1/2))*sin(3*Pi*n/5) -cos(4*Pi*n/5) +( -1/10*(3*(5-5^(1/2))^(1/2) +(5 +5^(1/2))^(1/2))*2^(1/2))*sin(4*Pi*n/5) -0.5*(-1)^n. - Richard Choulet, Dec 12 2008
a(n) = n^k mod 10; for k > 0 where k mod 4 = 3. - Doug Bell, Jun 15 2015
G.f.: x*(1+8*x+7*x^2+4*x^3+5*x^4+6*x^5+3*x^6+2*x^7+9*x^8) / ((1-x)*(1+x)*(1-x+x^2-x^3+x^4)*(1+x+x^2+x^3+x^4)). - Colin Barker, Nov 30 2015

A010077 a(n) = sum of digits of a(n-1) + sum of digits of a(n-2); a(0) = 0, a(1) = 1.

Original entry on oeis.org

0, 1, 1, 2, 3, 5, 8, 13, 12, 7, 10, 8, 9, 17, 17, 16, 15, 13, 10, 5, 6, 11, 8, 10, 9, 10, 10, 2, 3, 5, 8, 13, 12, 7, 10, 8, 9, 17, 17, 16, 15, 13, 10, 5, 6, 11, 8, 10, 9, 10, 10, 2, 3, 5, 8, 13, 12, 7, 10, 8, 9, 17, 17, 16, 15, 13, 10, 5, 6, 11, 8, 10, 9, 10, 10

Offset: 0

N. J. A. Sloane, Leonid Broukhis

The digital sum analog (in base 10) of the Fibonacci recurrence. - Hieronymus Fischer, Jun 27 2007
a(n) and Fibonacci(n) = A000045(n) are congruent modulo 9 which implies that (a(n) mod 9) is equal to (Fibonacci(n) mod 9) = A007887(n). Thus (a(n) mod 9) is periodic with the Pisano period A001175(9)=24. - Hieronymus Fischer, Jun 27 2007
a(n) == A004090(n) (mod 9) (A004090(n) = digital sum of Fibonacci(n)). - Hieronymus Fischer, Jun 27 2007
For general bases p > 2, we have the inequality 2 <= a(n) <= 2p-3 (for n > 2). Actually, a(n) <= 17 = A131319(10) for the base p=10. - Hieronymus Fischer, Jun 27 2007

David A. Corneth, Table of n, a(n) for n = 0..9999
Index entries for Colombian or self numbers and related sequences

Cf. A007953, A007612, A065076.
Cf. A000045, A010073, A010074, A010075, A010076.
Cf. A131294, A131295, A131296, A131297, A131318, A131319, A131320.

a[0] = 0; a[1] = 1; a[n_] := a[n] = Apply[ Plus, IntegerDigits[ a[n - 1] ]] + Apply[ Plus, IntegerDigits[ a[n - 2] ]]; Table[ a[n], {n, 0, 100} ]
nxt[{a_,b_}]:={b, Total[IntegerDigits[a]]+Total[IntegerDigits[b]]}; NestList[ nxt,{0,1},80][[All,1]] (* Harvey P. Dale, Apr 15 2018 *)

first(n) = {n = max(n, 2); my(res = vector(n)); res[2] = 1; for(i = 3, n, res[i] = sumdigits(res[i-1]) + sumdigits(res[i-2]) ); res } \\ David A. Corneth, May 26 2021

Periodic from n=3 with period 24. - Franklin T. Adams-Watters, Mar 13 2006
a(n) = A030132(n-4) + A030132(n-3) for n>3. - Reinhard Zumkeller, Jul 04 2007
a(n) = a(n-1) + a(n-2) - 9*(floor(a(n-1)/10) + floor(a(n-2)/10)). - Hieronymus Fischer, Jun 27 2007
a(n) = floor(a(n-1)/10) + floor(a(n-2)/10) + (a(n-1) mod 10) + (a(n-2) mod 10). - Hieronymus Fischer, Jun 27 2007
a(n) = A059995(a(n-1)) + A059995(a(n-2)) + A010879(a(n-1)) + A010879(a(n-2)). - Hieronymus Fischer, Jun 27 2007
a(n) = Fibonacci(n) - 9*Sum_{k=2..n-1} Fibonacci(n-k+1)*floor(a(k)/10) where Fibonacci(n) = A000045(n). - Hieronymus Fischer, Jun 27 2007

cp's OEIS Frontend

A008959 Final digit of squares: a(n) = n^2 mod 10.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Maple

Mathematica

PARI

Sage

Formula

A002418 4-dimensional figurate numbers: a(n) = (5*n-1)*binomial(n+2,3)/4.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Magma

Magma

Mathematica

PARI

Sage

Formula

A106621 a(n) = numerator of n/(n+20).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A023961 First digit after decimal point of square root of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A131242 Partial sums of A059995: a(n) = sum_{k=0..n} floor(k/10).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

A168184 Characteristic function of numbers that are not multiples of 10.

Views

Author

Keywords

Links

Crossrefs

A002418 4-dimensional figurate numbers: a(n) = (5n-1)binomial(n+2,3)/4.

A007584 9-gonal (or enneagonal) pyramidal numbers: a(n) = n(n+1)(7*n-4)/6.

A007585 10-gonal (or decagonal) pyramidal numbers: a(n) = n(n + 1)(8*n - 5)/6.