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.

Previous Showing 21-30 of 41 results. Next

A067726 a(n) = 6*n^2 + 12*n.

Original entry on oeis.org

18, 48, 90, 144, 210, 288, 378, 480, 594, 720, 858, 1008, 1170, 1344, 1530, 1728, 1938, 2160, 2394, 2640, 2898, 3168, 3450, 3744, 4050, 4368, 4698, 5040, 5394, 5760, 6138, 6528, 6930, 7344, 7770, 8208, 8658, 9120, 9594, 10080, 10578, 11088, 11610
Offset: 1

Views

Author

Robert G. Wilson v, Feb 05 2002

Keywords

Comments

Positive numbers k such that 6*(6 + k) is a perfect square.

Crossrefs

Cf. numbers k such that k*(k + m) is a perfect square: A028560 (k=9), A067728 (k=8), A067727 (k=7), A067724 (k=5), A028347 (k=4), A067725 (k=3), A054000 (k=2), A005563 (k=1).

Programs

Formula

G.f.: 6*x*(3 - x)/(1 - x)^3. - Vincenzo Librandi, Jul 08 2012
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Vincenzo Librandi, Jul 08 2012
E.g.f.: 6*x*(3 + x)*exp(x). - G. C. Greubel, Sep 01 2019
From Amiram Eldar, Feb 25 2022: (Start)
Sum_{n>=1} 1/a(n) = 1/8.
Sum_{n>=1} (-1)^(n+1)/a(n) = 1/24. (End)
a(n) = A003215(2*n) - A003154(n). - Leo Tavares, May 20 2023
a(n) = 6*A005563(n). - Hugo Pfoertner, May 24 2023

A098832 Square array read by antidiagonals: even-numbered rows of the table are of the form n*(n+m) and odd-numbered rows are of the form n*(n+m)/2.

Original entry on oeis.org

1, 3, 3, 6, 8, 2, 10, 15, 5, 5, 15, 24, 9, 12, 3, 21, 35, 14, 21, 7, 7, 28, 48, 20, 32, 12, 16, 4, 36, 63, 27, 45, 18, 27, 9, 9, 45, 80, 35, 60, 25, 40, 15, 20, 5, 55, 99, 44, 77, 33, 55, 22, 33, 11, 11, 66, 120, 54, 96, 42, 72, 30, 48, 18, 24, 6, 78, 143, 65, 117, 52, 91, 39, 65, 26, 39, 13, 13
Offset: 1

Views

Author

Eugene McDonnell (eemcd(AT)mac.com), Nov 02 2004

Keywords

Comments

The rows of this table and that in A098737 are related. Given a function f = n/( 1 + (1+n) mod(2) ), row n of A098737 can be derived from row n of T by multiplying the latter by f(n); row n of T can be derived from row n of A098737 by dividing the latter by f(n).

Examples

			Array begins as:
  1,  3,  6, 10, 15, 21,  28,  36,  45 ... A000217;
  3,  8, 15, 24, 35, 48,  63,  80,  99 ... A005563;
  2,  5,  9, 14, 20, 27,  35,  44,  54 ... A000096;
  5, 12, 21, 32, 45, 60,  77,  96, 117 ... A028347;
  3,  7, 12, 18, 25, 33,  42,  52,  63 ... A027379;
  7, 16, 27, 40, 55, 72,  91, 112, 135 ... A028560;
  4,  9, 15, 22, 30, 39,  49,  60,  72 ... A055999;
  9, 20, 33, 48, 65, 84, 105, 128, 153 ... A028566;
  5, 11, 18, 26, 35, 45,  56,  68,  81 ... A056000;
Antidiagonals begin as:
   1;
   3,  3;
   6,  8,  2;
  10, 15,  5,  5;
  15, 24,  9, 12,  3;
  21, 35, 14, 21,  7,  7;
  28, 48, 20, 32, 12, 16,  4;
  36, 63, 27, 45, 18, 27,  9,  9;
  45, 80, 35, 60, 25, 40, 15, 20,  5;
  55, 99, 44, 77, 33, 55, 22, 33, 11, 11;
		

Crossrefs

Row m of array: A000217 (m=1), A005563 (m=2), A000096 (m=3), A028347 (m=4), A027379 (m=5), A028560 (m=6), A055999 (m=7), A028566 (m=8), A056000 (m=9), A098603 (m=10), A056115 (m=11), A098847 (m=12), A056119 (m=13), A098848 (m=14), A056121 (m=15), A098849 (m=16), A056126 (m=17), A098850 (m=18), A051942 (m=19).
Column m of array: A026741 (m=1), A022998 (m=2), A165351 (m=3).

Programs

  • Magma
    A098832:= func< n,k | (1/4)*(3+(-1)^k)*(n+1)*(n-k+1) >;
    [A098832(n,k): k in [1..n], n in [1..15]]; // G. C. Greubel, Jul 31 2022
    
  • Mathematica
    A098832[n_, k_]:= (1/4)*(3+(-1)^k)*(n+1)*(n-k+1);
    Table[A098832[n,k], {n,15}, {k,n}]//Flatten (* G. C. Greubel, Jul 31 2022 *)
  • SageMath
    def A098832(n,k): return (1/4)*(3+(-1)^k)*(n+1)*(n-k+1)
    flatten([[A098832(n,k) for k in (1..n)] for n in (1..15)]) # G. C. Greubel, Jul 31 2022

Formula

Item m of row n of T is given (in infix form) by: n T m = n * (n + m) / (1 + m (mod 2)). E.g. Item 4 of row 3 of T: 3 T 4 = 14.
From G. C. Greubel, Jul 31 2022: (Start)
A(n, k) = (1/4)*(3 + (-1)^n)*k*(k+n) (array).
T(n, k) = (1/4)*(3 + (-1)^k)*(n+1)*(n-k+1) (antidiagonal triangle).
Sum_{k=1..n} T(n, k) = (1/8)*(n+1)*( (3*n-1)*(n+1) + (1+(-1)^n)/2 ).
T(2*n-1, n) = A181900(n).
T(2*n+1, n) = 2*A168509(n+1). (End)

Extensions

Missing terms added by G. C. Greubel, Jul 31 2022

A132768 a(n) = n*(n + 26).

Original entry on oeis.org

0, 27, 56, 87, 120, 155, 192, 231, 272, 315, 360, 407, 456, 507, 560, 615, 672, 731, 792, 855, 920, 987, 1056, 1127, 1200, 1275, 1352, 1431, 1512, 1595, 1680, 1767, 1856, 1947, 2040, 2135, 2232, 2331, 2432, 2535, 2640, 2747, 2856, 2967, 3080, 3195, 3312, 3431
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = n*(n + 26).
a(n) = 2*n + a(n-1) + 25, with a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(26)/26 = A001008(26)/A102928(26) = 34395742267/232016584800, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 18051406831/696049754400. (End)
From G. C. Greubel, Mar 13 2022: (Start)
G.f.: x*(27 - 25*x)/(1-x)^3.
E.g.f.: x*(27 + x)*exp(x). (End)

A132769 a(n) = n*(n + 27).

Original entry on oeis.org

0, 28, 58, 90, 124, 160, 198, 238, 280, 324, 370, 418, 468, 520, 574, 630, 688, 748, 810, 874, 940, 1008, 1078, 1150, 1224, 1300, 1378, 1458, 1540, 1624, 1710, 1798, 1888, 1980, 2074, 2170, 2268, 2368, 2470, 2574, 2680, 2788, 2898, 3010, 3124, 3240, 3358, 3478
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = 2*n + a(n-1) + 26, with a(0)=0. - Vincenzo Librandi, Aug 03 2010
a(0)=0, a(1)=28, a(2)=58; for n > 2, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Oct 14 2012
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(27)/27 = A001008(27)/A102928(27) = 312536252003/2168462696400, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/27 - 57128792093/2168462696400. (End)
From Elmo R. Oliveira, Nov 29 2024: (Start)
G.f.: 2*x*(14 - 13*x)/(1 - x)^3.
E.g.f.: exp(x)*x*(28 + x).
a(n) = 2*A132756(n). (End)

A020742 Pisot sequence T(7,9).

Original entry on oeis.org

7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29, 31, 33, 35, 37, 39, 41, 43, 45, 47, 49, 51, 53, 55, 57, 59, 61, 63, 65, 67, 69, 71, 73, 75, 77, 79, 81, 83, 85, 87, 89, 91, 93, 95, 97, 99, 101, 103, 105, 107, 109, 111, 113, 115, 117, 119, 121, 123, 125, 127, 129, 131, 133, 135, 137, 139
Offset: 0

Views

Author

Keywords

Crossrefs

Subsequence of A005408, A020735. See A008776 for definitions of Pisot sequences.
Essentially the same as A005818.

Programs

  • Mathematica
    T[x_, y_, z_] := Block[{a}, a[0] = x; a[1] = y; a[n_] := a[n] = Floor[a[n - 1]^2/a[n - 2]]; Table[a[n], {n, 0, z}]]; T[7, 9, 66] (* Michael De Vlieger, Aug 08 2016 *)
  • PARI
    pisotT(nmax, a1, a2) = {
      a=vector(nmax); a[1]=a1; a[2]=a2;
      for(n=3, nmax, a[n] = floor(a[n-1]^2/a[n-2]));
      a
    }
    pisotT(50, 7, 9) \\ Colin Barker, Aug 08 2016

Formula

a(n) = 2*n + 7.
a(n) = 2*a(n-1) - a(n-2).
From Elmo R. Oliveira, Oct 30 2024: (Start)
G.f.: (7 - 5*x)/(1 - x)^2.
E.g.f.: (7 + 2*x)*exp(x).
a(n) = A016825(n+3)/2 = A028560(n+1) - A028560(n). (End)

A105020 Array read by antidiagonals: row n (n >= 0) contains the numbers m^2 - n^2, m >= n+1.

Original entry on oeis.org

1, 3, 4, 5, 8, 9, 7, 12, 15, 16, 9, 16, 21, 24, 25, 11, 20, 27, 32, 35, 36, 13, 24, 33, 40, 45, 48, 49, 15, 28, 39, 48, 55, 60, 63, 64, 17, 32, 45, 56, 65, 72, 77, 80, 81, 19, 36, 51, 64, 75, 84, 91, 96, 99, 100, 21, 40, 57, 72, 85, 96, 105, 112, 117, 120, 121
Offset: 0

Views

Author

Keywords

Comments

A "Goldbach Conjecture" for this sequence: when there are n terms between consecutive odd integers (2n+1) and (2n+3) for n > 0, at least one will be the product of 2 primes (not necessarily distinct). Example: n=3 for consecutive odd integers a(7)=7 and a(11)=9 and of the 3 sequence entries a(8)=12, a(9)=15 and a(10)=16 between them, one is the product of 2 primes a(9)=15=3*5. - Michael Hiebl, Jul 15 2007
A024352 gives distinct values in the array, minus the first row (1, 4, 9, 16, etc.). a(n) gives all solutions to the equation x^2 + xy = n, with y mod 2 = 0, x > 0, y >= 0. - Andrew S. Plewe, Oct 19 2007
Alternatively, triangular sequence of coefficients of Dynkin diagram weights for the Cartan groups C_n: t(n,m) = m*(2*n - m). Row sums are A002412. - Roger L. Bagula, Aug 05 2008

Examples

			Array begins:
  1  4  9 16 25 36  49  64  81 100 ...
  3  8 15 24 35 48  63  80  99 120 ...
  5 12 21 32 45 60  77  96 117 140 ...
  7 16 27 40 55 72  91 112 135 160 ...
  9 20 33 48 65 84 105 128 153 180 ...
  ...
Triangle begins:
   1;
   3,  4;
   5,  8,  9;
   7, 12, 15, 16;
   9, 16, 21, 24, 25;
  11, 20, 27, 32, 35, 36;
  13, 24, 33, 40, 45, 48, 49;
  15, 28, 39, 48, 55, 60, 63, 64;
  17, 32, 45, 56, 65, 72, 77, 80, 81;
  19, 36, 51, 64, 75, 84, 91, 96, 99, 100;
		

References

  • R. N. Cahn, Semi-Simple Lie Algebras and Their Representations, Dover, NY, 2006, ISBN 0-486-44999-8, p. 139.

Crossrefs

Programs

  • Magma
    [(k+1)*(2*n-k+1): k in [0..n], n in [0..15]]; // G. C. Greubel, Mar 15 2023
    
  • Mathematica
    t[n_, m_]:= (n^2 - m^2); Flatten[Table[t[i, j], {i,12}, {j,i-1,0,-1}]]
    (* to view table *) Table[t[i, j], {j,0,6}, {i,j+1,10}]//TableForm (* Robert G. Wilson v, Jul 11 2005 *)
    Table[(k+1)*(2*n-k+1), {n,0,15}, {k,0,n}]//Flatten (* Roger L. Bagula, Aug 05 2008 *)
  • SageMath
    def A105020(n,k): return (k+1)*(2*n-k+1)
    flatten([[A105020(n,k) for k in range(n+1)] for n in range(16)]) # G. C. Greubel, Mar 15 2023

Formula

a(n) = r^2 - (r^2 + r - m)^2/4, where r = round(sqrt(m)) and m = 2*n+2. - Wesley Ivan Hurt, Sep 04 2021
a(n) = A128076(n+1) * A105020(n+1). - Wesley Ivan Hurt, Jan 07 2022
From G. C. Greubel, Mar 15 2023: (Start)
Sum_{k=0..n} T(n, k) = A002412(n+1).
Sum_{k=0..n} (-1)^k*T(n, k) = (1/2)*((1+(-1)^n)*A000384((n+2)/2) - (1- (-1)^n)*A000384((n+1)/2)). (End)

Extensions

More terms from Robert G. Wilson v, Jul 11 2005

A164007 Zero together with row 7 of the array in A163280.

Original entry on oeis.org

0, 13, 26, 33, 52, 55, 78, 91, 112, 135, 160, 187, 216, 247, 280, 315, 352, 391, 432, 475, 520, 567, 616, 667, 720, 775, 832, 891, 952, 1015, 1080, 1147, 1216, 1287, 1360, 1435, 1512, 1591, 1672, 1755, 1840, 1927, 2016, 2107, 2200, 2295, 2392, 2491, 2592, 2695
Offset: 0

Views

Author

Omar E. Pol, Aug 08 2009

Keywords

Crossrefs

Programs

  • Maple
    A033676 := proc(n) local a,d; a := 0 ; for d in numtheory[divisors](n) do if d^2 <= n then a := max(a,d) ; fi; od: a; end: A163280 := proc(n,k) local r,T ; r := 0 ; for T from k^2 by k do if A033676(T) = k then r := r+1 ; if r = n then RETURN(T) ; fi; fi; od: end: A164007 := proc(n) if n = 0 then 0; else A163280(7,n) ; fi; end: seq(A164007(n),n=0..80) ;  # R. J. Mathar, Aug 09 2009
  • Mathematica
    Join[{0, 13, 26, 33, 52, 55, 78}, Table[n*(n + 6), {n, 7, 50}]] (* G. C. Greubel, Aug 28 2017 *)
    LinearRecurrence[{3,-3,1},{0,13,26,33,52,55,78,91,112,135},50] (* Harvey P. Dale, Jul 03 2020 *)
  • PARI
    my(x='x+O('x^50)); concat([0], Vec(x*(13 - 13*x - 6*x^2 + 18*x^3 - 28*x^4 + 36*x^5 - 30*x^6 + 18*x^7 - 6*x^8)/(1 - x)^3)) \\ G. C. Greubel, Aug 28 2017

Formula

From G. C. Greubel, Aug 28 2017: (Start)
a(n) = n*(n+6), n >= 7.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3), n >= 7.
G.f.: x*(13 - 13*x - 6*x^2 + 18*x^3 - 28*x^4 + 36*x^5 - 30*x^6 + 18*x^7 - 6*x^8)/(1 - x)^3.
E.g.f.: (7*x + x^2)*exp(x) + 6*x +5*x^2 + x^3 + x^4/2 + x^6/120. (End)

Extensions

Extended by R. J. Mathar, Aug 09 2009

A132770 a(n) = n*(n + 28).

Original entry on oeis.org

0, 29, 60, 93, 128, 165, 204, 245, 288, 333, 380, 429, 480, 533, 588, 645, 704, 765, 828, 893, 960, 1029, 1100, 1173, 1248, 1325, 1404, 1485, 1568, 1653, 1740, 1829, 1920, 2013, 2108, 2205, 2304, 2405, 2508, 2613, 2720, 2829, 2940, 3053, 3168, 3285, 3404, 3525
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = 2*n + a(n-1) + 27, with a(0)=0. - Vincenzo Librandi, Aug 03 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(28)/28 = A001008(28)/A102928(28) = 315404588903/2248776129600, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 7751493599/321253732800. (End)
G.f.: x*(29 - 27*x)/(1-x)^3. - Harvey P. Dale, Aug 03 2021
E.g.f.: x*(29 + x)*exp(x). - G. C. Greubel, Mar 13 2022

A132771 a(n) = n*(n + 29).

Original entry on oeis.org

0, 30, 62, 96, 132, 170, 210, 252, 296, 342, 390, 440, 492, 546, 602, 660, 720, 782, 846, 912, 980, 1050, 1122, 1196, 1272, 1350, 1430, 1512, 1596, 1682, 1770, 1860, 1952, 2046, 2142, 2240, 2340, 2442, 2546, 2652, 2760, 2870, 2982, 3096, 3212, 3330, 3450, 3572
Offset: 0

Views

Author

Omar E. Pol, Aug 28 2007

Keywords

Crossrefs

Programs

Formula

a(n) = 2*n + a(n-1) + 28 (with a(0)=0). - Vincenzo Librandi, Aug 03 2010
From Amiram Eldar, Jan 16 2021: (Start)
Sum_{n>=1} 1/a(n) = H(29)/29 = A001008(29)/A102928(29) = 9227046511387/67543597321200, where H(k) is the k-th harmonic number.
Sum_{n>=1} (-1)^(n+1)/a(n) = 2*log(2)/29 - 236266661971/9649085331600. (End)
From G. C. Greubel, Mar 13 2022: (Start)
G.f.: 2*(15*x - 14*x^2)/(1-x)^3.
E.g.f.: x*(30 + x)*exp(x). (End)

A001533 a(n) = (8*n+1)*(8*n+7).

Original entry on oeis.org

7, 135, 391, 775, 1287, 1927, 2695, 3591, 4615, 5767, 7047, 8455, 9991, 11655, 13447, 15367, 17415, 19591, 21895, 24327, 26887, 29575, 32391, 35335, 38407, 41607, 44935, 48391, 51975, 55687, 59527, 63495, 67591, 71815, 76167, 80647, 85255, 89991, 94855, 99847, 104967
Offset: 0

Views

Author

Keywords

Comments

From Klaus Purath, Aug 18 2022: (Start)
This is A028560(8*n+1), and thus a(n) + 9 is a square. (See formulas)
7 is the only prime number of this sequence in which all odd prime factors occur.
Each prime factor p appears exactly twice in any interval of p consecutive terms. If a(m) and a(n) are within such an interval containing p, then m + n == -1 (mod p). (End)

Crossrefs

Programs

Formula

a(n) = 4*A001539(n) - 5.
a(n) = 128*n + a(n-1) with a(0)=7. - Vincenzo Librandi, Nov 12 2010
Sum_{n>=0} 1/a(n) = (Psi(7/8)-Psi(1/8))/48 = 0.1580099..., see A250129. - R. J. Mathar, May 30 2022 [ = (sqrt(2)+1)*Pi/48. - Amiram Eldar, Sep 08 2022]
From Klaus Purath, Aug 18 2022: (Start)
a(n) = A028560(8*n+1).
a(n) + 9 = ((a(n+1) - a(n-1))/32)^2 = A017113(n)^2.
a(2*n) = (a(n+1) - a(n-1))*n + 7. (End)
From Amiram Eldar, Feb 19 2023: (Start)
a(n) = A017077(n)*A004771(n).
Sum_{n>=0} (-1)^n/a(n) = (cos(Pi/8) * log(cot(Pi/16)) + sin(Pi/8) * log(cot(3*Pi/16)))/12.
Product_{n>=0} (1 - 1/a(n)) = cosec(Pi/8)*cos(sqrt(5/2)*Pi/4).
Product_{n>=0} (1 + 1/a(n)) = cosec(Pi/8)*cos(sqrt(2)*Pi/4). (End)
G.f.: -(7+114*x+7*x^2)/(x-1)^3. - R. J. Mathar, Apr 23 2024
From Elmo R. Oliveira, Oct 25 2024: (Start)
E.g.f.: exp(x)*(7 + 64*x*(2 + x)).
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n > 2. (End)
Previous Showing 21-30 of 41 results. Next