%I M4133 #42 Jan 12 2023 10:16:09
%S 6,20,28,88,104,272,304,350,368,464,490,496,550,572,650,748,770,910,
%T 945,1184,1190,1312,1330,1376,1430,1504,1575,1610,1696,1870,1888,1952,
%U 2002,2030,2090,2170,2205,2210,2470,2530,2584,2590,2870,2990,3010,3128,3190,3230,3290,3410,3465,3496,3710,3770,3944,4070,4095,4130,4216,4270,4288,4408,4510,4544,4672,4690,4712,4730,4970
%N Primitive pseudoperfect numbers.
%C A primitive pseudoperfect number is a pseudoperfect number that is not a multiple of any other pseudoperfect number.
%C The odd entries so far are identical to the odd primitive abundant A006038. - _Walter Kehowski_, Aug 12 2005
%C Zachariou and Zachariou (1972) called these numbers "irreducible semiperfect numbers". - _Amiram Eldar_, Dec 04 2020
%D Richard K. Guy, Unsolved Problems in Number Theory, 3rd edition, Springer, 2004, Section B2, pp. 74-75.
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Donovan Johnson, <a href="/A006036/b006036.txt">Table of n, a(n) for n = 1..10000</a>
%H Richard K. Guy, <a href="/A001599/a001599_1.pdf">Letter to N. J. A. Sloane with attachment, Jun. 1991</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PrimitivePseudoperfectNumber.html">Primitive Pseudoperfect Number</a>.
%H Andreas Zachariou and Eleni Zachariou, <a href="http://www.hms.gr/apothema/?s=sap&i=261">Perfect, Semi-Perfect and Ore Numbers</a>, Bull. Soc. Math. Grèce (New Ser.), Vol. 13, No. 13A (1972), pp. 12-22; <a href="https://eudml.org/doc/238923">alternative link</a>.
%p with(numtheory): with(combinat): issemiperfect := proc(n) local b, S;
%p b:=false; S:=subsets(divisors(n) minus {n}); while not S[finished] do if
%p convert(S[nextvalue](),`+`)=n then b:=true; break fi od; return b end:
%p L:=remove(proc(z) isprime(z) end,[$1..5000]): PP:=[]: for zz from 1 to 1 do
%p for n in L do if issemiperfect(n) then PP:=[op(PP),n] fi od od;
%p sr := proc(l::list) local x, R, S, P, L; S:=sort(l); R:=[]; P:=S;
%p for x in S do
%p if not(x in R) then
%p L:=selectremove(proc(z) z>x and z mod x = 0 end, P);
%p R:=[op(R),op(L[1])]; P:=L[2];
%p fi; od; return P; end:
%p PPP:=sr(PP); # primitive pseudoperfect numbers less than 5000 # _Walter Kehowski_, Aug 12 2005
%t (* First run one of the programs for A005835 *) A006036 = A005835; curr = 1; max = A005835[[-1]]; While[curr < Length[A006036], currMult = A006036[[curr]]; A006036 = Complement[A006036, Range[2currMult, Ceiling[max/currMult] currMult, currMult]]; curr++]; A006036 (* _Alonso del Arte_, Sep 08 2012 *)
%o (Haskell)
%o a006036 n = a006036_list !! (n-1)
%o a006036_list = filter (all (== 0) . map a210455 . a027751_row) a005835_list
%o -- _Reinhard Zumkeller_, Jan 21 2013
%Y Cf. A005835.
%Y Cf. A210455, A027751.
%K nonn,nice
%O 1,1
%A _N. J. A. Sloane_, _R. K. Guy_
%E More terms from _Walter Kehowski_, Aug 12 2005