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

A114207 Smallest solution to 10^m == 1 (mod m) having the prime divisor A066364(n).

Original entry on oeis.org

3, 111, 13203, 20439, 1997001, 22494039, 116226009, 761157, 278522253, 206613747, 17677747557, 835525881, 12933400720959, 228717562653, 5465090439, 13095850041, 431138536893, 4734551277, 58199580096201, 59875330325409, 228520359, 3003003, 257494085001, 1029221499627, 136635497220969
Offset: 1

Views

Author

Max Alekseyev, Nov 17 2005

Keywords

Examples

			a(6)=m(5477)=22494039 since it is the smallest m such that 10^m == 1 (mod m) and 5477|m.
		

Crossrefs

Cf. A066364.

Programs

  • PARI
    { m(p) = my(f,l,q); f=factorint(p)[,1]; l=p; for(i=1,length(f),q=znorder(Mod(10,f[i])); l=lcm(l,q); l=lcm(l,m(q)) ); l }

Formula

a(n)=m(p), where p=A066364(n) and m(p)=lcm(p, ord_p(10), m(q)) with q going over all prime divisors of ord_p(10).

A014950 Numbers m such that m divides 10^m - 1.

Original entry on oeis.org

1, 3, 9, 27, 81, 111, 243, 333, 729, 999, 2187, 2997, 4107, 6561, 8991, 12321, 13203, 19683, 20439, 26973, 36963, 39609, 59049, 61317, 80919, 110889, 118827, 151959, 177147, 183951, 242757, 332667, 356481, 455877, 488511, 531441, 551853, 728271
Offset: 1

Views

Author

Keywords

Comments

Also, m such that m | R(m) = A002275(m). - Lekraj Beedassy, Mar 25 2005
For n > 1, 3 divides a(n). If m is in the sequence and d divides m then for each positive integer k, d^k*m is in the sequence. So if m is in the sequence then m^k is in the sequence for each positive integer k. In particular, 3^k is in this sequence for all k. - Farideh Firoozbakht, Apr 14 2010
Numbers m such that m divides s(m), where s(1) = 1, s(k) = s(k-1) + k*10^(k-1).
Number of terms <= 10^k, beginning with k = 0: 1, 3, 5, 10, 15, 25, 41, 68, 108, 178, 291, ... - Robert G. Wilson v, Nov 30 2013
Numbers m such that m divides A033713(m). - Hans Havermann, Jan 25 2014

References

  • J. D. E. Konhauser et al., Which Way Did The Bicycle Go? Problem 80 pp. 26; 133, Dolciani Math. Exp., No. 18, MAA, Washington DC, 1996.

Crossrefs

Programs

Formula

Solutions to 10^m == 1 (mod m). - Vladeta Jovovic

Extensions

More terms from Vladeta Jovovic, Dec 18 2001
More terms from Larry Reeves (larryr(AT)acm.org), Jan 06 2005
Edited by Max Alekseyev, May 20 2011

A057719 Prime factors of numbers in A006521 (numbers k that divide 2^k + 1).

Original entry on oeis.org

3, 19, 163, 571, 1459, 8803, 9137, 17497, 41113, 52489, 78787, 87211, 135433, 139483, 144667, 164617, 174763, 196579, 274081, 370009, 370387, 478243, 760267, 941489, 944803, 1041619, 1220347, 1236787, 1319323, 1465129, 1663579, 1994659
Offset: 1

Views

Author

Keywords

Comments

A prime p is in this sequence iff all prime divisors of ord_p(2)/2 are in this sequence, where ord_p(2) is the order of 2 modulo p. - Max Alekseyev, Jul 30 2006

Examples

			2^171 + 1 == 0 (mod 171), 171 = 3^2*19, 2^13203+1 == 0 (mod 13203), 13203 = 3^4*163.
		

Crossrefs

Programs

  • Mathematica
    S = {2}; Reap[For[p = 3, p < 2 10^6, p = NextPrime[p], f = FactorInteger[ MultiplicativeOrder[2, p]]; If[f[[1, 1]] != 2 || f[[1, 2]] != 1, Continue[]]; f = f[[All, 1]]; If[Length[Intersection[S, f]] == Length[f], S = Union[S, {p}]; Print[p]; Sow[p]]]][[2, 1]] (* Jean-François Alcover, Nov 11 2018, from PARI *)
  • PARI
    { A057719() = local(S,f); S=Set([2]); forprime(p=3,10^7, f=factorint(znorder(Mod(2,p))); if(f[1,1]!=2||f[1,2]!=1,next); f=f[,1]; if(length(setintersect(S,Set(f)))==length(f), S=setunion(S,[p]); print1(p,", "))) }

Extensions

Edited by Max Alekseyev, Jul 30 2006

A171980 Prime divisors of elements of A129066.

Original entry on oeis.org

5, 3001, 120041, 532501, 720241, 2160721, 3937501, 9375001, 16505501, 120040001, 158453021, 165055001, 202567501, 289312501, 562500061, 900307501, 985937501, 1500512501, 1512504701, 3169060421, 3301100021, 3908604433, 3993757501
Offset: 1

Views

Author

Max Alekseyev, Jan 20 2010

Keywords

Comments

Corresponding smallest multiples from A129066 are given in A171981.
Prime p>5 is in this sequence if the multiplicative order of (sqrt(5)-3)/2 modulo p is the product of smaller terms of this sequence.

Crossrefs

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

A342810 Numbers k that divide the smallest number whose sum of digits is k.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 21, 27, 81, 191, 243, 729, 999, 2187, 2997, 6561, 8991, 19683, 26973, 33321, 36963, 39049, 59049, 80919, 100389, 110889, 118827, 177147, 177897, 183951, 242757, 332667, 356481, 531441, 551853, 728271, 998001, 1069443, 1367631, 1594323, 1655559, 2184813
Offset: 1

Views

Author

Ruediger Jehn, Mar 22 2021

Keywords

Comments

By definition, if k divides A051885(k), then k is a term of this sequence.
From Ruediger Jehn, Jun 17 2021: (Start)
None of the terms is divisible by 2*5*11*13.
If a term x has the form 3^m * y where m > 1 (which is the case for the overwhelming number of terms of this sequence), then all prime factors of y are terms of A066364.
If a term x has the form 3^m * p * q where m > 1, where p is a term of A066364 and where q is the product of all other factors of the prime factorization of x, then all numbers 3^m * p^i * q are also terms for any integer i. (End)

Examples

			21 is a term because the smallest number with a digital sum of 21 is 399 (A051885(21) = 399) which is divisible by 21.
		

Crossrefs

Programs

  • PARI
    MAX=10000; for (e = 0, MAX, for (d = 1, 9, k =(d+1)*10^e - 1; x = d+9*e; if (k%x==0, print1(x, ", ");)))
    
  • Python
    A342810_list = [n for n in range(1,10**6) if n==1 or ((n % 9)+1)*pow(10,n//9,n) % n == 1] # Chai Wah Wu, Apr 04 2021

Extensions

Name clarified by Jon E. Schoenfield, Apr 27 2021

A122787 a(n) is the smallest prime p such that the multiplicative order of 10 modulo p is 3^n.

Original entry on oeis.org

3, 37, 333667, 757, 163, 411361786890737698932559, 313471, 2558791, 618846643, 2238862519, 396319276163359, 34720813
Offset: 0

Views

Author

Farideh Firoozbakht, Oct 06 2006

Keywords

Comments

For n>0, a(n) is the smallest prime p>3 such that 3^n*p but not 3^(n-1)*p is a solution to 10^x==1 (mod x). A014950 gives solutions of this equation. It's obvious that if n is a term of A014950 then 3n is also a term of A014950. So according to the definition of a(n), for each m>n-1, 3^m*a(n) is in the sequence A014950.
a(n) is the smallest prime divisor of \Phi_{3^n}(10)/3, where \Phi_k(x) is k-th cyclotomic polynomial. a(n) is congruent to 1 modulo 3^n and 1, 3, 9, 13, 27, 31, 37, or 39 modulo 40. - Max Alekseyev, Nov 18 2014
a(12)>10^17, a(13)=796884087799, a(14)=86093443, a(15)=70367039929, a(16)>8*10^18, a(17)=662489036191, a(18)>10^19, a(19)>10^19, a(20)=38180289190951, a(21)=28305715767319, a(22)>10^20, a(23)>10^20, a(24)=63829075244707. - Ray Chandler, Dec 25 2013

Examples

			p=333667 is the smallest prime such that multiplicative order of 10 modulo p is 3^2, so a(2)=333667.
		

Crossrefs

Programs

  • PARI
    a(n) = factor(polcyclo(3^n,10)/3)[1,1] \\ Max Alekseyev, Nov 18 2014

Extensions

Revised and extended by Max Alekseyev, Apr 25 2009
a(13), a(15), a(17), a(20), a(21), a(24) from Ray Chandler, Dec 25 2013

A164816 Prime factors in a divisibility sequence of the Lucas sequence v(P=3,Q=5) of the second kind.

Original entry on oeis.org

3, 17, 103, 163, 373, 487, 1733, 3469, 4373, 8803, 10259, 15607, 16069, 26237, 26297, 31193, 31517, 35153, 37987, 38047, 38149, 39367, 52817, 60427, 60589, 61553, 74357, 76837, 78713, 100733, 103979, 114377, 119891, 152189, 181277, 231131, 235891, 238307, 239783, 280927, 289243, 316903, 338581
Offset: 1

Views

Author

Jonathan Vos Post, Aug 26 2009

Keywords

Comments

This is the last sequence on p. 15 of Smyth. [WARNING: Smyth lists 2 as a possible prime factor, which, in fact, is not possible. - Max Alekseyev, Sep 17 2024]
The Lucas sequence with P = 3, Q = 5 is defined as v=2,3,-1,-18,-49,-57,.. where v(n) = P*v(n-1)-Q*v(n-2), with g.f. (2-3x)/(1-3x+5x^2).
The indices n such that n|v(n) define the sequence T = 1,3,9,27,81,153,243,459,... as listed by Smyth.
The OEIS sequence shows all distinct prime factors of elements of T.

Crossrefs

Extensions

More detailed definition, comments rephrased, non-ascii characters in URL's removed - R. J. Mathar, Sep 09 2009
a(8)-a(9), a(11), a(18) from Jean-François Alcover, Dec 08 2017
Incorrect codes (depending on a search limit) removed, prime 2 removed, terms a(10), (12)-a(17), and a(19) onward added by Max Alekseyev, Sep 17 2024

A354026 Primes that divide some k dividing 4^k + 3^k (A045584).

Original entry on oeis.org

7, 379, 14407, 689431, 4235659, 41647747, 137534083, 239900179, 242121643, 349909477, 1245283747, 1478065891, 1605314383, 2500276549, 2748751303, 5618210347, 7490947129, 11236420693, 11260421089, 16948514941, 29440659361, 74163546829, 75093609319, 82188727303
Offset: 1

Views

Author

Max Alekseyev, May 15 2022

Keywords

Comments

Prime p > 3 is in this sequence iff all prime factors of the multiplicative order of -3/4 modulo p belong to this sequence.

Crossrefs

Programs

  • PARI
    S=[]; forprime(p=5,oo, f=Set(factor(znorder(Mod(-3/4,p)))[,1]); if(#setintersect(S,f)==#f, S=setunion(S,[p]); print1(p,", ")));

Extensions

a(18)-a(24) from Jinyuan Wang, Jan 29 2025
Showing 1-9 of 9 results.