%I #21 Mar 14 2021 18:43:49

%S 9,19,41,85,155,235,341,457,691,929,1179,1555,1805,2065,2539,3181,

%T 3659,4149,4825,5255,5841,6637,7471,8723,9895,10505,11125,11771,12427,

%U 14465,16765,18079,19181,20851,22649,23859,25749,27385,29059,31141,32579,34753,37055,38215,39401,42189,47265,50845,52211

%N a(n) = (prime(1+n)*prime(n)) + prime(n) + 1.

%C 9 is the only perfect square in this sequence. - _Altug Alkan_, Jul 01 2017

%H Antti Karttunen, <a href="/A286624/b286624.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = (A000040(1+n)*A000040(n)) + A000040(n) + 1.

%F a(n) = 1 + A123134(n).

%F a(n) = A000040(n) + A023523(1+n).

%t Array[(#1 #2) + #1 + 1 & @@ Prime[# + {0, 1}] &, 49] (* _Michael De Vlieger_, Mar 14 2021 *)

%o (Scheme)

%o (define (A286624 n) (+ (* (A000040 (+ 1 n)) (A000040 n)) (A000040 n) 1))

%o (define (A286624 n) (+ 1 (* (+ 1 (A000040 (+ 1 n))) (A000040 n))))

%o (PARI) lista(nn) = forprime(p=2, nn, print1(p*(nextprime(p+1)+1)+1, ", ")); \\ _Altug Alkan_, Jul 01 2017

%Y Row 6 of A286625 (column 6 of A286623). Column 4 of A328464.

%Y One more than A123134.

%Y Cf. A000040, A023523, A180932 (primes in this sequence).

%K nonn

%O 1,1

%A _Antti Karttunen_, Jun 28 2017