cp's OEIS Frontend

A238048 Square array A(n,k), n>=1, k>=1, read by antidiagonals, where column k is the increasing list of all primes p such that (p+k)^2+k is also prime.

%I A238048 #25 Dec 15 2015 00:56:05
%S A238048 3,7,5,5,13,13,3,7,19,19,7,11,11,31,23,5,31,13,19,37,53,3,13,43,23,47,
%T A238048 43,73,7,5,19,67,29,59,79,83,11,13,11,29,73,31,61,97,89,3,23,43,19,59,
%U A238048 109,41,67,103,109,13,17,29,73,23,73,157,43,71,109,149
%N A238048 Square array A(n,k), n>=1, k>=1, read by antidiagonals, where column k is the increasing list of all primes p such that (p+k)^2+k is also prime.
%C A238048 Prime 2 is not contained in this array.
%H A238048 Alois P. Heinz, <a href="/A238048/b238048.txt">Antidiagonals n = 1..150, flattened</a>
%e A238048 Column k=3 contains prime 47 because (47+3)^2+3 = 2503 is prime.
%e A238048 Square array A(n,k) begins:
%e A238048    3,  7,  5,  3,   7,   5,  3,   7, ...
%e A238048    5, 13,  7, 11,  31,  13,  5,  13, ...
%e A238048   13, 19, 11, 13,  43,  19, 11,  43, ...
%e A238048   19, 31, 19, 23,  67,  29, 19,  73, ...
%e A238048   23, 37, 47, 29,  73,  59, 23,  79, ...
%e A238048   53, 43, 59, 31, 109,  73, 29, 103, ...
%e A238048   73, 79, 61, 41, 157,  83, 31, 109, ...
%e A238048   83, 97, 67, 43, 163, 103, 41, 127, ...
%p A238048 A:= proc(n, k) option remember; local p;
%p A238048       p:= `if`(n=1, 1, A(n-1, k));
%p A238048       do p:= nextprime(p);
%p A238048          if isprime((p+k)^2+k) then return p fi
%p A238048       od
%p A238048     end:
%p A238048 seq(seq(A(n, 1+d-n), n=1..d), d=1..11);
%t A238048 A[n_, k_] := A[n, k] = Module[{p}, For[p = If[n == 1, 1, A[n-1, k]] // NextPrime, True, p = NextPrime[p], If[PrimeQ[(p+k)^2+k], Return[p]]]]; Table[Table[A[n, 1+d-n], {n, 1, d}], {d, 1, 11}] // Flatten (* _Jean-François Alcover_, Jan 19 2015, after _Alois P. Heinz_ *)
%Y A238048 Column k=1 gives A157468.
%Y A238048 Cf. A238086.
%K A238048 nonn,tabl,look
%O A238048 1,1
%A A238048 _Alois P. Heinz_, Feb 17 2014