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.

User: Wouter Meeussen

Wouter Meeussen's wiki page.

Wouter Meeussen has authored 453 sequences. Here are the ten most recent ones:

A382247 Number of fixed points of solid partitions under twice the 'time-lapse' operation.

Original entry on oeis.org

1, 0, 2, 2, 3, 4, 7, 12, 16, 22, 32, 50, 68, 96, 134, 195, 261, 364, 497, 701, 941, 1288, 1738
Offset: 1

Author

Wouter Meeussen, Mar 19 2025

Keywords

Comments

Permutes the 4 axes of the 4D-Ferrers plot of the solid partitions as 2143.

Examples

			z[{{2},{2}}] -> z[{{1,1}},{{1,1}}] -> z[{{2},{2}}] under the 'lapse' operation.
		

Programs

  • Mathematica
    Tr/@Table[Count[solidformBTK[par], arg_z/; Nest[lapse,arg,2]===arg], {n, 20}, {par, IntegerPartitions[n]}]

A381332 a(n) is the number of different hooklength lists of the plane partitions of n.

Original entry on oeis.org

1, 1, 2, 4, 6, 11, 19, 31, 52, 86, 146, 231, 392, 615, 1006, 1594, 2612, 4062, 6518, 10116, 15958, 24557, 38565, 58548
Offset: 1

Author

Wouter Meeussen, Feb 20 2025

Keywords

Comments

The hooklength list of a plane partition is the sorted list of 3D hooklengths of its 3D Ferrers plot, analogous to the classic 2D case.

Examples

			The plane partition {{2,1},{2}} has hooklengths {{{4,2},{1}},{{2,1}}} and so hooklength list is {4,2,2,1,1}. So a(2) = 1.
The 24 plane partitions of n=5 generate only these 6 hooklength lists: {4,2,2,1,1}, {4,3,2,1,1}, {5,2,1,1,1}, {5,2,2,1,1}, {5,3,2,1,1}, {5,4,3,2,1}. So a(5) = 6.
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Union[planehooks/@planepartitions[n]]],{n,20}]

A381265 a(n) is the number of solid (3D) partitions of n with 2 layers and second layer a plane partition of 3.

Original entry on oeis.org

6, 21, 57, 138, 294, 606, 1170, 2208, 4008, 7176, 12492, 21510, 36348, 60801, 100281, 164019, 265263, 425853
Offset: 3

Author

Wouter Meeussen, Feb 18 2025

Keywords

Comments

Conjecture: equal to 3*(2*A000219 -A000990 -2*A000041 +1) tested up to n=20.

Examples

			a(3)=6 since the 6 solid partitions of {3,3} are:
  z[{{3}},{{3}}],
  z[{{2,1}},{{2,1}}],
  z[{{1,1,1}},{{1,1,1}}],z[{{2},{1}},{{2},{1}}],
  z[{{1,1},{1}},{{1,1},{1}}],
  z[{{1},{1},{1}},{{1},{1},{1}}].
		

Crossrefs

Programs

  • Mathematica
    Table[Length@solidformBTK[{n,3}],{n,3,20}] (* or *)
    g=20;3 CoefficientList[Series[2/Product[(1-x^m)^m,{m,g}]+ 1/(1-x)-(1-x)/Product[(1-x^m)^2,{m,g}]-2/Product[(1-x^m),{m,g}],{x,0,g}],x]

Formula

G.f.: 3*(2*Product_{k>0} 1/(1-x^k)^k -(1-x)*Product_{k>0} 1/(1-x^k)^2 - 2*Product_{k>0} 1/(1-x^k) + 1/(1 - x)) (conjectured).

A378854 a(n) is the number of solid (3D) partitions of n with 2 layers and second layer a plane partition of 2.

Original entry on oeis.org

3, 9, 24, 51, 111, 213, 414, 756, 1374, 2409, 4206, 7152, 12096, 20109, 33198, 54111, 87648, 140520, 223953, 354006, 556419, 868584, 1348857, 2082225, 3198927, 4888944, 7438548, 11265141, 16990077, 25516401, 38175240, 56894490, 84490935, 125028489, 184400952
Offset: 4

Author

Wouter Meeussen, Feb 10 2025

Keywords

Comments

Appears to equal 3*A097391(n-3).

Examples

			a(4)=3 since the 3 solid partitions of {2,2} are:
  [{{2}},{{2}}], [{{1,1}},{{1,1}}], [{{1},{1}},{{1},{1}}]
a(5)=9 since the 9 solid partitions of {3,2} are:
  [{{3}},{{2}}], [{{2,1}},{{2}}], [{{2,1}},{{1,1}}], [{{1,1,1}},{{1,1}}], [{{2},{1}},{{2}}], [{{2},{1}},{{1},{1}}],[{{1,1},{1}},{{1,1}}], [{{1,1},{1}},{{1},{1}}], [{{1},{1},{1}},{{1},{1}}]
		

Crossrefs

Programs

  • Mathematica
    Table[Length@solidformBTK[{n-2,2}],{n,4,22}] (* uses functions defined in link above *)
  • PARI
    A_x(N) = {my(x='x+O('x^N)); Vec(3*x^2*(prod(i=1,N, 1/(1-x^i)^i)-prod(i=1,N, 1/(1-x^i))))}
    A_x(40) \\ John Tyler Rascoe, Feb 20 2025

Formula

From John Tyler Rascoe, Feb 20 2025: (Start)
a(n) = 3*(A000219(n-2) - A000041(n-2)).
G.f.: 3*x^2 * (Product_{i>0} (1/(1-x^i)^i) - Product_{i>0} (1/(1-x^i)^i)). (End)

Extensions

a(23) onwards from John Tyler Rascoe, Feb 20 2025

A380893 Triangle read by rows: T(n,m) = number of solid partitions of n with shape of a plane partition of m.

Original entry on oeis.org

1, 1, 3, 1, 3, 6, 1, 6, 6, 13, 1, 6, 15, 13, 24, 1, 9, 21, 37, 24, 48, 1, 9, 30, 58, 75, 48, 86, 1, 12, 39, 95, 132, 159, 86, 160, 1, 12, 54, 128, 231, 297, 299, 160, 282, 1, 15, 63, 197, 345, 552, 593, 574, 282, 500, 1, 15, 81, 251, 546, 873, 1156, 1180, 1038, 500, 859, 1, 18, 96, 345, 771, 1452, 1933, 2390, 2208, 1874, 859, 1479, 1, 18, 114, 432, 1110, 2151, 3340, 4154, 4614, 4082, 3268, 1479, 2485, 1, 21, 132, 558, 1491, 3276, 5214, 7430, 8310, 8758, 7276, 5685, 2485, 4167
Offset: 1

Author

Wouter Meeussen, Feb 07 2025

Keywords

Comments

A solid (or 3D) partition of n describes a piling of boxes in a corner with heights nonincreasing away from the corner, and containing integers, similarly nonincreasing, that sum to n.
The shape of a solid partitions is defined as the plane partition containing the heights of the piling, irrespective of the numerical content of the boxes.
Row sums equal A000293, T(n,n) = T(n+1,n) equals A000219;
Equals number of solid partitions with total by layer equal to partitions of n with largest part m.

Examples

			Table starts as:
  1,
  1,3
  1,3,6
  1,6,6,13
  1,6,15,13,24
  1,9,21,37,24,48
T(4,2) = 6 since the solid partitions of 4 with shapes a plane partition of 2 are:
 z[{{2,2}}], z[{{3,1}}], z[{{2},{2}}], z[{{3},{1}}], z[{{3}},{{1}}], z[{{2}},{{2}}]
with shapes equal to these plane partitions:
  {{2}}, {{2}}, {{1,1}}, {{1,1}}, {{1},{1}}, {{1},{1}}
		

Crossrefs

Programs

  • Mathematica
    Table[Tr@(Count[First[lapse[#]]&/@Flatten[sols=Table[solidformBTK[par],{par,IntegerPartitions[n]}] ],#]&/@planepartitions[k]),{n,10},{k,n}] (* using functions from link above, or with the faster second program: *)
    Table[ Sum[Length[solidformBTK[TransposePartition@par]],{par,IntegerPartitions[n,{k}]} ],{n,16},{k,n}] (* with transposePartition[par:{Integer..}]:=Count[par,i/;i>=#]&/@Range[Max[par]] *)

A380611 Irregular triangle read by rows: T(r,c) is the product of the number of standard Young tableaux (A117506) and the number of semistandard Young tableaux (A262030) for partitions of r.

Original entry on oeis.org

1, 1, 3, 1, 10, 16, 1, 35, 135, 40, 45, 1, 126, 896, 875, 756, 375, 96, 1, 462, 5250, 10206, 8400, 2450, 14336, 2800, 875, 1701, 175, 1, 1716, 28512, 90552, 74250, 65856, 257250, 48000, 74088, 55566, 102900, 8100, 10976, 5488, 288, 1, 6435, 147147, 686400, 567567, 931392, 3244032, 606375, 194040, 2910600, 1448832, 2673000, 202125, 666792, 846720, 1029000, 491520, 19845, 24696, 65856, 14400, 441, 1
Offset: 0

Author

Wouter Meeussen, Jan 28 2025

Keywords

Comments

Partitions are generated in reverse lexicographic order.
Remark that A262030 uses Abramowitz-Stegun (A-St) order.
Sum of row r equals r^r for r > 0 (Robinson-Schensted correspondence).

Examples

			Triangle begins:
    1;
    1;
    3,    1;
   10,   16,     1;
   35,  135,    40,   45,    1;
  126,  896,   875,  756,  375,    96,    1;
  462, 5250, 10206, 8400, 2450, 14336, 2800, 875, 1701, 175, 1;
  ...
Fourth row is 1*35, 3*45, 2*20, 3*15, 1*1 with sum 256 = 4^4.
		

Crossrefs

Row sums give A000312.
Row lengths give A000041.
Leftmost column gives A088218.

Programs

  • Mathematica
    Needs["Combinatorica`"];
    hooklength[par_?PartitionQ]:=Table[Count[par,q_/;q>=j]+1-i+par[[i]]-j,{i,Length[par]},{j,par[[i]]}];
    countSYT[par_?PartitionQ]:=Tr[par]!/Times@@Flatten[hooklength[par]];
    content[par_?PartitionQ]:=Table[j-i,{i,Length[par]},{j,par[[i]]}];
    countSSYT[par_?PartitionQ,t_Integer_]:=Times@@((t+Flatten[content[par]])/Flatten[hooklength[par]]);
    Table[countSYT[par] countSSYT[par,n],{n,8},{par,IntegerPartitions[n]}]

A379149 Specialization of the Elementary Symmetric Functions e(n) at x_i -> Euler phi(i).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 4, 5, 2, 1, 6, 13, 12, 4, 1, 10, 37, 64, 52, 16, 1, 12, 57, 138, 180, 120, 32, 1, 18, 129, 480, 1008, 1200, 752, 192, 1, 22, 201, 996, 2928, 5232, 5552, 3200, 768, 1, 28, 333, 2202, 8904, 22800, 36944, 36512, 19968, 4608, 1, 32, 445, 3534, 17712, 58416, 128144, 184288, 166016, 84480, 18432
Offset: 0

Author

Wouter Meeussen, Dec 16 2024

Keywords

Comments

Triangular table with alternating signed sum equal to 0 for n>0,
1
1,-1
1,-2,1
1,-4,5,-2
1,-6,13,-12,4
..
and with alternating signed weighted sum (first moment) also equal to 0 for n>1,
0
0,-1
0,-2,2
0,-4,10,-6
0,-6,26,-36,16
..
also when shifting the weights to start at 1,
1
1,-2
1,-4,3
1,-8,15,-8
1,-12,39,-48,20

Examples

			Triangle begins:
  1;
  1,  1;
  1,  2,  1;
  1,  4,  5,   2;
  1,  6, 13,  12,   4;
  1, 10, 37,  64,  52,  16;
  1, 12, 57, 138, 180, 120, 32;
  ...
		

Crossrefs

Columns k=0-1 give: A000012, A002088.
Main diagonal gives A001088.
T(n,n-1) gives A067578.
Cf. A000010.

Programs

  • Maple
    b:= proc(n) option remember; `if`(n=0, 1,
           b(n-1)*(1+x*numtheory[phi](n)))
        end:
    T:= (n, k)-> coeff(b(n), x, k):
    seq(seq(T(n, k), k=0..n), n=0..10);  # Alois P. Heinz, Dec 16 2024
  • Mathematica
    Table[CoefficientList[Expand@Product[z EulerPhi[k]+1,{k,0,n}],z,n+1],{n,0,10}]
  • PARI
    row(n) = Vecrev(prod(k=1, n, 1 + 'x * eulerphi(k))) \\ Andrew Howroyd, Dec 16 2024

Formula

T(n,k) = [x^k] Product_{j=1..n} (1 + x*phi(j)). - Andrew Howroyd, Dec 16 2024

A370852 Irregular triangle read by rows: row n is the list of residues mod n that occur among the Markov numbers.

Original entry on oeis.org

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

Author

Keywords

Comments

Length of row n is A370164(n).

Examples

			The first rows are:
n
1:   0
2:   0  1
3:   1  2
4:   1  2
5:   0  1  2  3  4
6:   1  2  4  5
7:   1  2  5  6
8:   1  2  5
9:   1  2  4  5  7  8
10:  0  1  2  3  4  5  6  7  8  9
11:  1  2  4  5  6  7  9 10
12:  1  2  5 10
13:  0  1  2  3  4  5  6  7  8  9 10 11 12
14:  1  2  5  6  8  9 12 13
15:  1  2  4  5  7  8 10 11 13 14
16:  1  2  5  9 13
17:  0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16
18:  1  2  4  5  7  8 10 11 13 14 16 17
19:  1  2  3  4  5  6  8  9 10 11 13 14 15 16 17 18
20:  1  2  5  6  9 10 13 14 17 18
For n = 14 residues congruent to 0, 3, or 4 mod 7 are forbidden. (See comments to A370164 for explanation.) All other residues occur. For example, the Markov numbers 1, 2, 5, 34, 610, 1325, 194, and 13 produce the residues shown in row 14 of the triangle (mod 14).
		

References

  • Martin Aigner, Markov's theorem and 100 years of the uniqueness conjecture. A mathematical journey from irrational numbers to perfect matchings. Springer, 2013. x+257 pp. ISBN: 978-3-319-00887-5; 978-3-319-00888-2 MR3098784.

Crossrefs

Markov numbers: A002559.
Markov tree: A327345, A368546.
Cf. A370164.

Programs

  • SageMath
    def generateAllMarkovTreeResidues(n):
        row = [[1 % n,5 % n,2 % n]]
        residuesFound = []
        triplesFound = []
        while row != []:
            newRow = []
            for trpl in row:
                if trpl[1] not in residuesFound:
                    residuesFound.append(trpl[1])
                if trpl[2] < trpl[0]:
                    trpl.reverse()
                if trpl not in triplesFound:
                    triplesFound.append(trpl)
                    newRow.append([trpl[0],(3*trpl[0]*trpl[1]-trpl[2]) % n,trpl[1]])
                    newRow.append([trpl[1],(3*trpl[1]*trpl[2]-trpl[0]) % n,trpl[2]])
            row = newRow
        residuesFound.sort()
        return(residuesFound)
    [r for n in range(1,16) for r in generateAllMarkovTreeResidues(n)]

A370164 The number of residues mod n that occur among the Markov numbers.

Original entry on oeis.org

1, 2, 2, 2, 5, 4, 4, 3, 6, 10, 8, 4, 13, 8, 10, 5, 17, 12, 16, 10, 8, 16, 20, 6, 25, 26, 18, 8, 29, 20, 28, 9, 16, 34, 20, 12, 37, 32, 26, 15, 41, 16, 40, 16, 30, 40, 44, 10, 28, 50, 34, 26, 53, 36, 40, 12, 32, 58, 56, 20, 61, 56, 24, 18, 65, 32, 64, 34, 40
Offset: 1

Author

Keywords

Comments

No Markov number is divisible by any prime congruent to 3 (mod 4). Proof: Let m be a Markov number and let p be an odd prime dividing m. Then there exist positive integers b and c such that m^2 + b^2 + c^2 = 3 * m * b * c (Markov's equation). Furthermore it is known that b and c must be coprime to m. Evaluating Markov's equation (mod p) gives b^2 = -c^2 (mod p), which implies -1 is a square (mod p). This is only possible if p = 1 (mod 4).
Given a prime p congruent to 3 (mod 4) and greater than 3, no Markov number is congruent to 1/3 * (p - 2 * sigma(p)) or 2/3 * (p + sigma(p)) (mod p), where sigma(p) = -1 if p is congruent to 7 (mod 12) and sigma(p) = 1 if p is congruent to 11 (mod 12). Proof: Let m be a Markov number congruent to one of the forbidden residues (mod p). Then evaluating Markov's equation (mod p) using the forms of the forbidden residues above and doing a few simplifications implies that -4 is a square (mod p). This contradicts that p = 3 (mod 4).
The first comment above implies that an odd Markov number is congruent to 1 (mod 4). It has also been shown that any positive even integer satisfying all of the constraints of the first two comments above is congruent to 2 (mod 32).
It is conjectured that, modulo n, all residues not forbidden by the constraints in the three comments above are actually realized by some Markov number. Specifically, mod n the forbidden residues include the following: (1) if p is a prime congruent to 3 (mod 4) that divides n, then any residue congruent to 0 (mod p) or congruent to either of the forbidden residues listed in the second comment (mod p), (2) if 4 divides n, then any odd residue congruent to 3 (mod 4), (3) if 2^r is the highest power of 2 dividing n and 2 <= r <= 5, then any even residue not congruent to 2 (mod 2^r), (4) if 2^r is the highest power of 2 dividing n and r > 5, then any even residue not congruent to 2 (mod 32). It is conjectured that all other residues occur. This has been verified for all n <= 38000.
If the conjecture in the previous comment is correct, then it follows from the Chinese remainder theorem that a(n) may be found by writing n = 2^s * 3^t * u * p_1^r_1 * p_2^r_2 * ... * p_k^r_k, where p_1, ..., p_k are distinct primes greater than 3 and congruent to 3 (mod 4) and r_1, ..., r_k are positive and where the prime divisors of u are all congruent to 1 (mod 4). Then a(n) = u * C_s * D_t * Product_{j=1..k} (p_j - 3) * p_j^(r_j - 1), where C_s = 2^s if s < 2, 1 + 2^(s-2) if 2 <= s <= 5, and 2^(s - 5) + 2^(s - 2) if s > 5, and where D_t = 1 if t = 0 and 2 * 3^(t-1) if t > 0. If the conjecture holds then a(n) is a multiplicative function.

Examples

			If n = 56 = 7 * 8 then, since only the residues 1, 2, 5, 6 are allowed (mod 7) and only the residues 1, 2, 5 are allowed (mod 8), the number of potential residues (mod 56) is 4 * 3 = 12, and these residues are 1, 2, 5, 9, 13, 26, 29, 33, 34, 37, 41, 50. That these residues are realized by Markov numbers is witnessed by 1, 2, 5, 233, 13, 194, 29, 89, 34, 1325, 433, 610.
		

References

  • Martin Aigner, Markov's theorem and 100 years of the uniqueness conjecture. A mathematical journey from irrational numbers to perfect matchings. Springer, 2013. x+257 pp. ISBN: 978-3-319-00887-5; 978-3-319-00888-2 MR3098784.

Crossrefs

Markov numbers: A002559.
Markov tree: A327345, A368546.
Triangle giving list of residues mod n: A370852.

Formula

Conjectured: write n = 2^s * 3^t * u * p_1^r_1 * p_2^r_2 * ... * p_k^r_k, where p_1, ..., p_k are distinct primes greater than 3 and congruent to 3 (mod 4) and r_1, ..., r_k are positive and where the prime divisors of u are all congruent to 1 (mod 4). Then a(n) = u * C_s * D_t * Product_{j=1..k} (p_j - 3) * p_j^(r_j - 1), where C_s = 2^s if s < 2, 1 + 2^(s-2) if 2 <= s <= 5, and 2^(s - 5) + 2^(s - 2) if s > 5, and where D_t = 1 if t = 0 and 2 * 3^(t-1) if t > 0.

A362548 Number of partitions of n with at least three parts larger than 1.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 1, 2, 5, 9, 16, 25, 40, 58, 85, 119, 166, 224, 303, 399, 526, 681, 880, 1122, 1430, 1801, 2266, 2827, 3521, 4354, 5378, 6601, 8092, 9870, 12020, 14576, 17652, 21294, 25653, 30804, 36937, 44162, 52732, 62798, 74690, 88627, 105028, 124201, 146696, 172924, 203600, 239292, 280912
Offset: 0

Author

Wouter Meeussen, Apr 24 2023

Keywords

Comments

Both following comments are empirical observations:
1) also accumulant of A119907;
2) the characters of exactly these partitions do not occur in the decomposition of the count of parts 1<=k<=n into the characters of the symmetric group of n (Elders' Theorem).
3) the complement (partitions with no more than 2 parts >1) is counted by A033638.

Crossrefs

Programs

  • Mathematica
    Table[PartitionsP[n]-(1 + Floor[n^2/4]),{n,0,30}];
    Table[ Count[Partitions[n], pa_ /; Length[DeleteCases[pa, 1]] > 2] , {n,0,30}]
  • Python
    from sympy import npartitions
    def A362548(n): return npartitions(n)-1-(n**2>>2) # Chai Wah Wu, Apr 27 2023

Formula

a(n) = A000041(n) - A033638(n).