OFFSET
1,2
COMMENTS
Define b(0)=2, b(1)=5 and b(n)=6*b(n-1)-b(n-2)-2 for n>1. A prime number p is in the sequence iff (p^2-p-2)/2 is a square iff p=b(n) for some n. The next prime in the sequence is b(21)=8946229758127349, followed by b(n) for n=33, 51, 57 and 75.
a(21) > 2*10^11. - Donovan Johnson, Jul 09 2011
EXAMPLE
The sum of the non-divisors of 14 between 1 and 14 is 3 + 4 + 5 + 6 + 8 + 9 + 10 + 11 + 12 + 13 = 81 = 9^2. 1, 2, 7 & 14 are divisors. Hence 14 is a term of the sequence.
MATHEMATICA
Select[ Range[14*10^6], IntegerQ[Sqrt[(# (# + 1)/2) - DivisorSigma[1, # ]]] &]
CROSSREFS
KEYWORD
nonn
AUTHOR
Joseph L. Pe, Oct 22 2002
EXTENSIONS
Edited by Robert G. Wilson v and Dean Hickerson, Oct 25 2002
a(16)-a(17) from Donovan Johnson, Oct 14 2009
a(18)-a(20) from Donovan Johnson, Jul 09 2011
STATUS
approved