%I #6 Aug 17 2020 20:52:14
%S 4,9,16,18,20,32,36,44,45,48,49,64,68,72,80,81,90,98,99,100,112,116,
%T 124,144,153,160,162,164,169,176,180,192,196,198,208,220,225,236,240,
%U 244,245,252,256,261,279,284,288,292,304,306,320,324,336,338,340,352,356,360,361,369,392,396,400,404,405,428,432,441
%N Numbers k such that A000005(k) does not divide A003973(k); numbers k for which A336839(k) > 1.
%C Numbers k such that A003961(k) is not in A003601, but in A049642.
%H Antti Karttunen, <a href="/A336919/b336919.txt">Table of n, a(n) for n = 1..20000</a>
%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>
%o (PARI)
%o A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
%o isA336919(n) = !!(sigma(A003961(n))%numdiv(n));
%Y Cf. A000005, A003601, A003961, A003973, A049642, A336839.
%Y Cf. A336918 (complement).
%K nonn
%O 1,1
%A _Antti Karttunen_, Aug 12 2020