OFFSET
1,8
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
a(n) = Sum_{d|n} A061358(d). - Robert Israel, Oct 17 2023
MAPLE
N:= 100: # for a(1)..a(N)
P:= select(isprime, [2, seq(i, i=3..N, 2)]): nP:= nops(P):
V:= Vector(N):
for i from 1 to nP do
for j from i to nP do
v:= P[i]+P[j];
J:= [seq(t, t=v..N, v)];
V[J]:= V[J] +~ 1
od od:
convert(V, list); # Robert Israel, Oct 17 2023
PROG
(PARI) a(n)=sum(i=1, n, sum(j=1, i, if(n%(prime(i)+prime(j)), 0, 1)))
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Benoit Cloitre, Feb 02 2003
STATUS
approved