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

A033428 a(n) = 3*n^2.

Original entry on oeis.org

0, 3, 12, 27, 48, 75, 108, 147, 192, 243, 300, 363, 432, 507, 588, 675, 768, 867, 972, 1083, 1200, 1323, 1452, 1587, 1728, 1875, 2028, 2187, 2352, 2523, 2700, 2883, 3072, 3267, 3468, 3675, 3888, 4107, 4332, 4563, 4800, 5043, 5292, 5547, 5808, 6075, 6348
Offset: 0

Views

Author

Keywords

Comments

The number of edges of a complete tripartite graph of order 3n, K_n,n,n. - Roberto E. Martinez II, Oct 18 2001
From Floor van Lamoen, Jul 21 2001: (Start)
Write 1,2,3,4,... in a hexagonal spiral around 0; then a(n) is the sequence found by reading the line from 0 in the direction 0,3,.... The spiral begins:
.
33--32--31--30
/ \
34 16--15--14 29
/ / \ \
35 17 5---4 13 28
/ / / \ \ \
36 18 6 0---3--12--27--48-->
/ / / / / / / /
37 19 7 1---2 11 26 47
\ \ \ / / /
38 20 8---9--10 25 46
\ \ / /
39 21--22--23--24 45
\ /
40--41--42--43--44
(End)
Number of edges of the complete bipartite graph of order 4n, K_n,3n. - Roberto E. Martinez II, Jan 07 2002
Also the number of partitions of 6n + 3 into at most 3 parts. - R. K. Guy, Oct 23 2003
Also the number of partitions of 6n into exactly 3 parts. - Colin Barker, Mar 23 2015
Numbers n such that the imaginary quadratic field Q[sqrt(-n)] has six units. - Marc LeBrun, Apr 12 2006
The denominators of Hoehn's sequence (recalled by G. L. Honaker, Jr.) and the numerators of that sequence reversed. The sequence is 1/3, (1+3)/(5+7), (1+3+5)/(7+9+11), (1+3+5+7)/(9+11+13+15), ...; reduced to 1/3, 4/12, 9/27, 16/48, ... . For the reversal, the reduction is 3/1, 12/4, 27/9, 48/16, ... . - Enoch Haga, Oct 05 2007
Right edge of tables in A200737 and A200741: A200737(n, A000292(n)) = A200741(n, A100440(n)) = a(n). - Reinhard Zumkeller, Nov 21 2011
The Wiener index of the crown graph G(n) (n>=3). The crown graph G(n) is the graph with vertex set {x(1), x(2), ..., x(n), y(1), y(2), ..., y(n)} and edge set {(x(i), y(j)): 1<=i, j<=n, i/=j} (= the complete bipartite graph K(n,n) with horizontal edges removed). Example: a(3)=27 because G(3) is the cycle C(6) and 6*1 + 6*2 + 3*3 = 27. The Hosoya-Wiener polynomial of G(n) is n(n-1)(t+t^2)+nt^3. - Emeric Deutsch, Aug 29 2013
From Michel Lagneau, May 04 2015: (Start)
Integer area A of equilateral triangles whose side lengths are in the commutative ring Z[3^(1/4)] = {a + b*3^(1/4) + c*3^(1/2) + d*3^(3/4), a,b,c and d in Z}.
The area of an equilateral triangle of side length k is given by A = k^2*sqrt(3)/4. In the ring Z[3^(1/4)], if k = q*3^(1/4), then A = 3*q^2/4 is an integer if q is even. Example: 27 is in the sequence because the area of the triangle (6*3^(1/4), 6*3^(1/4), 6*3^(1/4)) is 27. (End)
a(n) is 2*sqrt(3) times the area of a 30-60-90 triangle with short side n. Also, 3 times the area of an n X n square. - Wesley Ivan Hurt, Apr 06 2016
Consider the hexagonal tiling of the plane. Extract any four hexagons adjacent by edge. This can be done in three ways. Fold the four hexagons so that all opposite faces occupy parallel planes. For all parallel projections of the resulting object, at least two correspond to area a(n) for side length of n of the original hexagons. - Torlach Rush, Aug 17 2022
The sequence terms are the exponents in the expansion of Product_{n >= 1} (1 - q^(3*n))/(1 + q^(3*n)) = ( Sum_{n in Z} q^(n*(3*n+1)/2) ) / ( Product_{n >= 1} 1 + q^n ) = 1 - 2*q^3 + 2*q^12 - 2*q^27 + 2*q^48 - 2*q^75 + - .... - Peter Bala, Dec 30 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  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=1      n=2            n=3                 n=4
(End)
		

Crossrefs

Programs

  • Haskell
    a033428 = (* 3) . (^ 2)
    a033428_list = 0 : 3 : 12 : zipWith (+) a033428_list
       (map (* 3) $ tail $ zipWith (-) (tail a033428_list) a033428_list)
    -- Reinhard Zumkeller, Jul 11 2013
    
  • Magma
    [3*n^2: n in [0..50]]; // Vincenzo Librandi, May 18 2015
    
  • Maple
    seq(3*n^2, n=0..46); # Nathaniel Johnston, Jun 26 2011
  • Mathematica
    3 Range[0, 50]^2
    LinearRecurrence[{3, -3, 1}, {0, 3, 12}, 50] (* Harvey P. Dale, Feb 16 2013 *)
  • Maxima
    makelist(3*n^2,n,0,30); /* Martin Ettl, Nov 12 2012 */
    
  • PARI
    a(n)=3*n^2
    
  • Python
    def a(n): return 3 * (n**2) # Torlach Rush, Aug 25 2022

Formula

a(n) = 3*a(n-1)-3*a(n-2)+a(n-3) for n>2.
G.f.: 3*x*(1+x)/(1-x)^3. - R. J. Mathar, Sep 09 2008
Main diagonal of triangle in A132111: a(n) = A132111(n,n). - Reinhard Zumkeller, Aug 10 2007
A214295(a(n)) = -1. - Reinhard Zumkeller, Jul 12 2012
a(n) = A215631(n,n) for n > 0. - Reinhard Zumkeller, Nov 11 2012
a(n) = A174709(6n+2). - Philippe Deléham, Mar 26 2013
a(n) = a(n-1) + 6*n - 3, with a(0)=0. - Jean-Bernard François, Oct 04 2013
E.g.f.: 3*x*(1 + x)*exp(x). - Ilya Gutkovskiy, Apr 13 2016
a(n) = t(3*n) - 3*t(n), where t(i) = i*(i+k)/2 for any k. Special case (k=1): A000217(3*n) - 3*A000217(n). - Bruno Berselli, Aug 31 2017
a(n) = A000326(n) + A005449(n). - Bruce J. Nicholson, Jan 10 2020
From Amiram Eldar, Jul 03 2020: (Start)
Sum_{n>=1} 1/a(n) = Pi^2/18 (A086463).
Sum_{n>=1} (-1)^(n+1)/a(n) = Pi^2/36. (End)
From Amiram Eldar, Feb 03 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = sqrt(3)*sinh(Pi/sqrt(3))/Pi.
Product_{n>=1} (1 - 1/a(n)) = sqrt(3)*sin(Pi/sqrt(3))/Pi. (End)
a(n) = A003215(n) - A016777(n). - Leo Tavares, Apr 29 2023

Extensions

Better description from N. J. A. Sloane, May 15 1998

A214295 a(n) = 1 if n is a square, -1 if n is three times a square, 0 otherwise.

Original entry on oeis.org

1, 0, -1, 1, 0, 0, 0, 0, 1, 0, 0, -1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0
Offset: 1

Views

Author

Michael Somos, Jul 10 2012

Keywords

Comments

Ramanujan theta functions: f(q) (see A121373), phi(q) (A000122), psi(q) (A010054), chi(q) (A000700).
a(A092206(n)) = 0; a(A000290(n)) = 1; a(A033428(n)) = -1.

Examples

			G.f. = q - q^3 + q^4 + q^9 - q^12 + q^16 + q^25 - q^27 + q^36 - q^48 + q^49 + ...
		

Crossrefs

Programs

  • Haskell
    a214295 n = a010052 n - a010052 (3*n)  -- Reinhard Zumkeller, Jul 12 2012
    
  • Magma
    Basis( ModularForms( Gamma1(12), 1/2), 50) [2] ; /* Michael Somos, Jun 10 2014 */
  • Mathematica
    a[ n_] := SeriesCoefficient[ (EllipticTheta[ 3, 0, q] - EllipticTheta[ 3, 0, q^3]) / 2, {q, 0, n}];
    a[ n_] := Boole[ IntegerQ[ Sqrt[ n]]] - Boole[ IntegerQ[ Sqrt[ 3 n]]]; (* Michael Somos, Jun 10 2014 *)
    Table[Which[IntegerQ[Sqrt[n]],1,IntegerQ[Sqrt[n/3]],-1,True,0],{n,120}] (* Harvey P. Dale, Apr 08 2013 *)
  • PARI
    {a(n) = issquare(n) - issquare(3*n)};
    
  • PARI
    {a(n) = if( n<1, 0, direuler( p=2, n, if( p==3, 1 - X, 1) / (1 - X^2 ))[n])};
    

Formula

Expansion of q * psi(q^3) * f(-q^2, -q^10) / f(-q^5, -q^7) in powers of q where psi(), f() are Ramanujan theta functions.
Multiplicative with a(3^e) = (-1)^e, a(p^e) = 1 if e even, 0 otherwise.
G.f.: (theta_3(q) - theta_3(q^3)) / 2 = Sum_{k>0} x^(k^2) - x^(3*k^2).
Dirichlet g.f.: zeta(2*s) * (1 - 3^(-s)). [corrected by Amiram Eldar, Oct 24 2023]
a(3*n) = - a(n). - Reinhard Zumkeller, Jul 12 2012
Expansion of (phi(q) - phi(q^3)) / 2 = q * chi(q) * f(-q, -q^11) in powers fof q where phi(), chi(), f() are Ramanujan theta functions. - Michael Somos, Jan 10 2015
Euler transform of period 12 sequence [ 0, -1, 1, 0, 1, -1, 1, 0, 1, -1, 0, -1, ...]. - Michael Somos, Jan 10 2015
Convolution product of A000700 and A247133. - Michael Somos, Jan 10 2015
Sum_{k=1..n} a(k) ~ c*sqrt(n), where c = 1 - 1/sqrt(3) = 0.42264973... . - Amiram Eldar, Oct 24 2023

A092205 Number of units in the imaginary quadratic field Q(sqrt(-n)).

Original entry on oeis.org

4, 2, 6, 4, 2, 2, 2, 2, 4, 2, 2, 6, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 6, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 6, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 6, 2, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2
Offset: 1

Views

Author

Eric W. Weisstein, Feb 24 2004

Keywords

Comments

Sequence of n such that a(n)=2 gives A092206; a(n)=4 gives A000290; a(n)=6 gives A033428. - Marc LeBrun, Apr 12 2006

Examples

			For n=1, the units are +/-1, +/-i, so a(1) = 4.
For n=3, the units are +/-1, +/-w, +/-w^2, where w is a cube root of unity, so a(3) = 6. [Corrected by _Jonathan Sondow_, Jan 29 2014]
		

Crossrefs

Programs

  • Maple
    A092205 := proc(n) if(type(sqrt(n),integer))then return 4: elif(n mod 3 = 0 and type(sqrt(n/3),integer))then return 6: else return 2: fi: end: seq(A092205(n),n=1..105); # Nathaniel Johnston, Jun 26 2011
  • Mathematica
    a[n_] := Which[ IntegerQ[ Sqrt[n] ], 4, Mod[n, 3] == 0 && IntegerQ[ Sqrt[n/3] ], 6, True, 2]; Table[a[n], {n, 1, 105}] (* Jean-François Alcover, Oct 30 2012, after Nathaniel Johnston *)
  • PARI
    a(n)=if(issquare(n),return(4));if(n%3==0&&issquare(n/3),6,2) \\ Charles R Greathouse IV, Oct 30 2012

Formula

a(A005117(n)) = A236213(n). - Jonathan Sondow, Jan 29 2014

A028737 Nonsquares mod 24.

Original entry on oeis.org

2, 3, 5, 6, 7, 8, 10, 11, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23
Offset: 1

Views

Author

Keywords

Comments

The first 18 terms of A092206. - R. J. Mathar, Mar 21 2023

Crossrefs

Cf. A010386.

Programs

  • Mathematica
    Complement[Range[0, 23], PowerMod[Range[0, 23], 2, 24]] (* Alonso del Arte, Nov 24 2019 *)
  • Scala
    val squaresMod24 = (0 to 23).map(n => n * n).map(_ % 24)
    (0 to 23).diff(squaresMod24) // Alonso del Arte, Nov 24 2019
Showing 1-4 of 4 results.