%I #15 Oct 20 2022 12:43:47
%S 7,8,9,12,15,16,22,25,27,34,53,83,85,88,95,107,108,144,149,187,196,
%T 223,234,238,249,255,268,274,315,324,350,355,358,367,386,410,411,416,
%U 424,436,440,445,450,462,469,471,481,494,501,509,511,517,522,549,554,564
%N Numbers m such that prime(m^2+1)-prime(m^2)=2, where prime(j) is the j-th prime.
%C Square roots of squares in A029707. - _Michel Marcus_, Oct 20 2022
%e 25 is here because 626th and 625th primes are twin: 4639-4637=2.
%t t=Table[0, {250}]; t1=Table[0, {250}]; s=0; k=0; Do[s=Prime[1+n^2]-Prime[n^2]; If[s==2, k=k+1; t[[k]]=n; t1[[k]]=Prime[n^2]; Print[{k, n, Prime[n^2]}]], {n, 1, 2500}] t t1
%o (PARI) isok(m) = my(p=prime(m^2)); nextprime(p+1) - p == 2; \\ _Michel Marcus_, Oct 20 2022
%Y Cf. A029707, A001359, A006512.
%K nonn
%O 1,1
%A _Labos Elemer_, Aug 21 2002