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.

A005232 Expansion of (1-x+x^2)/((1-x)^2*(1-x^2)*(1-x^4)).

Original entry on oeis.org

1, 1, 3, 4, 8, 10, 16, 20, 29, 35, 47, 56, 72, 84, 104, 120, 145, 165, 195, 220, 256, 286, 328, 364, 413, 455, 511, 560, 624, 680, 752, 816, 897, 969, 1059, 1140, 1240, 1330, 1440, 1540, 1661, 1771, 1903, 2024, 2168, 2300, 2456, 2600, 2769, 2925, 3107, 3276
Offset: 0

Views

Author

Keywords

Comments

Also number of n X 2 binary matrices under row and column permutations and column complementations (if offset is 0).
Also Molien series for certain 4-D representation of dihedral group of order 8.
With offset 4, number of bracelets (turnover necklaces) of n-bead of 2 colors with 4 red beads. - Washington Bomfim, Aug 27 2008
From Vladimir Shevelev, Apr 23 2011: (Start)
Also number of non-equivalent necklaces of 4 beads each of them painted by one of n colors.
The sequence solves the so-called Reis problem about convex k-gons in case k=4 (see our comment to A032279). (End)
Number of 2 X 2 matrices with nonnegative integer values totaling n under row and column permutations. - Gabriel Burns, Nov 08 2016
From Petros Hadjicostas, Jan 12 2019: (Start)
By "necklace", Vladimir Shevelev (above) means "turnover necklace", i.e., a bracelet. Zagaglia Salvi (1999) also uses this terminology: she calls a bracelet "necklace" and a necklace "cycle".
According to Cyvin et al. (1997), the sequence (a(n): n >= 0) consists of "the total numbers of isomers of polycyclic conjugated hydrocarbons with q + 1 rings and q internal carbons in one ring (class Q_q)", where q = 4 and n is the hydrogen content (i.e., we count certain isomers of C_{n+2*q} H_n with q = 4 and n >= 0). (End)

Examples

			G.f. = 1 + x + 3*x^2 + 4*x^3 + 8*x^4 + 10*x^5 + 16*x^6 + 20*x^7 + 29*x^8 + ...
There are 8 4 X 2 matrices up to row and column permutations and column complementations:
  [1 1] [1 0] [1 0] [0 1] [0 1] [0 1] [0 1] [0 0]
  [1 1] [1 1] [1 0] [1 0] [1 0] [1 0] [0 1] [0 1]
  [1 1] [1 1] [1 1] [1 1] [1 0] [1 0] [1 0] [1 0]
  [1 1] [1 1] [1 1] [1 1] [1 1] [1 0] [1 0] [1 1].
There are 8 2 X 2 matrices of nonnegative integers totaling 4 up to row and column permutations:
  [4 0] [3 1] [2 2] [2 1] [2 1] [3 0] [2 0] [1 1]
  [0 0] [0 0] [0 0] [0 1] [1 0] [1 0] [2 0] [1 1].
		

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • N. Zagaglia Salvi, Ordered partitions and colourings of cycles and necklaces, Bull. Inst. Combin. Appl., 27 (1999), 37-40.

Crossrefs

Row n=2 of A343875.
Column k=4 of A052307.

Programs

  • Maple
    A005232:=-(-1-z-2*z**3+2*z**2+z**7-2*z**6+2*z**4)/(z**2+1)/(1+z)**2/(z-1)**4; # conjectured by Simon Plouffe in his 1992 dissertation; gives sequence apart from an initial 1
  • Mathematica
    k = 4; Table[(Apply[Plus, Map[EulerPhi[ # ]Binomial[n/#, k/# ] &, Divisors[GCD[n, k]]]]/n + Binomial[If[OddQ[n], n - 1, n - If[OddQ[k], 2, 0]]/2, If[OddQ[k], k - 1, k]/2])/2, {n, k, 50}] (* Robert A. Russell, Sep 27 2004 *)
    CoefficientList[ Series[(1 - x + x^2)/((1 - x)^2(1 - x^2)(1 - x^4)), {x, 0, 51}], x] (* Robert G. Wilson v, Mar 29 2006 *)
    LinearRecurrence[{2,0,-2,2,-2,0,2,-1},{1,1,3,4,8,10,16,20},60] (* Harvey P. Dale, Oct 24 2012 *)
    k=4 (* Number of red beads in bracelet problem *); CoefficientList[Series[(1/k Plus@@(EulerPhi[#] (1-x^#)^(-(k/#))&/@Divisors[k])+(1+x)/(1-x^2)^Floor[(k+2)/2])/2,{x,0,50}],x] (* Herbert Kociemba, Nov 04 2016 *)
  • PARI
    {a(n) = (n^3 + 9*n^2 + (32-9*(n%2))*n + [48, 15, 36, 15][n%4+1]) / 48}; \\ Michael Somos, Feb 01 2007
    
  • PARI
    {a(n) = my(s=1); if( n<-5, n = -6-n; s=-1); if( n<0, 0, s * polcoeff( (1 - x + x^2) / ((1 - x)^2 * (1 - x^2) * (1 - x^4)) + x * O(x^n), n))}; \\ Michael Somos, Feb 01 2007
    
  • PARI
    a(n) = round((n^3 +9*n^2 +(32-9*(n%2))*n)/48 +0.6) \\ Washington Bomfim, Jul 17 2008
    
  • PARI
    a(n) = ceil((n+1)*(2*n^2+16*n+39+9*(-1)^n)/96) \\ Tani Akinari, Aug 23 2013
    
  • Python
    a=lambda n: sum((k//2+1)*((n-k)//2+1) for k in range((n-1)//2+1))+(n+1)%2*(((n//4+1)*(n//4+2))//2)  # Gabriel Burns, Nov 08 2016

Formula

G.f.: (1+x^3)/((1-x)*(1-x^2)^2*(1-x^4)).
G.f.: (1/8)*(1/(1-x)^4+3/(1-x^2)^2+2/(1-x)^2/(1-x^2)+2/(1-x^4)). - Vladeta Jovovic, Aug 05 2000
Euler transform of length 6 sequence [ 1, 2, 1, 1, 0, -1 ]. - Michael Somos, Feb 01 2007
a(2n+1) = A006918(2n+2)/2;
a(2n) = (A006918(2n+1) + A008619(n))/2.
a(n) = -a(-6 - n) for all n in Z. - Michael Somos, Feb 05 2011
From Vladimir Shevelev, Apr 22 2011: (Start)
if n == 0 (mod 4), then a(n) = n*(n^2-3*n+8)/48;
if n == 1, 3 (mod 4), then a(n) = (n^2-1)*(n-3)/48;
if n == 2 (mod 4), then a(n) = (n-2)*(n^2-n+6)/48. (End)
a(n) = 2*a(n-1) - 2*a(n-3) + 2*a(n-4) - 2*a(n-5) + 2*a(n-7) - a(n-8) with a(0) = 1, a(1) = 1, a(2) = 3, a(3) = 4, a(4) = 8, a(5) = 10, a(6) = 16, a(7) = 20. - Harvey P. Dale, Oct 24 2012
a(n) = ((n+3)*(2*n^2+12*n+19+9*(-1)^n) + 6*(-1)^((2*n-1+(-1)^n)/4)*(1+(-1)^n))/96. - Luce ETIENNE, Mar 16 2015
a(n) = |A128498(n)| + |A128498(n-3)|. - R. J. Mathar, Jun 11 2019

Extensions

Sequence extended by Christian G. Bower

A343097 Array read by antidiagonals: T(n,k) is the number of k-colorings of an n X n grid, up to rotations and reflections.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, 3, 6, 1, 0, 1, 4, 21, 102, 1, 0, 1, 5, 55, 2862, 8548, 1, 0, 1, 6, 120, 34960, 5398083, 4211744, 1, 0, 1, 7, 231, 252375, 537157696, 105918450471, 8590557312, 1, 0, 1, 8, 406, 1284066, 19076074375, 140738033618944, 18761832172500795, 70368882591744, 1, 0
Offset: 0

Views

Author

Andrew Howroyd, Apr 14 2021

Keywords

Examples

			Array begins:
====================================================================
n\k | 0 1       2            3               4                 5
----+---------------------------------------------------------------
  0 | 1 1       1            1               1                 1 ...
  1 | 0 1       2            3               4                 5 ...
  2 | 0 1       6           21              55               120 ...
  3 | 0 1     102         2862           34960            252375 ...
  4 | 0 1    8548      5398083       537157696       19076074375 ...
  5 | 0 1 4211744 105918450471 140738033618944 37252918396015625 ...
  ...
		

Crossrefs

Programs

  • PARI
    T(n,k) = {(k^(n^2) + 2*k^((n^2 + 3*(n%2))/4) + k^((n^2 + (n%2))/2) + 2*k^(n*(n+1)/2) + 2*k^(n*(n+n%2)/2) )/8}

Formula

T(n,k) = (k^(n^2) + 2*k^((n^2 + 3*(n mod 2))/4) + k^((n^2 + (n mod 2))/2) + 2*k^(n*(n+1)/2) + 2*k^(n*(n + n mod 2)/2) )/8.

A343874 Array read by antidiagonals: T(n,k) is the number of n X n nonnegative integer matrices with sum of elements equal to k, up to rotational symmetry.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 1, 1, 1, 0, 1, 3, 3, 1, 0, 1, 5, 13, 4, 1, 0, 1, 10, 43, 36, 7, 1, 0, 1, 14, 129, 204, 85, 9, 1, 0, 1, 22, 327, 980, 735, 171, 13, 1, 0, 1, 30, 761, 3876, 5145, 2109, 313, 16, 1, 0, 1, 43, 1619, 13596, 29715, 20610, 5213, 528, 21, 1
Offset: 0

Views

Author

Andrew Howroyd, May 06 2021

Keywords

Examples

			Array begins:
=====================================================
n\k | 0  1   2    3     4      5       6        7
----+------------------------------------------------
  0 | 1  0   0    0     0      0       0        0 ...
  1 | 1  1   1    1     1      1       1        1 ...
  2 | 1  1   3    5    10     14      22       30 ...
  3 | 1  3  13   43   129    327     761     1619 ...
  4 | 1  4  36  204   980   3876   13596    42636 ...
  5 | 1  7  85  735  5145  29715  148561   657511 ...
  6 | 1  9 171 2109 20610 164502 1124382  6744582 ...
  7 | 1 13 313 5213 67769 717509 6457529 50732669 ...
  ...
		

Crossrefs

Rows n=0..4 are A000007, A000012, A008610, A054771, A054773.
Columns k=0..1 are A000012, A004652.
Cf. A054772 (binary case), A318795, A343095, A343875.

Programs

  • PARI
    U(n,s)={(s(1)^(n^2) + s(1)^(n%2)*(2*s(4)^(n^2\4) + s(2)^(n^2\2)))/4}
    T(n,k)={polcoef(U(n,i->1/(1-x^i) + O(x*x^k)), k)}

A054343 Number of nonnegative integer 3 X 3 matrices with sum of elements equal to n, under action of dihedral group of the square D_4.

Original entry on oeis.org

1, 3, 11, 31, 84, 198, 440, 904, 1766, 3266, 5802, 9906, 16384, 26284, 41104, 62752, 93831, 137589, 198309, 281249, 393148, 542154, 738480, 994320, 1324668, 1747220, 2283396, 2958228, 3801600, 4848120, 6138624, 7720032, 9647133, 11982423, 14798223, 18176499
Offset: 0

Views

Author

Vladeta Jovovic, May 05 2000

Keywords

Examples

			There are 11 nonisomorphic nonnegative integer 3 X 3 matrices with sum of elements equal to 2, under action of D_4:
[0 0 0] [0 0 0] [0 0 0] [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 0 0] [0 0 0] [0 0 1] [0 0 1] [0 1 0] [0 1 0] [1 0 1] [0 0 0] [0 0 0] [0 0 0] [0 2 0]
[0 1 1] [1 0 1] [0 1 0] [1 0 0] [0 0 1] [0 1 0] [0 0 0] [1 0 0] [0 0 2] [0 2 0] [0 0 0].
		

Crossrefs

Row n=3 of A343875.

Programs

  • PARI
    Vec((2*x^6+2*x^5+x^4+4*x^2-2*x+1)/((1-x^4)^2*(1-x^2)^2*(1-x)^5) + O(x^40)) \\ Colin Barker, Apr 26 2019

Formula

G.f.: (2*x^6+2*x^5+x^4+4*x^2-2*x+1)/((1-x^4)^2*(1-x^2)^2*(1-x)^5).
a(n) = 5*a(n-1) - 8*a(n-2) + 16*a(n-4) - 24*a(n-5) + 16*a(n-6) + 8*a(n-7) - 34*a(n-8) + 34*a(n-9) - 8*a(n-10) - 16*a(n-11) + 24*a(n-12) - 16*a(n-13) + 8*a(n-15) - 5*a(n-16) + a(n-17) for n>16. - Colin Barker, Apr 26 2019
Showing 1-4 of 4 results.