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.

A284044 Largest positive k among all primes p < n such that n^(p-1) == 1 (mod p^k).

Original entry on oeis.org

1, 1, 2, 1, 2, 2, 3, 2, 1, 1, 2, 1, 1, 1, 4, 3, 3, 1, 2, 2, 2, 2, 3, 3, 2, 3, 2, 2, 2, 2, 5, 1, 2, 1, 2, 2, 1, 2, 3, 3, 2, 2, 2, 2, 1, 2, 4, 2, 2, 1, 3, 2, 3, 1, 3, 1, 1, 2, 2, 2, 2, 2, 6, 1, 2, 3, 2, 2, 2, 1, 3, 2, 2, 2, 2, 2, 2, 4, 4, 4, 1, 1, 2, 1, 1, 1, 3
Offset: 3

Views

Author

Felix Fröhlich, Apr 02 2017

Keywords

Comments

a(n) > 1 iff A255920(n) > 0, i.e., iff n is a term of A273786.

Examples

			For n = 7: the maximal exponents k in the congruence 7^(p-1) == 1 (mod p^k) for p = 2, 3, 5 are 1, 1, 2, respectively. Since 2 is the largest exponent among that list, a(7) = 2.
		

Crossrefs

Programs

  • PARI
    a(n) = my(r=1); forprime(p=1, n-1, my(k=1); while(1, if(Mod(n, p^k)^(p-1)!=1, k--; break, k++)); if(k > r, r=k)); r