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 24 results. Next

A333973 Numbers k such that A019320(k) is greater than A064078(k) and the latter is a prime or a prime power.

Original entry on oeis.org

18, 20, 21, 54, 147, 342, 602, 889, 258121
Offset: 1

Views

Author

Jeppe Stig Nielsen, Sep 22 2020

Keywords

Comments

The unique prime factor of A064078(k) is then a unique prime to base 2 (see A161509), but not a cyclotomic number.
Subsequence of A161508. In fact, subsequence of the set difference A161508 \ A072226.
In all known examples, A064078(k) is a prime. If A064078(k) was a prime power p^j with j>1, then p would be both a Wieferich prime (A001220) and a unique prime to base 2.
Subsequence of A093106 (the characterization of A093106 can be useful when searching for more terms).
Should this sequence be infinite?

Crossrefs

Programs

  • PARI
    for(n=1,+oo,c=polcyclo(n,2); c % n < 2 && next(); c/=(c%n); ispseudoprime(if(ispower(c,,&b),b,c))&&print1(n, ", "))

A112927 a(n) is the least prime such that the multiplicative order of 2 mod a(n) equals n, or a(n)=1 if no such prime exists.

Original entry on oeis.org

1, 3, 7, 5, 31, 1, 127, 17, 73, 11, 23, 13, 8191, 43, 151, 257, 131071, 19, 524287, 41, 337, 683, 47, 241, 601, 2731, 262657, 29, 233, 331, 2147483647, 65537, 599479, 43691, 71, 37, 223, 174763, 79, 61681, 13367, 5419, 431, 397, 631, 2796203, 2351, 97, 4432676798593, 251, 103, 53, 6361, 87211
Offset: 1

Views

Author

Vladimir Shevelev, Aug 25 2008

Keywords

Comments

If a(n) differs from 1, then a(n) is the minimal prime divisor of A064078(n);
a(n)=n+1 iff n+1 is prime from A001122; a(n)=2n+1 iff 2n+1 is prime from A115591.
If a(n) > 1 then a(n) is the index where n occurs first in A014664. - M. F. Hasler, Feb 21 2016
Bang's theorem (special case of Zsigmondy's theorem, see links): a(n)>1 for all n>6. - Jeppe Stig Nielsen, Aug 31 2020

Crossrefs

Cf. A112927 (base 2), A143663 (base 3), A112092 (base 4), A143665 (base 5), A379639 (base 6), A379640 (base 7), A379641 (base 8), A379642 (base 9), A007138 (base 10), A379644 (base 11), A252170 (base 12).

Programs

  • PARI
    A112927(n,f=factor(2^n-1)[,1])=!for(i=1,#f,znorder(Mod(2,f[i]))==n&&return(f[i])) \\ Use the optional 2nd arg to give a list of pseudoprimes to try when factoring of 2^n-1 is too slow. You may try factor(2^n-1,0)[,1]. - M. F. Hasler, Feb 21 2016

A085021 Number of prime factors of cyclotomic(n,2), which is A019320(n), the value of the n-th cyclotomic polynomial evaluated at x=2.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 3, 1, 1, 1, 1, 1, 2, 2, 2, 1, 2, 1, 2, 1, 3, 2, 2, 1, 3, 2, 1, 2, 3, 3, 3, 2, 3, 1, 2, 2, 2, 2, 1, 1, 2, 2, 1, 2, 2, 3, 1, 2, 3, 2, 3, 2, 2, 3, 1, 1, 3, 1, 3, 2, 2, 2, 1, 1, 2, 2, 1, 1, 3, 4, 1, 2, 3, 2, 2, 1, 3, 4
Offset: 1

Views

Author

T. D. Noe, Jun 19 2003

Keywords

Comments

The Mobius transform of this sequence yields A046051, the number of prime factors of Mersenne number 2^n-1.
The number of prime factors in the primitive part of 2^n-1. - T. D. Noe, Jul 19 2008

Examples

			a(11) = 2 because cyclotomic(11,2) = 2047, which has two factors: 23 and 89.
		

Crossrefs

omega(Phi(n,x)): this sequence (x=2), A085028 (x=3), A085029 (x=4), A085030 (x=5), A085031 (x=6), A085032 (x=7), A085033 (x=8), A085034 (x=9), A085035 (x=10).

Programs

  • Mathematica
    Join[{0}, Table[Plus@@Transpose[FactorInteger[Cyclotomic[n, 2]]][[2]], {n, 2, 100}]]
  • PARI
    a(n) = #factor(polcyclo(n, 2))~; \\ Michel Marcus, Mar 06 2015

A064080 Zsigmondy numbers for a = 4, b = 1: Zs(n, 4, 1) is the greatest divisor of 4^n - 1^n (A024036) that is relatively prime to 4^m - 1^m for all positive integers m < n.

Original entry on oeis.org

3, 5, 7, 17, 341, 13, 5461, 257, 1387, 41, 1398101, 241, 22369621, 3277, 49981, 65537, 5726623061, 4033, 91625968981, 61681, 1826203, 838861, 23456248059221, 65281, 1100586419201, 13421773, 22906579627, 15790321, 96076792050570581
Offset: 1

Views

Author

Jens Voß, Sep 04 2001

Keywords

Comments

By Zsigmondy's theorem, the n-th Zsigmondy number for bases a and b is not 1 except in the three cases (1) a = 2, b = 1, n = 1, (2) a = 2, b = 1, n = 6, (3) n = 2 and a+b is a power of 2.

Crossrefs

Formula

For even n, a(n) = A064078(2*n); for odd n, a(n) = A064078(n) * A064078(2*n). - Max Alekseyev, Apr 28 2022

Extensions

Corrected and extended by Vladeta Jovovic, Sep 05 2001
Definition corrected by Jerry Metzger, Nov 04 2009

A086251 Number of primitive prime factors of 2^n - 1.

Original entry on oeis.org

0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 2, 3, 1, 1, 1, 1, 1, 2, 2, 2, 1, 2, 1, 2, 1, 3, 2, 2, 1, 3, 2, 1, 2, 3, 3, 3, 1, 3, 1, 2, 2, 2, 2, 1, 1, 2, 2, 1, 2, 2, 3, 1, 2, 3, 2, 3, 2, 2, 3, 1, 1, 3, 1, 3, 2, 2, 2, 1, 1, 2, 2, 1, 1, 3, 4, 1, 2, 3, 2, 2, 1, 3, 3, 2, 3, 2, 2, 3
Offset: 1

Views

Author

T. D. Noe, Jul 14 2003

Keywords

Comments

A prime factor of 2^n - 1 is called primitive if it does not divide 2^r - 1 for any r < n. Equivalently, p is a primitive prime factor of 2^n - 1 if ord(2,p) = n. Zsigmondy's theorem says that there is at least one primitive prime factor for n > 1, except for n=6. See A086252 for those n that have a record number of primitive prime factors.
Number of odd primes p such that A002326((p-1)/2) = n. Number of occurrences of number n in A014664. - Thomas Ordowski, Sep 12 2017
The prime factors are not counted with multiplicity, which matters for a(364)=4 and a(1755)=6. - Jeppe Stig Nielsen, Sep 01 2020

Examples

			a(11) = 2 because 2^11 - 1 = 23*89 and both 23 and 89 have order 11.
		

Crossrefs

Cf. A046800, A046051 (number of prime factors, with repetition, of 2^n-1), A086252, A002588, A005420, A002184, A046801, A049093, A049094, A059499, A085021, A097406, A112927, A237043.

Programs

  • Mathematica
    Join[{0}, Table[cnt=0; f=Transpose[FactorInteger[2^n-1]][[1]]; Do[If[MultiplicativeOrder[2, f[[i]]]==n, cnt++ ], {i, Length[f]}]; cnt, {n, 2, 200}]]
  • PARI
    a(n) = sumdiv(n, d, moebius(n/d)*omega(2^d-1)); \\ Michel Marcus, Sep 12 2017
    
  • PARI
    a(n) = my(m=polcyclo(n, 2)); omega(m/gcd(m,n)) \\ Jeppe Stig Nielsen, Sep 01 2020

Formula

a(n) = Sum{d|n} mu(n/d) A046800(d), inverse Mobius transform of A046800.
a(n) <= A182590(n). - Thomas Ordowski, Sep 14 2017
a(n) = A001221(A064078(n)). - Thomas Ordowski, Oct 26 2017

Extensions

Terms to a(500) in b-file from T. D. Noe, Nov 11 2010
Terms a(501)-a(1200) in b-file from Charles R Greathouse IV, Sep 14 2017
Terms a(1201)-a(1206) in b-file from Max Alekseyev, Sep 11 2022

A060443 Table T(n,k) in which n-th row lists prime factors of 2^n - 1 (n >= 2), without repetition.

Original entry on oeis.org

0, 1, 3, 7, 3, 5, 31, 3, 7, 127, 3, 5, 17, 7, 73, 3, 11, 31, 23, 89, 3, 5, 7, 13, 8191, 3, 43, 127, 7, 31, 151, 3, 5, 17, 257, 131071, 3, 7, 19, 73, 524287, 3, 5, 11, 31, 41, 7, 127, 337, 3, 23, 89, 683, 47, 178481, 3, 5, 7, 13, 17, 241
Offset: 0

Views

Author

Keywords

Comments

For n > 1, the length of row n is A046800(n). - T. D. Noe, Aug 06 2007

Examples

			From _Wolfdieter Lang_, Sep 23 2017: (Start)
The irregular triangle T(n,k) begins for n >= 2:
n\k      1   2    3    4   5
2:       3
3:       7
4:       3   5
5:      31
6:       3   7
7:     127
8:       3   5   17
9:       7  73
10:      3  11   31
11:     23  89
12:      3   5    7   13
13:   8191
14:      3  43  127
15:      7  31  151
16:      3   5   17  257
17: 131071
18:      3   7   19   73
19: 524287
20:      3   5   11   31  41
... (End)
		

References

  • J. Brillhart et al., Factorizations of b^n +- 1. Contemporary Mathematics, Vol. 22, Amer. Math. Soc., Providence, RI, 2nd edition, 1985; and later supplements.

Crossrefs

Programs

  • Mathematica
    Array[FactorInteger[2^# - 1][[All, 1]] &, 25, 0] (* Paolo Xausa, Apr 18 2024 *)

A064081 Zsigmondy numbers for a = 5, b = 1: Zs(n, 5, 1) is the greatest divisor of 5^n - 1^n (A024049) that is relatively prime to 5^m - 1^m for all positive integers m < n.

Original entry on oeis.org

4, 3, 31, 13, 781, 7, 19531, 313, 15751, 521, 12207031, 601, 305175781, 13021, 315121, 195313, 190734863281, 5167, 4768371582031, 375601, 196890121, 8138021, 2980232238769531, 390001, 95397958987501, 203450521, 3814699218751, 234750601, 46566128730773925781, 464881, 1164153218269348144531
Offset: 1

Views

Author

Jens Voß, Sep 04 2001

Keywords

Comments

By Zsigmondy's theorem, the n-th Zsigmondy number for bases a and b is not 1 except in the three cases (1) a = 2, b = 1, n = 1, (2) a = 2, b = 1, n = 6, (3) n = 2 and a+b is a power of 2.

Crossrefs

Extensions

More terms from Vladeta Jovovic, Sep 06 2001
Definition corrected by Jerry Metzger, Nov 04 2009
More terms from Robert Israel, Feb 21 2025

A064083 Zsigmondy numbers for a = 7, b = 1: Zs(n, 7, 1) is the greatest divisor of 7^n - 1^n (A024075) that is relatively prime to 7^m - 1^m for all positive integers m < n.

Original entry on oeis.org

6, 1, 19, 25, 2801, 43, 137257, 1201, 39331, 2101, 329554457, 2353, 16148168401, 102943, 4956001, 2882401, 38771752331201, 117307, 1899815864228857, 1129901, 11898664849, 247165843, 4561457890013486057, 5762401, 79797014141614001
Offset: 1

Views

Author

Jens Voß, Sep 04 2001

Keywords

Comments

By Zsigmondy's theorem, the n-th Zsigmondy number for bases a and b is not 1 except in the three cases (1) a = 2, b = 1, n = 1, (2) a = 2, b = 1, n = 6, (3) n = 2 and a+b is a power of 2.

Crossrefs

Extensions

More terms from Vladeta Jovovic, Sep 06 2001
Definition corrected by Jerry Metzger, Nov 04 2009

A063982 Number of divisors of 2^n - 1 that are relatively prime to 2^m - 1 for all 0 < m < n.

Original entry on oeis.org

1, 2, 2, 2, 2, 1, 2, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 4, 2, 4, 2, 2, 4, 8, 2, 2, 2, 2, 2, 4, 4, 4, 2, 4, 2, 4, 2, 8, 4, 4, 2, 8, 4, 2, 4, 8, 8, 8, 2, 8, 2, 4, 4, 4, 4, 2, 2, 4, 4, 2, 4, 4, 8, 2, 4, 8, 4, 8, 4, 4, 8, 2, 2, 8, 2, 8, 4, 4, 4, 2, 2, 4, 4, 2, 2, 8, 16, 2, 4, 8, 4, 4, 2, 8, 8
Offset: 1

Views

Author

Vladeta Jovovic, Sep 06 2001

Keywords

Comments

a(364) = 24 is the first term not a power of 2. - Jianing Song, Apr 29 2018
a(n) is the number of divisors of A064078(n). - Jianing Song, Apr 20 2019

Examples

			Divisors of 2^8-1 are {1, 3, 5, 15, 17, 51, 85, 255}, but only 1 and 17 are relatively prime to 2^m - 1 for all m < 8, thus a(8)=2.
		

Crossrefs

Programs

  • Haskell
    a063982 n = a063982_list !! (n-1)
    a063982_list = f [] $ tail a000225_list where
       f us (v:vs) = (length ds) : f (v:us) vs where
         ds = [d | d <- a027750_row v, all ((== 1). (gcd d)) us]
    -- Reinhard Zumkeller, Jan 04 2013
    
  • Mathematica
    a = {1}; Do[ d = Divisors[2^n - 1]; l = Length[d]; c = 0; k = 1; While[ k < l + 1, If[ Union[ GCD[a, d[[k]] ]] == {1}, c++ ]; k++ ]; Print[c]; a = Union[ Flatten[ Append[a, Transpose[ FactorInteger[2^n - 1]][[ 1]] ]]], {n, 1, 100} ]
  • PARI
    a(n) = {my(v = vector(n-1, k, 2^k-1), na = 0, nb); fordiv(2^n-1, d, nb = 0; for (k=1, n-1, if (gcd(d, v[k]) == 1, nb++, break);); if (nb == n-1, na++);); return (na);} \\ Michel Marcus, Apr 30 2018

Extensions

More terms from Robert G. Wilson v, Sep 10 2001

A064079 Zsigmondy numbers for a = 3, b = 1: Zs(n, 3, 1) is the greatest divisor of 3^n - 1^n (A024023) that is relatively prime to 3^m - 1^m for all positive integers m < n.

Original entry on oeis.org

2, 1, 13, 5, 121, 7, 1093, 41, 757, 61, 88573, 73, 797161, 547, 4561, 3281, 64570081, 703, 581130733, 1181, 368089, 44287, 47071589413, 6481, 3501192601, 398581, 387440173, 478297, 34315188682441, 8401, 308836698141973, 21523361
Offset: 1

Views

Author

Jens Voß, Sep 04 2001

Keywords

Comments

By Zsigmondy's theorem, the n-th Zsigmondy number for bases a and b is not 1 except in the three cases (1) a = 2, b = 1, n = 1, (2) a = 2, b = 1, n = 6, (3) n = 2 and a+b is a power of 2.

Crossrefs

Extensions

More terms from Vladeta Jovovic, Sep 06 2001
Definition corrected by Jerry Metzger, Nov 04 2009
Showing 1-10 of 24 results. Next