OFFSET
1,4
COMMENTS
See also graph of A045917 ("Goldbach's comet"). - Altug Alkan, Jan 30 2017
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
Altug Alkan, Illustration Of Residue Classes Modulo 6
EXAMPLE
a(6) = 3 because 2 * 6 = 12 = 2 + 10 = 4 + 8 = 6 + 6 and 2, 4, 6, 8, 10 are in A002202.
MAPLE
N:= 1000: V:= Vector(2*N):
V[1]:= 1:
for n from 2 to 2*N by 2 do
if nops(numtheory:-invphi(n))>1 then V[n]:= 1 fi
od:
C:= map(round, SignalProcessing:-Convolution(V, V)):
seq((C[2*i-1]+V[i])/2, i=1..N); # Robert Israel, Jan 27 2017
PROG
(PARI) a(n) = sum(k=1, n, istotient(k) && istotient(2*n-k));
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Altug Alkan, Jan 27 2017
STATUS
approved