OFFSET
1,2
COMMENTS
Contains 36*k^2 unless 6*k-1 is in A005384.
Are 2, 8, 18, 24, 74, 170, 614, 704, 1010, 1164, 15054, 21774, 24476 the only terms not == 0 or 36 (mod 144)?
LINKS
Robert Israel, Table of n, a(n) for n = 1..652
A. V. Kumchev and J. Y. Liu, Sums of primes and squares of primes in short intervals, Monatsh. Math. 157 (2009), 335-363.
EXAMPLE
24 is in the sequence because 24 == 0 (mod 6) and 24 can't be written as p+q^2 where p and q are primes.
MAPLE
N:= 50000: # to get all terms <= N
P:= select(isprime, [2, seq(i, i=3..N, 2)]):
P2:= select(`<=`, map(`^`, P, 2), N):
PP2:= {seq(seq(s+t, s=P), t=P2)}:
sort(convert({seq(seq(6*i+j, i=1=0..(N-j)/6), j=[0, 2])} minus PP2, list));
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Israel, Feb 12 2019
STATUS
approved