A084732 -id:A084732 - cp's OEIS frontend

A066180 a(n) = smallest base b so that repunit (b^prime(n) - 1) / (b - 1) is prime, where prime(n) = n-th prime; or 0 if no such base exists.

2, 2, 2, 2, 5, 2, 2, 2, 10, 6, 2, 61, 14, 15, 5, 24, 19, 2, 46, 3, 11, 22, 41, 2, 12, 22, 3, 2, 12, 86, 2, 7, 13, 11, 5, 29, 56, 30, 44, 60, 304, 5, 74, 118, 33, 156, 46, 183, 72, 606, 602, 223, 115, 37, 52, 104, 41, 6, 338, 217, 13, 136, 220, 162, 35, 10, 218, 19, 26, 39

Offset: 1

Frank Ellermann, Dec 15 2001

nonn

Is a(n) = 0 possible?
Let p be the n-th prime; Cp(x) be the p-th cyclotomic polynomial (x^p - 1)/(x - 1); a(n) is the least k > 1 such that Cp(k) is prime.
The values associated with a(5) and a(8) through a(70) have been certified prime with Primo. (a(1) through a(4), a(6) and a(7) give prime(2), prime(4), prime(11), prime(31), prime(1028) and prime(12251), respectively.)

			a(5) = 5 because 11 is the 5th prime; (b^5 - 1)/(b - 1) is composite for b = 2,3,4 and prime ((5^11 - 1)/4 = 12207031) for b = 5.
b = 61 for prime(12) = 37 because (61^37 - 1)/60 is prime and 61 is the least base b that makes (b^37 - 1)/(b - 1) a prime.

Paulo Ribenboim, "The New Book of Prime Numbers Records", Springer, 1996, p. 353.

Robert G. Wilson v, Table of n, a(n) for n = 1..300 (terms 1..200 from Charles R Greathouse IV).
H. Dubner, Generalized repunit primes, Math. Comp., 61 (1993), 927-930.
Andy Steward, Titanic Prime Generalized Repunits.
Eric Weisstein's World of Mathematics, Repunit.
H. C. Williams and E. Seah, Some primes of the form: (a^n - 1)/(a - 1), Mathematics of Computation 23, 1979.

Cf. A004023 (prime repunits in base 10), A000043 (prime repunits in base 2, Mersenne primes), A055129 (table of repunits).

Cf. A084732, A085398.

Table[p = Prime[n]; b = 1; While[b++; ! PrimeQ[(b^p - 1)/(b - 1)]]; b, {n, 1, 70}] (* Lei Zhou, Oct 07 2011 *)

/* This program assumes (probable) primes exist for each n. */
/* All 70 (probable) primes found by this program have been proved prime. */
gen_repunit(b,n) = (b^prime(n)-1)/(b-1);
for(n=1,70, b=1; until(isprime(p), b++; p=gen_repunit(b,n)); print1(b,","));

a(n) = A085398(prime(n)).

Sequence extended to 16 terms by Don Reble, Dec 18 2001
More terms from Rick L. Shepherd, Sep 14 2002
Entry revised by N. J. A. Sloane, Jul 23 2006

A128164 Least k > 2 such that (n^k - 1)/(n-1) is prime, or 0 if no such prime exists.

Original entry on oeis.org

3, 3, 0, 3, 3, 5, 3, 0, 19, 17, 3, 5, 3, 3, 0, 3, 25667, 19, 3, 3, 5, 5, 3, 0, 7, 3, 5, 5, 5, 7, 0, 3, 13, 313, 0, 13, 3, 349, 5, 3, 1319, 5, 5, 19, 7, 127, 19, 0, 3, 4229, 103, 11, 3, 17, 7, 3, 41, 3, 7, 7, 3, 5, 0, 19, 3, 19, 5, 3, 29, 3, 7, 5, 5, 3, 41, 3, 3, 5, 3, 0, 23, 5, 17, 5, 11, 7, 61, 3, 3

Offset: 2

Alexander Adamchuk, Feb 20 2007

nonn

a(n) = A084740(n) for all n except n = p-1, where p is an odd prime, for which A084740(n) = 2.
All nonzero terms are odd primes.
a(n) = 0 for n = {4,9,16,25,32,36,49,64,81,100,121,125,144,...}, which are the perfect powers with exceptions of the form n^(p^m) where p>2 and (n^(p^(m+1))-1)/(n^(p^m)-1) are prime and m>=1 (in which case a(n^(p^m))=p). - Max Alekseyev, Jan 24 2009
a(n) = 3 for n in A002384, i.e., for n such that n^2 + n + 1 is prime.
a(152) > 20000. - Eric Chen, Jun 01 2015
a(n) is the least number k such that (n^k - 1)/(n-1) is a Brazilian prime, or 0 if no such Brazilian prime exists. - Bernard Schott, Apr 23 2017
These corresponding Brazilian primes are in A285642. - Bernard Schott, Aug 10 2017
a(152) = 270217, see the top PRP link. - Eric Chen, Jun 04 2018
a(184) = 16703, a(200) = 17807, a(210) = 19819, a(306) = 26407, a(311) = 36497, a(326) = 26713, a(331) = 25033; a(185) > 66337, a(269) > 63659, a(281) > 63421, and there are 48 unknown a(n) for n <= 1024. - Eric Chen, Jun 04 2018
Six more terms found: a(522)=20183, a(570)=12907, a(684)=22573, a(731)=15427, a(820)=12043, a(996)=14629. - Michael Stocker, Apr 09 2020

			a(7) = 5 because (7^5 - 1)/6 = 2801 = 11111_7 is prime and (7^k - 1)/6 = 1, 8, 57, 400 for k = 1, 2, 3, 4. - _Bernard Schott_, Apr 23 2017

Max Alekseyev and Eric Chen, Table of n, a(n) for n = 2..184 (terms 2..151 from Max Alekseyev)
Eric Chen, Table of n, a(n) for n = 2..1024 status (updated by Jinyuan Wang)
H. Dubner, Generalized repunit primes, Math. Comp., 61 (1993), 927-930.
Richard Fischer, Generalized repunit primes of the form (B^N-1)/(B-1)
Top PRPs, Search by (152^n-1)/(152-1)
Top PRPs, Search by (b^n-1)/a
Eric Weisstein's World of Mathematics, Repunit

Cf. A084738, A065854, A084740, A084741, A065507, A084742, A066180, A084732, A285642, A085104.
Cf. A002384, A049409, A100330, A162862, A217070-A217089. (numbers b such that (b^p-1)/(b-1) is prime for prime p = 3 to 97)
Cf. A000043, A028491, A004061, A004062, A004063, A004023, A005808, A004064, A016054, A006032, A006033, A006034, A133857, A006035, A127995, A127996, A127997, A204940, A127998, A127999, A128000, A181979, A098438, A128002, A209120, A185073, A128003, A128004, A181987, A128005, A239637, A240765, A294722, A242797, A243279, A267375, A245237, A245442, A173767. (numbers n such that (b^n-1)/(b-1) is prime for b = 2 to 53)
A126589 gives locations of zeros.

Table[Function[m, If[m > 0, k = 3; While[! PrimeQ[(m^k - 1)/(m - 1)], k++]; k, 0]]@ If[Set[e, GCD @@ #[[All, -1]]] > 1, {#, IntegerExponent[n, #]} &@ Power[n, 1/e] /. {{k_, m_} /; Or[Not[PrimePowerQ@ m], Prime@ m, FactorInteger[m][[1, 1]] == 2] :> 0, {k_, m_} /; m > 1 :> n}, n] &@ FactorInteger@ n, {n, 2, 17}] (* Michael De Vlieger, Apr 24 2017 *)

a052409(n) = my(k=ispower(n)); if(k, k, n>1)
a052410(n) = if (ispower(n, , &r), r, n)
is(n) = issquare(n) || (ispower(n) && !ispseudoprime((n^a052410(a052409(n))-1)/(n-1)))
a(n) = if(is(n), 0, forprime(p=3, 2^16, if(ispseudoprime((n^p-1)/(n-1)), return(p)))) \\ Eric Chen, Jun 01 2015, corrected by Eric Chen, Jun 04 2018, after Charles R Greathouse IV in A052409 and Michel Marcus in A052410

a(18) = 25667 found by Henri Lifchitz, Sep 26 2007

A123487 Smallest prime q such that (q^p-1)/(q-1) is prime, where p = prime(n); or 0 if no such prime q exists.

Original entry on oeis.org

2, 2, 2, 2, 5, 2, 2, 2, 113, 151, 2, 61, 53, 89, 5, 307, 19, 2, 491, 3, 11, 271, 41, 2, 271, 359, 3, 2, 79, 233, 2, 7, 13, 11, 5, 29, 71, 139, 127, 139, 2003, 5, 743, 673, 593, 383, 653, 661, 251, 6389, 2833, 223, 163, 37, 709, 131, 41, 2203, 941, 2707, 13, 1283, 383

Offset: 1

Alexander Adamchuk, Sep 30 2006, Oct 02 2006

nonn

Corresponding primes (q^p-1)/(q-1) are listed in A123488.
a(n) coincides with A066180(n) when A066180(n) is prime or 0.
From Robert G. Wilson v, Dec 28 2016: (Start)
Conjecture: Never is a(n) equal to 0.
Records: 2, 5, 113, 151, 307, 491, 2003, 6389, 7883, 11813, 18587, 31721, 40763, ... ;
First occurrence of the k_th prime: 1, 20, 5, 32, 21, 33, 81, 17, ... ;
Positions where two occurs: 1, 2, 3, 4, 6, 7, 8, 11, 18, 24, 28, 31, 98, 111, ... ;
Positions where three occurs: 20, 27, 100, 182, ... ;
Positions where five occurs: 5, 15, 35, 42, 114, 158, ... ; etc. (End)
Jones & Zvonkin conjecture (as did Robert G. Wilson v above) that a(n) > 0 for all n. - Charles R Greathouse IV, Jul 23 2021

Robert G. Wilson v, Table of n, a(n) for n = 1..205
Gareth A. Jones and Alexander K. Zvonkin, Groups of prime degree and the Bateman-Horn Conjecture, arXiv:2106.00346 [math.GR], 2021.

Cf. A123488, A066180, A084732.

f[n_] := NestWhile[NextPrime, 2, ! PrimeQ[Cyclotomic[Prime[n], #]] &]; Array[f, 63](* Davin Park, Dec 28 2016 and Robert G. Wilson v, Dec 28 2016 *)

a(n) = {my(x = 2); while (!isprime(polcyclo(prime(n), x)), x= nextprime(x+1)); x;} \\ Michel Marcus, Dec 10 2016

A085399 Let Cn(x) be the n-th cyclotomic polynomial; a(n) is the first prime Cn(x) after Cn(1).

Original entry on oeis.org

2, 3, 7, 5, 31, 3, 127, 17, 73, 11, 12207031, 13, 8191, 43, 151, 257, 131071, 46441, 524287, 61681, 368089, 683, 11111111111111111111111, 241, 705429635566498619547944801, 2731, 262657, 15790321, 7369130657357778596659, 331, 2147483647, 65537, 599479, 43691

Offset: 1

Don Reble, Jun 28 2003

nonn

			a(11)=12207031 because C11(x) is composite for x=2,3,4 and C11(5)=12207031 is prime.

Jinyuan Wang, Table of n, a(n) for n = 1..396

Cf. A084732, A085398, A307687.

a(n) = Phi(n,A085398(n)) where Phi(n,k) is the n-th cyclotomic polynomial evaluated at k. - Jinyuan Wang, Sep 04 2022

More terms from Jinyuan Wang, Sep 04 2022

cp's OEIS Frontend

A066180 a(n) = smallest base b so that repunit (b^prime(n) - 1) / (b - 1) is prime, where prime(n) = n-th prime; or 0 if no such base exists.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A128164 Least k > 2 such that (n^k - 1)/(n-1) is prime, or 0 if no such prime exists.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A123487 Smallest prime q such that (q^p-1)/(q-1) is prime, where p = prime(n); or 0 if no such prime q exists.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A085399 Let Cn(x) be the n-th cyclotomic polynomial; a(n) is the first prime Cn(x) after Cn(1).

Views

Author

Keywords

Examples

Links

Crossrefs

Formula

Extensions

A123488 Smallest prime of the form (q^p-1)/(q-1), where p = prime(n) and q is also prime (q = A123487(n)).

Views

Author

Keywords

Comments

Crossrefs

Formula