OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = 1: p1 = 17 and p2 = 19 are the first such pair, with p1 + p2 = 36 = 6^2, (17 + 19)/36 = 1;
a(2) = 4: p1 = 71, p2 = 73; p1 + p2 = 144 = 12^2, (71 + 73)/36 = 4.
MAPLE
isa := n -> isprime(n) and isprime(n+2) and iperfpow(2*n+2) <> FAIL:
select(isa, [$4..1000000]): map(n -> (n+1)/18, %); # Peter Luschny, Jan 05 2020
PROG
(PARI) my(pp=5); forprime(p=7, 130000, if(p-pp==2, if(ispower(p+pp), print1((p+pp)/36, ", "))); pp=p)
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Jan 05 2020
STATUS
approved