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

A187023 a(n) is the smallest prime factor of n^n-1 having the form k*n+1.

Original entry on oeis.org

3, 13, 5, 11, 7, 29, 17, 19, 11, 15797, 13, 53, 7027567, 61, 17, 10949, 19, 109912203092239643840221, 41, 43, 23, 461, 73, 101, 937, 109, 29, 59, 31, 568972471024107865287021434301977158534824481, 257, 67, 103, 281, 37, 149, 191, 157, 41
Offset: 2

Views

Author

Michel Lagneau, Mar 02 2011

Keywords

Comments

The values of k are in A187025.

Examples

			7^7-1 = 2*3*29*4733; the smallest prime divisor of the form k*n+1 is 29 = 4*7+1, hence a(7) = 29.

Links

Crossrefs

Cf. A187022, A187025.

Programs

  • Magma
    A187023:=function(n); for d in PrimeDivisors(n^n-1) do if d mod n eq 1 then return d; end if; end for; return 0; end function; [ A187023(n): n in [2..50] ]; // Klaus Brockhaus, Mar 02 2011
  • Mathematica
    Table[p=First/@FactorInteger[n^n-1]; Select[p, Mod[#1, n] == 1 &, 1][[1]], {n, 2, 40}]

A187024 a(n) is the least number k such that k*n+1 is a prime dividing n^n+1.

Original entry on oeis.org

2, 2, 64, 104, 2, 16, 12, 8, 10, 2, 16, 1032, 2, 2, 17136, 2703399548648159360, 2, 5700, 7436, 16, 437174342164, 2, 1392, 9568, 2, 6, 16, 8, 2, 12, 20, 20, 176764, 2, 623673825204293256537467494040862720, 16, 5340, 2, 16, 440, 16, 22, 8, 990520
Offset: 2

Views

Author

Michel Lagneau, Mar 02 2011

Keywords

Comments

The smallest prime factor of n^n+1 of the form k*n+1 is A187022(n).

Examples

			12^12+1 = 89*193*233*2227777; the smallest prime divisor of the form k*n+1 is 193 = 16*12+1, hence a(12)=16.

Links

Crossrefs

Cf. A187022.

Programs

  • Magma
    A187024:=function(n); for d in PrimeDivisors(n^n+1) do if d mod n eq 1 then return (d-1)/n; end if; end for; return 0; end function; [ A187024(n): n in [2..40] ]; // Klaus Brockhaus, Mar 02 2011
  • Mathematica
    Table[p=First/@FactorInteger[n^n+1]; (Select[p, Mod[#1, n] == 1 &, 1][[1]]
  - 1)/n, {n, 2, 40}]

A216487 Smallest prime factor of n^(2n) - 1 having the form k*n+1.

Original entry on oeis.org

3, 7, 5, 11, 7, 29, 17, 19, 11, 23, 13, 53, 29, 31, 17, 10949, 19, 108301, 41, 43, 23, 47, 73, 101, 53, 109, 29, 59, 31, 373, 257, 67, 103, 71, 37, 149, 191, 79, 41, 83, 43, 173, 89, 181, 47, 659, 97, 197, 101, 103, 53, 107, 109, 881, 113, 229, 59, 709, 61, 977
Offset: 2

Views

Author

Michel Lagneau, Sep 11 2012

Keywords

Comments

The corresponding values of k are in A216506.

Examples

			a(7) = 29 because 7^14 - 1 = 2 ^ 4 * 3 * 29 * 113 * 911 * 4733 and the smallest prime divisor of the form k*n+1 is 29 = 4*7+1.

Links

Crossrefs

Cf. A006486, A007571, A187022, A187023, A216506.

Programs

  • Mathematica
    Table[p=First/@FactorInteger[n^(2*n)-1]; Select[p, Mod[#1, n] == 1 &, 1][[1]], {n, 2, 41}]
a[n_] := Module[{m = n + 1}, While[!PrimeQ[m] || PowerMod[n, 2*n, m] != 1, m += n]; m]; Array[a, 100, 2] (* Amiram Eldar, May 17 2024 *)
  • PARI
    a(n) = {my(m = n + 1); while(!isprime(m) || Mod(n, m)^(2*n) != 1, m += n); m;} \\ Amiram Eldar, May 17 2024

Formula

a(n) = Min{A187022(n), A187023(n)}.

Extensions

Data corrected by Amiram Eldar, May 17 2024

A216506 Least number k such that k*n+1 is a prime dividing n^(2n) - 1.

Original entry on oeis.org

1, 2, 1, 2, 1, 4, 2, 2, 1, 2, 1, 4, 2, 2, 1, 644, 1, 5700, 2, 2, 1, 2, 3, 4, 2, 4, 1, 2, 1, 12, 8, 2, 3, 2, 1, 4, 5, 2, 1, 2, 1, 4, 2, 4, 1, 14, 2, 4, 2, 2, 1, 2, 2, 16, 2, 4, 1, 12, 1, 16, 273, 2, 3, 2, 1, 4, 2, 2, 1, 246, 1, 4, 2, 2, 16, 8, 1, 4, 15, 2, 1, 2, 4, 12
Offset: 2

Views

Author

Michel Lagneau, Sep 11 2012

Keywords

Comments

The corresponding prime factors of n^(2n)-1 of the form k*n+1 is in A216487.

Examples

			a(7) = 4 because 7^14 - 1 = 2 ^ 4 * 3 * 29 * 113 * 911 * 4733 and the smallest prime divisor of the form k*n+1 is 29 = 4*7+1 => k = 4.

Links

Crossrefs

Cf. A006486, A007571, A187022, A187023, A216487.

Programs

  • Mathematica
    Table[p=First/@FactorInteger[n^(2*n)-1]; (Select[p, Mod[#1, n] == 1 &, 1][[1]]-1)/n, {n, 2, 50}]
a[n_] := Module[{m = n + 1}, While[!PrimeQ[m] || PowerMod[n, 2*n, m] != 1, m += n]; (m - 1)/n]; Array[a, 100, 2] (* Amiram Eldar, May 17 2024 *)
  • PARI
    a(n) = {my(m = n + 1); while(!isprime(m) || Mod(n, m)^(2*n) != 1, m += n); (m - 1)/n;} \\ Amiram Eldar, May 17 2024

Extensions

Data corrected and extended by Amiram Eldar, May 17 2024
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.