cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 19 results. Next

A194046 Natural interspersion of A052905, a rectangular array, by antidiagonals.

Original entry on oeis.org

1, 5, 2, 10, 6, 3, 16, 11, 7, 4, 23, 17, 12, 8, 9, 31, 24, 18, 13, 14, 15, 40, 32, 25, 19, 20, 21, 22, 50, 41, 33, 26, 27, 28, 29, 30, 61, 51, 42, 34, 35, 36, 37, 38, 39, 73, 62, 52, 43, 44, 45, 46, 47, 48, 49, 86, 74, 63, 53, 54, 55, 56, 57, 58, 59, 60, 100, 87, 75
Offset: 1

Views

Author

Clark Kimberling, Aug 13 2011

Keywords

Comments

See A194029 for definitions of natural fractal sequence and natural interspersion. Every positive integer occurs exactly once (and every pair of rows intersperse), so that as a sequence, A194046 is a permutation of the positive integers; its inverse is A194047.

Examples

			 Northwest corner:
1...5...10...16...23
2...6...11...17...24
3...7...12...18...25
4...8...13...19...26
9...14..20...27...35
		

Crossrefs

Programs

  • Mathematica
    z = 30;
    c[k_] := (k^2 + 5 k - 4)/2;
    c = Table[c[k], {k, 1, z}]  (* A052905 *)
    f[n_] := If[MemberQ[c, n], 1, 1 + f[n - 1]]
    f = Table[f[n], {n, 1, 255}]  (* fractal sequence [A002260] *)
    r[n_] := Flatten[Position[f, n]]
    t[n_, k_] := r[n][[k]]
    TableForm[Table[t[n, k], {n, 1, 7}, {k, 1, 7}]]
    p = Flatten[Table[t[k, n - k + 1], {n, 1, 13}, {k, 1, n}]]  (* A194046 *)
    q[n_] := Position[p, n]; Flatten[Table[q[n], {n, 1, 70}]] (* A194047 *)

A001477 The nonnegative integers.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77
Offset: 0

Views

Author

Keywords

Comments

Although this is a list, and lists normally have offset 1, it seems better to make an exception in this case. - N. J. A. Sloane, Mar 13 2010
The subsequence 0,1,2,3,4 gives the known values of n such that 2^(2^n)+1 is a prime (see A019434, the Fermat primes). - N. J. A. Sloane, Jun 16 2010
Also: The identity map, defined on the set of nonnegative integers. The restriction to the positive integers yields the sequence A000027. - M. F. Hasler, Nov 20 2013
The number of partitions of 2n into exactly 2 parts. - Colin Barker, Mar 22 2015
The number of orbits of Aut(Z^7) as function of the infinity norm n of the representative lattice point of the orbit, when the cardinality of the orbit is equal to 8960 or 168.- Philippe A.J.G. Chevalier, Dec 29 2015
Partial sums give A000217. - Omar E. Pol, Jul 26 2018
First differences are A000012 (the "all 1's" sequence). - M. F. Hasler, May 30 2020
See A061579 for the transposed infinite square matrix, or triangle with rows reversed. - M. F. Hasler, Nov 09 2021
This is the unique sequence (a(n)) that satisfies the inequality a(n+1) > a(a(n)) for all n in N. This simple and surprising result comes from the 6th problem proposed by Bulgaria during the second day of the 19th IMO (1977) in Belgrade (see link and reference). - Bernard Schott, Jan 25 2023

Examples

			Triangular view:
   0
   1   2
   3   4   5
   6   7   8   9
  10  11  12  13  14
  15  16  17  18  19  20
  21  22  23  24  25  26  27
  28  29  30  31  32  33  34  35
  36  37  38  39  40  41  42  43  44
  45  46  47  48  49  50  51  52  53  54
		

References

  • Maurice Protat, Des Olympiades à l'Agrégation, suite vérifiant f(n+1) > f(f(n)), Problème 7, pp. 31-32, Ellipses, Paris 1997.

Crossrefs

Cf. A000027 (n>=1).
Cf. A000012 (first differences).
Partial sums of A057427. - Jeremy Gardiner, Sep 08 2002
Cf. A038608 (alternating signs), A001787 (binomial transform).
Cf. A055112.
Cf. Boustrophedon transforms: A231179, A000737.
Cf. A245422.
Number of orbits of Aut(Z^7) as function of the infinity norm A000579, A154286, A102860, A002412, A045943, A115067, A008586, A008585, A005843, A000217.
When written as an array, the rows/columns are A000217, A000124, A152948, A152950, A145018, A167499, A166136, A167487... and A000096, A034856, A055998, A046691, A052905, A055999... (with appropriate offsets); cf. analogous lists for A000027 in A185787.
Cf. A000290.
Cf. A061579 (transposed matrix / reversed triangle).

Programs

Formula

a(n) = n.
a(0) = 0, a(n) = a(n-1) + 1.
G.f.: x/(1-x)^2.
Multiplicative with a(p^e) = p^e. - David W. Wilson, Aug 01 2001
When seen as array: T(k, n) = n + (k+n)*(k+n+1)/2. Main diagonal is 2*n*(n+1) (A046092), antidiagonal sums are n*(n+1)*(n+2)/2 (A027480). - Ralf Stephan, Oct 17 2004
Dirichlet generating function: zeta(s-1). - Franklin T. Adams-Watters, Sep 11 2005
E.g.f.: x*e^x. - Franklin T. Adams-Watters, Sep 11 2005
a(0)=0, a(1)=1, a(n) = 2*a(n-1) - a(n-2). - Jaume Oliver Lafont, May 07 2008
Alternating partial sums give A001057 = A000217 - 2*(A008794). - Eric Desbiaux, Oct 28 2008
a(n) = 2*A080425(n) + 3*A008611(n-3), n>1. - Eric Desbiaux, Nov 15 2009
a(n) = A007966(n)*A007967(n). - Reinhard Zumkeller, Jun 18 2011
a(n) = Sum_{k>=0} A030308(n,k)*2^k. - Philippe Deléham, Oct 20 2011
a(n) = 2*A028242(n-1) + (-1)^n*A000034(n-1). - R. J. Mathar, Jul 20 2012
a(n+1) = det(C(i+1,j), 1 <= i, j <= n), where C(n,k) are binomial coefficients. - Mircea Merca, Apr 06 2013
a(n-1) = floor(n/e^(1/n)) for n > 0. - Richard R. Forberg, Jun 22 2013
a(n) = A000027(n) for all n>0.
a(n) = floor(cot(1/(n+1))). - Clark Kimberling, Oct 08 2014
a(0)=0, a(n>0) = 2*z(-1)^[( |z|/z + 3 )/2] + ( |z|/z - 1 )/2 for z = A130472(n>0); a 1 to 1 correspondence between integers and naturals. - Adriano Caroli, Mar 29 2015
G.f. as triangle: x*(1 + (x^2 - 5*x + 2)*y + x*(2*x - 1)*y^2)/((1 - x)^3*(1 - x*y)^3). - Stefano Spezia, Jul 22 2025

A028387 a(n) = n + (n+1)^2.

Original entry on oeis.org

1, 5, 11, 19, 29, 41, 55, 71, 89, 109, 131, 155, 181, 209, 239, 271, 305, 341, 379, 419, 461, 505, 551, 599, 649, 701, 755, 811, 869, 929, 991, 1055, 1121, 1189, 1259, 1331, 1405, 1481, 1559, 1639, 1721, 1805, 1891, 1979, 2069, 2161, 2255, 2351, 2449, 2549, 2651
Offset: 0

Views

Author

Keywords

Comments

a(n+1) is the least k > a(n) + 1 such that A000217(a(n)) + A000217(k) is a square. - David Wasserman, Jun 30 2005
Values of Fibonacci polynomial n^2 - n - 1 for n = 2, 3, 4, 5, ... - Artur Jasinski, Nov 19 2006
A127701 * [1, 2, 3, ...]. - Gary W. Adamson, Jan 24 2007
Row sums of triangle A135223. - Gary W. Adamson, Nov 23 2007
Equals row sums of triangle A143596. - Gary W. Adamson, Aug 26 2008
a(n-1) gives the number of n X k rectangles on an n X n chessboard (for k = 1, 2, 3, ..., n). - Aaron Dunigan AtLee, Feb 13 2009
sqrt(a(0) + sqrt(a(1) + sqrt(a(2) + sqrt(a(3) + ...)))) = sqrt(1 + sqrt(5 + sqrt(11 + sqrt(19 + ...)))) = 2. - Miklos Kristof, Dec 24 2009
When n + 1 is prime, a(n) gives the number of irreducible representations of any nonabelian group of order (n+1)^3. - Andrew Rupinski, Mar 17 2010
a(n) = A176271(n+1, n+1). - Reinhard Zumkeller, Apr 13 2010
The product of any 4 consecutive integers plus 1 is a square (see A062938); the terms of this sequence are the square roots. - Harvey P. Dale, Oct 19 2011
Or numbers not expressed in the form m + floor(sqrt(m)) with integer m. - Vladimir Shevelev, Apr 09 2012
Left edge of the triangle in A214604: a(n) = A214604(n+1,1). - Reinhard Zumkeller, Jul 25 2012
Another expression involving phi = (1 + sqrt(5))/2 is a(n) = (n + phi)(n + 1 - phi). Therefore the numbers in this sequence, even if they are prime in Z, are not prime in Z[phi]. - Alonso del Arte, Aug 03 2013
a(n-1) = n*(n+1) - 1, n>=0, with a(-1) = -1, gives the values for a*c of indefinite binary quadratic forms [a, b, c] of discriminant D = 5 for b = 2*n+1. In general D = b^2 - 4ac > 0 and the form [a, b, c] is a*x^2 + b*x*y + c*y^2. - Wolfdieter Lang, Aug 15 2013
a(n) has prime factors given by A038872. - Richard R. Forberg, Dec 10 2014
A253607(a(n)) = -1. - Reinhard Zumkeller, Jan 05 2015
An example of a quadratic sequence for which the continued square root map (see A257574) produces the number 2. There are infinitely many sequences with this property - another example is A028387. See Popular Computing link. - N. J. A. Sloane, May 03 2015
Left edge of the triangle in A260910: a(n) = A260910(n+2,1). - Reinhard Zumkeller, Aug 04 2015
Numbers m such that 4m+5 is a square. - Bruce J. Nicholson, Jul 19 2017
The numbers represented as 131 in base n: 131_4 = 29, 131_5 = 41, ... . If 'digits' larger than the base are allowed then 131_2 = 11 and 131_1 = 5 also. - Ron Knott, Nov 14 2017
From Klaus Purath, Mar 18 2019: (Start)
Let m be a(n) or a prime factor of a(n). Then, except for 1 and 5, there are, if m is a prime, exactly two squares y^2 such that the difference y^2 - m contains exactly one pair of factors {x,z} such that the following applies: x*z = y^2 - m, x + y = z with
x < y, where {x,y,z} are relatively prime numbers. {x,y,z} are the initial values of a sequence of the Fibonacci type. Thus each a(n) > 5, if it is a prime, and each prime factor p > 5 of an a(n) can be assigned to exactly two sequences of the Fibonacci type. a(0) = 1 belongs to the original Fibonacci sequence and a(1) = 5 to the Lucas sequence.
But also the reverse assignment applies. From any sequence (f(i)) of the Fibonacci type we get from its 3 initial values by f(i)^2 - f(i-1)*f(i+1) with f(i-1) < f(i) a term a(n) or a prime factor p of a(n). This relation is also valid for any i. In this case we get the absolute value |a(n)| or |p|. (End)
a(n-1) = 2*T(n) - 1, for n>=1, with T = A000217, is a proper subsequence of A089270, and the terms are 0,-1,+1 (mod 5). - Wolfdieter Lang, Jul 05 2019
a(n+1) is the number of wedged n-dimensional spheres in the homotopy of the neighborhood complex of Kneser graph KG_{2,n}. Here, KG_{2,n} is a graph whose vertex set is the collection of subsets of cardinality 2 of set {1,2,...,n+3,n+4} and two vertices are adjacent if and only if they are disjoint. - Anurag Singh, Mar 22 2021
Also the number of squares between (n+2)^2 and (n+2)^4. - Karl-Heinz Hofmann, Dec 07 2021
(x, y, z) = (A001105(n+1), -a(n-1), -a(n)) are solutions of the Diophantine equation x^3 + 4*y^3 + 4*z^3 = 8. - XU Pingya, Apr 25 2022
The least significant digit of terms of this sequence cycles through 1, 5, 1, 9, 9. - Torlach Rush, Jun 05 2024

Examples

			From _Ilya Gutkovskiy_, Apr 13 2016: (Start)
Illustration of initial terms:
                                        o               o
                        o           o   o o           o o
            o       o   o o       o o   o o o       o o o
    o   o   o o   o o   o o o   o o o   o o o o   o o o o
o   o o o   o o o o o   o o o o o o o   o o o o o o o o o
n=0  n=1       n=2           n=3               n=4
(End)
From _Klaus Purath_, Mar 18 2019: (Start)
Examples:
a(0) = 1: 1^1-0*1 = 1, 0+1 = 1 (Fibonacci A000045).
a(1) = 5: 3^2-1*4 = 5, 1+3 = 4 (Lucas A000032).
a(2) = 11: 4^2-1*5 = 11, 1+4 = 5 (A000285); 5^2-2*7 = 11, 2+5 = 7 (A001060).
a(3) = 19: 5^2-1*6 = 19, 1+5 = 6 (A022095); 7^2-3*10 = 19, 3+7 = 10 (A022120).
a(4) = 29: 6^2-1*7 = 29, 1+6 = 7 (A022096); 9^2-4*13 = 29, 4+9 = 13 (A022130).
a(11)/5 = 31: 7^2-2*9 = 31, 2+7 = 9 (A022113); 8^2-3*11 = 31, 3+8 = 11 (A022121).
a(24)/11 = 59: 9^2-2*11 = 59, 2+9 = 11 (A022114); 12^2-5*17 = 59, 5+12 = 17 (A022137).
(End)
		

Crossrefs

Complement of A028392. Third column of array A094954.
Cf. A000217, A002522, A062392, A062786, A127701, A135223, A143596, A052905, A162997, A062938 (squares of this sequence).
A110331 and A165900 are signed versions.
Cf. A002327 (primes), A094210.
Frobenius number for k successive numbers: this sequence (k=2), A079326 (k=3), A138984 (k=4), A138985 (k=5), A138986 (k=6), A138987 (k=7), A138988 (k=8).

Programs

Formula

a(n) = sqrt(A062938(n)). - Floor van Lamoen, Oct 08 2001
a(0) = 1, a(1) = 5, a(n) = (n+1)*a(n-1) - (n+2)*a(n-2) for n > 1. - Gerald McGarvey, Sep 24 2004
a(n) = A105728(n+2, n+1). - Reinhard Zumkeller, Apr 18 2005
a(n) = A109128(n+2, 2). - Reinhard Zumkeller, Jun 20 2005
a(n) = 2*T(n+1) - 1, where T(n) = A000217(n). - Gary W. Adamson, Aug 15 2007
a(n) = A005408(n) + A002378(n); A084990(n+1) = Sum_{k=0..n} a(k). - Reinhard Zumkeller, Aug 20 2007
Binomial transform of [1, 4, 2, 0, 0, 0, ...] = (1, 5, 11, 19, ...). - Gary W. Adamson, Sep 20 2007
G.f.: (1+2*x-x^2)/(1-x)^3. a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - R. J. Mathar, Jul 11 2009
a(n) = (n + 2 + 1/phi) * (n + 2 - phi); where phi = 1.618033989... Example: a(3) = 19 = (5 + .6180339...) * (3.381966...). Cf. next to leftmost column in A162997 array. - Gary W. Adamson, Jul 23 2009
a(n) = a(n-1) + 2*(n+1), with n > 0, a(0) = 1. - Vincenzo Librandi, Nov 18 2010
For k < n, a(n) = (k+1)*a(n-k) - k*a(n-k-1) + k*(k+1); e.g., a(5) = 41 = 4*11 - 3*5 + 3*4. - Charlie Marion, Jan 13 2011
a(n) = lower right term in M^2, M = the 2 X 2 matrix [1, n; 1, (n+1)]. - Gary W. Adamson, Jun 29 2011
G.f.: (x^2-2*x-1)/(x-1)^3 = G(0) where G(k) = 1 + x*(k+1)*(k+4)/(1 - 1/(1 + (k+1)*(k+4)/G(k+1))); (continued fraction, 3-step). - Sergei N. Gladkovskii, Oct 16 2012
Sum_{n>0} 1/a(n) = 1 + Pi*tan(sqrt(5)*Pi/2)/sqrt(5). - Enrique Pérez Herrero, Oct 11 2013
E.g.f.: exp(x) (1+4*x+x^2). - Tom Copeland, Dec 02 2013
a(n) = A005408(A000217(n)). - Tony Foster III, May 31 2016
From Amiram Eldar, Jan 29 2021: (Start)
Product_{n>=0} (1 + 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2).
Product_{n>=1} (1 - 1/a(n)) = -Pi*sec(sqrt(5)*Pi/2)/6. (End)
a(5*n+1)/5 = A062786(n+1). - Torlach Rush, Jun 05 2024

Extensions

Minor edits by N. J. A. Sloane, Jul 04 2010, following suggestions from the Sequence Fans Mailing List

A055998 a(n) = n*(n+5)/2.

Original entry on oeis.org

0, 3, 7, 12, 18, 25, 33, 42, 52, 63, 75, 88, 102, 117, 133, 150, 168, 187, 207, 228, 250, 273, 297, 322, 348, 375, 403, 432, 462, 493, 525, 558, 592, 627, 663, 700, 738, 777, 817, 858, 900, 943, 987, 1032, 1078, 1125, 1173, 1222, 1272
Offset: 0

Views

Author

Barry E. Williams, Jun 14 2000

Keywords

Comments

If X is an n-set and Y a fixed (n-3)-subset of X then a(n-3) is equal to the number of 2-subsets of X intersecting Y. - Milan Janjic, Aug 15 2007
Bisection of A165157. - Jaroslav Krizek, Sep 05 2009
a(n) is the number of (w,x,y) having all terms in {0,...,n} and w=x+y-1. - Clark Kimberling, Jun 02 2012
Numbers m >= 0 such that 8m+25 is a square. - Bruce J. Nicholson, Jul 26 2017
a(n-1) = 3*(n-1) + (n-1)*(n-2)/2 is the number of connected, loopless, non-oriented, multi-edge vertex-labeled graphs with n edges and 3 vertices. Labeled multigraph analog of A253186. There are 3*(n-1) graphs with the 3 vertices on a chain (3 ways to label the middle graph, n-1 ways to pack edges on one of connections) and binomial(n-1,2) triangular graphs (one way to label the graphs, pack 1 or 2 or ...n-2 on the 1-2 edge, ...). - R. J. Mathar, Aug 10 2017
a(n) is also the number of vertices of the quiver for PGL_{n+1} (see Shen). - Stefano Spezia, Mar 24 2020
Starting from a(2) = 7, this is the 4th column of the array: natural numbers written by antidiagonals downwards. See the illustration by Kival Ngaokrajang and the cross-references. - Andrey Zabolotskiy, Dec 21 2021

References

  • Albert H. Beiler, Recreations in the Theory of Numbers, Dover, N.Y., 1964, p. 193.

Crossrefs

a(n) = A095660(n+1, 2): third column of (1, 3)-Pascal triangle.
Row n=2 of A255961.

Programs

Formula

G.f.: x*(3-2*x)/(1-x)^3.
a(n) = A027379(n), n > 0.
a(n) = A126890(n,2) for n > 1. - Reinhard Zumkeller, Dec 30 2006
a(n) = A000217(n) + A005843(n). - Reinhard Zumkeller, Sep 24 2008
If we define f(n,i,m) = Sum_{k=0..n-i} binomial(n,k)*Stirling1(n-k,i)*Product_{j=0..k-1} (-m-j), then a(n) = -f(n,n-1,3), for n >= 1. - Milan Janjic, Dec 20 2008
a(n) = A167544(n+8). - Philippe Deléham, Nov 25 2009
a(n) = a(n-1) + n + 2 with a(0)=0. - Vincenzo Librandi, Aug 07 2010
a(n) = Sum_{k=1..n} (k+2). - Gary Detlefs, Aug 10 2010
a(n) = A034856(n+1) - 1 = A000217(n+2) - 3. - Jaroslav Krizek, Sep 05 2009
Sum_{n>=1} 1/a(n) = 137/150. - R. J. Mathar, Jul 14 2012
a(n) = 3*n + A000217(n-1) = 3*n - floor(n/2) + floor(n^2/2). - Wesley Ivan Hurt, Jun 15 2013
a(n) = Sum_{i=3..n+2} i. - Wesley Ivan Hurt, Jun 28 2013
a(n) = 3*A000217(n) - 2*A000217(n-1). - Bruno Berselli, Dec 17 2014
a(n) = A046691(n) + 1. Also, a(n) = A052905(n-1) + 2 = A055999(n-1) + 3 for n>0. - Andrey Zabolotskiy, May 18 2016
E.g.f.: x*(6+x)*exp(x)/2. - G. C. Greubel, Apr 05 2019
Sum_{n>=1} (-1)^(n+1)/a(n) = 4*log(2)/5 - 47/150. - Amiram Eldar, Jan 10 2021
From Amiram Eldar, Feb 12 2024: (Start)
Product_{n>=1} (1 - 1/a(n)) = -5*cos(sqrt(33)*Pi/2)/(4*Pi).
Product_{n>=1} (1 + 1/a(n)) = 15*cos(sqrt(17)*Pi/2)/(2*Pi). (End)

A185787 Sum of first k numbers in column k of the natural number array A000027; by antidiagonals.

Original entry on oeis.org

1, 7, 25, 62, 125, 221, 357, 540, 777, 1075, 1441, 1882, 2405, 3017, 3725, 4536, 5457, 6495, 7657, 8950, 10381, 11957, 13685, 15572, 17625, 19851, 22257, 24850, 27637, 30625, 33821, 37232, 40865, 44727, 48825, 53166, 57757, 62605, 67717, 73100, 78761, 84707, 90945, 97482, 104325, 111481, 118957, 126760, 134897, 143375
Offset: 1

Views

Author

Clark Kimberling, Feb 03 2011

Keywords

Comments

This is one of many interesting sequences and arrays that stem from the natural number array A000027, of which a northwest corner is as follows:
1....2.....4.....7...11...16...22...29...
3....5.....8....12...17...23...30...38...
6....9....13....18...24...31...39...48...
10...14...19....25...32...40...49...59...
15...20...26....33...41...50...60...71...
21...27...34....42...51...61...72...84...
28...35...43....52...62...73...85...98...
Blocking out all terms below the main diagonal leaves columns whose sums comprise A185787. Deleting the main diagonal and then summing give A185787. Analogous treatments to the left of the main diagonal give A100182 and A101165. Further sequences obtained directly from this array are easily obtained using the following formula for the array: T(n,k)=n+(n+k-2)(n+k-1)/2.
Examples:
row 1: A000124
row 2: A022856
row 3: A016028
row 4: A145018
row 5: A077169
col 1: A000217
col 2: A000096
col 3: A034856
col 4: A055998
col 5: A046691
col 6: A052905
col 7: A055999
diag. (1,5,...) ...... A001844
diag. (2,8,...) ...... A001105
diag. (4,12,...)...... A046092
diag. (7,17,...)...... A056220
diag. (11,23,...) .... A132209
diag. (16,30,...) .... A054000
diag. (22,38,...) .... A090288
diag. (3,9,...) ...... A058331
diag. (6,14,...) ..... A051890
diag. (10,20,...) .... A005893
diag. (15,27,...) .... A097080
diag. (21,35,...) .... A093328
antidiagonal sums: (1,5,15,34,...)=A006003=partial sums of A002817.
Let S(n,k) denote the n-th partial sum of column k. Then
S(n,k)=n*(n^2+3k*n+3*k^2-6*k+5)/6.
S(n,1)=n(n+1)(n+2)/6
S(n,2)=n(n+1)(n+5)/6
S(n,3)=n(n+2)(n+7)/6
S(n,4)=n(n^2+12n+29)/6
S(n,5)=n(n+5)(n+10)/6
S(n,6)=n(n+7)(n+11)/6
S(n,7)=n(n+10)(n+11)/6
Weight array of T: A144112
Accumulation array of T: A185506
Second rectangular sum array of T: A185507
Third rectangular sum array of T: A185508
Fourth rectangular sum array of T: A185509

Crossrefs

Programs

  • Magma
    [n*(7*n^2-6*n+5)/6: n in [1..50]]; // Vincenzo Librandi, Jul 04 2012
  • Mathematica
    f[n_,k_]:=n+(n+k-2)(n+k-1)/2;
    s[k_]:=Sum[f[n,k],{n,1,k}];
    Factor[s[k]]
    Table[s[k],{k,1,70}]  (* A185787 *)
    CoefficientList[Series[(3*x^2+3*x+1)/(1-x)^4,{x,0,50}],x] (* Vincenzo Librandi, Jul 04 2012 *)

Formula

a(n)=n*(7*n^2-6*n+5)/6.
G.f.: x*(3*x^2+3*x+1)/(1-x)^4. - Vincenzo Librandi, Jul 04 2012

Extensions

Edited by Clark Kimberling, Feb 25 2023

A059605 a(n) = (1/3!)*(n^3 + 24*n^2 + 107*n + 90), compare A059604.

Original entry on oeis.org

15, 37, 68, 109, 161, 225, 302, 393, 499, 621, 760, 917, 1093, 1289, 1506, 1745, 2007, 2293, 2604, 2941, 3305, 3697, 4118, 4569, 5051, 5565, 6112, 6693, 7309, 7961, 8650, 9377, 10143, 10949, 11796, 12685, 13617, 14593, 15614, 16681, 17795, 18957
Offset: 0

Views

Author

Vladeta Jovovic, Jan 29 2001

Keywords

Crossrefs

Programs

  • Magma
    [(1/6)*(n^3+24*n^2+107*n+90) : n in [0..50]]; // Vincenzo Librandi, Nov 13 2011

Formula

G.f.: (15 - 23*x + 10*x^2 - x^3)/(1-x)^4, compare A059340.

Extensions

More terms from James Sellers, Feb 01 2001

A101882 Write three numbers, skip one, write three, skip two, write three, skip three... and so on.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 12, 16, 17, 18, 23, 24, 25, 31, 32, 33, 40, 41, 42, 50, 51, 52, 61, 62, 63, 73, 74, 75, 86, 87, 88, 100, 101, 102, 115, 116, 117, 131, 132, 133, 148, 149, 150, 166, 167, 168, 185, 186, 187, 205, 206, 207, 226, 227, 228, 248, 249, 250, 271
Offset: 1

Views

Author

Candace Mills (scorpiocand(AT)yahoo.com), Dec 19 2004

Keywords

Comments

Union of A052905, A052905+1, and A052905+2. - Ivan Neretin, Aug 03 2016
First terms of the 3 repeated terms belong to A052905. - Michael De Vlieger, Aug 03 2016

Crossrefs

Programs

  • Mathematica
    Flatten@Table[(n^2 + 5 n - 4)/2 + {0, 1, 2}, {n, 20}] (* Ivan Neretin, Aug 03 2016 *)
    Table[Range[#, # + 2] &[(n^2 + 7 n + 2)/2], {n, 0, 20}] // Flatten (* or *)
    Rest@ CoefficientList[Series[x (1 + x + x^2 - x^4 - x^5)/((1 + x + x^2)^2 (1 - x)^3), {x, 0, 61}], x] (* Michael De Vlieger, Aug 03 2016 *)
    LinearRecurrence[{1,0,2,-2,0,-1,1},{1,2,3,5,6,7,10},70] (* Harvey P. Dale, Dec 26 2019 *)
  • PARI
    a(n)=my(k=n%3); if(k==2, n^2+17*n-2, k==1, n^2+19*n-2, n^2+15*n)/18 \\ Charles R Greathouse IV, Aug 03 2016

Formula

G.f.: x*(1+x+x^2-x^4-x^5)/ ((1+x+x^2)^2 * (1-x)^3). [Maksym Voznyy (voznyy(AT)mail.ru), Aug 11 2009; checked and corrected by R. J. Mathar, Sep 16 2009]
a(n) = n + k * (k+1) / 2 where k = floor((n-1) / 3). - Ziad Ahmed, Jun 17 2025

A164981 A triangle with Pell numbers in the first column.

Original entry on oeis.org

1, 2, 1, 5, 3, 1, 12, 10, 4, 1, 29, 30, 16, 5, 1, 70, 87, 56, 23, 6, 1, 169, 245, 185, 91, 31, 7, 1, 408, 676, 584, 334, 136, 40, 8, 1, 985, 1836, 1784, 1158, 546, 192, 50, 9, 1, 2378, 4925, 5312, 3850, 2052, 834, 260, 61, 10, 1, 5741, 13079, 15497, 12386, 7342, 3366, 1212, 341, 73, 11, 1
Offset: 1

Views

Author

Mark Dols, Sep 03 2009

Keywords

Comments

Rows sum up to A000244 (powers of 3), diagonals to A001654 (golden rectangles).
Up to reflection at the vertical axis, the triangle of numbers given here coincides with the triangle given in A210557, i.e. the numbers are the same just read row-wise in the opposite direction. [Christine Bessenrodt, Jul 20 2012]
Subtriangle of (0, 2, 1/2, -1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, 0, -1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938. - Philippe Deléham, Oct 10 2013

Examples

			Triangle begins
  1
  2,1
  5,3,1
  12,10,4,1
  29,30,16,5,1
  70,87,56,23,6,1
  169,245,185,91,31,7,1
  ...
From _Philippe Deléham_, Oct 10 2013: (Start)
Triangle (0, 2, 1/2, -1/2, 0, 0, ...) DELTA (1, 0, -1/2, 1/2, 0, 0, ...):
  1
  0, 1
  0, 2, 1
  0, 5, 3, 1
  0, 12, 10, 4, 1
  0, 29, 30, 16, 5, 1
  0, 70, 87, 56, 23, 6, 1
  0, 169, 245, 185, 91, 31, 7, 1
  ... (End)
		

Crossrefs

Programs

  • Maple
    A164981 := proc(n,k) option remember; if n <1 or k<1 or k>n then 0; elif n = 1 then 1; else 2*procname(n-1,k)+procname(n-1,k-1)+procname(n-2,k)-procname(n-2,k-1) ; end if; end proc:
  • Mathematica
    T[n_, k_] := T[n, k] = Which[n < 1 || k < 1 || k > n, 0, n == 1, 1, True, 2*T[n-1, k] + T[n-1, k-1] + T[n-2, k] - T[n-2, k-1]];
    Table[T[n, k], {n, 1, 11}, {k, 1, n}] // Flatten (* Jean-François Alcover, Aug 06 2023 *)
  • PARI
    T(n,k) = if ((n==1) && (k==1), return(1)); if ((n<=0) || (k<=0) || (nMichel Marcus, Feb 01 2023

Formula

From R. J. Mathar, Jan 27 2011: (Start)
T(1,1) =1. T(n,k)=0 if n<1 or k<1 or k>n. T(n,k) = 2*T(n-1,k)+T(n-1,k-1)+T(n-2,k)-T(n-2,k-1) otherwise.
T(n,1) = A000129(n).
T(n,n-1) = n.
T(n,n-2) = A052905(n-2).
T(n,2) = A026937(n-2). (End)
G.f. x*y/(1-2*x-x^2+x^2*y-x*y). - R. J. Mathar, Aug 11 2015

Extensions

Rows 10-11 from Michel Marcus, Feb 01 2023

A210552 Triangle of coefficients of polynomials u(n,x) jointly generated with A210553; see the Formula section.

Original entry on oeis.org

1, 1, 2, 1, 3, 3, 1, 4, 5, 5, 1, 5, 7, 10, 8, 1, 6, 9, 16, 18, 13, 1, 7, 11, 23, 31, 33, 21, 1, 8, 13, 31, 47, 62, 59, 34, 1, 9, 15, 40, 66, 101, 119, 105, 55, 1, 10, 17, 50, 88, 151, 205, 227, 185, 89, 1, 11, 19, 61, 113, 213, 321, 414, 426, 324, 144, 1, 12, 21, 73
Offset: 1

Views

Author

Clark Kimberling, Mar 22 2012

Keywords

Comments

Let T(n,k) denote the term in row n, column k.
T(n,n): A000045 (Fibonacci numbers)
T(n,n-1): A010049 (second-order Fibonacci numbers)
T(n,1): 1,1,1,1,1,1,1,1,1,1,1,,...
T(n,2): 2,3,4,5,6,7,8,9,10,11,...
T(n,3): 3,5,7,9,11,13,15,17,19,...
T(n,4): A052905
Row sums: A000225
Alternating row sums: A094024 (signed)
For a discussion and guide to related arrays, see A208510.

Examples

			First five rows:
1
1...2
1...3...3
1...4...5...5
1...5...7...10...8
First three polynomials u(n,x): 1, 1 + 2x, 1 + 3x + 3x^2.
		

Crossrefs

Programs

  • Mathematica
    u[1, x_] := 1; v[1, x_] := 1; z = 16;
    u[n_, x_] := x*u[n - 1, x] + x*v[n - 1, x] + 1;
    v[n_, x_] := x*u[n - 1, x] + v[n - 1, x] + 1;
    Table[Expand[u[n, x]], {n, 1, z/2}]
    Table[Expand[v[n, x]], {n, 1, z/2}]
    cu = Table[CoefficientList[u[n, x], x], {n, 1, z}];
    TableForm[cu]
    Flatten[%]   (* A210552 *)
    Table[Expand[v[n, x]], {n, 1, z}]
    cv = Table[CoefficientList[v[n, x], x], {n, 1, z}];
    TableForm[cv]
    Flatten[%]   (* A210553 *)
    Table[u[n, x] /. x -> 1, {n, 1, z}]  (* A000225 *)
    Table[v[n, x] /. x -> 1, {n, 1, z}]  (* A000225 *)
    Table[u[n, x] /. x -> -1, {n, 1, z}] (* A094024 *)
    Table[v[n, x] /. x -> -1, {n, 1, z}] (* A052551 *)

Formula

u(n,x)=x*u(n-1,x)+x*v(n-1,x)+1,
v(n,x)=x*u(n-1,x)+v(n-1,x)+1,
where u(1,x)=1, v(1,x)=1.

A302537 a(n) = (n^2 + 13*n + 2)/2.

Original entry on oeis.org

1, 8, 16, 25, 35, 46, 58, 71, 85, 100, 116, 133, 151, 170, 190, 211, 233, 256, 280, 305, 331, 358, 386, 415, 445, 476, 508, 541, 575, 610, 646, 683, 721, 760, 800, 841, 883, 926, 970, 1015, 1061, 1108, 1156, 1205, 1255, 1306, 1358, 1411, 1465, 1520, 1576
Offset: 0

Views

Author

Keywords

Comments

Binomial transform of [1, 7, 1, 0, 0, 0, ...].
Numbers m > 0 such that 8*m + 161 is a square.

Examples

			Illustration of initial terms (by the formula a(n) = A052905(n) + 3*n):
.                                                                    o
.                                                                  o o
.                                                    o           o o o
.                                                  o o         o o o o
.                                      o         o o o       o o o o o
.                                    o o       o o o o     o o o o o o
.                          o       o o o     o o o o o   o . . . . . o
.                        o o     o o o o   o . . . . o   o . . . . . o
.                o     o o o   o . . . o   o . . . . o   o . . . . . o
.              o o   o . . o   o . . . o   o . . . . o   o . . . . . o
.        o   o . o   o . . o   o . . . o   o . . . . o   o . . . . . o
.      o o   o . o   o . . o   o . . . o   o . . . . o   o . . . . . o
.  o   o o   o o o   o o o o   o o o o o   o o o o o o   o o o o o o o
.        o     o o     o o o     o o o o     o o o o o     o o o o o o
.        o     o o     o o o     o o o o     o o o o o     o o o o o o
.        o     o o     o o o     o o o o     o o o o o     o o o o o o
----------------------------------------------------------------------
.  1     8      16        25          35            46              58
		

References

  • R. L. Graham, D. E. Knuth and O. Patashnik, Concrete Mathematics: A Foundation for Computer Science, Addison-Wesley, 1994.

Crossrefs

Sequences whose n-th terms are of the form binomial(n, 2) + n*k + 1:
A152947 (k = 0); A000124 (k = 1); A000217 (k = 2); A034856 (k = 3);
A052905 (k = 4); A051936 (k = 5); A246172 (k = 6).

Programs

  • Magma
    A302537:= func< n | ((n+1)^2 +12*n +1)/2 >;
    [A302537(n): n in [0..50]]; // G. C. Greubel, Jan 21 2025
    
  • Maple
    a := n -> (n^2 + 13*n + 2)/2;
    seq(a(n), n = 0 .. 100);
  • Mathematica
    Table[(n^2 + 13 n + 2)/2, {n, 0, 100}]
    CoefficientList[ Series[(5x^2 - 5x - 1)/(x - 1)^3, {x, 0, 50}], x] (* or *)
    LinearRecurrence[{3, -3, 1}, {1, 8, 16}, 51] (* Robert G. Wilson v, May 19 2018 *)
  • Maxima
    makelist((n^2 + 13*n + 2)/2, n, 0, 100);
    
  • PARI
    a(n) = (n^2 + 13*n + 2)/2; \\ Altug Alkan, Apr 12 2018
    
  • Python
    def A302537(n): return (n**2 + 13*n + 2)//2
    print([A302537(n) for n in range(51)]) # G. C. Greubel, Jan 21 2025

Formula

a(n) = binomial(n + 1, 2) + 6*n + 1 = binomial(n, 2) + 7*n + 1.
a(n) = a(n-1) + n + 6.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n >= 3, where a(0) = 1, a(1) = 8 and a(2) = 16.
a(n) = 2*a(n-1) - a(n-2) + 1.
a(n) = A004120(n+1) for n > 1.
a(n) = A056119(n) + 1.
a(n) = A152947(n+1) + A008589(n).
a(n) = A060544(n+1) - A002939(n).
a(n) = A000578(n+1) - A162261(n) for n > 0.
G.f.: (1 + 5*x - 5*x^2)/(1 - x)^3.
E.g.f.: (1/2)*(2 + 14*x + x^2)*exp(x).
Sum_{n>=0} 1/a(n) = 24097/45220 + 2*Pi*tan(sqrt(161)*Pi/2) / sqrt(161) = 1.4630922534498496... - Vaclav Kotesovec, Apr 11 2018
Showing 1-10 of 19 results. Next