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.

Previous Showing 11-14 of 14 results.

A014959 Integers k such that k divides 22^k - 1.

Original entry on oeis.org

1, 3, 7, 9, 21, 27, 39, 49, 63, 81, 117, 147, 189, 243, 273, 343, 351, 441, 507, 567, 729, 819, 1029, 1053, 1143, 1323, 1521, 1701, 1911, 2187, 2401, 2457, 2943, 3081, 3087, 3159, 3429, 3549, 3969, 4401, 4563, 5103, 5733, 6561, 6591, 7203, 7371
Offset: 1

Views

Author

Keywords

Comments

Also, numbers n such that n divides s(n), where s(1)=1, s(k)=s(k-1)+k*22^(k-1) (cf. A014940).

Crossrefs

Integers n such that n divides b^n - 1: A067945 (b=3), A014945 (b=4), A067946 (b=5), A014946 (b=6), A067947 (b=7), A014949 (b=8), A068382 (b=9), A014950 (b=10), A068383 (b=11), A014951 (b=12), A116621 (b=13), A014956 (b=14), A177805 (b=15), A014957 (b=16), A177807 (b=17), A128358 (b=18), A125000 (b=19), A128360 (b=20), A014960 (b=24).

Programs

  • Mathematica
    nxt[{n_,s_}]:={n+1,s+(n+1)*22^n}; Transpose[Select[NestList[nxt,{1,1},7500], Divisible[ Last[#],First[#]]&]][[1]] (* Harvey P. Dale, Jan 27 2015 *)

Extensions

Edited by Max Alekseyev, Nov 16 2019

A333433 a(n) is the n-th number m that divides n^m - 1 (or 0 if m does not exist).

Original entry on oeis.org

1, 0, 4, 21, 8, 1555, 9, 6223, 40, 999, 20, 130801, 24, 4484077, 128, 117, 60, 118285781329, 42, 1432001198261, 104, 819, 72, 302508121, 81, 75625, 200, 61731, 78, 14507145975869, 72, 21958351241, 820, 12321, 289, 4375, 144
Offset: 1

Views

Author

Seiichi Manyama, Mar 21 2020

Keywords

Comments

From Jinyuan Wang, Mar 25 2020: (Start)
For n > 2, n < a(n) < q^(n-1), where q is a prime factor of n - 1.
If p is a prime, then a(p^e+1) is divisible by p. Proof: we can prove that p | m for m > 1 and n = p^e + 1. If n^m == 1 (mod m) and m > 1 is the minimum value that cannot be divisible by p, then gcd(m, eulerphi(m)) = 1. Thus, m must be of the form q*p_2*...*p_k, where q < p_2 < ... < p_k. Note that q | (n^m - 1) = (n^q - 1)*(Sum_{i=0..(m/q)-1} n)^(i*q)) and n^q - 1 can never be divisible by q. Therefore, Sum_{i=0..(m/q)-1} n^(i*q) == n^(m/q) - 1 == 0 (mod q). Because n^(q-1) == 1 (mod q) and gcd(m/q, q - 1) = 1, then n == 1 (mod q), a contradiction! (End)
a(38) <= 14948925257859919. - Giovanni Resta, Apr 15 2020

Crossrefs

Main diagonal of A333432.

Programs

  • PARI
    {a(n) = if(n==2, 0, my(cnt=0, k=0); while(cnt
    				

Formula

a(n) = A333432(n,n).

Extensions

a(30)-a(37) from Giovanni Resta, Apr 15 2020

A333506 Numbers k that divide 19^k-1.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 10, 12, 16, 18, 20, 24, 27, 30, 32, 36, 40, 42, 48, 50, 54, 60, 64, 72, 80, 81, 84, 90, 96, 100, 108, 110, 120, 126, 128, 136, 144, 150, 156, 160, 162, 168, 180, 192, 200, 210, 216, 220, 240, 243, 250, 252, 256, 270, 272, 288, 294, 300, 312, 320, 324, 330, 336
Offset: 1

Views

Author

Seiichi Manyama, Mar 25 2020

Keywords

Crossrefs

Columns k=19 of A333432.

Programs

  • Mathematica
    Join[{1}, Select[Range[2, 336], PowerMod[19, #, #] == 1 &]] (* Amiram Eldar, May 05 2021 *)
  • PARI
    for(k=1, 1e3, if(Mod(19, k)^k==1, print1(k", ")))

A014962 Odd numbers k that divide 25^k - 1.

Original entry on oeis.org

1, 3, 9, 21, 27, 63, 81, 93, 147, 171, 189, 243, 279, 441, 513, 567, 609, 651, 729, 837, 903, 1029, 1197, 1323, 1539, 1701, 1827, 1953, 2187, 2511, 2667, 2709, 2883, 2943, 3087, 3249, 3591, 3969, 4263, 4401, 4557, 4617, 5103, 5301, 5481, 5859, 6321
Offset: 1

Views

Author

Keywords

Comments

Also, numbers k such that k divides s(k), where s(1)=1, s(j) = s(j-1) + j*25^(j-1).
Equivalently, numbers k that divide ((24*k - 1)*25^k + 1) / 24^2 (cf. A014943).

Crossrefs

Programs

  • Maple
    select(t -> 25 &^ t - 1 mod t = 0, [seq(i,i=1..10^4,2)]); # Robert Israel, Oct 04 2020

Extensions

Edited by Max Alekseyev, Nov 16 2019
Previous Showing 11-14 of 14 results.