OFFSET
1,1
LINKS
Robert Israel, Table of n, a(n) for n = 1..1849
EXAMPLE
We have sigma(68)=sigma(82)=126, and gcd(68, 82) = 2, hence 126 is in the sequence.
On the other hand, for x in [20, 26, 41], sigma(x)=42, and gcd(20, 26, 41)=1, hence 42 is not here, although gcd(20, 26) is 2.
MAPLE
N:= 10^5: # for terms <= N
V:= Vector(N):
for x from 1 to N do
s:= numtheory:-sigma(x);
if s <= N then
if V[s] = 0 then V[s]:= x
else V[s]:= igcd(V[s], x)
fi
fi
od: select(t -> V[t]=2, [$1..N]); # Robert Israel, Aug 18 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Apr 26 2014
STATUS
approved