%I #25 Jun 21 2021 06:05:40
%S -1,4,-6,30,-18,42,-30,-22,128,-112,98,90,-78,-70,36,248,-232,158,150,
%T -280,182,-142,-130,420,210,-198,222,-210,-1074,1326,-238,560,-1092,
%U 1212,-592,36,350,-310,36,728,-1428,1548,-378,402,-1966,144,1832,462,-450,-442
%N Image of primes (A000040) under "little Hankel" transform that sends [c_0, c_1, ...] to [d_0, d_1, ...] where d_n = c_n^2 - c_{n+1}*c_{n-1}.
%C a(n) > 0 if and only if n+1 is in A046868. a(n) < 0 if and only if n+1 is in A233671. - _Chai Wah Wu_, Sep 10 2019
%H L. Panaitopol, <a href="http://jipam.vu.edu.au/v3n4/053_02.html">On the sequence p(n)^2=p(n-1)*p(n+1)</a>, J. Inequal. Pure Appl. Math., Volume 3, Issue 4, Article 53, 2002.
%F a(n) = determinant of matrix
%F | prime(n+1) prime(n)|
%F | prime(n+2) prime(n+1)|. - _Zak Seidov_, Jul 23 2008, indices corrected by _Gary Detlefs_, Dec 09 2011
%F a(n) = 2*A342567(n+1) for n >= 2. - _Hugo Pfoertner_, Jun 20 2021
%p A056221 := proc(n)
%p ithprime(n+1)^2-ithprime(n)*ithprime(n+2) ;
%p end proc:
%p seq(A056221(n),n=1..10) ; # _R. J. Mathar_, Dec 10 2011
%Y Cf. A046868, A233671, A342567.
%K sign
%O 1,2
%A _N. J. A. Sloane_, Aug 06 2000
|