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”).

A348525
Noninfinitary weird numbers: noninfinitary abundant numbers (A348274) that are not equal to the sum of any subset of their noninfinitary divisors.
3
3344, 12636, 88900, 95900, 109900, 116900, 121100, 181424, 271472, 365552, 476272, 504016, 975568, 1016048, 1354288, 1375504, 1407824, 1552304, 1628528, 1641904, 1862608, 1882672, 1902736, 1909424, 1929488, 1962928, 1982992, 2003056, 2009744, 2029808, 2049872
OFFSET
1,1
LINKS
EXAMPLE
3344 is a term since the sum of its noninfinitary divisors, {2, 4, 8, 22, 38, 44, 76, 88, 152, 418, 836, 1672}, is 3360 > 3344, and no subset of these divisors sums to 3344.
MATHEMATICA
q[n_] := !IntegerQ@ Log2@ DivisorSigma[0, n]; nidiv[1] = {}; nidiv[n_] := Complement[Divisors[n], Sort@ Flatten@ Outer[Times, Sequence @@ (FactorInteger[n] /. {p_, m_Integer} :> p^Select[Range[0, m], BitOr[m, #] == m &])]]; s = {}; Do[If[! q[n], Continue[]]; d = nidiv[n]; If[Total[d] <= n, Continue[]]; c = SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n]; If[c == 0, AppendTo[s, n]], {n, 1, 13000}]; s
CROSSREFS
Similar sequences: A006037, A064114, A292986, A306984, A321146, A327948.
Sequence in context: A270338 A054356 A358312 * A255731 A204884 A204880
KEYWORD
nonn
AUTHOR
Amiram Eldar, Oct 21 2021
EXTENSIONS
More terms from Amiram Eldar, Mar 25 2023
STATUS
approved