|
%I #10 May 13 2013 01:54:10
%S 199720105201,18162880851661706161,108242541698965770601,
%T 170461841497718037541,314304417137765143321,522983530625286783781,
%U 8464133962630480939309,11418010921003118271841
%N Values of the Legendre polynomial P_6 that are the greater of a twin prime pair.
%H Charles R Greathouse IV, <a href="/A144173/b144173.txt">Table of n, a(n) for n = 1..10000</a>
%t TwinPrimeQ[n_]:=If[PrimeQ[n],If[PrimeQ[n-2]||PrimeQ[n+2],True,False],False](*TwinPrimeQ*) lst={};Do[p=LegendreP[6,n];If[TwinPrimeQ[p],AppendTo[lst,p]],{n,1,22^3,2}];lst
%o (PARI) forstep(n=1,1e4,2,if(isprime(k=pollegendre(6,n))&isprime(k-2),print1(k", ")))
%K nonn
%O 1,1
%A _Vladimir Joseph Stephan Orlovsky_, Sep 12 2008
%E New definition from _N. J. A. Sloane_ and _Charles R Greathouse IV_, Feb 03 2011 to Feb 05 2011.
|
