Helmut Richter - cp's OEIS frontend

User: Helmut Richter

Helmut Richter has authored 2 sequences.

A092131 Distance from 2^n to the next prime.

0, 1, 3, 1, 5, 3, 3, 1, 9, 7, 5, 3, 17, 27, 3, 1, 29, 3, 21, 7, 17, 15, 9, 43, 35, 15, 29, 3, 11, 3, 11, 15, 17, 25, 53, 31, 9, 7, 23, 15, 27, 15, 29, 7, 59, 15, 5, 21, 69, 55, 21, 21, 5, 159, 3, 81, 9, 69, 131, 33, 15, 135, 29, 13, 131, 9, 3, 33, 29, 25, 11, 15, 29, 37, 33, 15, 11, 7, 23

Offset: 1

Helmut Richter (richter(AT)lrz.de), Mar 30 2004

Essentially the same as A013597. - T. D. Noe, Jul 17 2007
From Jianing Song, May 28 2024: (Start)
Not every odd number is present, as no term can be equal to a Sierpiński number (for example 78557); cf. A076336. See also A067760.
Conjecture: Every odd number which is not a Sierpiński number is a term. In other words, for every odd k which is not a Sierpiński number, there exists some n >= 1 such that 2^n + 1, 2^n + 3, ..., 2^n + (k-2) are all composite while 2^n + k is prime. (End)

			a(13)=17 because 2^13=8192 and the next prime is 8209=8192+17.

T. D. Noe, Table of n, a(n) for n=1..5000

Cf. A013597.
Equivalent sequence for previous prime: A013603.
Cf. A000079, A007920, A067760, A076336.

Join[{0},NextPrime[#]-#&/@(2^Range[2,80])] (* Harvey P. Dale, Jun 06 2012 *)

for(i=1,100,x=2^i;print1(nextprime(x)-x,","))

a(n) = nextprime(2^n) - 2^n.

a(n) = A007920(A000079(n)). - Michel Marcus, Oct 19 2022

A045616 Primes p such that 10^(p-1) == 1 (mod p^2).

Original entry on oeis.org

3, 487, 56598313

Offset: 1

Helmut Richter, Dec 11 1999

Primes p such that the decimal fraction 1/p has same period length as 1/p^2, i.e., the multiplicative order of 10 modulo p is the same as the multiplicative order of 10 modulo p^2. [extended by Felix Fröhlich, Feb 05 2017]
No further terms below 1.172*10^14 (as of Feb 2020, cf. Fischer's table).
56598313 was announced in the paper by Brillhart et al. - Helmut Richter, May 17 2004
A265012(A049084(a(n))) = 1. - Reinhard Zumkeller, Nov 30 2015
From Jianing Song, Jun 21 2025: (Start)
The ring of integers of Q(10^(1/k)) is Z[10^(1/k)] if and only if k does not have a prime factor in this sequence. See Theorem 5.3 of the paper of Keith Conrad. For example, we have:
  (1 + 10^(1/3) + 10^(2/3))/3 is an algebraic integer, but it is not in Z[10^(1/3)];
  (1 + 10^(486/487) + 10^(2*486/487) + ... + 10^(486*486/487))/487 is an algebraic integer, but it is not in Z[10^(1/487)];
  (1 + 10^(56598312/56598313) + 10^(2*56598312/56598313) + ... + 10^(56598312*56598312/56598313))/56598313 is an algebraic integer, but it is not in Z[10^(1/56598313)]. (End)

J. Brillhart, J. Tonascia, and P. Weinberger, On the Fermat quotient, pp. 213-222 of A. O. L. Atkin and B. J. Birch, editors, Computers in Number Theory. Academic Press, NY, 1971.
Richard K. Guy, Unsolved Problems in Number Theory, Springer, 2004, A3.

Amir Akbary and Sahar Siavashi, The Largest Known Wieferich Numbers, INTEGERS, 18(2018), A3. See Table 1 p. 5.
Keith Conrad, The ring of integers in a radical extension.
Richard Fischer, Fermat quotients B^(P-1) == 1 (mod P^2).
Wilfrid Keller and Jörg Richstein, Solutions of the congruence a^(p-1) == 1 (mod p^r), Math. Comp. 74 (2005), 927-936.
Peter L. Montgomery, New solutions of a^(p-1) == 1 (mod p^2), Math. Comp. 61 (1993), 361-363.
Math Overflow, Is the smallest primitive root modulo p a primitive root modulo p^2?, Jun 09 2010.
Helmut Richter, The period length of the decimal expansion of a fraction.
Helmut Richter, The Prime Factors Of 10^486-1.
Siqiong Yao and Akira Toyohara, The length of the repeating decimal, arXiv:2507.01295 [math.NT], 2025. See p. 19.
Samuel Yates, The Mystique of Repunits, Math. Mag. 51 (1978), 22-28.

Cf. A001220, A014127, A123692, A212583, A123693, A111027, A128667, A234810, A242741, A128668, A244260, A090968, A242982, A128669, A039951.

Cf. A265012, A049084, A000040.

import Math.NumberTheory.Moduli (powerMod)
a045616 n = a045616_list !! (n-1)
a045616_list = filter
               (\p -> powerMod 10 (p - 1) (p ^ 2) == 1) a000040_list'
-- Reinhard Zumkeller, Nov 30 2015

A045616Q = PrimeQ@# && PowerMod[10, # - 1, #^2] == 1 &; Select[Range[1000000], A045616Q] (* JungHwan Min, Feb 04 2017 *)
Select[Prime[Range[34*10^5]],PowerMod[10,#-1,#^2]==1&] (* Harvey P. Dale, Apr 10 2018 *)

lista(nn) = forprime(p=2, nn, if (Mod(10, p^2)^(p-1)==1, print1(p, ", "))); \\ Michel Marcus, Aug 16 2015

cp's OEIS Frontend

User: Helmut Richter

A092131 Distance from 2^n to the next prime.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula