%I #10 Feb 29 2020 13:02:28
%S 81,46225,151321,588289,769129,1108809,1692601,2283121,3367225,
%T 6195121,9715689,12341169,14814801,18344089,18983449,20657025,
%U 35271721,38551681,38775529,48511225,50993881,54597321,60730849,63059481,80156209
%N Squares that are a sum of twin primes and their indices.
%C Add both primes of the pair as well as their prime indices. Example: prime(11) - prime(10) = 31 - 29 = 2; prime(10) + prime(11) + 10 + 11 = 81 = 9^2 = a(1).
%t Select[#[[1]]+#[[2]]+PrimePi[#[[1]]]+PrimePi[#[[2]]]&/@Select[ Partition[ Prime[Range[3*10^6]],2,1],#[[2]]-#[[1]]==2&],IntegerQ[Sqrt[#]]&] (* _Harvey P. Dale_, Feb 29 2020 *)
%Y Squares in A088186. Cf. A088187.
%K easy,nonn
%O 1,1
%A _Zak Seidov_, Sep 22 2003
%E More terms from _Ray Chandler_, Sep 25 2003