|
%I #8 Sep 24 2023 04:30:00
%S 3,5,7,7,11,13,13,11,19,33,19,19,17,23,23,31,34,51,31,31,37,29,37,43,
%T 43,91,95,47,58,43,47,119,53,41,143,74,59,67,67,161,59,82,73,67,71,
%U 209,71,221,59,71,89,185,106,89,221,79,71,122,83,222,103,177,103,101,107,109
%N A064623(n)/s, where s = sum of the prime divisors of A064623(n).
%o (PARI) x = vector(80); B = vector(80); i = 3; count = 0; while (count < 80, i++; if (omega(i) > 1, f = factor(i); s = sum(j = 1, matsize(f)[1], f[j, 1]); if (!(i%s), works = 1; j = 1; while (works && j <= count, if (i%x[j], j++, works = 0)); if (works, count++; x[count] = i; B[count] = i/s)))); B \\ _David Wasserman_, Jan 04 2005
%Y Cf. A064623.
%K nonn
%O 0,1
%A _Amarnath Murthy_, Jun 16 2003
%E More terms from _David Wasserman_, Jan 04 2005
|
