cp's OEIS Frontend

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.

Showing 1-2 of 2 results.

A250206 Least base b > 1 such that b^A000010(n) = 1 (mod n^2).

Original entry on oeis.org

2, 5, 8, 7, 7, 17, 18, 15, 26, 7, 3, 17, 19, 19, 26, 31, 38, 53, 28, 7, 19, 3, 28, 17, 57, 19, 80, 19, 14, 107, 115, 63, 118, 65, 18, 53, 18, 69, 19, 7, 51, 19, 19, 3, 26, 63, 53, 17, 18, 57, 134, 19, 338, 161, 3, 31, 28, 41, 53, 107, 264, 115, 19, 127, 99, 161, 143, 65, 28, 99, 11, 55
Offset: 1

Views

Author

Eric Chen, Feb 21 2015

Keywords

Comments

a(n) = least base b > 1 such that n is a Wieferich number (see A077816).
At least, b = n^2+1 can satisfy this equation, so a(n) is defined for all n.
Least Wieferich number (>1) to base n: 2, 1093, 11, 1093, 2, 66161, 4, 3, 2, 3, 71, 2693, 2, 29, 4, 1093, 2, 5, 3, 281, 2, 13, 4, 5, 2, ...; each is a prime or 4. It is 4 if and only if n mod 72 is in the set {7, 15, 23, 31, 39, 47, 63}.
Does every natural number (>1) appear in this sequence? If yes, do they appear infinitely many times?
For prime n, a(n) = A185103(n), does there exist any composite n such that a(n) = A185103(n)?

Examples

			a(30) = 107 since A000010(30) = 8, 30^2 = 900, and 107 is the least base b > 1 such that b^8 = 1 (mod 900).
		

Crossrefs

Programs

  • Mathematica
    f[n_] := Block[{b = 2, m = EulerPhi[n]}, While[ PowerMod[b, m, n^2] != 1, b++]; b]; f[1] = 2; Array[f, 72] (* Robert G. Wilson v, Feb 28 2015 *)
  • PARI
    a(n)=for(k=2,2^24,if((k^eulerphi(n))%(n^2)==1, return(k)))

Formula

a(prime(n)) = A039678(n) = A185103(prime(n)).
a(A077816(n)) = 2.
a(A242958(n)) <= 3.

A257660 Numbers n such that 13^phi(n) == 1 (mod n^2), where phi(n) = A000010(n).

Original entry on oeis.org

2, 863, 1726, 3452, 371953, 743906, 1487812, 1747591, 1859765, 2975624, 3495182, 3719530, 5242773, 6990364, 7439060, 8737955, 10485546, 14878120, 15993979, 17475910, 20971092, 26213865, 29756240, 31987958, 34951820, 41942184, 47981937, 52427730, 59512480
Offset: 1

Views

Author

Felix Fröhlich, Jul 26 2015

Keywords

Comments

The subsequence of primes in this sequence is A128667.

Crossrefs

Programs

  • Mathematica
    Select[Range@ 1000000, Mod[13^EulerPhi[#], #^2] == 1 &] (* Michael De Vlieger, Jul 27 2015 *)
  • PARI
    for(n=2, 1e9, if(Mod(13, n^2)^(eulerphi(n))==1, print1(n, ", ")))
Showing 1-2 of 2 results.