# Greetings from The On-Line Encyclopedia of Integer Sequences! http://oeis.org/ Search: id:a045616 Showing 1-1 of 1 %I A045616 #54 Feb 29 2020 06:46:31 %S A045616 3,487,56598313 %N A045616 Primes p such that 10^(p-1) == 1 (mod p^2). %C A045616 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] %C A045616 No further terms below 1.172*10^14 (as of Feb 2020, cf. Fischer's table). %C A045616 56598313 was announced in the paper by Brillhart et al. - _Helmut Richter_, May 17 2004 %C A045616 A265012(A049084(a(n))) = 1. - _Reinhard Zumkeller_, Nov 30 2015 %D A045616 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. %D A045616 Richard K. Guy, Unsolved Problems in Number Theory, Springer, 2004, A3. %H A045616 Amir Akbary and Sahar Siavashi, The Largest Known Wieferich Numbers, INTEGERS, 18(2018), A3. See Table 1 p. 5. %H A045616 Richard Fischer, Fermat quotients B^(P-1) == 1 (mod P^2). %H A045616 Wilfrid Keller and Jörg Richstein, Solutions of the congruence a^(p-1) == 1 (mod p^r), Math. Comp. 74 (2005), 927-936. %H A045616 Peter L. Montgomery, New solutions of a^(p-1) == 1 (mod p^2), Math. Comp. 61 (1993), 361-363. %H A045616 Math Overflow, Is the smallest primitive root modulo p a primitive root modulo p^2?, Jun 09 2010. %H A045616 Helmut Richter, The period length of the decimal expansion of a fraction. %H A045616 Helmut Richter, The Prime Factors Of 10^486-1. %H A045616 Samuel Yates, The Mystique of Repunits, Math. Mag. 51 (1978), 22-28. %t A045616 A045616Q = PrimeQ@# && PowerMod[10, # - 1, #^2] == 1 &; Select[Range[1000000], A045616Q] (* _JungHwan Min_, Feb 04 2017 *) %t A045616 Select[Prime[Range[34*10^5]],PowerMod[10,#-1,#^2]==1&] (* _Harvey P. Dale_, Apr 10 2018 *) %o A045616 (PARI) lista(nn) = forprime(p=2, nn, if (Mod(10, p^2)^(p-1)==1, print1(p, ", "))); \\ _Michel Marcus_, Aug 16 2015 %o A045616 (Haskell) %o A045616 import Math.NumberTheory.Moduli (powerMod) %o A045616 a045616 n = a045616_list !! (n-1) %o A045616 a045616_list = filter %o A045616 (\p -> powerMod 10 (p - 1) (p ^ 2) == 1) a000040_list' %o A045616 -- _Reinhard Zumkeller_, Nov 30 2015 %Y A045616 Cf. A001220, A014127, A123692, A212583, A123693, A111027, A128667, A234810, A242741, A128668, A244260, A090968, A242982, A128669, A039951. %Y A045616 Cf. A265012, A049084, A000040. %K A045616 bref,hard,nonn,nice,more %O A045616 1,1 %A A045616 _Helmut Richter_, Dec 11 1999 # Content is available under The OEIS End-User License Agreement: http://oeis.org/LICENSE