|
%I #9 Sep 13 2015 21:11:13
%S 325,2133,10693,16513,19521,51301,159841,176661,214273,306181,1433701,
%T 1570153,1950625,2469601,2924101,5199013,9398593,10445221,15407173,
%U 23548753,28600321,39147301,60110701,62722153,88347781,112803841
%N Nonsemiprime hyperperfect numbers.
%C The other 25 of the first 35 values of A007592 are all semiprimes A001358.
%H Donovan Johnson, <a href="/A133447/b133447.txt">Table of n, a(n) for n = 1..192</a> (terms < 10^12)
%F A100959 INTERSECTION A007592.
%e a(1) = 325 = 5^2 * 13.
%e a(2) = 2133 = 3^3 * 79.
%e a(3) = 10693 = 17^2 * 37.
%e a(4) = 16513 = 7^2 * 337.
%e a(5) = 19521 = 3^4 * 341.
%e a(6) = 51301 = 29^2 * 61.
%e a(7) = 159841 = 11^2 * 1321.
%e a(8) = 176661 = 3^5 * 727.
%e a(9) = 214273 = 47^2 * 97.
%e a(10) = 306181 = 53^2 * 109.
%p a034897 := [] : fd := fopen("b034897.txt",READ) : bf := fscanf(fd,"%d %d") : while nops(bf) <> 0 do a034897 := [op(a034897), op(2,bf) ] ; bf := fscanf(fd,"%d %d") ; od: a007592 := [] : for n in a034897 do m := (n-1)/( numtheory[sigma](n)-n-1) ; if m > 1 then a007592 := [op(a007592),n] ; fi ; od: isA100959 := proc(n) if numtheory[bigomega](n) <> 2 then true ; else false ; end: end: for n in a007592 do if isA100959(n) then printf("%d, ",n) ; fi ; od: # _R. J. Mathar_, Jan 08 2008
%Y Cf. A001358, A007592, A100959.
%K nonn
%O 1,1
%A _Jonathan Vos Post_, Dec 22 2007
%E More terms from _R. J. Mathar_, Jan 08 2008
| 