OFFSET
1,1
COMMENTS
This seems to have a disproportionately high probability of generating a prime number.
FORMULA
a(n) = A063533(n) - 1.
EXAMPLE
a(1) = prime(1)^2 + (prime(1)+2)^2 - 1 = 3^2 + (3+2)^2 - 1 = 33;
a(2) = prime(2)^2 + (prime(2)+2)^2 - 1 = 5^2 + (5+2)^2 - 1 = 73;
a(3) = prime(3)^2 + (prime(3)+2)^2 - 1 = 11^2 + (11+2)^2 - 1 = 289.
MATHEMATICA
Total/@(Select[Partition[Prime[Range[500]], 2, 1], #[[2]]-#[[1]]==2&]^2)-1 (* Harvey P. Dale, Feb 24 2011 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert Mohr, Feb 13 2011
STATUS
approved