%I #36 Nov 15 2017 04:29:14
%S 12,18,20,30,42,56,66,70,78,88,102,104,114,138,174,186,196,222,246,
%T 258,272,282,304,308,318,354,364,366,368,402,426,438,464,474,476,498,
%U 532,534,550,572,582,606,618,642,644,650,654,678,748,762,786,812,822
%N Primitive abundant numbers: abundant numbers (A005101) having no abundant proper divisor.
%C A080224(a(n)) = 1.
%C This is a supersequence of the primitive abundant number sequence A071395, since many of these numbers will be positive integer multiples of the perfect numbers (A000396). - _Timothy L. Tiffin_, Jul 15 2016
%C If the terms of A071395 are removed from this sequence, then the resulting sequence will be A275082. - _Timothy L. Tiffin_, Jul 16 2016
%H T. D. Noe, <a href="/A091191/b091191.txt">Table of n, a(n) for n = 1..10000</a>
%H P. Erdős, <a href="http://www.renyi.hu/~p_erdos/1934-04.pdf">On the density of the abundant numbers</a>, J. London Math. Soc. 9 (1934), pp. 278-282.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/AbundantNumber.html">Abundant Number</a>
%F Erdős shows that a(n) >> n log^2 n. - _Charles R Greathouse IV_, Dec 05 2012
%e 12 is a term since 1, 2, 3, 4, and 6 (the proper divisors of 12) are either deficient or perfect numbers, and thus not abundant. - _Timothy L. Tiffin_, Jul 15 2016
%p isA005101 := proc(n) is(numtheory[sigma](n) > 2*n ); end proc:
%p isA091191 := proc(n) local d; if isA005101(n) then for d in numtheory[divisors](n) minus {1,n} do if isA005101(d) then return false; end if; end do: return true; else false; end if; end proc:
%p for n from 1 to 200 do if isA091191(n) then printf("%d\n",n) ; end if;end do: # _R. J. Mathar_, Mar 28 2011
%t t = {}; n = 1; While[Length[t] < 100, n++; If[DivisorSigma[1, n] > 2*n && Intersection[t, Divisors[n]] == {}, AppendTo[t, n]]]; t (* _T. D. Noe_, Mar 28 2011 *)
%t Select[Range@ 840, DivisorSigma[1, #] > 2 # && Times @@ Boole@ Map[DivisorSigma[1, #] <= 2 # &, Most@ Divisors@ #] == 1 &] (* _Michael De Vlieger_, Jul 16 2016 *)
%o (PARI) is(n)=sumdiv(n,d,sigma(d,-1)>2)==1 \\ _Charles R Greathouse IV_, Dec 05 2012
%o (Haskell)
%o a091191 n = a091191_list !! (n-1)
%o a091191_list = filter f [1..] where
%o f x = sum pdivs > x && all (<= 0) (map (\d -> a000203 d - 2 * d) pdivs)
%o where pdivs = a027751_row x
%o -- _Reinhard Zumkeller_, Jan 31 2014
%Y Cf. A006038 (odd terms), A005101 (abundant numbers), A091192.
%Y Cf. A027751, A071395 (subsequence), supersequence of A275082.
%Y Cf. A294930 (characteristic function), A294890.
%K nonn
%O 1,1
%A _Reinhard Zumkeller_, Dec 27 2003