login
Primes of the form p^k - p^(k-1) + 1 for some prime p and integer k > 1.
3

%I #8 Jul 16 2012 12:10:23

%S 3,5,7,17,19,43,101,157,163,257,487,1459,2029,4423,6163,14407,19183,

%T 22651,23549,26407,37057,39367,62501,65537,77659,113233,121453,143263,

%U 208393,292141,342733,375157,412807,527803,564899,590593,697049,843643

%N Primes of the form p^k - p^(k-1) + 1 for some prime p and integer k > 1.

%C It is usually the case that, for prime p and k > 1, the first time the totient function phi(n) has value p^k - p^(k-1) is for n = p^k. However, this is not true when p^k - p^(k-1) + 1 is prime.

%H Vincenzo Librandi, <a href="/A087126/b087126.txt">Table of n, a(n) for n = 1..1000</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/TotientFunction.html">Totient Function</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/TotientValenceFunction.html">Totient Valence Function</a>

%t lst={}; maxNum=10^6; n=1; While[p=Prime[n]; p^2-p+1<maxNum, i=2; While[m=p^i-p^(i-1)+1; m<maxNum, If[PrimeQ[m], AppendTo[lst, m]]; i++ ]; n++ ]; Sort[lst]

%Y Cf. A002383 (primes of the form n^2 + n + 1, which is the same as n^2 - n + 1).

%Y Cf. A019434 (Fermat primes), A003306 (2*3^n + 1 is prime), A056799 (8*9^n + 1 is prime), A056797 (9*10^n + 1 is prime), A087139 (least k such that p^k - p^(k-1) + 1 is prime for p = prime(n)).

%K nonn

%O 1,1

%A _T. D. Noe_, Aug 15 2003