|
%I #16 Aug 22 2014 10:30:28
%S 2,7,11,17,19,23,31,43,47,53,67,71,79,97,103,107,127,151,163,179,191,
%T 199,211,223,241,251,271,283,331,337,367,383,431,449,463,467,487,499,
%U 523,547,577,587,607,631,647,683,691,719,751,787,823,827,881,907,967
%N Primes p such that (p + 1)/q is a perfect power for some prime q.
%H Jens Kruse Andersen, <a href="/A243009/b243009.txt">Table of n, a(n) for n = 1..10000</a>
%e 2 is in this sequence because (2 + 1)/3 = 1 = A001597(1), and 2 and 3 are both primes.
%o (PARI) isOK(p) = forprime(q=2, p, if(ispower((p+1)/q), return(1))); return(0)
%o concat(2, select(p->isOK(p), primes(1000))) \\ _Colin Barker_, Aug 18 2014
%Y Cf. A001597, A245303.
%K nonn
%O 1,1
%A _Juri-Stepan Gerasimov_, Aug 18 2014
%E Several terms corrected by _Colin Barker_, Aug 18 2014
|
