OFFSET
1,1
COMMENTS
It appears that this sequence is the union of 3 sets.
First term is 8, and is the only even known value.
Then we get terms that are equal to 25 * b with b a squarefree product of primes congruent to 1 modulo 4 (A002144), that is, terms of A231754.
And we get the following terms 77, 133, 209, 301, 469, 473, 737, 817, 1141, 1273, 1793, 2881, 3097, 7009, 10921. These numbers are the products of 2 distinct primes from this list: 7, 11, 19, 43, 67, 163 (a subsequence of A003173).
LINKS
Hugo Pfoertner, Table of n, a(n) for n = 1..10000
John P. Robertson, Shirali’s Questions About Sums of Residues of Squares
S. A. Shirali, A family portrait of primes-a case study in discrimination, Math. Mag. Vol. 70, No. 4 (Oct., 1997), pp. 263-272.
PROG
(PARI) isok(n) = A231589(n) == n*(n-1)/4 - n;
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Nov 13 2013
STATUS
approved