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

A128360 Numbers k such that k divides 20^k - 1.

Original entry on oeis.org

1, 19, 361, 6859, 130321, 2476099, 47045881, 148305659, 893871739, 2817807521, 4234136149, 10350100679, 16983563041, 53538342899, 80448586831, 196651912901, 322687697779, 815211156289, 1017228515081, 1432001198261, 1528523149789
Offset: 1

Views

Author

Alexander Adamchuk, Mar 02 2007

Keywords

Comments

19 divides a(n) for n > 1. All powers of 19 are terms. a(n) = 19^(n-1) for all to n < 8, while a(8) = A128356(8) = 148305659 = 410819*19^2.
Prime divisors of a(n) in the order of appearance are {19, 410819, 617311, 1508981, ...}. - Alexander Adamchuk, May 16 2010

Crossrefs

Programs

Extensions

a(9)-a(11) from Stefan Steinerberger, May 09 2007
a(12)-a(15) from Alexander Adamchuk, May 16 2010
Edited and a(16)-a(21) added by Max Alekseyev, Oct 02 2010

A128356 Least number k > 1 (that is not the power of prime p) such that k divides (p+1)^k-1, where p = prime(n).

Original entry on oeis.org

20, 21, 1555, 889, 253, 2041, 5846759, 148305659, 1081, 279241, 9641, 950123, 33661, 63213709997, 583223, 3775349, 72707647, 149070763, 196932497, 5091481, 25760459, 14307947980741, 13861, 9362711, 376457, 132766545553, 63757
Offset: 1

Views

Author

Alexander Adamchuk, Mar 02 2007

Keywords

Comments

All listed terms have 2 distinct prime divisors. Most listed terms are semiprimes, except a(7) = 20231*17^2 and a(8) = 410819*19^2. p = prime(n) divides a(n). Quotients a(n)/prime(n) are listed in A128357 = {10, 7, 311, 127, 23, 157, 343927, ...}. a(15) = 583223 = 47*12409. a(16) = 3775349 = 53*71233.

Crossrefs

Programs

  • Mathematica
    (* This program is not suitable to compute a large number of terms *) a[n_] := For[p = Prime[n]; k = 2, True, k++, If[Length[FactorInteger[k]] == 2, If[Mod[PowerMod[p + 1, k, k] - 1, k] == 0, Print[k]; Return[k]]]]; Table[a[n], {n, 1, 13}] (* Jean-François Alcover, Oct 07 2013 *)

Extensions

Terms a(14) onwards from Max Alekseyev, Feb 08 2010

A128357 Quotients A128356(n)/prime(n).

Original entry on oeis.org

10, 7, 311, 127, 23, 157, 343927, 7805561, 47, 9629, 311, 25679, 821, 1470086279, 12409, 71233, 1232333, 2443783, 2939291, 71711, 352883, 181113265579, 167, 105199, 3881, 1314520253, 619, 20759, 117503, 1162660843, 1880415721, 263
Offset: 1

Views

Author

Alexander Adamchuk, Mar 02 2007, Mar 09 2007

Keywords

Comments

A128356 = {20, 21, 1555, 889, 253, 2041, 5846759, ...} = Least number k>1 (that is not the power of prime p) such that k divides (p+1)^k-1, where p = prime(n). Most listed terms are primes, except a(7) = 20231*17 and a(8) = 410819*19. a(15) = 12409. a(16) = 71233.
Note that all prime listed terms of {a(n)} coincide with terms of A128456 = {2, 7, 311, 127, 23, 157, 7563707819165039903, 75368484119, 47, 9629, 311, 25679, 821, ...} = least prime factor of ((p+1)^p - 1)/p^2, where p = prime(n).

Crossrefs

Cf. A128356 (least number k > 1 (that is not a power of prime p) such that k divides (p+1)^k-1, where p = prime(n)).
Cf. A128456 (least prime factor of ((p+1)^p - 1)/p^2, where p = prime(n)).

Extensions

Terms a(14) onwards from Max Alekseyev, Feb 08 2010

A177805 Numbers k such that k divides 15^k - 1.

Original entry on oeis.org

1, 2, 4, 7, 8, 14, 16, 28, 32, 49, 56, 64, 98, 112, 128, 136, 196, 224, 256, 272, 343, 392, 448, 452, 512, 544, 686, 784, 812, 896, 904, 952, 1024, 1088, 1372, 1568, 1624, 1792, 1808, 1904, 2048, 2176, 2312, 2401, 2744, 3136, 3164, 3248, 3584, 3616, 3808, 4096
Offset: 1

Views

Author

Alexander Adamchuk, May 17 2010

Keywords

Comments

A000420 are the only odd terms of the sequence. - Robert Israel, Feb 25 2020

Crossrefs

Programs

A177807 Numbers k that divide 17^k - 1.

Original entry on oeis.org

1, 2, 4, 6, 8, 12, 16, 18, 20, 24, 32, 36, 40, 42, 48, 54, 60, 64, 72, 78, 80, 84, 96, 100, 108, 116, 120, 126, 128, 144, 156, 160, 162, 168, 180, 192, 200, 216, 220, 232, 234, 240, 252, 256, 288, 294, 300, 312, 320, 324, 336, 342, 348, 360, 378, 384, 400, 420
Offset: 1

Views

Author

Alexander Adamchuk, May 17 2010

Keywords

Crossrefs

Programs

  • Mathematica
    {1}~Join~Select[Range[420], PowerMod[17, #, #] == 1 &] (* Giovanni Resta, Jan 30 2020 *)

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

A087807 Prime factors of solutions to 24^n == 1 (mod n).

Original entry on oeis.org

23, 47, 14759, 49727, 124799, 304751, 497261, 609503, 1828507, 2685259, 10741037, 12872687, 13877879, 23462213, 23652649, 27755759, 29134267, 31908959, 53753807, 65205263, 132771091, 218148653, 341965703, 551361983, 734951759
Offset: 1

Views

Author

Thomas Baruchel, Oct 14 2003

Keywords

Comments

Primes that divide at least one term of A014960.
Prime p is in this sequence iff the multiplicative order of 24 modulo p is the product of smaller terms of this sequence. - Max Alekseyev, May 26 2010

Examples

			A014960(12) = 2870377 = 23 * 124799
		

Crossrefs

Extensions

Corrected and extended by Max Alekseyev, May 26 2010
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

A014942 a(n) = (1 + 24^n*(23*n - 1))/529.

Original entry on oeis.org

1, 49, 1777, 57073, 1715953, 49491697, 1387212529, 38078983921, 1028756811505, 27446832213745, 724884022832881, 18985057740861169, 493749574409596657, 12764586312924606193, 328300388160453422833, 8406016915457191128817, 214387788361524002631409, 5448748098049810036109041
Offset: 1

Views

Author

Keywords

Crossrefs

Cf. A014960.

Programs

  • Mathematica
    a[1] = 1; a[n_Integer] := a[n - 1] + n*24^(n - 1); Table[ a[n], {n, 1, 20}]
  • PARI
    Vec(-x/((x-1)*(24*x-1)^2) + O(x^19)) \\ Elmo R. Oliveira, May 18 2025

Formula

a(n) = a(n-1) + n*24^(n-1).
a(n)= 49*a(n-1) - 624*a(n-2) + 576*a(n-3) for n > 3.
G.f.: x/((1-x)*(1-24*x)^2). - Joerg Arndt, May 23 2010
From Elmo R. Oliveira, May 18 2025: (Start)
E.g.f.: exp(x)*(1 + exp(23*x)*(552*x - 1))/529.
a(n) = 48*a(n-1) - 576*a(n-2) + 1 for n > 2. (End)

Extensions

More terms from Robert G. Wilson v, Sep 15 2000
Exchanged Lava formula of Jan 2009 with the definition [R. J. Mathar, May 23 2010]

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