A245598 Smallest k > 1 such that prime(n)*k^prime(n)+1 is prime.

3, 4, 8, 4, 42, 60, 8, 198, 8, 54, 130, 778, 108, 178, 924, 44, 180, 706, 4, 170, 474, 30, 480, 1578, 214, 416, 34, 132, 2940, 60, 834, 666, 336, 168, 408, 216, 538, 114, 60, 266, 188, 58, 36, 1504, 4868, 2398, 430, 4, 1940, 408, 2036, 3038, 1146, 1902

Offset

1,1

Comments

From Robert G. Wilson v, Aug 05 2014: (Start)

For primes, p < 25000, for which p*k^p-1 is a prime:

k=1: only the first prime, 2;

k=2: none;

k=3: only the first prime, 2;

k=4: 3, 7, 67, 223, ..., ; indices: 2, 4, 19, 48, ..., ;

k=6: 2, 9901, 12043, ..., ; indices: 1, 1221, 1443, ..., ;

k=8: 5, 17, 23, ..., ; indices: 3, 7, 9, ..., ;

k=9: only the first prime, 2;

k=10: 3, 2161, ..., ; indices: 2, 326, ..., ; etc.

(End)

Links

Pierre CAMI, Table of n, a(n) for n = 1..300

Example

2*2^2+1=9 composite.
2*3^2+1=19 prime so a(1)=3.

Mathematica

a[n_] := Block[{k = 2, p = Prime@ n}, While[ !PrimeQ[p*k^p - 1], k += 2]; k];
Array[a, 60] (* Robert G. Wilson v, Aug 27 2014 *)

Prog

(PARI) a(n) = k=2; while(!isprime(prime(n)*k^prime(n)+1), k++); k
vector(40, n, a(n)) \\ Colin Barker, Jul 30 2014

Crossrefs

Cf. A245597.

Sequence in context: A198125 A127122 A086850 * A340533 A382557 A050274

Adjacent sequences: A245595 A245596 A245597 * A245599 A245600 A245601

Keyword

nonn

Author

Pierre CAMI, Jul 27 2014

Extensions

Definition corrected by Zak Seidov, Jul 27 2014

Status

approved