|
%I #13 Jun 22 2019 09:22:15
%S 2,12,12,16,16,18,20,24,24,24,24,28,24,24,28,32,36,36,32,32,40,36,32,
%T 32,48,48,36,36,32,48,36,48,48,32,48,30,36,48,48,40,48,36,36,48,48,40,
%U 60,40,48,48,64,40,48,48,48,64,40,48,48,36,64,48,64,60,48,48,64,60,36,48
%N Values associated with sequence A078511.
%H Amiram Eldar, <a href="/A088240/b088240.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A000005(A078511(n)) = A001414(A078511(n)).
%t s = {}; Do[If[(d = DivisorSigma[0, n]) == Total[Times @@@ FactorInteger[n]], AppendTo[s, d]], {n, 2, 100000}]; s (* _Amiram Eldar_, Jun 22 2019 *)
%o (PARI) isA078511(n) = {my(f = factor(n)); sum(i=1, #f~, f[i,1]*f[i,2]) == numdiv(n);}
%o lista(nn) = for (n=1, nn, if (isA078511(n), print1(numdiv(n), ", "))); \\ _Michel Marcus_, Feb 11 2015
%Y Cf. A000005, A001414, A078511.
%K nonn
%O 1,1
%A _Lekraj Beedassy_, Nov 03 2003
%E More terms from _Ray Chandler_, Nov 05 2003
|
