OFFSET
1,2
COMMENTS
The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
From Peter Munn, Feb 04 2022: (Start)
For every odd squarefree number, s, exactly one of s and 2s is a term.
Closed under the commutative operation A350066(.,.).
Closed under the commutative operation A059897(.,.) forming a subgroup of the positive integers considered as a group under A059897. As subgroups, this sequence and A028982 are each a transversal of the other.
(End)
FORMULA
{a(n) : n >= 1} = {A019565(A158704(n)) : n >= 1} = {A073675(A319242(n)) : n >= 1}. - Peter Munn, Feb 04 2022
EXAMPLE
30 is the Heinz number of (3,2,1), which is strict and has even weight, so 30 belongs to the sequence.
The sequence of all even-weight strict partitions begins: (), (2), (4), (3,1), (6), (8), (4,2), (5,1), (10), (3,2,1), (7,1), (12), (6,2), (14), (9,1), (16), (5,3), (8,2), (18), (11,1), (5,2,1), (4,3,1).
MATHEMATICA
Select[Range[100], And[SquareFreeQ[#], EvenQ[Total[Cases[FactorInteger[#], {p_, k_}:>k*PrimePi[p]]]]]&]
PROG
(PARI) isok(m) = issquarefree(m) && !(vecsum(apply(primepi, factor(m)[, 1])) % 2); \\ Michel Marcus, Feb 08 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Gus Wiseman, Sep 15 2018
STATUS
approved