A097986 Number of strict integer partitions of n with a part dividing all the other parts.
1, 1, 2, 2, 2, 4, 3, 5, 5, 7, 6, 12, 9, 13, 15, 20, 18, 28, 26, 37, 39, 47, 49, 71, 68, 85, 94, 117, 120, 159, 160, 201, 216, 257, 277, 348, 357, 430, 470, 562, 592, 720, 758, 901, 981, 1134, 1220, 1457, 1542, 1798, 1952, 2250, 2419, 2819, 3023, 3482, 3773, 4291
Offset: 1
Examples
From _Gus Wiseman_, Dec 01 2023: (Start) The a(1) = 1 through a(8) = 5 strict partitions with a part dividing all the other parts: (1) (2) (3) (4) (5) (6) (7) (8) (2,1) (3,1) (4,1) (4,2) (6,1) (6,2) (5,1) (4,2,1) (7,1) (3,2,1) (4,3,1) (5,2,1) The a(1) = 1 through a(8) = 5 uniform partitions containing 1: (1) (11) (21) (31) (41) (51) (61) (71) (111) (1111) (11111) (321) (421) (431) (2211) (1111111) (521) (111111) (3311) (11111111) (End)
Links
- Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 500 terms from John Tyler Rascoe)
Crossrefs
The non-strict version is A083710.
The case with no 1's is A098965.
The Heinz numbers of these partitions are A339563.
The strict complement is counted by A341450.
The version for "divisible by" instead of "dividing" is A343347.
The case where there is also a part divisible by all the others is A343378.
The case where there is no part divisible by all the others is A343381.
A000005 counts divisors.
A000009 counts strict partitions.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A015723 counts strict partitions with a selected part.
A167865 counts strict chains of divisors > 1 summing to n.
Programs
-
Mathematica
Take[ CoefficientList[ Expand[ Sum[x^k*Product[1 + x^(k*i), {i, 2, 62}], {k, 62}]], x], {2, 60}] (* Robert G. Wilson v, Nov 01 2004 *) Table[Length[Select[IntegerPartitions[n], UnsameQ@@#&&Or@@Table[And@@IntegerQ/@(#/x), {x,#}]&]], {n,0,30}] (* Gus Wiseman, Apr 23 2021 *)
-
PARI
A_x(N) = {my(x='x+O('x^N)); Vec(sum(k=1,N,x^k*prod(i=2,N-k, (1+x^(k*i)))))} A_x(50) \\ John Tyler Rascoe, Nov 19 2024
Formula
a(n) = Sum_{d|n} A025147(d-1).
G.f.: Sum_{k>=1} (x^k*Product_{i>=2} (1+x^(k*i))).
a(n) ~ exp(Pi*sqrt(n/3)) / (8*3^(1/4)*n^(3/4)). - Vaclav Kotesovec, Jul 06 2025
Extensions
More terms from Robert G. Wilson v, Nov 01 2004
Name shortened by Gus Wiseman, Apr 23 2021
Comments