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

A047971 Triangle of coefficients of Gaussian polynomials [ n+3,3 ].

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 1, 1, 1, 1, 2, 3, 3, 3, 3, 2, 1, 1, 1, 1, 2, 3, 4, 4, 5, 4, 4, 3, 2, 1, 1, 1, 1, 2, 3, 4, 5, 6, 6, 6, 6, 5, 4, 3, 2, 1, 1, 1, 1, 2, 3, 4, 5, 7, 7, 8, 8, 8, 7, 7, 5, 4, 3, 2, 1, 1
Offset: 0

Views

Author

Keywords

Comments

a(n) as illustrated is related to the following sequences: The row sum values are A001400. The column sums are A000292. The row lengths are the stuttering sequence A037915 (stutter values in A016777). The column lengths are the sequence A016777. The max values in each column are A001971. - Alford Arnold, Aug 16 2004
The entry a(p,w), p >= 0, w = 0,1,...,3*p, of this irregular triangle is the number of nonnegative solutions of m_0 + m_1 + m_2 + m_3 = p and 1*m_1 + 2*m_2 + 3*m_3 = w. See the Hawkins reference given in A008967, p. 264, (4,7),(4.8), concerning Cayley's counting problem. N(p,3,w) there equals a(p,w). The o.g.f. has been given in the formula section by Peter Bala. See also the Cayley reference given in A008967, p. 110, 35. with m = 3, Theta = p and q = w. - Wolfdieter Lang, Dec 02 2012
The entry a(p,w) p >= 0, w = 0,1,...,3*p, of this array gives the number of partitions of w into at most p parts, each at most 3. This follows from the preceding comment with the two Diophantine equations. From Andrews, p. 33 and p. 35, a(p,w) (called there p(N,M,n) with N=p, M=3, n=w) gives also the number of partitions of w into at most 3 parts, each at most p. This is in accordance with the symmetry of the q-binomials [p+3,p] = [p+3,3]. - Wolfdieter Lang, Dec 04 2012

Examples

			The table a(p,w) = [q^w][p+3,3]_q starts:
p\w 0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 ...
0:  1
1:  1  1  1  1
2:  1  1  2  2  2  1  1
3:  1  1  2  3  3  3  3  2  1  1
4:  1  1  2  3  4  4  5  4  4  3  2  1  1
5:  1  1  2  3  4  5  6  6  6  6  5  4  3  2  1  1
6:  1  1  2  3  4  5  7  7  8  8  8  7  7  5  4  3  2  1  1
... Reformatted and extended by _Wolfdieter Lang_, Dec 04 2012
Partition example: Row p=2 is 1 1 2 2 2 1 1 because there are ten solution for (m_0, m_1, m_2, m_3) of the first equation given in a comment above, namely (2,0,0,0), (0,2,0,0), (0,0,2,0), (0,0,0,2), (1,1,0,0), (1,0,1,0), (1,0,0,1), (0,1,1,0), (0,1,0,1) and (0,0,1,1) which have the w = 1*m_1 + 2*m_2 + 3*m_3 values 0, 2, 4, 6, 1, 2, 3, 3, 4 and 5, respectively. Therefore there are 1, 1, 2, 2, 2, 1, 1 solutions for w = 0, 1, 2, 3, 4, 5, 6, respectively. - _Wolfdieter Lang_, Dec 03 2012
a(4,5) = 4 because there are 4 partitions of 5 with 1, 2, 3  or 4 parts, each being <= 3, namely all partitions of 5 excluding 5, 14 and 11111. There are also 4 partitions of 5 with 1, 2, or 3 parts, each being <= 4, namely all partitions of 5 excluding 5, 1112 and 11111. - _Wolfdieter Lang_, Dec 04 2012
The table may also be arranged as follows (see the Alford Arnold comment above):
1
..1
..1..1
..1..1..1
..1..2..1..1
.....2..2..1..1
.....2..3..2..1..1
.....1..3..3..2..1..1
		

References

  • G. E. Andrews, The Theory of Partitions, Addison-Wesley, 1976, p. 242.

Crossrefs

Cf. A008967.
Cf. A001400.

Programs

  • Mathematica
    nmax = 6;
    se = Series[ 1/Product[1 - q^k*x, {k, 0, 3}], {x, 0, nmax}];
    row[n_] := CoefficientList[ SeriesCoefficient[se, n], q];
    Flatten[ Table[ row[n], {n, 0, nmax}]] (* Jean-François Alcover, Dec 19 2011 *)

Formula

O.g.f.: 1/((1-x)(1-qx)(1-q^2x)(1-q^3x)) = 1 + x(1 + q + q^2 + q^3) + x^2(1 + q + 2q^2 + 2q^3 + 2q^4 + q^5 + q^6) + .... - Peter Bala, Sep 23 2007

A195097 Fractalization of (1+[3n/4]), where [ ] = floor.

Original entry on oeis.org

1, 1, 2, 1, 2, 3, 1, 2, 3, 4, 1, 2, 3, 5, 4, 1, 2, 3, 5, 6, 4, 1, 2, 3, 5, 6, 7, 4, 1, 2, 3, 5, 6, 7, 8, 4, 1, 2, 3, 5, 6, 7, 9, 8, 4, 1, 2, 3, 5, 6, 7, 9, 10, 8, 4, 1, 2, 3, 5, 6, 7, 9, 10, 11, 8, 4, 1, 2, 3, 5, 6, 7, 9, 10, 11, 12, 8, 4, 1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 12, 8, 4, 1, 2, 3
Offset: 1

Views

Author

Clark Kimberling, Sep 08 2011

Keywords

Comments

See A194959 for a discussion of fractalization and the interspersion fractally induced by a sequence. The sequence (1+[3n/4]) is a subsequence ofy A037915.

Crossrefs

Programs

  • Mathematica
    r = 3/4; p[n_] := 1 + Floor[n*r] (* A037915 *)
    Table[p[n], {n, 1, 90}]
    g[1] = {1}; g[n_] := Insert[g[n - 1], n, p[n]]
    f[1] = g[1]; f[n_] := Join[f[n - 1], g[n]]
    f[20]    (* A195097 *)
    row[n_] := Position[f[30], n];
    u = TableForm[Table[row[n], {n, 1, 5}]]
    v[n_, k_] := Part[row[n], k];
    w = Flatten[Table[v[k, n - k + 1], {n, 1, 13},
    {k, 1, n}]](* A195098 *)
    q[n_] := Position[w, n]; Flatten[Table[q[n],
    {n, 1, 80}]](* A195099 *)

A195098 Interspersion fractally induced by (1+[3n/4]), where [ ] = floor; a rectangular array, by antidiagonals.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 15, 14, 16, 17, 18, 21, 19, 20, 22, 23, 24, 28, 25, 26, 27, 29, 30, 31, 36, 32, 33, 34, 35, 37, 38, 39, 45, 40, 41, 42, 44, 43, 46, 47, 48, 55, 49, 50, 51, 54, 52, 53, 56, 57, 58, 66, 59, 60, 61, 65, 62, 63, 64, 67, 68, 69
Offset: 1

Views

Author

Clark Kimberling, Sep 08 2011

Keywords

Comments

See A194959 for a discussion of fractalization and the interspersion fractally induced by a sequence. Every pair of rows eventually intersperse. As a sequence, A194998 is a permutation of the positive integers, with inverse A195099.

Examples

			Northwest corner:
1...2...4...7...11..16..22
3...5...8...12..17..23..30
6...9...13..18..24..31..39
10..15..21..28..36..45..55
14..19..25..32..40..49..59
		

Crossrefs

Programs

  • Mathematica
    r = 3/4; p[n_] := 1 + Floor[n*r] (* A037915 *)
    Table[p[n], {n, 1, 90}]
    g[1] = {1}; g[n_] := Insert[g[n - 1], n, p[n]]
    f[1] = g[1]; f[n_] := Join[f[n - 1], g[n]]
    f[20]    (* A195097 *)
    row[n_] := Position[f[30], n];
    u = TableForm[Table[row[n], {n, 1, 5}]]
    v[n_, k_] := Part[row[n], k];
    w = Flatten[Table[v[k, n - k + 1], {n, 1, 13},
    {k, 1, n}]](* A195098 *)
    q[n_] := Position[w, n]; Flatten[Table[q[n],
    {n, 1, 80}]](* A195099 *)

A090221 Array used for numerators of g.f.s for column sequences of array A090214 ((4,4)-Stirling2).

Original entry on oeis.org

1, 96, 72, 14400, 16, 38400, 3456000, 1, 27000, 22104000, 1270080000, 7200, 34905600, 16111872000, 682795008000, 856, 21154176, 48248363520, 15279164006400, 516193026048000, 48, 6064128, 54644474880, 78083415244800
Offset: 4

Views

Author

Wolfdieter Lang, Dec 01 2003

Keywords

Comments

The row length sequence for this array is A037915(k-4)= floor(3*(k-4)/4)+1, k>=4: [1, 1, 2, 3, 4, 4, 5, 6, 7, 7, 8, 9, 10, 10, 11, ...].
The g.f. G(k,x) for the k-th column (with leading zeros) of array A090214 is given there. The recurrence is G(k,x) = x*sum(binomial(k-r,4-r)*fallfac(4,4-r)*G(k-r,x),r=1..4))/(1-fallfac(k,4)*x), k>=4, with inputs G(k,x)=0 for k=1,2,3 and G(4,x)=x/(1-4!*x); where fallfac(n,m) := A008279(n,m) (falling factorials with fallfac(n,0) := 1). Computed from the Blasiak et al. reference, eqs. (20) and (21) with r=4: recurrence for S_{4,4}(n,k).

Examples

			[1]; [96]; [72,14400]; [16,38400,3456000]; [1,27000,22104000,1270080000]; ...
G(5,x)/x^2 = 96/((1-4!*x)*(1-5*4*3*2*x)). kmax(5)=0, hence P(5,x)=a(5,0)=96; x^2 from x^ceiling(5/4).
		

Formula

a(k, n) from: sum(a(k, n)*x^n, n=0..kmax(k)) = G(k, x)* product(1-fallfac(p, 4)*x, p=4..k)/x^ceiling(k/4), k>=4, with G(k, x) defined from the recurrence given above and kmax(k) := A057353(k-4)= floor(3*(k-4)/4)= A037915(k-4)-1.

A353212 Hadwiger number of the n-path complement graph.

Original entry on oeis.org

1, 1, 2, 2, 3, 4, 5, 6, 7, 7, 8, 9, 10, 10, 11, 12, 13, 13, 14, 15, 16, 16, 17, 18, 19, 19, 20, 21, 22, 22, 23, 24, 25, 25, 26, 27, 28, 28, 29, 30, 31, 31, 32, 33, 34, 34, 35, 36, 37, 37, 38, 39, 40, 40, 41, 42, 43, 43, 44, 45, 46, 46, 47, 48, 49, 49, 50, 51, 52, 52
Offset: 1

Views

Author

Eric W. Weisstein, Apr 30 2022

Keywords

Comments

A contraction in the complement of any set of paths reduces the total number of edges in the complement by at most 4. This gives an upper bound for the Hadwiger number which is obtainable for all path lengths except 4 and 5. In particular, for n >= 6, the complement of a P_n reduces to the complement of a P_{n-4} union 3 universal nodes by contracting the second and second to last nodes of the path. With P_8 and P_9 the 2nd and 6th nodes should be contracted (instead of reducing to P_4 or P_5 respectively). - Andrew Howroyd, Jun 18 2025

Crossrefs

Cf. A037915.

Programs

  • PARI
    a(n) = (3*n + 1)\4 - (n==4||n==5) \\ Andrew Howroyd, Jun 18 2025

Formula

a(n) = floor((3*n + 1)/4) = A037915(n+1) for n >= 6. - Andrew Howroyd, Jun 18 2025

Extensions

a(16) onwards from Andrew Howroyd, Jun 18 2025

A066530 Expansion of (1+x+x^3)/((1-x)*(1-x^4)).

Original entry on oeis.org

1, 2, 2, 3, 4, 5, 5, 6, 7, 8, 8, 9, 10, 11, 11, 12, 13, 14, 14, 15, 16, 17, 17, 18, 19, 20, 20, 21, 22, 23, 23, 24, 25, 26, 26, 27, 28, 29, 29, 30, 31, 32, 32, 33, 34, 35, 35, 36, 37, 38, 38, 39, 40, 41, 41, 42, 43, 44, 44, 45, 46, 47, 47, 48, 49, 50, 50, 51, 52, 53, 53, 54, 55
Offset: 0

Views

Author

Robert G. Wilson v, Jan 06 2002

Keywords

Crossrefs

Cf. A037915.

Programs

  • Magma
    [Floor((3*n + 5)/4): n in [0..100]]; // Wesley Ivan Hurt, Jan 02 2017
  • Maple
    A066530:=n->floor((3*n+5)/4): seq(A066530(n), n=0..100); # Wesley Ivan Hurt, Jan 02 2017
  • Mathematica
    CoefficientList[ Series[ (1 + x + x^3)/((1 - x)*(1 - x^4)), {x, 0, 75} ], x] (* or *) Table[Floor[(3 n + 5)/4], {n, 0, 75}]
    LinearRecurrence[{1, 0, 0, 1, -1}, {1, 2, 2, 3, 4}, 80] (* Harvey P. Dale, Apr 01 2013 *)

Formula

G.f.: (1 + x + x^3)/((1 - x)*(1 - x^4)).
a(n) = floor((3*n + 5)/4).
a(n) = a(n-1) + a(n-4) - a(n-5) for n>4, with a(0)=1, a(1)=2, a(2)=2, a(3)=3, a(4)=4. - Harvey P. Dale, Apr 01 2013
a(n) = (6*n+7+2*cos(n*Pi/2)+cos((n+1)*Pi)+2*sin(n*Pi/2))/8. - Wesley Ivan Hurt, Oct 01 2017
Sum_{n>=0} (-1)^n/a(n) = log(3)/2 + Pi/(6*sqrt(3)). - Amiram Eldar, Jan 31 2023

A137178 a(n) = sum_(1..n) [S2(n)mod 2 - floor(5*S2(n)/7)mod 2], where S2(n) = binary weight of n.

Original entry on oeis.org

0, 1, 2, 1, 2, 1, 0, 1, 2, 1, 0, 1, 0, 1, 2, 2, 3, 2, 1, 2, 1, 2, 3, 3, 2, 3, 4, 4, 5, 5, 5, 5, 6, 5, 4, 5, 4, 5, 6, 6, 5, 6, 7, 7, 8, 8, 8, 8, 7, 8, 9, 9, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 12, 11, 10, 11, 10, 11, 12, 12, 11, 12, 13, 13, 14, 14, 14, 14, 13, 14, 15, 15, 16, 16, 16
Offset: 0

Views

Author

Ctibor O. Zizka, Apr 04 2008, Apr 15 2008

Keywords

Comments

The graph of this sequence is a special case of de Rham's fractal curve. In general, the graph of any sequence of the form a(n)=sum_(1..n) [Sk(n)mod m - floor(p*Sk(n)/q)mod m], where Sk(n) is the digit sum of n, n in k-ary notation, p,q,m integers, gives a de Rham fractal curve. The self-symmetries of all of de Rham curves are given by the monoid that describes the symmetries of the infinite binary tree or Cantor set. This so-called period-doubling monoid is a subset of the modular group.

Crossrefs

Programs

  • Mathematica
    Accumulate@ Array[Mod[#2, 2] - Mod[Floor[5 #2/7], 2] & @@ {#, DigitCount[#, 2, 1]} &, 85, 0] (* Michael De Vlieger, Jan 23 2019 *)

Extensions

Converted references to links - R. J. Mathar, Oct 30 2009

A278703 Divide a full period sine wave into n equally spaced points along the x-axis, labeled 1 through n, from left to right. Rank the points according to their sine value and read by rows.

Original entry on oeis.org

1, 2, 1, 3, 2, 1, 4, 3, 2, 1, 4, 5, 3, 1, 2, 5, 6, 4, 3, 1, 2, 6, 5, 7, 4, 1, 3, 2, 7, 6, 8, 5, 4, 1, 3, 2, 7, 8, 6, 9, 5, 1, 4, 2, 3, 8, 9, 7, 10, 6, 5, 1, 4, 2, 3, 9, 8, 10, 7, 11, 6, 1, 5, 2, 4, 3, 10, 9, 11, 8, 12, 7, 6, 1, 5, 2, 4, 3, 10, 11, 9, 12, 8, 13, 7, 1, 6, 2, 5, 3, 4
Offset: 1

Views

Author

Robert G. Wilson v, Nov 26 2016

Keywords

Comments

Inspired by A276669.
Two other sequences in the same vein could be constructed, one with points from 0 to n-1 and the other with points from 0 to n. The latter would only insert a zero before n.
Column 1 appears to be A037915 and the last column appears to be A002265.

Examples

			Row  1:  1;
Row  2:  2,  1;
Row  3:  3,  2,  1;
Row  4:  4,  3,  2,  1;
Row  5:  4,  5,  3,  1,  2;
Row  6:  5,  6,  4,  3,  1,  2;
Row  7:  6,  5,  7,  4,  1,  3,  2;
Row  8:  7,  6,  8,  5,  4,  1,  3,  2;
Row  9:  7,  8,  6,  9,  5,  1,  4,  2,  3;
Row 10:  8,  9,  7, 10,  6,  5,  1,  4,  2,  3;
etc.
Row 3: The first point is (Pi/2, 1), the second point is (Pi,0) and the third point is (3*Pi/2, -1). Sorting by the Y value and reading the points by their index, we have 3, 2, 1.
Row 4: The first point is at (2*Pi/5, sqrt(5/8 + sqrt(5)/8)), point two is at (4*Pi/5, sqrt(5/8 - sqrt(5)/8)), point three is at (6*Pi/5, -sqrt(5/8 - sqrt(5)/8)) and point four is at (8*Pi/5, -sqrt(5/8 + sqrt(5)/8)). Sorting the point labels by their Y values in increasing order, we have 4, 3, 2, 1.
		

Crossrefs

Cf. A276669.

Programs

  • Mathematica
    f[n_] := Transpose[ Sort[ Table[{N[ Sin[ 2i*Pi/(n +1)], 9], i}, {i, n}]]][[2]]; Array[f, 13] // Flatten
Showing 1-8 of 8 results.