%I #13 Oct 01 2018 03:10:08
%S 2,10,16,22,26,34,46,50,52,58,64,66,70,76,82,86,88,92,94,100,106,116,
%T 118,122,130,134,136,142,146,148,154,166,170,172,178,184,188,190,196,
%U 202,206,208,214,218,220,226,232,236,238,244,246,250,254,262,268,274
%N Even numbers n such that sigma(x) = n has no solution.
%C Contains 2*p if p is a prime not in A005382. - _Robert Israel_, Sep 30 2018
%H Robert Israel, <a href="/A060658/b060658.txt">Table of n, a(n) for n = 1..10000</a>
%e a(2) = 10 because there is no x < 10 whose sigma(x) = 10, an even number.
%p N:= 1000: # to get all terms <= N
%p V:= Vector(N):
%p for x from 1 to N-1 do
%p v:= numtheory:-sigma(x);
%p if v <= N then V[v]:= 1 fi
%p od:
%p select(t -> V[t]=0, [seq(i,i=2..N,2)]); # _Robert Israel_, Sep 30 2018
%t a = {}; Do[ s = DivisorSigma[ 1, n ]; If[ EvenQ[ s ], a = Append[ a, s ] ], {n, 1, 400} ]; Complement[ Table[ 2n, {n, 1, 200} ], Union[ a ] ]
%o (PARI) isok(n) = !(n%2) && !sum(k=1, n, sigma(k)==n); \\ _Michel Marcus_, Mar 08 2017
%Y Cf. A005382, A007369.
%K nonn
%O 1,1
%A _Robert G. Wilson v_, Apr 18 2001