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 11 results. Next

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

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

A229445 T(n,k)=Number of nXk 0..2 arrays with horizontal differences mod 3 never 1, vertical differences mod 3 never -1, and rows and columns lexicographically nondecreasing.

Original entry on oeis.org

3, 4, 5, 5, 7, 8, 6, 10, 13, 12, 7, 14, 22, 25, 17, 8, 19, 37, 53, 47, 23, 9, 25, 60, 109, 128, 84, 30, 10, 32, 93, 212, 324, 293, 142, 38, 11, 40, 138, 387, 753, 915, 625, 228, 47, 12, 49, 197, 665, 1609, 2546, 2402, 1244, 350, 57, 13, 59, 272, 1083, 3184, 6374, 8024
Offset: 1

Views

Author

R. H. Hardin Sep 23 2013

Keywords

Comments

Table starts
..3...4....5....6.....7.....8......9.....10......11......12......13.......14
..5...7...10...14....19....25.....32.....40......49......59......70.......82
..8..13...22...37....60....93....138....197.....272.....365.....478......613
.12..25...53..109...212...387....665...1083....1684....2517....3637.....5105
.17..47..128..324...753..1609...3184...5890...10281...17075...27176....41696
.23..84..293..915..2546..6374..14536..30571...59969..110816..194535...326723
.30.142..625.2402..8024.23610..62205.149031..329106..677706.1314145..2419348
.38.228.1244.5843.23428.81177.247607.676983.1685570.3873314.8307126.16784531

Examples

			Some solutions for n=4 k=4
..0..2..2..2....0..2..2..2....0..0..2..2....0..0..2..2....0..2..2..2
..1..0..0..2....1..0..0..0....0..0..2..2....1..1..0..0....0..2..2..2
..2..1..1..0....2..1..1..1....1..1..0..0....1..1..1..1....0..2..2..2
..2..1..1..1....2..2..2..2....1..1..1..1....2..2..1..1....1..0..0..2
		

Crossrefs

Column 1 is A022856(n+4)
Row 2 is A145018(n+1)

Formula

Empirical for column k:
k=1: a(n) = (1/2)*n^2 + (1/2)*n + 2
k=2: a(n) = (1/24)*n^4 + (1/12)*n^3 - (1/24)*n^2 + (23/12)*n + 2
k=3: [polynomial of degree 6]
k=4: [polynomial of degree 8]
k=5: [polynomial of degree 10]
k=6: [polynomial of degree 12]
k=7: [polynomial of degree 14]
Empirical for row n:
n=1: a(n) = n + 2
n=2: a(n) = (1/2)*n^2 + (1/2)*n + 4
n=3: a(n) = (1/3)*n^3 + (8/3)*n + 5
n=4: a(n) = (1/4)*n^4 - (1/3)*n^3 + (13/4)*n^2 + (11/6)*n + 7
n=5: a(n) = (11/60)*n^5 - (1/2)*n^4 + (15/4)*n^3 - n^2 + (257/30)*n + 6
n=6: [polynomial of degree 6]
n=7: [polynomial of degree 7]

A027693 a(n) = n^2 + n + 8.

Original entry on oeis.org

8, 10, 14, 20, 28, 38, 50, 64, 80, 98, 118, 140, 164, 190, 218, 248, 280, 314, 350, 388, 428, 470, 514, 560, 608, 658, 710, 764, 820, 878, 938, 1000, 1064, 1130, 1198, 1268, 1340, 1414, 1490, 1568, 1648, 1730, 1814, 1900, 1988, 2078, 2170, 2264, 2360, 2458, 2558
Offset: 0

Views

Author

Keywords

Crossrefs

Programs

Formula

a(n) = 2*n + a(n-1) (with a(0)=8). - Vincenzo Librandi, Aug 05 2010
From Harvey P. Dale, Dec 13 2011: (Start)
a(0)=8, a(1)=10, a(2)=14, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
G.f.: (2*(7-4*x)*x-8)/(x-1)^3. (End)
Sum_{n>=0} 1/a(n) = Pi*tanh(Pi*sqrt(31)/2)/sqrt(31). - Amiram Eldar, Jan 17 2021
From Elmo R. Oliveira, Oct 31 2024: (Start)
E.g.f.: exp(x)*(8 + 2*x + x^2).
a(n) = 2*A145018(n+1). (End)

A279967 Square array read by antidiagonals upwards in which each term is the sum of prior elements in the same row, column, diagonal, or antidiagonal that divide n; the array is seeded with an initial value a(1)=1.

Original entry on oeis.org

1, 1, 2, 2, 2, 7, 2, 9, 10, 15, 2, 10, 1, 13, 17, 8, 0, 13, 1, 14, 9, 8, 0, 13, 3, 30, 13, 10, 2, 16, 1, 23, 5, 7, 14, 15, 2, 8, 28, 32, 2, 23, 2, 9, 49, 12, 0, 48, 2, 11, 1, 20, 3, 18, 13, 28, 0, 4, 1, 56, 5, 8, 16, 35, 46, 4, 2, 6, 2, 10
Offset: 1

Views

Author

Alec Jones, Dec 24 2016

Keywords

Comments

From Hartmut F. W. Hoft, Jan 23 2017: (Start)
Shown by induction and direct (modular) computations for
column 1: Every number is even, except for the first two 1's; in addition to row 3, value 2 occurs in rows 4*k and 4*k+1, and every value in rows 4*k+2 and 4*k+3 is divisible by 4, for all k>=1.
column 2: The first four entries, 2, 2, 9 and 10, contain the only odd number; no nonzero entry in row k>3 has 9 as a factor, and value 0 occurs in rows 4*k+1 and 4*k+2, for all k>=1.
Conjecture:
a({1, 6, 8, 9, 10, 15, 26, 45, 48, 84, 96, 112, 115, 252, 336, 343}) =
{1, 7, 9,10, 15, 17, 30, 49, 48,104,117, 115, 122, 257, 343, 395} are the only numbers in the sequence with the property a(n) >= n (verified through n=500500, i.e., the triangle with 1000 antidiagonals).
This conjecture together with Bouniakowsky's conjecture that certain quadratic integer polynomials generate infinitely many primes (e.g. see A002496 for n^2+1 and A188382 for 2*n^2+n+1) implies that in every column in the triangle infinitely many prime sequence indices occur and therefore infinitely many 0's whenever the column contains no 1's. The proof is based on the fact that for a large enough prime sequence index p in whose prior column no 1 occurs then a(p)=0; therefore infinitely many 0's occur in that column. Obviously, once value 1 occurs in a column no 0 value can occur in a subsequent row.
Conjecture:
Every row in the triangle contains exactly two 1's.
(End)

Examples

			After 6 terms, the array looks like:
.
1   2   7
1   2
2
We have a(6) = 7 because a(1) = 1, a(3) = 2, a(4) = 2, and a(5) = 2 divide 6; 1 + 2 + 2 + 2 = 7.
From _Hartmut F. W. Hoft_, Jan 23 2017: (Start)
1   2   7  15  17   9  10  15  49  13   4  31  22
1   2  10  13  14  13  14   9  18  46  12  66
2   9   1   1  30   7   2   3  35  12   3
2  10  13   3   5  23  20  16  14  17
2   0  13  23   2   1   8  11   2
8   0   1  32  11   5   3   6
8  16  28   2  56  42   8
2   8  48   1   2 104
2   0   4  10   1
12   0   2  10
28   6   2
2  42
2
.
Expanded the triangle to the first 13 antidiagonals of the array, i.e. a(1) ... a(91), to show the start of the 2- and 0-value patterns in columns 1 and 2. The first 0 beyond column 2 is a(677) in row 27, column 11 of the triangle.
A188382(n)=2*n^2+n+1 for n>=0 are the alternate sequence indices for column 1 starting in row 1, 2*n^2+n+2 for n>=1 are the alternate sequence indices for column 2 starting in row 2, and 2*n^2+n+11 for n>=5 are the alternate sequence indices for column 11 starting in row 1.
The sequence indices in the triangle for row positions k>=1 in columns 1,..., 5 are given in sequences A000124(k), A152948(k+3), A152950(k+3), A145018(k+4) and A167499(k+4).
(End)
		

Crossrefs

Cf. A279966 for the related sequence which counts prior terms.
Cf. A269347 for a one-dimensional version of this sequence.
Cf. also A279211, A279212.

Programs

  • Mathematica
    (*  printing of the triangle is commented out of function a279967[]  *)
    pCol[{i_, j_}] := Map[{#, j}&, Range[1, i-1]]
    pDiag[{i_, j_}] := If[j>=i, Map[{#, j-i+#}&, Range[1, i-1]], Map[{i-j+#, #}&, Range[1, j-1]]]
    pRow[{i_, j_}] := Map[{i, #}&, Range[1, j-1]]
    pAdiag[{i_, j_}] := Map[{i+j-#, #}&, Range[1, j-1]]
    priorPos[{i_, j_}] := Join[pCol[{i, j}], pDiag[{i, j}], pRow[{i, j}], pAdiag[{i, j}]]
    seqPos[{i_, j_}] := (i+j-2)(i+j-1)/2+j
    antiDiag[k_] := Map[{k+1-#, #}&, Range[1, k]]
    upperTriangle[k_] := Flatten[Map[antiDiag, Range[1, k]], 1]
    a279967[k_] := Module[{ut=upperTriangle[k], ms=Table[" ", {i, 1, k}, {j, 1, k}], h, pos, val, seqL={1}}, ms[[1, 1]]=1; For[h=2, h<=Length[ut], h++, pos=ut[[h]]; val=Apply[Plus, Select[Map[ms[[Apply[Sequence, #]]]&, priorPos[pos]], #!=0 && Mod[seqPos[pos], #]==0&]]; AppendTo[seqL, val]; ms[[Apply[Sequence, pos]]]=val]; (* Print[TableForm[ms]]; *) seqL]
    a279967[13] (* values in first 13 antidiagonals *)
    (* Hartmut F. W. Hoft, Jan 23 2017 *)

A166136 a(n) = n*(n+3)/2 + 7.

Original entry on oeis.org

9, 12, 16, 21, 27, 34, 42, 51, 61, 72, 84, 97, 111, 126, 142, 159, 177, 196, 216, 237, 259, 282, 306, 331, 357, 384, 412, 441, 471, 502, 534, 567, 601, 636, 672, 709, 747, 786, 826, 867, 909, 952, 996, 1041, 1087, 1134, 1182, 1231, 1281, 1332, 1384, 1437
Offset: 1

Views

Author

Vincenzo Librandi, Oct 08 2009

Keywords

Comments

Numbers m >= 9 such that 8*m - 47 is a square. - Bruce J. Nicholson, Jul 25 2017

Crossrefs

Cf. A145018.

Programs

  • Magma
    I:=[9, 12, 16]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2)+Self(n-3): n in [1..40]]; // Vincenzo Librandi, Mar 15 2012
    
  • Mathematica
    Table[n*(n+3)/2+7, {n, 1, 40}] (* or *) LinearRecurrence[{3,-3,1}, {9, 12, 16}, 40] (* Vincenzo Librandi, Mar 15 2012 *)
  • PARI
    for(n=1, 40, print1(n*(n+3)/2+7, ", ")); \\ Vincenzo Librandi, Mar 15 2012

Formula

a(n) = a(n-1) + n = 3*a(n-1) - 3*a(n-2) + a(n-3) = A145018(n+2) + 2.
G.f.: -x*(9 - 15*x + 7*x^2)/(x-1)^3.
E.g.f.: (1/2)*(14 + 4*x + x^2)*exp(x) - 7. - G. C. Greubel, Apr 26 2016
Sum_{n>=1} 1/a(n) = -13/42 + 2*Pi*tanh(sqrt(47)*Pi/2)/sqrt(47). - Amiram Eldar, Dec 13 2022

Extensions

Definition replaced by polynomial from R. J. Mathar, Oct 12 2009

A245300 Triangle T(n,k) = (n+k)*(n+k+1)/2 + k, 0 <= k <= n, read by rows.

Original entry on oeis.org

0, 1, 4, 3, 7, 12, 6, 11, 17, 24, 10, 16, 23, 31, 40, 15, 22, 30, 39, 49, 60, 21, 29, 38, 48, 59, 71, 84, 28, 37, 47, 58, 70, 83, 97, 112, 36, 46, 57, 69, 82, 96, 111, 127, 144, 45, 56, 68, 81, 95, 110, 126, 143, 161, 180, 55, 67, 80, 94, 109, 125, 142, 160, 179, 199, 220
Offset: 0

Views

Author

Reinhard Zumkeller, Jul 17 2014

Keywords

Examples

			First rows and their row sums (A245301):
   0                                                                  0;
   1,  4                                                              5;
   3,  7,  12                                                        22;
   6, 11,  17,  24                                                   58;
  10, 16,  23,  31,  40                                             120;
  15, 22,  30,  39,  49,  60                                        215;
  21, 29,  38,  48,  59,  71,  84                                   350;
  28, 37,  47,  58,  70,  83,  97, 112                              532;
  36, 46,  57,  69,  82,  96, 111, 127, 144                         768;
  45, 56,  68,  81,  95, 110, 126, 143, 161, 180                   1065;
  55, 67,  80,  94, 109, 125, 142, 160, 179, 199, 220              1430;
  66, 79,  93, 108, 124, 141, 159, 178, 198, 219, 241, 264         1870;
  78, 92, 107, 123, 140, 158, 177, 197, 218, 240, 263, 287, 312    2392.
		

Crossrefs

Programs

  • Haskell
    a245300 n k = (n + k) * (n + k + 1) `div` 2 + k
    a245300_row n = map (a245300 n) [0..n]
    a245300_tabl = map a245300_row [0..]
    a245300_list = concat a245300_tabl
    
  • Magma
    [k + Binomial(n+k+1,2): k in [0..n], n in [0..15]]; // G. C. Greubel, Apr 01 2021
    
  • Mathematica
    Table[k + Binomial[n+k+1,2], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Apr 01 2021 *)
  • Sage
    flatten([[k + binomial(n+k+1,2) for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Apr 01 2021

Formula

T(n, 0) = A000217(n).
T(n, n) = A046092(n).
T(2*n, n) = A062725(n) (central terms).
Sum_{k=0..n} T(n, k) = A245301(n).
From G. C. Greubel, Apr 01 2021: (Start)
T(n, 1) = A000124(n+1) = A134869(n+1), n >= 1.
T(n, 2) = A152948(n+4), n >= 2.
T(n, 3) = A152950(n+4), n >= 3.
T(n, 4) = A145018(n+5), n >= 4.
T(n, 5) = A167499(n+4), n >= 5.
T(n, 6) = A166136(n+5), n >= 6.
T(n, 7) = A167487(n+6), n >= 7.
T(n, n-1) = A056220(n), n >= 1.
T(n, n-2) = A142463(n-1), n >= 2.
T(n, n-3) = A054000(n-1), n >= 3.
T(n, n-4) = A090288(n-3), n >= 4.
T(n, n-5) = A268581(n-4), n >= 5.
T(n, n-6) = A059993(n-4), n >= 6.
T(n, n-7) = (-1)*A147973(n), n >= 7.
T(n, n-8) = A139570(n-5), n >= 8.
T(n, n-9) = A271625(n-5), n >= 9.
T(n, n-10) = A222182(n-4), n >= 10.
T(2*n, n-1) = A081270(n-1), n >= 1.
T(2*n, n+1) = A117625(n+1), n >= 1. (End)

A278436 Numbers such that A279966(n) = 0.

Original entry on oeis.org

17, 23, 47, 57, 93, 107, 173, 233, 353, 437, 467, 563, 677, 743, 817, 829, 851, 863, 955, 1037, 1187, 1213, 1277, 1387, 1433, 1487, 1549, 2089, 2147, 2213, 2287, 2293, 2417, 2473, 2689, 2777, 2911, 3083, 3323, 3391, 6691, 9337, 22969, 38557, 47347, 75391, 104999, 130927, 146719
Offset: 1

Views

Author

Alec Jones, Dec 24 2016

Keywords

Comments

Not all numbers in this list are prime; nonprime elements include 57, 93, 437, 817, 851, 955, 1037, 1387, 2147.
This sequence is neither a subset nor a superset of sequence A281533 since 155 and 817 are the first numbers in one, but not the other, respectively. - Hartmut F. W. Hoft, Jan 23 2017

Examples

			Number 817 = 19*43, equivalent to array position (4, 37), is in the sequence since none of the numbers in the prior column, diagonal, row and antidiagonal contain the counts of 1, 19, 43 and 817. - _Hartmut F. W. Hoft_, Jan 23 2017
		

Crossrefs

Programs

  • Mathematica
    (* support functions are in A279967 *)
    a278436[k_] := Module[{ut=upperTriangle[k], ms=Table[" ", {i, 1, k}, {j, 1, k}], h, pos, val, seqL={}}, ms[[1, 1]]=1; For[h=2, h<=Length[ut], h++, pos=ut[[h]]; val=Length[Select[Map[ms[[Apply[Sequence, #]]]&, priorPos[pos]], #!=0 && Mod[seqPos[pos], #]==0&]]; If[val==0, AppendTo[seqL, h]]; ms[[Apply[Sequence, pos]]]=val]; seqL]
    a278436[100] (* data through 3391. - Hartmut F. W. Hoft, Jan 23 2017 *)

A328000 a(n) = Sum_{k=0..n}(k!*(n - k)!)/(floor(k/2)!*floor((n - k)/2)!)^2.

Original entry on oeis.org

1, 2, 5, 16, 28, 96, 160, 512, 896, 2560, 4864, 12288, 25600, 57344, 131072, 262144, 655360, 1179648, 3211264, 5242880, 15466496, 23068672, 73400320, 100663296, 343932928, 436207616, 1593835520, 1879048192, 7314866176, 8053063680, 33285996544, 34359738368
Offset: 0

Views

Author

Peter Luschny, Oct 01 2019

Keywords

Crossrefs

Programs

  • Magma
    [IsOdd(n) select 2^(n - 1)*(n + 1) else 2^(n - 5)*(n*(n + 2) + 32):n in [0..30]]; // Marius A. Burtea, Feb 05 2020
  • Maple
    swing := n -> n!/iquo(n,2)!^2: a := n -> add(swing(k)*swing(n-k), k=0..n):
    seq(`if`(irem(n, 2) = 0, 2 + n*(n + 2)/16, n + 1)*2^(n - 1), n=0..31);
  • Mathematica
    A328000List[len_] := CoefficientList[Series[(4 x^2 - x - 1)^2 / (1 - 4 x^2)^3 , {x, 0, len}], x]; A328000List[31]
    LinearRecurrence[{0,12,0,-48,0,64},{1,2,5,16,28,96},40] (* Harvey P. Dale, Jun 19 2022 *)
  • PARI
    x='x + O('x^32);
    Vec(serlaplace(((3*x + 8)*sinh(2*x) + (2*x^2 + 16*(x + 1))*cosh(2*x))/16))
    
  • PARI
    Vec((1 + x - 4*x^2)^2 / ((1 - 2*x)^3*(1 + 2*x)^3) + O(x^30)) \\ Colin Barker, Feb 05 2020
    

Formula

a(n) = Sum_{k=0..n} s(k)*s(n-k) where s(n) = A056040(n).
a(n) = [x^n] (4*x^2 - x - 1)^2 / (1 - 4*x^2)^3.
a(n) = 2^(n - 5)*(n*(n + 2) + 32) if n even else 2^(n - 1)*(n + 1).
a(2*n) = A327999(n).
a(2*n-1) = A002699(n), (with a(-1) = 0).
a(2^n-1) = 2^(2^n - 2 + n) for n >= 1.
2*a(2*n)/2^n = A081908(n+1).
4*a(2*n)/4^n = A145018(n+1).
2*a(2*n-1)/4^n = A001477(n).
From Stefano Spezia, Oct 19 2019: (Start)
a(n) = n! [x^n] (1/32)*exp(-2*x)*(8 + exp(4*x)*(8 + x)*(3 + 2*x) + x*(13 + 2*x)).
a(n) = 12*a(n-2) - 48*a(n-4) + 64*a(n-6) for n > 5. (End)
Showing 1-10 of 11 results. Next