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-7 of 7 results.

A195040 Square array read by antidiagonals with T(n,k) = k*n^2/4+(k-4)*((-1)^n-1)/8, n>=0, k>=0.

Original entry on oeis.org

0, 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 3, 2, 1, 0, 1, 4, 5, 3, 1, 0, 0, 7, 8, 7, 4, 1, 0, 1, 9, 13, 12, 9, 5, 1, 0, 0, 13, 18, 19, 16, 11, 6, 1, 0, 1, 16, 25, 27, 25, 20, 13, 7, 1, 0, 0, 21, 32, 37, 36, 31, 24, 15, 8, 1, 0, 1, 25, 41, 48, 49, 45, 37, 28, 17, 9, 1, 0
Offset: 0

Views

Author

Omar E. Pol, Sep 27 2011

Keywords

Comments

Also, if k >= 2 and m = 2*k, then column k lists the numbers of the form k*n^2 and the centered m-gonal numbers interleaved.
For k >= 3, this is also a table of concentric polygonal numbers. Column k lists the concentric k-gonal numbers.
It appears that the first differences of column k are the numbers that are congruent to {1, k-1} mod k, if k >= 3.

Examples

			Array begins:
  0,   0,   0,   0,   0,   0,   0,   0,   0,   0, ...
  1,   1,   1,   1,   1,   1,   1,   1,   1,   1, ...
  0,   1,   2,   3,   4,   5,   6,   7,   8,   9, ...
  1,   3,   5,   7,   9,  11,  13,  15,  17,  19, ...
  0,   4,   8,  12,  16,  20,  24,  28,  32,  36, ...
  1,   7,  13,  19,  25,  31,  37,  43,  49,  55, ...
  0,   9,  18,  27,  36,  45,  54,  63,  72,  81, ...
  1,  13,  25,  37,  49,  61,  73,  85,  97, 109, ...
  0,  16,  32,  48,  64,  80,  96, 112, 128, 144, ...
  1,  21,  41,  61,  81, 101, 121, 141, 161, 181, ...
  0,  25,  50,  75, 100, 125, 150, 175, 200, 225, ...
  ...
		

Crossrefs

Rows n: A000004 (n=0), A000012 (n=1), A001477 (n=2), A005408 (n=3), A008586 (n=4), A016921 (n=5), A008591 (n=6), A017533 (n=7), A008598 (n=8), A215145 (n=9), A008607 (n=10).
Columns k: A000035 (k=0), A004652 (k=1), A000982 (k=2), A077043 (k=3), A000290 (k=4), A032527 (k=5), A032528 (k=6), A195041 (k=7), A077221 (k=8), A195042 (k=9), A195142 (k=10), A195043 (k=11), A195143 (k=12), A195045 (k=13), A195145 (k=14), A195046 (k=15), A195146 (k=16), A195047 (k=17), A195147 (k=18), A195048 (k=19), A195148 (k=20), A195049 (k=21), A195149 (k=22), A195058 (k=23), A195158 (k=24).

Programs

  • GAP
    nmax:=13;; T:=List([0..nmax],n->List([0..nmax],k->k*n^2/4+(k-4)*((-1)^n-1)/8));; b:=List([2..nmax],n->OrderedPartitions(n,2));;
    a:=Flat(List([1..Length(b)],i->List([1..Length(b[i])],j->T[b[i][j][2]][b[i][j][1]]))); # Muniru A Asiru, Jul 19 2018
  • Maple
    A195040 := proc(n,k)
            k*n^2/4+((-1)^n-1)*(k-4)/8 ;
    end proc:
    for d from 0 to 12 do
            for k from 0 to d do
                    printf("%d,",A195040(d-k,k)) ;
            end do:
    end do; # R. J. Mathar, Sep 28 2011
  • Mathematica
    t[n_, k_] := k*n^2/4+(k-4)*((-1)^n-1)/8; Flatten[ Table[ t[n-k, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, Dec 14 2011 *)

A195143 a(n) = n-th concentric 12-gonal number.

Original entry on oeis.org

0, 1, 12, 25, 48, 73, 108, 145, 192, 241, 300, 361, 432, 505, 588, 673, 768, 865, 972, 1081, 1200, 1321, 1452, 1585, 1728, 1873, 2028, 2185, 2352, 2521, 2700, 2881, 3072, 3265, 3468, 3673, 3888, 4105, 4332, 4561, 4800, 5041, 5292, 5545, 5808, 6073, 6348
Offset: 0

Views

Author

Omar E. Pol, Sep 17 2011

Keywords

Comments

Concentric dodecagonal numbers. [corrected by Ivan Panchenko, Nov 09 2013]
Sequence found by reading the line from 0, in the direction 0, 12,..., and the same line from 1, in the direction 1, 25,..., in the square spiral whose vertices are the generalized octagonal numbers A001082. Main axis, perpendicular to A028896 in the same spiral.
Partial sums of A091998. - Reinhard Zumkeller, Jan 07 2012
Column 12 of A195040. - Omar E. Pol, Sep 28 2011

Crossrefs

A135453 and A069190 interleaved.
Cf. A016921 (6n+1), A016969 (6n+5), A091998 (positive integers of the form 12*k +- 1), A092242 (positive integers of the form 12*k +- 5).

Programs

  • Haskell
    a195143 n = a195143_list !! n
    a195143_list = scanl (+) 0 a091998_list
    -- Reinhard Zumkeller, Jan 07 2012
  • Magma
    [(3*n^2+(-1)^n-1): n in [0..50]]; // Vincenzo Librandi, Sep 27 2011
    
  • Mathematica
    Table[Sum[2*(-1)^(n - k + 1) + 6*k - 3, {k, n}], {n, 0, 47}] (* L. Edson Jeffery, Sep 14 2014 *)

Formula

From Vincenzo Librandi, Sep 27 2011: (Start)
a(n) = 3*n^2+(-1)^n-1.
a(n) = -a(n-1) + 6*n^2 - 6*n + 1. (End)
G.f.: -x*(1+10*x+x^2) / ( (1+x)*(x-1)^3 ). - R. J. Mathar, Sep 18 2011
a(n) = Sum_{k=1..n} (2*(-1)^(n-k+1) + 3*(2*k-1)), n>0, a(0) = 0. - L. Edson Jeffery, Sep 14 2014
Sum_{n>=1} 1/a(n) = Pi^2/72 + tan(Pi/sqrt(6))*Pi/(4*sqrt(6)). - Amiram Eldar, Jan 16 2023

A195145 Concentric 14-gonal numbers.

Original entry on oeis.org

0, 1, 14, 29, 56, 85, 126, 169, 224, 281, 350, 421, 504, 589, 686, 785, 896, 1009, 1134, 1261, 1400, 1541, 1694, 1849, 2016, 2185, 2366, 2549, 2744, 2941, 3150, 3361, 3584, 3809, 4046, 4285, 4536, 4789, 5054, 5321, 5600, 5881, 6174, 6469, 6776, 7085, 7406
Offset: 0

Views

Author

Omar E. Pol, Sep 17 2011

Keywords

Comments

Also concentric tetradecagonal numbers or concentric tetrakaidecagonal numbers. Also sequence found by reading the line from 0, in the direction 0, 14, ..., and the same line from 1, in the direction 1, 29, ..., in the square spiral whose vertices are the generalized enneagonal numbers A118277. Main axis, perpendicular to A024966 in the same spiral.
Partial sums of A113801. - Reinhard Zumkeller, Jan 07 2012

Crossrefs

Programs

  • Haskell
    a195145 n = a195145_list !! n
    a195145_list = scanl (+) 0 a113801_list
    -- Reinhard Zumkeller, Jan 07 2012
  • Magma
    [(14*n^2+5*(-1)^n-5)/4: n in [0..50]]; // Vincenzo Librandi, Sep 27 2011
    
  • Mathematica
    LinearRecurrence[{2, 0, -2, 1}, {0, 1, 14, 29}, 50] (* Amiram Eldar, Jan 16 2023 *)

Formula

G.f.: -x*(1+12*x+x^2) / ( (1+x)*(x-1)^3 ). - R. J. Mathar, Sep 18 2011
From Vincenzo Librandi, Sep 27 2011: (Start)
a(n) = (14*n^2 + 5*(-1)^n - 5)/4;
a(n) = a(-n) = -a(n-1) + 7*n^2 - 7*n + 1. (End)
Sum_{n>=1} 1/a(n) = Pi^2/84 + tan(sqrt(5/7)*Pi/2)*Pi/(2*sqrt(35)). - Amiram Eldar, Jan 16 2023
E.g.f.: (7*x*(x + 1)*cosh(x) + (7*x^2 + 7*x - 5)*sinh(x))/2. - Stefano Spezia, Nov 30 2024

A175886 Numbers that are congruent to {1, 12} mod 13.

Original entry on oeis.org

1, 12, 14, 25, 27, 38, 40, 51, 53, 64, 66, 77, 79, 90, 92, 103, 105, 116, 118, 129, 131, 142, 144, 155, 157, 168, 170, 181, 183, 194, 196, 207, 209, 220, 222, 233, 235, 246, 248, 259, 261, 272, 274, 285, 287, 298, 300, 311, 313, 324, 326, 337, 339, 350
Offset: 1

Views

Author

Bruno Berselli, Oct 08 2010 - Nov 17 2010

Keywords

Comments

Cf. property described by Gary Detlefs in A113801: more generally, these numbers are of the form (2*h*n+(h-4)*(-1)^n-h)/4 (h, n natural numbers), therefore ((2*h*n+(h-4)*(-1)^n-h)/4)^2-1 == 0 (mod h); in this case, a(n)^2-1 == 0 (mod 13).

Crossrefs

Programs

  • Haskell
    a175886 n = a175886_list !! (n-1)
    a175886_list = 1 : 12 : map (+ 13) a175886_list
    -- Reinhard Zumkeller, Jan 07 2012
    
  • Magma
    [n: n in [1..350] | n mod 13 in [1, 12]]; // Bruno Berselli, Feb 29 2012
    
  • Magma
    [(26*n+9*(-1)^n-13)/4: n in [1..55]]; // Vincenzo Librandi, Aug 19 2013
    
  • Mathematica
    Select[Range[1, 350], MemberQ[{1, 12}, Mod[#, 13]]&] (* Bruno Berselli, Feb 29 2012 *)
    CoefficientList[Series[(1 + 11 x + x^2) / ((1 + x) (1 - x)^2), {x, 0, 55}], x] (* Vincenzo Librandi, Aug 19 2013 *)
    LinearRecurrence[{1,1,-1},{1,12,14},60] (* Harvey P. Dale, Oct 23 2015 *)
  • PARI
    a(n)=(26*n+9*(-1)^n-13)/4 \\ Charles R Greathouse IV, Sep 24 2015

Formula

G.f.: x*(1+11*x+x^2)/((1+x)*(1-x)^2).
a(n) = (26*n+9*(-1)^n-13)/4.
a(n) = -a(-n+1) = a(n-1)+a(n-2)-a(n-3).
a(n) = a(n-2)+13.
a(n) = 13*A000217(n-1)+1 - 2*Sum_{i=1..n-1} a(i) for n>1.
Sum_{n>=1} (-1)^(n+1)/a(n) = (Pi/13)*cot(Pi/13). - Amiram Eldar, Dec 04 2021
E.g.f.: 1 + ((26*x - 13)*exp(x) + 9*exp(-x))/4. - David Lovler, Sep 04 2022
From Amiram Eldar, Nov 25 2024: (Start)
Product_{n>=1} (1 - (-1)^n/a(n)) = 2*cos(Pi/13).
Product_{n>=2} (1 + (-1)^n/a(n)) = (Pi/13)*cosec(Pi/13). (End)

A195043 Concentric 11-gonal numbers.

Original entry on oeis.org

0, 1, 11, 23, 44, 67, 99, 133, 176, 221, 275, 331, 396, 463, 539, 617, 704, 793, 891, 991, 1100, 1211, 1331, 1453, 1584, 1717, 1859, 2003, 2156, 2311, 2475, 2641, 2816, 2993, 3179, 3367, 3564, 3763, 3971, 4181, 4400, 4621, 4851, 5083, 5324, 5567
Offset: 0

Views

Author

Omar E. Pol, Sep 27 2011

Keywords

Comments

Also concentric hendecagonal numbers. A033584 and A069173 interleaved.
Partial sums of A175885. - Reinhard Zumkeller, Jan 07 2012

Crossrefs

Programs

  • Haskell
    a195043 n = a195043_list !! n
    a195043_list = scanl (+) 0 a175885_list
    -- Reinhard Zumkeller, Jan 07 2012
    
  • Magma
    [11*n^2/4+7*((-1)^n-1)/8: n in [0..50]]; // Vincenzo Librandi, Sep 30 2011
    
  • Mathematica
    LinearRecurrence[{2,0,-2,1},{0,1,11,23},50] (* Harvey P. Dale, May 20 2019 *)
  • PARI
    Vec(-x*(x^2+9*x+1)/((x-1)^3*(x+1)) + O(x^100)) \\ Colin Barker, Sep 15 2013

Formula

a(n) = 11*n^2/4 + 7*((-1)^n - 1)/8.
a(n) = -a(n-1) + A069125(n). - Vincenzo Librandi, Sep 30 2011
From Colin Barker, Sep 15 2013: (Start)
a(n) = 2*a(n-1) - 2*a(n-3) + a(n-4).
G.f.: -x*(x^2+9*x+1) / ((x-1)^3*(x+1)). (End)
Sum_{n>=1} 1/a(n) = Pi^2/66 + tan(sqrt(7/11)*Pi/2)*Pi/sqrt(77). - Amiram Eldar, Jan 16 2023

A195046 Concentric 15-gonal numbers.

Original entry on oeis.org

0, 1, 15, 31, 60, 91, 135, 181, 240, 301, 375, 451, 540, 631, 735, 841, 960, 1081, 1215, 1351, 1500, 1651, 1815, 1981, 2160, 2341, 2535, 2731, 2940, 3151, 3375, 3601, 3840, 4081, 4335, 4591, 4860, 5131, 5415, 5701, 6000, 6301, 6615, 6931, 7260, 7591
Offset: 0

Views

Author

Omar E. Pol, Sep 27 2011

Keywords

Crossrefs

Programs

  • Mathematica
    Table[15n^2/4+11((-1)^n-1)/8,{n,0,50}] (* or *) LinearRecurrence[ {2,0,-2,1},{0,1,15,31},50] (* Harvey P. Dale, Feb 23 2012 *)
  • PARI
    a(n)=15*n^2/4+11*((-1)^n-1)/8 \\ Charles R Greathouse IV, Oct 07 2015

Formula

a(n) = 15*n^2/4+11*((-1)^n-1)/8.
From Harvey P. Dale, Feb 23 2012: (Start)
a(0)=0, a(1)=1, a(2)=15, a(3)=31, a(n)=2*a(n-1)-2*a(n-3)+a(n-4).
G.f.: -((x*(1+x*(13+x)))/((-1+x)^3*(1+x))). (End)
Sum_{n>=1} 1/a(n) = Pi^2/90 + tan(sqrt(11/15)*Pi/2)*Pi/sqrt(165). - Amiram Eldar, Jan 16 2023

Extensions

a(1)=1 added by Harvey P. Dale, Feb 23 2012

A270693 Alternating sum of centered 25-gonal numbers.

Original entry on oeis.org

1, -25, 51, -100, 151, -225, 301, -400, 501, -625, 751, -900, 1051, -1225, 1401, -1600, 1801, -2025, 2251, -2500, 2751, -3025, 3301, -3600, 3901, -4225, 4551, -4900, 5251, -5625, 6001, -6400, 6801, -7225, 7651, -8100, 8551, -9025, 9501, -10000, 10501
Offset: 0

Views

Author

Ilya Gutkovskiy, Mar 21 2016

Keywords

Comments

The absolute value alternating sum of centered k-gonal numbers gives concentric k-gonal numbers.
More generally, the ordinary generating function for the alternating sum of centered k-gonal numbers is (1 - (k - 2)*x + x^2)/((1 - x)*(1 + x)^3).

Crossrefs

Programs

  • Magma
    [((-1)^n*(50*n^2 + 100*n + 29) - 21)/8 : n in [0..40]]; // Wesley Ivan Hurt, Mar 21 2016
  • Maple
    A270693:=n->((-1)^n*(50*n^2 + 100*n + 29) - 21)/8: seq(A270693(n), n=0..100); # Wesley Ivan Hurt, Sep 18 2017
  • Mathematica
    LinearRecurrence[{-2, 0, 2, 1}, {1, -25, 51, -100}, 41]
    Table[((-1)^n (50 n^2 + 100 n + 29) - 21)/8, {n, 0, 40}]
  • PARI
    x='x+O('x^100); Vec((1-23*x+x^2)/((1-x)*(1+x)^3)) \\ Altug Alkan, Mar 21 2016
    

Formula

G.f.: (1 - 23*x + x^2)/((1 - x)*(1 + x)^3).
E.g.f.: (1/8)*(-21*exp(x) + (29 - 150*x + 50*x^2)*exp(-x)).
a(n) = -2*a(n-1) + 2*a(n-3) + a(n-4).
a(n) = ((-1)^n*(50*n^2 + 100*n + 29) - 21)/8.
Showing 1-7 of 7 results.