A299926 -id:A299926 - cp's OEIS frontend

A000898 a(n) = 2(a(n-1) + (n-1)a(n-2)) for n >= 2 with a(0) = 1.

1, 2, 6, 20, 76, 312, 1384, 6512, 32400, 168992, 921184, 5222208, 30710464, 186753920, 1171979904, 7573069568, 50305536256, 342949298688, 2396286830080, 17138748412928, 125336396368896, 936222729254912, 7136574106003456, 55466948299223040, 439216305474605056, 3540846129311916032

Offset: 0

N. J. A. Sloane

Number of solutions to the rook problem on a 2n X 2n board having a certain symmetry group (see Robinson for details).
Also the value of the n-th derivative of exp(x^2) evaluated at 1. - N. Calkin, Apr 22 2010
For n >= 1, a(n) is also the sum of the degrees of the irreducible representations of the group of n X n signed permutation matrices (described in sequence A066051). The similar sum for the "ordinary" symmetric group S_n is in sequence A000085. - Sharon Sela (sharonsela(AT)hotmail.com), Jan 12 2002
It appears that this is also the number of permutations of 1, 2, ..., n+1 such that each term (after the first) is within 2 of some preceding term. Verified for n+1 <= 6. E.g., a(4) = 20 because of the 24 permutations of 1, 2, 3, 4, the only ones not permitted are 1, 4, 2, 3; 1, 4, 3, 2; 4, 1, 2, 3; and 4, 1, 3, 2. - Gerry Myerson, Aug 06 2003
Hankel transform is A108400. - Paul Barry, Feb 11 2008
From Emeric Deutsch, Jun 19 2010: (Start)
Number of symmetric involutions of [2n]. Example: a(2)=6 because we have 1234, 2143, 1324, 3412, 4231, and 4321. See the Egge reference, pp. 419-420.
Number of symmetric involutions of [2n+1]. Example: a(2)=6 because we have 12345, 14325, 21354, 45312, 52341, and 54321. See the Egge reference, pp. 419-420.
(End)
Binomial convolution of sequence A000085: a(n) = Sum_{k=0..n} binomial(n,k)*A000085(k)*A000085(n-k). - Emanuele Munarini, Mar 02 2016
The sequence can be obtained from the infinite product of 2 X 2 matrices [(1,N); (1,1)] by extracting the upper left terms, where N = (1, 3, 5, ...), the odd integers. - Gary W. Adamson, Jul 28 2016
Apparently a(n) is the number of standard domino tableaux of size 2n, where a domino tableau is a generalized Young tableau in which all rows and columns are weakly increasing and all regions are dominos. - Gus Wiseman, Feb 25 2018

			G.f. = 1 + 2*x + 6*x^2 + 20*x^3 + 76*x^4 + 312*x^5 + 1384*x^6 + 6512*x^7 + ...
The a(3) = 20 domino tableaux:
1 1 2 2 3 3
.
1 2 2 3 3
1
.
1 2 3 3   1 1 3 3   1 1 2 2
1 2       2 2       3 3
.
1 1 3 3   1 1 2 2
2         3
2         3
.
1 2 3   1 2 2   1 1 3
1 2 3   1 3 3   2 2 3
.
1 3 3   1 2 2
1       1
2       3
2       3
.
1 2   1 1   1 1
1 2   2 3   2 2
3 3   2 3   3 3
.
1 3   1 2   1 1
1 3   1 2   2 2
2     3     3
2     3     3
.
1 1
2
2
3
3
.
1
1
2
2
3
3 - _Gus Wiseman_, Feb 25 2018

D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 3, Sect 5.1.4 Exer. 31.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181.
R. W. Robinson, Counting arrangements of bishops, pp. 198-214 of Combinatorial Mathematics IV (Adelaide 1975), Lect. Notes Math., 560 (1976).
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n=0..200
Y. Alp and E. G. Kocer, Exponential Almost-Riordan Arrays, Results Math 79, 173 (2024). See page 10.
Arvind Ayyer, Hiranya Kishore Dey, and Digjoy Paul, How large is the character degree sum compared to the character table sum for a finite group?, arXiv:2406.06036 [math.RT], 2024. See p. 10.
C. Banderier, M. Bousquet-Mélou, A. Denise, P. Flajolet, D. Gardy and D. Gouyou-Beauchamps, Generating functions for generating trees, Discrete Mathematics 246(1-3), March 2002, pp. 29-55.
Paul Barry, The Gamma-Vectors of Pascal-like Triangles Defined by Riordan Arrays, arXiv:1804.05027 [math.CO], 2018.
R. A. Brualdi, Shi-Mie Ma, Enumeration of involutions by descents and symmetric matrices, Eur. J. Combin. 43 (2015) 220-228
P. J. Cameron, Sequences realized by oligomorphic permutation groups, J. Integ. Seqs. Vol. 3 (2000), #00.1.5.
C.-O. Chow, Counting involutory, unimodal, and alternating signed permutations, Discr. Math. 306 (2006), 2222-2228.
R. Donaghey, Binomial self-inverse sequences and tangent coefficients, J. Combin. Theory, Series A, 21 (1976), 155-163.
Eric S. Egge, Restricted symmetric permutations, Ann. Combin., 11 (2007), 405-434.
Adam M. Goyt and Lara K. Pudwell, Avoiding colored partitions of two elements in the pattern sense, arXiv preprint arXiv:1203.3786 [math.CO], 2012. - From _N. J. A. Sloane_, Sep 17 2012
T. Halverson and M. Reeks, Gelfand Models for Diagram Algebras, arXiv:1302.6150 [math.RT], 2013-2014.
Guo-Niu Han, Enumeration of Standard Puzzles, 2011. [Cached copy]
Guo-Niu Han, Enumeration of Standard Puzzles, arXiv:2006.14070 [math.CO], 2020.
L. C. Larson, The number of essentially different nonattacking rook arrangements, J. Recreat. Math., 7 (No. 3, 1974), circa pages 180-181. [Annotated scan of pages 180 and 181 only]
Yen-chi R. Lin, Asymptotic Formula for Symmetric Involutions, arXiv:1310.0988 [math.CO], 2013.
E. Lucas, Théorie des Nombres, Gauthier-Villars, Paris, 1891, Vol. 1, p. 221.
E. Lucas, Théorie des nombres (annotated scans of a few selected pages)
J. Riordan, Letter to N. J. A. Sloane, Feb 03 1975 (with notes by njas)
D. P. Roberts and A. Venkatesh, Hurwitz monodromy and full number fields, 2014.
Index entries for sequences related to Hermite polynomials

Column k=2 of A294042 and A376826.
Row sums of A122832 and A185296.
Cf. A000085, A000712, A004003, A066051, A099390, A100861, A135401, A138178, A153452, A297388, A299699, A299926, A300056, A300060.

a000898 n = a000898_list !! n
a000898_list = 1 : 2 : (map (* 2) $
   zipWith (+) (tail a000898_list) (zipWith (*) [1..] a000898_list))
-- Reinhard Zumkeller, Oct 10 2011

# For Maple program see A000903.
seq(simplify((-I)^n*HermiteH(n, I)), n=0..25); # Peter Luschny, Oct 23 2015

a[n_] := Sum[ 2^k*StirlingS1[n, k]*BellB[k], {k, 0, n}]; Table[a[n], {n, 0, 21}] (* Jean-François Alcover, Nov 17 2011, after Vladeta Jovovic *)
RecurrenceTable[{a[0]==1,a[1]==2,a[n]==2(a[n-1]+(n-1)a[n-2])},a,{n,30}] (* Harvey P. Dale, Aug 04 2012 *)
Table[Abs[HermiteH[n, I]], {n, 0, 20}] (* Vladimir Reshetnikov, Oct 22 2015 *)
a[ n_] := Sum[ 2^(n - 2 k) n! / (k! (n - 2 k)!), {k, 0, n/2}]; (* Michael Somos, Oct 23 2015 *)

makelist((%i)^n*hermite(n,-%i),n,0,12); /* Emanuele Munarini, Mar 02 2016 */

{a(n) = if( n<0, 0, n! * polcoeff( exp(2*x + x^2 + x * O(x^n)), n))}; /* Michael Somos, Feb 08 2004 */

{a(n) = if( n<2, max(0, n+1), 2*a(n-1) + (2*n - 2) * a(n-2))}; /* Michael Somos, Feb 08 2004 */

my(x='x+O('x^66)); Vec(serlaplace(exp(2*x+x^2))) \\ Joerg Arndt, Oct 04 2013

{a(n) = sum(k=0, n\2, 2^(n - 2*k) * n! / (k! * (n - 2*k)!))}; /* Michael Somos, Oct 23 2015 */

a(n) = Sum_{m=0..n} |A060821(n,m)| = H(n,-i)*i^n, with the Hermite polynomials H(n,x); i.e., these are row sums of the unsigned triangle A060821.
E.g.f.: exp(x*(x + 2)).
a(n) = 2 * A000902(n) for n >= 1.
a(n) = Sum_{k=0..n} binomial(n,2k)*binomial(2k,k)*k!*2^(n-2k). - N. Calkin, Apr 22 2010
Binomial transform of A047974. - Paul Barry, May 09 2003
a(n) = Sum_{k=0..n} Stirling1(n, k)*2^k*Bell(k). - Vladeta Jovovic, Oct 01 2003
From Paul Barry, Aug 29 2005: (Start)
a(n) = Sum_{k=0..floor(n/2)} A001498(n-k, k) * 2^(n-k).
a(n) = Sum_{k=0..n} A001498((n+k)/2, (n-k)/2) * 2^((n+k)/2) * (1+(-1)^(n-k))/2. (End)
For asymptotics, see the Robinson paper. [This is disputed by Yen-chi R. Lin. See below, Sep 30 2013.]
a(n) = Sum_{k=0..floor(n/2)} 2^(n-2*k) * C(n,2*k) * (2*k)!/k!. - Paul Barry, Feb 11 2008
G.f.: 1/(1 - 2*x - 2*x^2/(1 - 2*x - 4*x^2/(1 - 2*x - 6*x^2/(1 - 2*x - 8*x^2/(1 - ... (continued fraction). - Paul Barry, Feb 25 2010
E.g.f.: exp(x^2 + 2*x) = Q(0); Q(k) = 1 + (x^2 + 2*x)/(2*k + 1 - (x^2 + 2*x)*(2*k + 1)/((x^2 + 2*x) + (2*k + 2)/Q(k+1))); (continued fraction). - Sergei N. Gladkovskii, Nov 24 2011
G.f.: 1/Q(0), where Q(k) = 1 + 2*x*k - x - x/(1 - 2*x*(k + 1)/Q(k+1) ); (continued fraction). - Sergei N. Gladkovskii, Mar 07 2013
a(n) = (2*n/e)^(n/2) * exp(sqrt(2*n)) / sqrt(2*e) * (1 + sqrt(2/n)/3 + O(n^(-1))). - Yen-chi R. Lin, Sep 30 2013
0 = a(n)*(2*a(n+1) + 2*a(n+2) - a(n+3)) + a(n+1)*(-2*a(n+1) + a(n+2)) for all n >= 0. - Michael Somos, Oct 23 2015
a(n) = Sum_{k=0..floor(n/2)} 2^(n-k)*B(n, k), where B are the Bessel numbers A100861. - Peter Luschny, Jun 04 2021

More terms from Larry Reeves (larryr(AT)acm.org), Feb 21 2001
Initial condition a(0)=1 added to definition by Jon E. Schoenfield, Oct 01 2013
More terms from Joerg Arndt, Oct 04 2013

A300121 Number of normal generalized Young tableaux, of shape the integer partition with Heinz number n, with all rows and columns weakly increasing and all regions connected skew partitions.

Original entry on oeis.org

1, 1, 2, 2, 4, 5, 8, 4, 11, 12, 16, 12, 32, 28, 31, 8, 64, 31, 128, 33, 82, 64, 256, 28, 69, 144, 69, 86, 512, 105, 1024, 16, 208, 320, 209, 82, 2048, 704, 512, 86, 4096, 318, 8192, 216, 262, 1536, 16384, 64, 465, 262, 1232, 528, 32768, 209, 588, 245, 2912, 3328

Offset: 1

Gus Wiseman, Feb 25 2018

nonn

The diagram of a connected skew partition is required to be connected as a polyomino but can have empty rows or columns. A generalized Young tableau of shape y is an array obtained by replacing the dots in the Ferrers diagram of y with positive integers. A tableau is normal if its entries span an initial interval of positive integers. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			The a(9) = 11 tableaux:
1 1
1 1
.
2 1   1 1   1 1   1 2
1 1   1 2   2 2   1 2
.
1 1   1 2   1 2   1 3
2 3   1 3   3 3   2 3
.
1 2   1 3
3 4   2 4

Cf. A000085, A000898, A056239, A006958, A138178, A153452, A238690, A259479, A259480, A296150, A296561, A297388, A299699, A299925, A299926, A300056, A300060, A300118, A300120, A300122, A300123, A300124.

undcon[y_]:=Select[Tuples[Range[0,#]&/@y],Function[v,GreaterEqual@@v&&With[{r=Select[Range[Length[y]],y[[#]]=!=v[[#]]&]},Or[Length[r]<=1,And@@Table[v[[i]]Table[PrimePi[p],{k}]]]];
Table[Length[cos[Reverse[primeMS[n]]]],{n,50}]

A238690 Let each integer m (1 <= m <= n) be factorized as m = prime_m(1)prime_m(2)...*prime_m(bigomega(m)), with the primes sorted in nonincreasing order. Then a(n) is the number of values of m such that each prime_m(i) <= prime_n(i).

Original entry on oeis.org

1, 2, 3, 3, 4, 5, 5, 4, 6, 7, 6, 7, 7, 9, 9, 5, 8, 9, 9, 10, 12, 11, 10, 9, 10, 13, 10, 13, 11, 14, 12, 6, 15, 15, 14, 12, 13, 17, 18, 13, 14, 19, 15, 16, 16, 19, 16, 11, 15, 16, 21, 19, 17, 14, 18, 17, 24, 21, 18, 19, 19, 23, 22, 7, 22, 24, 20, 22, 27, 23, 21

Offset: 1

Matthew Vandermast, Apr 28 2014

nonn

Equivalently, a(n) equals the number of values of m such that each value of A238689 T(m,k) <= A238689 T(n,k). (Since the prime factorization of 1 is the empty factorization, we consider each prime_1(i) not to be greater than prime_n(i) for all positive integers n.)
Suppose we say that n "covers" m iff both m and n are factorized as described in the sequence definition and each prime_m(i) <= prime_n(i). At least three sequences (A037019, A108951 and A181821) have the property that a(m) divides a(n) iff n "covers" m.  These sequences are also divisibility sequences (i.e., sequences with the property that a(m) divides a(n) if m divides n), since any positive integer "covers" each of its divisors.
For any positive integers m and k, the following integer sequences (with n >= 0) are arithmetic progressions:
1. The sequence b(n) = a(m*(2^n)).
2. The sequence b(n) = a(m*(prime(n+k))) if prime(k) >= A006530(m).
Also, a(n) = the number of distinct prime signatures that occur among the divisors of any integer m such that A181819(m) = n and/or A238745(m) = n.
Number of skew partitions whose numerator has Heinz number n, where a skew partition is a pair y/v of integer partitions such that the diagram of v fits inside the diagram of y. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Feb 24 2018

			The prime factorizations of integers 1 through 9, with prime factors sorted from largest to smallest:
1 - the empty factorization (no prime factors)
2 = 2
3 = 3
4 = 2*2
5 = 5
6 = 3*2
7 = 7
8 = 2*2*2
9 = 3*3
To find a(9), we consider 9 = 3*3. There are 6 positive integers (1, 2, 3, 4, 6 and 9) which satisfy the following criteria:
1) The largest prime factor, if one exists, is not greater than 3;
2) The second-largest prime factor, if one exists, is not greater than 3;
3) The total number of prime factors (counting repeated factors) does not exceed 2.
Therefore, a(9) = 6.
From _Gus Wiseman_, Feb 24 2018: (Start)
Heinz numbers of the a(15) = 9 partitions contained within the partition (32) are 1, 2, 3, 4, 5, 6, 9, 10, 15. The a(15) = 9 skew partitions are (32)/(), (32)/(1), (32)/(11), (32)/(2), (32)/(21), (32)/(22), (32)/(3), (32)/(31), (32)/(32).
Corresponding diagrams are:
  o o o   . o o   . o o   . . o   . . o   . . o   . . .   . . .   . . .
  o o     o o     . o     o o     . o     . .     o o     . o     . .    (End)

Amiram Eldar, Table of n, a(n) for n = 1..10000

Rearrangement of A115728, A115729 and A238746. A116473(n) is the number of times n appears in the sequence.

Cf. A000041, A000085, A000720, A056239, A063834, A112798, A122111, A153452, A215366, A238689, A259478, A259480, A296150, A296188, A296561, A297388, A299925, A299926, A299966, A299967.

undptns[y_]:=Select[Tuples[Range[0,#]&/@y],OrderedQ[#,GreaterEqual]&];
primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length[undptns[Reverse[primeMS[n]]]],{n,100}] (* Gus Wiseman, Feb 24 2018 *)

a(n) = A085082(A108951(n)) = A085082(A181821(n)).
a(n) = a(A122111(n)).
a(prime(n)) = a(2^n) = n+1.
a((prime(n))^m) = a((prime(m))^n) = binomial(n+m, n).
a(A002110(n)) = A000108(n+1).
A000005(n) <= a(n) <= n.

A299968 Number of normal generalized Young tableaux of size n with all rows and columns strictly increasing.

Original entry on oeis.org

1, 1, 2, 5, 15, 51, 189, 753, 3248, 14738, 70658, 354178, 1857703, 10121033, 57224955, 334321008, 2017234773, 12530668585, 80083779383, 525284893144, 3533663143981, 24336720018666, 171484380988738, 1234596183001927, 9075879776056533, 68052896425955296

Offset: 0

Gus Wiseman, Feb 26 2018

nonn

			The a(4) = 15 tableaux:
1 2 3 4
.
1 2 3   1 2 4   1 3 4   1 2 3   1 2 3
4       3       2       2       3
.
1 2   1 3   1 2
3 4   2 4   2 3
.
1 2   1 3   1 2   1 4   1 3
3     2     2     2     2
4     4     3     3     3
.
1
2
3
4

Alois P. Heinz, Table of n, a(n) for n = 0..50 (first 46 terms from Ludovic Schwob)
D. E. Knuth, Permutations, matrices, and generalized Young tableaux, Pacific Journal of Mathematics, Vol. 34, No. 3 (1970), 709-727.
Wikipedia, Young tableau

Cf. A000085, A000898, A001221, A005117, A006958, A015128, A138178, A238121, A238690, A285175, A296561, A297388, A299926, A300120, A300122.

unddis[y_]:=DeleteCases[y-#,0]&/@Tuples[Table[If[y[[i]]>Append[y,0][[i+1]],{0,1},{0}],{i,Length[y]}]];
dos[y_]:=With[{sam=Rest[unddis[y]]},If[Length[sam]===0,If[Total[y]===0,{{}},{}],Join@@Table[Prepend[#,y]&/@dos[sam[[k]]],{k,1,Length[sam]}]]];
Table[Sum[Length[dos[y]],{y,IntegerPartitions[n]}],{n,1,8}]

a(n) = Sum_{k=0..n} 2^k * A238121(n,k). - Ludovic Schwob, Sep 23 2023

More terms from Ludovic Schwob, Sep 23 2023

A026010 a(n) = number of (s(0), s(1), ..., s(n)) such that s(i) is a nonnegative integer and |s(i) - s(i-1)| = 1 for i = 1,2,...,n and s(0) = 2. Also a(n) = sum of numbers in row n+1 of array T defined in A026009.

Original entry on oeis.org

1, 2, 4, 7, 14, 25, 50, 91, 182, 336, 672, 1254, 2508, 4719, 9438, 17875, 35750, 68068, 136136, 260338, 520676, 999362, 1998724, 3848222, 7696444, 14858000, 29716000, 57500460, 115000920, 222981435, 445962870, 866262915, 1732525830, 3370764540

Offset: 0

Clark Kimberling

nonn

Conjecture: a(n) is the number of integer compositions of n + 2 in which the even parts appear as often at even positions as at odd positions (confirmed up to n = 19). - Gus Wiseman, Mar 17 2018

			The a(3) = 7 compositions of 5 in which the even parts appear as often at even positions as at odd positions are (5), (311), (131), (113), (221), (122), (11111). Missing are (41), (14), (32), (23), (212), (2111), (1211), (1121), (1112). - _Gus Wiseman_, Mar 17 2018

Michael De Vlieger, Table of n, a(n) for n = 0..3326
Christian Krattenthaler, Daniel Yaqubi, Some determinants of path generating functions, II, arXiv:1802.05990 [math.CO], 2018; Adv. Appl. Math. 101 (2018), 232-265.

First differences of A050168. Pairwise sums of A037952.

Cf. A000712, A001405, A005774, A045931, A063886, A097613, A130780, A171966, A239241, A299926, A300061, A300787, A300788, A300789.

[(&+[Binomial(Floor((n+k)/2), Floor(k/2)): k in [0..n]]): n in [0..40]]; // G. C. Greubel, Nov 08 2018

Array[Sum[Binomial[Floor[(# + k)/2], Floor[k/2]], {k, 0, #}] &, 34, 0] (* Michael De Vlieger, May 16 2018 *)
Table[2^(-1 + n)*(((2 + 3*#)*Gamma[(1 + #)/2])/(Sqrt[Pi]*Gamma[2 + #/2]) &[n + Mod[n, 2]]), {n,0,40}] (* Peter Pein, Nov 08 2018 *)
Table[(1/2)^((5 - (-1)^n)/2)*(6*n + 7 - 3*(-1)^n)*CatalanNumber[(2*n + 1 - (-1)^n)/4], {n, 0, 40}] (* G. C. Greubel, Nov 08 2018 *)

vector(40, n, n--; sum(k=0,n, binomial(floor((n+k)/2), floor(k/2)))) \\ G. C. Greubel, Nov 08 2018

a(2*n) = ((3*n + 1)/(2*n + 1))*C(2*n + 1, n)= A051924(1+n), n>=0, a(2*n-1) = a(2*n)/2 = A097613(1+n), n >= 1. - Herbert Kociemba, May 08 2004
a(n) = Sum_{k=0..n} binomial(floor((n+k)/2), floor(k/2)). - Paul Barry, Jul 15 2004
Inverse binomial transform of A005774: (1, 3, 9, 26, 75, 216, ...). - Gary W. Adamson, Oct 22 2007
Conjecture: (n+3)*a(n) - 2*a(n-1) + (-5*n-3)*a(n-2) + 2*a(n-3) + 4*(n-3)*a(n-4) = 0. - R. J. Mathar, Jun 20 2013
a(n) = (1/2)^((5 - (-1)^n)/2)*(6*n + 7 - 3*(-1)^n)*Catalan((2*n + 1 - (-1)^n)/4), where Catalan is the Catalan number = A000108. - G. C. Greubel, Nov 08 2018

A300060 Number of domino tilings of the diagram of the integer partition with Heinz number n.

Original entry on oeis.org

1, 0, 1, 1, 0, 0, 1, 0, 2, 1, 0, 1, 1, 0, 0, 1, 0, 0, 1, 0, 2, 1, 0, 0, 3, 0, 3, 1, 1, 0, 0, 0, 0, 1, 0, 2, 1, 0, 2, 1, 0, 0, 1, 0, 0, 1, 0, 1, 5, 0, 0, 1, 1, 0, 3, 0, 2, 0, 0, 0, 1, 1, 3, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 4, 1, 0, 0, 1, 0, 5, 1, 0, 2, 3, 0, 2, 1, 1, 1, 5, 0, 0, 1, 0, 0, 0, 0, 0, 3, 1, 0, 0, 0, 0, 0

Offset: 1

Gus Wiseman, Feb 23 2018

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Alois P. Heinz, Table of n, a(n) for n = 1..100000
Wikipedia, Domino tiling
Gus Wiseman, The a(91) = 5 domino tilings of (6,4).

Cf. A000085, A000720, A000712, A000898, A001222, A004003, A056239, A099390, A138178, A153452, A238690, A296150, A296188, A299925, A299926, A300056, A300061, A304662.

h:= proc(l, f) option remember; local k; if min(l[])>0 then
     `if`(nops(f)=0, 1, h(map(x-> x-1, l[1..f[1]]), subsop(1=[][], f)))
    else for k from nops(l) while l[k]>0 by -1 do od;
        `if`(nops(f)>0 and f[1]>=k, h(subsop(k=2, l), f), 0)+
        `if`(k>1 and l[k-1]=0, h(subsop(k=1, k-1=1, l), f), 0)
      fi
    end:
g:= l-> `if`(add(`if`(l[i]::odd, (-1)^i, 0), i=1..nops(l))=0,
        `if`(l=[], 1, h([0$l[1]], subsop(1=[][], l))), 0):
a:= n-> g(sort(map(i-> numtheory[pi](i[1])$i[2], ifactors(n)[2]), `>`)):
seq(a(n), n=1..120);  # Alois P. Heinz, May 22 2018

h[l_, f_] := h[l, f] = Module[{k}, If[Min[l] > 0, If[Length[f] == 0, 1, h[Map[Function[x, x-1], l[[Range @ f[[1]]]]], ReplacePart[f, 1 -> Nothing]]], For[k = Length[l], l[[k]] > 0, k-- ]; If[Length[f] > 0 && f[[1]] >= k, h[ReplacePart[l, k -> 2], f], 0] + If[k > 1 && l[[k-1]] == 0, h[ReplacePart[l, {k -> 1, k - 1 -> 1}], f], 0]]];
g[l_] := If[Sum[If[OddQ @ l[[i]], (-1)^i, 0], {i, 1, Length[l]}] == 0, If[l == {}, 1, h[Table[0, l[[1]]], ReplacePart[l, 1 -> Nothing]]], 0];
a[n_] := g[Reverse @ Sort[ Flatten[ Map[ Function[i, Table[PrimePi[i[[1]]], i[[2]]]], FactorInteger[n]]]]];
Array[a, 120] (* Jean-François Alcover, May 28 2018, after Alois P. Heinz *)

A generalized Young tableau of shape y is an array obtained by replacing the dots in the Ferrers diagram of y with positive integers. A tableau is normal if its entries span an initial interval of positive integers. A standard domino tableau is a generalized Young tableau in which all rows and columns are weakly increasing and all regions are dominos. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

The diagram of a connected skew partition is required to be connected as a polyomino but can have empty rows or columns. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

cp's OEIS Frontend

A000898 a(n) = 2*(a(n-1) + (n-1)*a(n-2)) for n >= 2 with a(0) = 1.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

Maxima

PARI

PARI

PARI

PARI

Formula

Extensions

A300121 Number of normal generalized Young tableaux, of shape the integer partition with Heinz number n, with all rows and columns weakly increasing and all regions connected skew partitions.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A238690 Let each integer m (1 <= m <= n) be factorized as m = prime_m(1)*prime_m(2)*...*prime_m(bigomega(m)), with the primes sorted in nonincreasing order. Then a(n) is the number of values of m such that each prime_m(i) <= prime_n(i).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A299968 Number of normal generalized Young tableaux of size n with all rows and columns strictly increasing.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A026010 a(n) = number of (s(0), s(1), ..., s(n)) such that s(i) is a nonnegative integer and |s(i) - s(i-1)| = 1 for i = 1,2,...,n and s(0) = 2. Also a(n) = sum of numbers in row n+1 of array T defined in A026009.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A300060 Number of domino tilings of the diagram of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

A300056 Number of normal standard domino tableaux whose shape is the integer partition with Heinz number n.

Views

Author

Keywords

A000898 a(n) = 2(a(n-1) + (n-1)a(n-2)) for n >= 2 with a(0) = 1.

A238690 Let each integer m (1 <= m <= n) be factorized as m = prime_m(1)prime_m(2)...*prime_m(bigomega(m)), with the primes sorted in nonincreasing order. Then a(n) is the number of values of m such that each prime_m(i) <= prime_n(i).