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-1 of 1 results.

A202036 Smallest prime residue of x^n (mod n) for x=0..n-1, or 0 if no such prime exists.

Original entry on oeis.org

0, 0, 2, 0, 2, 3, 2, 0, 0, 5, 2, 0, 2, 2, 2, 0, 2, 0, 2, 5, 7, 3, 2, 0, 7, 3, 0, 0, 2, 19, 2, 0, 2, 2, 2, 0, 2, 5, 5, 0, 2, 7, 2, 5, 17, 2, 2, 0, 19, 0, 2, 13, 2, 0, 11, 0, 7, 5, 2, 0, 2, 2, 0, 0, 2, 3, 2, 13, 2, 11, 2, 0, 2, 3, 7, 5, 2, 13, 2, 0, 0, 2, 2, 0
Offset: 1

Views

Author

Michel Lagneau, Dec 09 2011

Keywords

Examples

			a(7) = 2 because k^7 == 0, 1, 2, 3, 4, 5, 6 (mod 7) => 2 is the smallest prime.
		

Crossrefs

Programs

  • Maple
    for n from 1 to 100 do: W:={}:for k from 0 to n-1 do:z:= irem(k^n,n): if type(z,prime)=true then W:=W union {z}:else fi:od: x:=nops(W): if x<>0 then printf(`%d, `,W[1]): else printf(`%d, `,0):fi: od:
  • Mathematica
    Table[SelectFirst[Sort[PowerMod[Range[n-1],n,n]],PrimeQ],{n,90}]/.Missing["NotFound"]->0 (* Harvey P. Dale, May 01 2023 *)
  • PARI
    A202036(n) = { my(z,y=n); for(x=1,n-1,z = lift(Mod(x,n)^n); if(isprime(z), y = min(z,y))); if(y==n,0,y); }; \\ - Antti Karttunen, May 19 2021
Showing 1-1 of 1 results.