OFFSET
1,2
COMMENTS
Numbers n such that either A005722(n)-+1 is prime.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
EXAMPLE
1 is in the sequence because (1st prime-1)^2-1=0 is nonprime and (1st prime-1)^2+1=2 is prime;
3 is in the sequence because (3rd prime-1)^2-1=15 is nonprime and (3rd prime-1)^2+1=17 is prime.
MAPLE
A005722 := proc(n) (ithprime(n)-1)^2 ; end proc: for n from 1 to 800 do a := A005722(n) ; if isprime(a-1) <> isprime(a+1) then printf("%d, ", n) ; end if; end do: # R. J. Mathar, Apr 24 2010
MATHEMATICA
ppQ[n_]:=Module[{c=(Prime[n]-1)^2}, Sort[PrimeQ[{c+1, c-1}]]== {False, True}]; Select[Range[400], ppQ] (* Harvey P. Dale, Jun 24 2011 *)
Select[Range[400], Total[Boole[PrimeQ[(Prime[#]-1)^2+{1, -1}]]]==1&] (* Harvey P. Dale, Feb 01 2023 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Juri-Stepan Gerasimov, Feb 18 2010, Mar 27 2010
EXTENSIONS
More terms from R. J. Mathar, Apr 24 2010
Definition clarified by Harvey P. Dale, Jun 24 2011
STATUS
approved
