A061167 -id:A061167 - cp's OEIS frontend

A052530 a(n) = 4*a(n-1) - a(n-2), with a(0) = 0, a(1) = 2.

0, 2, 8, 30, 112, 418, 1560, 5822, 21728, 81090, 302632, 1129438, 4215120, 15731042, 58709048, 219105150, 817711552, 3051741058, 11389252680, 42505269662, 158631825968, 592022034210, 2209456310872, 8245803209278, 30773756526240

Offset: 0

encyclopedia(AT)pommard.inria.fr, Jan 25 2000

a(n-1) and a(n+1) are the solutions for c if b = a(n) in (b^2 + c^2)/(b*c + 1) = 4 and there are no other pairs of solutions apart from consecutive pairs of terms in this sequence. Cf. A061167. - Henry Bottomley, Apr 18 2001
a(n)^2 for n >= 1 gives solutions to A007913(3*x+4) = A007913(x). - Benoit Cloitre, Apr 07 2002
For all terms k of the sequence, 3*k^2 + 4 is a perfect square. Limit_{n->oo} a(n)/a(n-1) = 2 + sqrt(3). - Gregory V. Richardson, Oct 06 2002
a(n) = the number of compositions of the integer 2*n into even parts, where each part 2*i comes in 2*i colors. (Dedrickson, Theorem 3.2.6) An example is given below. Cf. A052529, A095263. - Peter Bala, Sep 17 2013
Except for an initial 1, this is the p-INVERT of (1, 1, 1, 1, 1, ...) for p(S) = 1 - 2*S - 2*S^2; see A291000.  - Clark Kimberling, Aug 24 2017
a(n+1) is the number of spanning trees of the graph P_n, where P_n is a 2 X n grid with two additional vertices, u and v, where u is adjacent to (1,1) and (2,1), and v is adjacent to (1,n) and (2,n). - Kevin Long, May 04 2018
a(n) is also the output of Tesler's formula for the number of perfect matchings of an m X n Mobius band where m is even and n is odd, specialized to m=2. (The twist is on the length-n side.) - Sarah-Marie Belcastro, Feb 15 2022
In general, values of x and y which satisfy (x^2 + y^2)/(x*y + 1) = k^2 are any two adjacent terms of a second-order recurrence with initial terms 0 and k and signature (k^2,-1). This can also be expressed as a first-order recurrence a(n+1) = (k^2*a(n) + sqrt((k^4-4)*a(n)^2 + 4*k^2))/2, n > 1. - Gary Detlefs, Feb 27 2024

			Colored compositions. a(2) = 8: There are two compositions of 4 into even parts, namely 4 and 2 + 2. Using primes to indicate the coloring of parts, the 8 colored compositions are 4, 4', 4'', 4''', 2 + 2, 2 + 2', 2' + 2 and 2' + 2'. - _Peter Bala_, Sep 17 2013

T. D. Noe, Table of n, a(n) for n = 0..200
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
Hacène Belbachir, Soumeya Merwa Tebtoub, and László Németh, Ellipse Chains and Associated Sequences, J. Int. Seq., Vol. 23 (2020), Article 20.8.5.
Daniel Birmajer, Juan B. Gil, and Michael D. Weiner, (an + b)-color compositions, arXiv:1707.07798 [math.CO], 2017.
Niccolò Castronuovo, On the number of fixed points of the map gamma, arXiv:2102.02739 [math.NT], 2021. Mentions this sequence.
C. R. Dedrickson III, Compositions, Bijections, and Enumerations Thesis, 2012, Jack N. Averitt College of Graduate Studies, Georgia Southern University.
J.-P. Ehrmann et al., Problem POLYA002, Integer pairs (x,y) for which (x^2+y^2)/(1+pxy) is an integer.
F. Goebel and A. A. Jagers, On a conjecture of Tutte concerning minimal tree numbers, J. Combin. Theory Ser. B 26 (1979), no. 3, 346-348. MR0535948 (80m:05064). [From _N. J. A. Sloane_, Feb 20 2012]
A. F. Horadam, Basic Properties of a Certain Generalized Sequence of Numbers, Fibonacci Quarterly, Vol. 3, No. 3, 1965, pp. 161-176.
INRIA Algorithms Project, Encyclopedia of Combinatorial Structures 460
N. J. A. Sloane, Transforms
G. Tesler, Matchings in graphs on non-orientable surfaces, Journal of Combinatorial Theory B, 78(2000), 198-231.
Index entries for linear recurrences with constant coefficients, signature (4,-1).

Cf. A001353, A007913, A003699, A217233.

Cf. A052529, A095263.

a052530 n = a052530_list !! n
a052530_list =
   0 : 2 : zipWith (-) (map (* 4) $ tail a052530_list) a052530_list
-- Reinhard Zumkeller, Sep 29 2011

I:=[0,2]; [n le 2 select I[n] else 4*Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Feb 25 2019

spec := [S,{S=Sequence(Prod(Union(Z,Z),Sequence(Z),Sequence(Z)))},unlabeled]: seq(combstruct[count](spec, size=n), n=0..20);
s := sqrt(3): a := n -> ((2-s)^n-(s+2)^n)/(s*(s-2)*(s+2)):
seq(simplify(a(n)), n=0..24); # Peter Luschny, Apr 28 2020

p=1; c=2; a[0]=0; a[1]=c; a[n_]:=a[n]=p*c^2*a[n-1]-a[n-2]; Table[a[n], {n, 0, 20}]
NestList[2 # + Sqrt[4 + 3 #^2]&, 0, 200] (* Zak Seidov, Mar 31 2011 *)
LinearRecurrence[{4, -1}, {0, 2}, 25] (* T. D. Noe, Jan 09 2012 *)
CoefficientList[Series[2x/(1-4x+x^2),{x,0,30}],x] (* Harvey P. Dale, May 31 2023 *)

{ polya002(p,c,m) = local(v,w,j,a); w=0; print1(w,", "); v=c; print1(v,", "); j=1; while(j<=m,a=p*c^2*v-w; print1(a,", "); w=v; v=a; j++) };
polya002(1,2,25)

my(x='x+O('x^30)); concat([0], Vec(2*x/(1-4*x+x^2))) \\ G. C. Greubel, Feb 25 2019

first(n) = n = max(n, 2); my(res = vector(n)); res[1] = 0; res[2] = 2; for(i = 3, n, res[i] = 4 * res[i-1] - res[i-2]); res \\ David A. Corneth, Apr 28 2020

(2*x/(1-4*x+x^2)).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Feb 25 2019

G.f.: 2*x/(1 - 4*x + x^2).
Invert transform of even numbers: a(n) = 2*Sum_{k=1..n} k*a(n-k). - Vladeta Jovovic, Apr 27 2001
From Gregory V. Richardson, Oct 06 2002: (Start)
a(n) = Sum_{alpha} -(1/3)*(-1 + 2*alpha)*alpha^(-1 - n), alpha = root of (1 - 4*Z + Z^2).
a(n) = (((2+sqrt(3))^(n+1) - (2-sqrt(3))^(n+1)) - ((2+sqrt(3))^n - (2-sqrt(3))^n) + ((2+sqrt(3))^(n-1) - (2-sqrt(3))^(n-1)))/(3*sqrt(3)). (End)
a(n) = A071954(n) - 2. - N. J. A. Sloane, Feb 20 2005
a(n) = (2*sinh(2n*arcsinh(1/sqrt(2))))/sqrt(3). - Herbert Kociemba, Apr 24 2008
a(n) = 2*A001353(n). - R. J. Mathar, Oct 26 2009
a(n) = ((3 - 2*sqrt(3))/3)*(2 - sqrt(3))^(n - 1) + ((3 + 2*sqrt(3))/3)*(2 + sqrt(3))^(n - 1). - Vincenzo Librandi, Nov 20 2010
a(n) = floor((2 + sqrt(3))^n/sqrt(3)). - Zak Seidov, Mar 31 2011
a(n) = ((2 + sqrt(3))^n - (2 - sqrt(3))^n)/sqrt(3). (See Horadam for construction.) - Johannes Boot, Jan 08 2012
a(n) = A217233(n) + A217233(n-1) with A217233(-1) = -1. - Bruno Berselli, Oct 01 2012
a(n) = A001835(n+1) - A001835(n). - Kevin Long, May 04 2018
E.g.f.: (exp((2 + sqrt(3))*x) - exp((2 - sqrt(3))*x))/sqrt(3). - Franck Maminirina Ramaharo, Nov 12 2018
a(n+1) = 2*a(n) + sqrt(3*a(n)^2 + 4), n > 1. - Gary Detlefs, Feb 27 2024

More terms from James Sellers, Jun 06 2000
Edited by N. J. A. Sloane, Nov 11 2006
a(0) changed to 0 and entry revised accordingly by Max Alekseyev, Nov 15 2007
Signs in definition corrected by John W. Layman, Nov 20 2007

A143324 Table T(n,k) by antidiagonals. T(n,k) is the number of length n primitive (=aperiodic or period n) k-ary words (n,k >= 1).

Original entry on oeis.org

1, 2, 0, 3, 2, 0, 4, 6, 6, 0, 5, 12, 24, 12, 0, 6, 20, 60, 72, 30, 0, 7, 30, 120, 240, 240, 54, 0, 8, 42, 210, 600, 1020, 696, 126, 0, 9, 56, 336, 1260, 3120, 4020, 2184, 240, 0, 10, 72, 504, 2352, 7770, 15480, 16380, 6480, 504, 0, 11, 90, 720, 4032, 16800, 46410, 78120, 65280, 19656, 990, 0

Offset: 1

Alois P. Heinz, Aug 07 2008

Column k is Dirichlet convolution of mu(n) with k^n.

The coefficients of the polynomial of row n are given by the n-th row of triangle A054525; for example row 4 has polynomial -k^2+k^4.

			T(2,3)=6, because there are 6 primitive words of length 2 over 3-letter alphabet {a,b,c}: ab, ac, ba, bc, ca, cb; note that the non-primitive words aa, bb and cc don't belong to the list; secondly note that the words in the list need not be Lyndon words, for example ba can be derived from ab by a cyclic rotation of the positions.
Table begins:
  1,  2,   3,    4,    5, ...
  0,  2,   6,   12,   20, ...
  0,  6,  24,   60,  120, ...
  0, 12,  72,  240,  600, ...
  0, 30, 240, 1020, 3120, ...

Alois P. Heinz, Antidiagonals n = 1..141, flattened
C. J. Smyth, A coloring proof of a generalisation of Fermat's little theorem, Amer. Math. Monthly 93, No. 6 (1986), pp. 469-471.
Lucas Vieira Barbosa, Nonclassical Temporal Correlations Under Finite-Memory Constraints, Doctoral Thesis, Univ. Wien (Austria 2024). See p. 35.
Index entries for sequences related to Lyndon words

Columns k=1-10 give: A000007(n-1), A027375, A054718, A054719, A054720, A054721, A218124, A218125, A218126, A218127.
Rows n=1-10 give: A000027, A002378(k-1), A007531(k+1), A047928(k+1), A061167, A218130, A133499, A218131, A218132, A218133.
Main diagonal gives A252764.
Cf. A074650, A143325, A008683, A054525.

with(numtheory): f0:= proc(n) option remember; unapply(k^n-add(f0(d)(k), d=divisors(n)minus{n}), k) end; T:= (n,k)-> f0(n)(k); seq(seq(T(n, 1+d-n), n=1..d), d=1..12);

f0[n_] := f0[n] = Function [k, k^n - Sum[f0[d][k], {d, Complement[Divisors[n], {n}]}]]; t[n_, k_] := f0[n][k]; Table[Table[t[n, 1 + d - n], {n, 1, d}], {d, 1, 12}] // Flatten (* Jean-François Alcover, Dec 12 2013, translated from Maple *)

T(n,k) = Sum_{d|n} k^d * mu(n/d).

T(n,k) = k^n - Sum_{d

T(n,k) = A143325(n,k) * k.

T(n,k) = A074650(n,k) * n.

So Sum_{d|n} k^d * mu(n/d) == 0 (mod n), this is a generalization of Fermat's little theorem k^p - k == 0 (mod p) for primes p to an arbitrary modulus n (see the Smyth link). - Franz Vrabec, Feb 09 2021

A033455 Convolution of nonzero squares A000290 with themselves.

Original entry on oeis.org

1, 8, 34, 104, 259, 560, 1092, 1968, 3333, 5368, 8294, 12376, 17927, 25312, 34952, 47328, 62985, 82536, 106666, 136136, 171787, 214544, 265420, 325520, 396045, 478296, 573678, 683704, 809999, 954304, 1118480, 1304512, 1514513, 1750728, 2015538, 2311464

Offset: 1

N. J. A. Sloane

Total area of all square regions from an n X n grid. E.g., at n = 3, there are nine individual squares, four 2 X 2's and one 3 X 3, total area 9 + 16 + 9 = 34, hence a(3) = 34. - Jon Perry, Jul 29 2003
If X is an n-set and Y and Z disjoint 2-subsets of X then a(n) is equal to the number of 7-subsets of X intersecting both Y and Z. - Milan Janjic, Aug 26 2007
Every fourth term is odd. However, there are no primes in the sequence. - Zak Seidov, Feb 28 2011
-120*a(n) is the real part of (n + n*i)*(n + 2 + n*i)*(n + (n + 2)i)*(n + 2+(n + 2)*i)*(n + 1 + (n + 1)*i), where i = sqrt(-1). - Jon Perry, Feb 05 2014
The previous formula rephrases the factorization of the 5th-order polynomial a(n) = (n+1)*((n+1)^4-1) = (n+1)*A123864(n+1) based on the factorization in A123865. - R. J. Mathar, Feb 08 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Abderrahim Arabi, Hacène Belbachir, Jean-Philippe Dubernard, Plateau Polycubes and Lateral Area, arXiv:1811.05707 [math.CO], 2018. See Column 1 Table 1 p. 8.
Milan Janjic, Two Enumerative Functions
Milan Janjić, On Restricted Ternary Words and Insets, arXiv:1905.04465 [math.CO], 2019.
C. P. Neuman and D. I. Schonbach, Evaluation of sums of convolved powers using Bernoulli numbers, SIAM Rev. 19 (1977), no. 1, 90--99. MR0428678 (55 #1698). See Table 2. - _N. J. A. Sloane_, Mar 23 2014
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).

Cf. A000290, A001249, A123865, A219086.

List([1..40], n-> ((n+1)^5 -(n+1))/30); # G. C. Greubel, Jul 05 2019

[(n^5-n)/30: n in [2..41]]; // Vincenzo Librandi, Mar 24 2014

Table[(n^5 - n)/30, {n,2,41}] (* Vladimir Joseph Stephan Orlovsky, Feb 28 2011 *)
CoefficientList[Series[(1+x)^2/(1-x)^6, {x,0,40}], x] (* Vincenzo Librandi, Mar 24 2014 *)

vector(40, n, ((n+1)^5 -n-1)/30) \\ G. C. Greubel, Jul 05 2019

[((n+1)^5 -n-1)/30 for n in (1..40)] # G. C. Greubel, Jul 05 2019

a(n-1) = n*(n^4 - 1)/30 = A061167(n)/30. - Henry Bottomley, Apr 18 2001
G.f.: x*(1+x)^2/(1-x)^6. - Philippe Deléham, Feb 21 2012
a(n) = Sum_{k=1..n+1} k^2*(n+1-k)^2. - Kolosov Petro, Feb 07 2019
E.g.f.: x*(30 +90*x +65*x^2 +15*x^3 +x^4)*exp(x)/30. - G. C. Greubel, Jul 05 2019

More terms from Vincenzo Librandi, Mar 24 2014

A138404 a(n) = prime(n)^5 - prime(n).

Original entry on oeis.org

30, 240, 3120, 16800, 161040, 371280, 1419840, 2476080, 6436320, 20511120, 28629120, 69343920, 115856160, 147008400, 229344960, 418195440, 714924240, 844596240, 1350125040, 1804229280, 2073071520, 3077056320, 3939040560

Offset: 1

Artur Jasinski, Mar 19 2008

Subsequence of A061167. - Bernard Schott, Feb 06 2023

Vincenzo Librandi, Table of n, a(n) for n = 1..200
Index to sequences related to prime powers

Cf. A000040, A006093, A008864, A050997, A061167, A066872, A138430.

[NthPrime((n))^5 - NthPrime((n)): n in [1..30] ]; // Vincenzo Librandi, Jun 17 2011

a = {}; Do[p = Prime[n]; AppendTo[a, p^5 - p], {n, 1, 50}]; a
#^5-#&/@Prime[Range[30]] (* Harvey P. Dale, Dec 25 2022 *)

forprime(p=2,1e3,print1(p^5-p", ")) \\ Charles R Greathouse IV, Jun 16 2011

a(n) = A050997(n) - A000040(n). - Elmo R. Oliveira, Jan 27 2023
From Bernard Schott, Feb 09 2023: (Start)
a(n) = A061167(A000040(n)).
a(n) = 30 * A138430(n).
a(n) = A000040(n) * A006093(n) * A008864(n) * A066872(n). (End)

A106512 Array read by antidiagonals: a(n,k) = number of k-colorings of a circle of n nodes (n >= 1, k >= 1).

Original entry on oeis.org

0, 0, 0, 0, 2, 0, 0, 6, 0, 0, 0, 12, 6, 2, 0, 0, 20, 24, 18, 0, 0, 0, 30, 60, 84, 30, 2, 0, 0, 42, 120, 260, 240, 66, 0, 0, 0, 56, 210, 630, 1020, 732, 126, 2, 0, 0, 72, 336, 1302, 3120, 4100, 2184, 258, 0, 0, 0, 90, 504, 2408, 7770, 15630, 16380, 6564, 510, 2, 0, 0, 110

Offset: 1

Joshua Zucker, May 29 2005

Note that we keep one edge in the circular graph even when there's only one node (so there are 0 colorings of one node with k colors).

Number of closed walks of length n on the complete graph K_{k}. - Andrew Howroyd, Mar 12 2017

			From _Andrew Howroyd_, Mar 12 2017: (Start)
Table begins:
  0 0   0     0      0       0        0        0         0 ...
  0 2   6    12     20      30       42       56        72 ...
  0 0   6    24     60     120      210      336       504 ...
  0 2  18    84    260     630     1302     2408      4104 ...
  0 0  30   240   1020    3120     7770    16800     32760 ...
  0 2  66   732   4100   15630    46662   117656    262152 ...
  0 0 126  2184  16380   78120   279930   823536   2097144 ...
  0 2 258  6564  65540  390630  1679622  5764808  16777224 ...
  0 0 510 19680 262140 1953120 10077690 40353600 134217720 ...
(End)
a(4,3) = 18 because there are three choices for the first node's color (call it 1) and then two choices for the second node's color (call it 2) and then the remaining two nodes can be 12, 13, or 32. So in total there are 3*2*3 = 18 ways. a(3,4) = 4*3*2 = 24 because the three nodes must be three distinct colors.

Andrew Howroyd, Table of n, a(n) for n = 1..1274

Columns include A092297, A226493. Main diagonal is A118537.
Rows 2-7 are A002378, A007531, A091940, A061167, A131472, A133499.
Cf. A090860, A208535.

a(n, k) = (k-1)^n + (-1)^n * (k-1).

a(67) corrected by Andrew Howroyd, Mar 12 2017

A133499 a(n) = n^7 - n.

Original entry on oeis.org

0, 0, 126, 2184, 16380, 78120, 279930, 823536, 2097144, 4782960, 9999990, 19487160, 35831796, 62748504, 105413490, 170859360, 268435440, 410338656, 612220014, 893871720, 1279999980, 1801088520, 2494357866, 3404825424, 4586471400, 6103515600, 8031810150

Offset: 0

Parthasarathy Nambi, Dec 01 2007

All the terms are divisible by 3.

			If n=2, then n^7 - n = 126.
If n=13, then n^7 - n = 62748504.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (8,-28,56,-70,56,-28,8,-1).

Cf. A030180, A061167.

[n^7-n: n in [0..40]]; // Vincenzo Librandi, May 02 2011

Table[n^7 - n, {n, 0, 40}] (* and *) LinearRecurrence[{8, -28, 56, -70, 56, -28, 8, -1}, {0, 0, 126, 2184, 16380, 78120, 279930, 823536}, 40] (* Vladimir Joseph Stephan Orlovsky, Feb 20 2012 *)

G.f.: 42*x^2*(3+28*x+58*x^2+28*x^3+3*x^4) / (x-1)^8 . - R. J. Mathar, Mar 13 2015

a(n) = 42*A030180(n). - Philippe Deléham, Mar 17 2016

Corrected and extended by Jonathan Vos Post, Dec 02 2007

A164938 a(n) = (n^5 - n)/10, which is always an integer.

Original entry on oeis.org

0, 3, 24, 102, 312, 777, 1680, 3276, 5904, 9999, 16104, 24882, 37128, 53781, 75936, 104856, 141984, 188955, 247608, 319998, 408408, 515361, 643632, 796260, 976560, 1188135, 1434888, 1721034, 2051112, 2429997, 2862912, 3355440, 3913536

Offset: 1

Bill Welsh (bill.welsh.75(AT)gmail.com), Aug 31 2009

Seiichi Manyama, Table of n, a(n) for n = 1..10000
Index entries for linear recurrences with constant coefficients, signature (6,-15,20,-15,6,-1).

Cf. A061167.

Table[(n^5 - n)/10, {n, 1, 50}] (* Stefan Steinerberger, Sep 03 2009 *)
LinearRecurrence[{6,-15,20,-15,6,-1},{0,3,24,102,312,777},50] (* Harvey P. Dale, Jan 14 2012 *)

a(n) = 6*a(n-1) - 15*a(n-2) + 20*a(n-3) - 15*a(n-4) + 6*a(n-5) - a(n-6) for n > 6. - Harvey P. Dale, Jan 14 2012
G.f.: (3*(x^3+2*x^2+x))/(x-1)^6. - Harvey P. Dale, Jan 14 2012
a(n) = A061167(n)/10. - Michel Marcus, Sep 04 2013
E.g.f.: exp(x)*x^2*(15 + 25*x + 10*x^2 + x^3)/10. - Stefano Spezia, Dec 27 2021

More terms from Stefan Steinerberger, Sep 03 2009

A262706 Triangle: Newton expansion of C(n,m)^5, read by rows.

Original entry on oeis.org

1, 0, 1, 0, 30, 1, 0, 150, 240, 1, 0, 240, 6810, 1020, 1, 0, 120, 63540, 94890, 3120, 1, 0, 0, 271170, 2615340, 740640, 7770, 1, 0, 0, 604800, 32186070, 47271840, 4029690, 16800, 1, 0, 0, 730800, 214628400, 1281612570, 518276640, 17075940, 32760, 1, 0, 0, 453600, 859992000, 18459063000, 26947757970, 4027831080, 60171300, 59040, 1

Offset: 0

Giuliano Cabrele, Sep 30 2015

Triangle here T_5(n,m) is such that C(n,m)^5 = Sum_{j=0..n} C(n,j)*T_5(j,m).
Equivalently, lower triangular matrix T_5 such that
|| C(n,m)^5 ||  =  P * T_5 = A007318 * T_5.
T_5(n,m) = 0 for n < m and for 5*m < n.
Refer to comment to A262704.
Example:
C(x,2)^5 = x^5*(x-1)^5/32 = 1*C(x,2) + 240*C(x,3) + 6810*C(x,4) + 63540*C(x,5) + 271170*C(x,6) + 604800*C(x,7) + 730800*C(x,8) + 453600*C(x,9) + 113400*C(x,10);
C(5,2)^5 = C(5,3)^5 = 100000 = 1*C(5,2) + 240*C(5,3) + 6810*C(5,4) + 63540*C(5,5) = 1*C(5,3) + 1020*C(5,4) + 94890*C(5,5).

			Triangle starts:
[1];
[0,   1];
[0,  30,      1];
[0, 150,    240,       1];
[0, 240,   6810,    1020,      1];
[0, 120,  63540,   94890,   3120,    1];
[0,   0, 271170, 2615340, 740640, 7770, 1];

P. Blasiak, K. A. Penson and A. I. Solomon, The general boson normal ordering problem, arXiv:quant-ph/0402027, 2004.

Second diagonal (T_5(n+1,n)) is A061167(n+1).
Column T_5(n,2) is A122193(5,n).
Cf. A109983 (transpose of), A262704, A262705.
Cf. A078739, A078741.

[&+[(-1)^(n-j)*Binomial(n,j)*Binomial(j,m)^5: j in [0..n]]: m in [0..n], n in [0..10]]; // Bruno Berselli, Oct 01 2015

T5[n_, m_] := Sum[(-1)^(n - j) * Binomial[n, j] * Binomial[j, m]^5, {j, 0, n}]; Table[T5[n, m], {n, 0, 9}, {m, 0, n}] // Flatten (* Jean-François Alcover, Oct 01 2015 *)

// as a function
T_5:=(n,m)->_plus((-1)^(n-j)*binomial(n,j)*binomial(j,m)^5 $ j=0..n):
// as a matrix h x h
_P:=h->matrix([[binomial(n,m) $m=0..h]$n=0..h]):
_P_5:=h->matrix([[binomial(n,m)^5 $m=0..h]$n=0..h]):
_T_5:=h->_P(h)^-1*_P_5(h):

T_5(nmax) = {for(n=0, nmax, for(m=0, n, print1(sum(j=0, n, (-1)^(n-j)*binomial(n,j)*binomial(j,m)^5), ", ")); print())} \\ Colin Barker, Oct 01 2015

T_5(n,m) = Sum_{j=0..n} (-1)^(n-j)*C(n,j)*C(j,m)^5.

Also, let S(r,s)(n,m) denote the Generalized Stirling2 numbers as defined in the link above, then T_5(n,m) = n! / (m!)^5 * S(m,m)(5,n).

A363916 Array read by descending antidiagonals. A(n, k) = Sum_{d=0..k} A363914(k, d) * n^d.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 0, 2, 1, 0, 0, 2, 3, 1, 0, 0, 6, 6, 4, 1, 0, 0, 12, 24, 12, 5, 1, 0, 0, 30, 72, 60, 20, 6, 1, 0, 0, 54, 240, 240, 120, 30, 7, 1, 0, 0, 126, 696, 1020, 600, 210, 42, 8, 1, 0, 0, 240, 2184, 4020, 3120, 1260, 336, 56, 9, 1

Offset: 0

Peter Luschny, Jul 04 2023

Row n gives the number of n-ary sequences with primitive period k.

See A074650 and A143324 for combinatorial interpretations.

			Array A(n, k) starts:
[0] 1, 0,  0,   0,    0,     0,      0,       0,        0, ... A000007
[1] 1, 1,  0,   0,    0,     0,      0,       0,        0, ... A019590
[2] 1, 2,  2,   6,   12,    30,     54,     126,      240, ... A027375
[3] 1, 3,  6,  24,   72,   240,    696,    2184,     6480, ... A054718
[4] 1, 4, 12,  60,  240,  1020,   4020,   16380,    65280, ... A054719
[5] 1, 5, 20, 120,  600,  3120,  15480,   78120,   390000, ... A054720
[6] 1, 6, 30, 210, 1260,  7770,  46410,  279930,  1678320, ... A054721
[7] 1, 7, 42, 336, 2352, 16800, 117264,  823536,  5762400, ... A218124
[8] 1, 8, 56, 504, 4032, 32760, 261576, 2097144, 16773120, ... A218125
A000012|A002378| A047928   |   A218130     |      A218131
    A001477,A007531,    A061167,        A133499,   (diagonal A252764)
.
Triangle T(n, k) starts:
[0] 1;
[1] 0, 1;
[2] 0, 1,  1;
[3] 0, 0,  2,   1;
[4] 0, 0,  2,   3,   1;
[5] 0, 0,  6,   6,   4,   1;
[6] 0, 0, 12,  24,  12,   5,  1;
[7] 0, 0, 30,  72,  60,  20,  6, 1;
[8] 0, 0, 54, 240, 240, 120, 30, 7, 1;

E. N. Gilbert and J. Riordan, Symmetry types of periodic sequences, Illinois J. Math., 5 (1961), 657-665.

Variant: A143324.
Rows: A000007 (n=0), A019590 (n=1), A027375 (n=2), A054718 (n=3), A054719 (n=4), A054720, A054721, A218124, A218125.
Columns: A000012 (k=0), A001477 (k=1), A002378 (k=2), A007531(k=3), A047928, A061167, A218130, A133499, A218131.
Cf. A252764 (main diagonal), A074650, A363914.

A363916 := (n, k) -> local d; add(A363914(k, d) * n^d, d = 0 ..k):
for n from 0 to 9 do seq(A363916(n, k), k = 0..8) od;

def A363916(n, k): return sum(A363914(k, d) * n^d for d in range(k + 1))
for n in range(9): print([A363916(n, k) for k in srange(9)])
def T(n, k): return A363916(k, n - k)

If k > 0 then k divides A(n, k), see the transposed array of A074650.

If k > 0 then n divides A(n, k), see the transposed array of A143325.

A342112 Drop the final digit of n^5.

Original entry on oeis.org

0, 0, 3, 24, 102, 312, 777, 1680, 3276, 5904, 10000, 16105, 24883, 37129, 53782, 75937, 104857, 141985, 188956, 247609, 320000, 408410, 515363, 643634, 796262, 976562, 1188137, 1434890, 1721036, 2051114, 2430000, 2862915, 3355443, 3913539, 4543542, 5252187, 6046617

Offset: 0

Stefano Spezia, Feb 28 2021

Why exponent 5? Because it is the smallest nontrivial exponent e such that for an integer k not ending in 0, 1, 5 and 6, k^e has the same unit digit of k in base 10.

Index entries for linear recurrences with constant coefficients, signature (5,-10, 10,-5,1,0,0,0,0,1,-5,10,-10,5,-1).
Index entries for sequences related to final digits of numbers.

Cf. A000584, A010879, A016813, A056865, A059995, A061167, A164938.

Table[(n^5-Last[IntegerDigits[n]])/10,{n,0,36}]

G.f.: x^2*(3 + 9*x + 12*x^2 + 12*x^3 + 12*x^4 + 12*x^5 + 12*x^6 + 12*x^7 + 13*x^8 + 8*x^9 + 15*x^10 - x^11 + x^12)/((1 - x)^6*(1 + x + x^2 + x^3 + x^4 + x^5 + x^6 + x^7 + x^8 + x^9)).
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5) + a(n-10) - 5*a(n-11) + 10*a(n-12) - 10*a(n-13) + 5*a(n-14) - a(n-15) for n > 14.
a(n) = floor(n^5/10).
a(n) = (A000584(n) - A010879(n))/10.
a(n) = A164938(n) + A059995(n).

Showing 1-10 of 10 results.

cp's OEIS Frontend

A052530 a(n) = 4*a(n-1) - a(n-2), with a(0) = 0, a(1) = 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Sage

Formula

Extensions

A143324 Table T(n,k) by antidiagonals. T(n,k) is the number of length n primitive (=aperiodic or period n) k-ary words (n,k >= 1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A033455 Convolution of nonzero squares A000290 with themselves.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Mathematica

PARI

Sage

Formula

Extensions

A138404 a(n) = prime(n)^5 - prime(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A106512 Array read by antidiagonals: a(n,k) = number of k-colorings of a circle of n nodes (n >= 1, k >= 1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Formula

Extensions

A133499 a(n) = n^7 - n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica