This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A247072 #73 Sep 24 2018 02:40:33 %S A247072 2,1093,11,1093,20771,66161,5,3,11,487,71,2693,863,29,29131,1093, %T A247072 46021,5,7,281 %N A247072 Smallest Wieferich prime (> sqrt(n)) in base n. %C A247072 a(n) = Smallest prime such that n appears in A143548. - _Eric Chen_, Nov 26 2014 %C A247072 The square of a(n) is the smallest squared prime that is a pseudoprime (> n) in base n; for example, a(2) = 1093, and 1093^2 = 1194649 is the smallest squared prime that is pseudoprime in base 2. - _Eric Chen_, Nov 26 2014 %C A247072 Is a(n) defined for all n? - _Eric Chen_, Nov 26 2014 %C A247072 Does every prime appear in this sequence? - _Eric Chen_, Nov 26 2014 %C A247072 a(22)..a(28) = {13, 13, 5, 20771, 71, 11, 19}, a(30)..a(46) = {7, 7, 1093, 233, 46145917691, 1613, 66161, 77867, 17, 8039, 11, 29, 23, 103, 229, 1283, 829}, a(48)..a(49) = {7, 491531}, a(51)..a(60) = {41, 461, 47, 19, 30109, 647, 47699, 131, 2777, 29}, a(62)..a(71) = {19, 23, 1093, 17, 89351671, 47, 19, 19, 13, 47}, a(74)..a(81) = {1251922253819, 17, 37, 32687, 43, 263, 13, 11}, a(83)..a(100) = {4871, 163, 11779, 68239, 1999, 2535619637, 13, 6590291053, 293, 727, 509, 11, 2137, 109, 2914393, 28627, 13, 487}; a(n) is currently unknown for n = {21, 29, 47, 50, 61, 72, 73, 82, 126, 132, 154, 186, 187, 188, 200, 203, 222, 231, 237, 301, 304, 309, 311, 327, 335, 347, 351, 355, 357, 435, 441, 454, 458, 496, 541, 542, 546, 554, 570, 593, 609, 610, 639, 640, 654, 662, 668, 674, 692, 697, 698, 701, 718, 724, 725, 727, 733, 743, 760, 772, 775, 777, 784, 798, 807, 808, 812, 829, 841, 858, 860, 871, 883, 912, 919, 944, 980, 983, 986, ...}. - _Eric Chen_, Nov 26 2014 %C A247072 a(21) > 3.4 * 10^13. - _Eric Chen_, Nov 26 2014 %H A247072 Eric Chen, <a href="/A247072/a247072_1.txt">Table of known a(n) up to a(1000)</a> %H A247072 Richard Fischer, <a href="http://www.fermatquotient.com/FermatQuotienten/"> Fermat quotients B^(P-1) == 1 (mod P^2)</a> %H A247072 Wikipedia, <a href="https://en.wikipedia.org/wiki/Wieferich_prime"> Wieferich prime</a> %e A247072 a(12) = 2693 because the Wieferich primes to base 12 are 2693, 123653, ..., and 2693 is greater than sqrt(12), so a(12) = 2693. %e A247072 a(17) = 46021 because the Wieferich primes to base 17 are 2, 3, 46021, 48947, 478225523351, ..., but neither 2 nor 3 is greater than sqrt(17), so a(17) = 46021. %t A247072 a247072[n_] := Block[{p = Int[Sqrt[n]]+1}, While[!PrimeQ[p] || [p < 10^8 && PowerMod[n, p - 1, p^2] != 1], p++]; If[p == 10^8, 0, p]]; Table[ a247072[n], {n, 100}] (* _Eric Chen_, Nov 27 2014 *) %o A247072 (PARI) a(n)=forprime(p=sqrtint(n)+1,,if(Mod(n^(p-1),p^2)==1,return(p))) %o A247072 n=1; while(n<101, print1(a(n), ", "); n++) \\ _Charles R Greathouse IV_, Nov 16 2014 %Y A247072 Cf. A039951, A096082, A174422, A178871. %Y A247072 Cf. A001220, A014127, A123692, A212583, A123693, A045616, A111027, A128667, A234810, A242741, A128668, A244260, A090968, A242982, A128669. %K A247072 nonn,more %O A247072 1,1 %A A247072 _Eric Chen_, Nov 16 2014