login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A342399
Unitary pseudoperfect numbers k such that no subset of the nontrivial unitary divisors {d|k : 1 < d < k, gcd(d, k/d) = 1} adds up to k.
1
3510, 3770, 5670, 5810, 6790, 7630, 7910, 9590, 9730, 544310, 740870, 2070970, 4017310, 4095190, 5368510, 5569690, 5762330, 5838770, 5855290, 5856130, 5887630, 5902470, 5985770, 6006070, 6039530, 6075370, 6083630, 6181210, 6259610, 6471290, 7038710, 7065730, 7285390
OFFSET
1,1
COMMENTS
Numbers that are the sum of a proper subset of their aliquot unitary divisors but are not the sum of any subset of their nontrivial unitary divisors.
The unitary perfect numbers (A002827) which are a subset of the unitary pseudoperfect numbers (A293188) are excluded from this sequence since otherwise they would all be trivial terms: if k is a unitary perfect number then the sum of the divisors {d|k : 1 < d < k, gcd(d, k/d) = 1} is k-1, so any subset of them has a sum smaller than k.
The unitary pseudoperfect numbers are thus a disjoint union of the unitary perfect numbers, this sequence and A342398.
The unitary abundant numbers (A034683) are a disjoint union of the unitary weird numbers (A064114), this sequence and A342398.
LINKS
EXAMPLE
3510 is a term since it is a unitary pseudoperfect number, 3510 = 1 + 2 + 5 + 13 + 27 + 54 + 65 + 130 + 135 + 270 + 351 + 702 + 1755, and the set of nontrivial unitary divisors of 3510, {d|3510 : 1 < d < 3510, gcd(d, 3510/d) = 1} = {2, 5, 10, 13, 26, 27, 54, 65, 130, 135, 270, 351, 702, 1755}, has no subset that adds up to 3510.
MATHEMATICA
q[n_] := Module[{d = Most @ Select[Divisors[n], CoprimeQ[#, n/#] &], x}, Plus @@ d > n && SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] > 0 && SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, 2, Length[d]}], {x, 0, n}], n] == 0]; Select[Range[10^4], q]
CROSSREFS
The unitary version of A339343.
Subsequence of A034683 and A293188.
Sequence in context: A274237 A053804 A338392 * A043448 A281002 A273472
KEYWORD
nonn
AUTHOR
Amiram Eldar, Mar 10 2021
STATUS
approved