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

A161509 The unique primitive prime factor of 2^n-1 for the n in A161508.

Original entry on oeis.org

3, 7, 5, 31, 127, 17, 73, 11, 13, 8191, 43, 151, 257, 131071, 19, 524287, 41, 337, 683, 241, 2731, 262657, 331, 2147483647, 65537, 599479, 43691, 174763, 61681, 5419, 2796203, 4432676798593, 87211, 15790321, 2305843009213693951, 715827883
Offset: 1

Views

Author

T. D. Noe, Jun 17 2009

Keywords

Comments

For these primes p, the binary expansion of 1/p has a unique period length. The binary analog of A007615.

Crossrefs

Cf. A144755 (sorted).

Programs

  • Mathematica
    Reap[Do[c=Cyclotomic[n,2]; q=c/GCD[c,n]; If[PrimePowerQ[q], Sow[FactorInteger[q][[1,1]]]],{n,100}]][[2,1]]

A007498 Unique period lengths of primes mentioned in A007615.

Original entry on oeis.org

1, 2, 3, 4, 9, 10, 12, 14, 19, 23, 24, 36, 38, 39, 48, 62, 93, 106, 120, 134, 150, 196, 294, 317, 320, 385, 586, 597, 654, 738, 945, 1031, 1172, 1282, 1404, 1426, 1452, 1521, 1752, 1812, 1836, 1844, 1862, 2134, 2232, 2264, 2667, 3750, 3903, 3927, 4274, 4354
Offset: 1

Views

Author

Keywords

Comments

Let {Zs(m, 10, 1)} be the Zsigmondy numbers for a = 10, b = 1: Zs(m, 10, 1) is the greatest divisor of 10^m - 1^m that is coprime to 10^r - 1^r for all positive integers r < m. Then this sequence gives m such that Zs(m, 10, 1) is a prime power (e.g., Zs(1, 10, 1) = 9 = 3^2, Zs(2, 10, 1) = 11, Zs(3, 10, 1) = 37, Zs(4, 10, 1) = 101). It is very likely that Zs(m, 10, 1) is prime if m > 1 is in this sequence (note that the Mathematica and PARI programs below are based on this assumption). - Jianing Song, Aug 12 2020

References

  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Samuel Yates, Period Lengths of Exactly One or Two Prime Numbers, J. Rec. Math., 18 (1985), 22-24.

Crossrefs

Cf. A161508 (unique period lengths in base 2).

Programs

  • Mathematica
    lst={1}; Do[p=Cyclotomic[n, 10]/GCD[n, Cyclotomic[n, 10]]; If[PrimeQ[p], AppendTo[lst, n]], {n, 3000}]; lst (* T. D. Noe, Sep 08 2005 *)
  • PARI
    isok(n) = if (n==1, 1, my(p = polcyclo(n, 10)); isprime(p/gcd(p, n))); \\ Michel Marcus, Jun 20 2018

Extensions

More terms from T. D. Noe, Sep 08 2005
a(48)-a(52) from Ray Chandler, Jul 09 2008

A250197 Numbers k such that the left Aurifeuillian primitive part of 2^k+1 is prime.

Original entry on oeis.org

10, 14, 18, 22, 26, 30, 42, 54, 58, 66, 70, 86, 94, 98, 106, 110, 126, 130, 138, 146, 158, 174, 186, 210, 222, 226, 258, 302, 334, 434, 462, 478, 482, 522, 566, 602, 638, 706, 734, 750, 770, 782, 914, 1062, 1086, 1114, 1126, 1226, 1266, 1358, 1382, 1434, 1742, 1926
Offset: 1

Views

Author

Eric Chen, Jan 18 2015

Keywords

Comments

All terms are congruent to 2 modulo 4.
Phi_n(x) is the n-th cyclotomic polynomial.
Numbers n such that Phi_{2nL(n)}(2) is prime.
Let J(n) = 2^n+1, J*(n) = the primitive part of 2^n+1, this is Phi_{2n}(2).
Let L(n) = the Aurifeuillian L-part of 2^n+1, L(n) = 2^(n/2) - 2^((n+2)/4) + 1 for n congruent to 2 (mod 4).
Let L*(n) = GCD(L(n), J*(n)).
This sequence lists all n such that L*(n) is prime.

Examples

			14 is in this sequence because the left Aurifeuillian primitive part of 2^14+1 is 113, which is prime.
34 is not in this sequence because the left Aurifeuillian primitive part of 2^34+1 is 130561, which equals 137 * 953 and is not prime.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[2000], Mod[#, 4] == 2 && PrimeQ[GCD[2^(#/2) - 2^((#+2)/4) + 1, Cyclotomic[2*#, 2]]] &]
  • PARI
    isok(n) = isprime(gcd(2^(n/2) - 2^((n+2)/4) + 1, polcyclo(2*n, 2))); \\ Michel Marcus, Jan 27 2015

A143584 Integers that are equal to the multiplicative order of 2 modulo some overpseudoprime to base 2.

Original entry on oeis.org

11, 23, 25, 28, 29, 35, 36, 37, 39, 41, 43, 44, 45, 47, 48, 50, 51, 52, 53, 55, 57, 58, 59, 60, 63, 64, 66, 67, 68, 70, 71, 72, 73, 74, 75, 76, 79, 81, 82, 83, 84, 87, 88, 91, 92, 94, 95, 96, 97, 99, 100, 101, 102, 103, 104, 105, 106, 108, 109, 110, 111, 112
Offset: 1

Views

Author

Vladimir Shevelev, Aug 25 2008

Keywords

Comments

A064078(a(n)) is a composite number. The sequence has a positive density since it contains, in particular, numbers of the form 8n+20 for n >= 1 (C. Pomerance, private correspondence). Since, e.g., 38 is not in the sequence, there is not an overpseudoprime m such that ord_m(2)=38.
Phi_{a(n)}(2), the a(n)-th cyclotomic polynomial of x evaluated at x=2 has at least 2 distinct prime factors that are not prime factors of the Phi_k(2) for any positive integer k < a(n). For example, Phi_11(2) = 2^11 - 1 = 2047 = 23 * 89 and Phi_25(2) = 2^20 + 2^15 + 2^10 + 2^5 + 1 = 1082401 = 601 * 1801. Note that p = a(n) is prime if and only if Phi_p(2) = 2^p - 1 is composite. - David Terr, Sep 09 2018
It is easy to prove the statement above. We use the fact that Phi_j(n) and Phi_k(n) are coprime whenever j and k are coprime as well as the fact that an overpseudoprime has at least 2 distinct prime factors. - David Terr, Oct 10 2018
A number k is included iff either 2^k-1 has more than one primitive prime factor (cf. A086251, A161508) or the only primitive prime factor of 2^k-1 is a Wieferich prime (no examples known). - Jeppe Stig Nielsen, Sep 01 2020

Crossrefs

Cf. A131952 (for the corresponding maximal overpseudoprimes).

Programs

  • PARI
    isok(k) = my(m=polcyclo(k,2)); m/=gcd(m,k); m!=1&&!isprime(m) \\ Jeppe Stig Nielsen, Sep 01 2020

Extensions

Name edited by Michel Marcus, Oct 06 2018
More terms from Michel Marcus, Oct 11 2018
Data for terms >= 100 corrected by Jeppe Stig Nielsen, Sep 01 2020

A247071 Numbers n such that 2^n-1 has only one primitive prime factor, sorted according to the magnitude of the corresponding prime.

Original entry on oeis.org

2, 4, 3, 10, 12, 8, 18, 5, 20, 14, 9, 7, 15, 24, 16, 30, 21, 22, 26, 42, 13, 34, 40, 32, 54, 17, 38, 27, 19, 33, 46, 56, 90, 78, 62, 31, 80
Offset: 1

Views

Author

Eric Chen, Nov 16 2014

Keywords

Comments

Periods associated with A144755 in base 2. The binary analog of A051627.

Examples

			2^12 - 1 = 4095 = 3 * 3 * 5 * 7 * 13, but none of 3, 5, 7 is a primitive prime factor, so the only primitive prime factor of 2^12 - 1 is 13.
		

Crossrefs

Programs

  • Mathematica
    nmax = 65536; primesPeriods = Reap[Do[p = Cyclotomic[n, 2]/GCD[n, Cyclotomic[n, 2]]; If[PrimeQ[p], Print[n]; Sow[{p, n}]], {n, 1, nmax}]][[2, 1]]; Sort[primesPeriods][[All, 2]]

Formula

a(n) = A002326((A144755(n+1)-1)/2). - Max Alekseyev, Feb 11 2024

Extensions

Sequence trimmed to the established terms of A144755 by Max Alekseyev, Feb 11 2024

A250198 Numbers k such that the right Aurifeuillian primitive part of 2^k+1 is prime.

Original entry on oeis.org

2, 6, 10, 14, 18, 22, 30, 34, 38, 42, 54, 58, 66, 70, 90, 102, 110, 114, 126, 138, 170, 178, 242, 294, 314, 326, 350, 378, 462, 566, 646, 726, 758, 1150, 1242, 1302, 1482, 1558, 1638, 1710, 1770, 1970, 1994
Offset: 1

Views

Author

Eric Chen, Jan 18 2015

Keywords

Comments

All terms are congruent to 2 modulo 4.
Let Phi_n(x) denote the n-th cyclotomic polynomial.
Numbers n such that Phi_{2nM(n)}(2) is prime.
Let J(n) = 2^n+1, J*(n) = the primitive part of 2^n+1, and this is Phi_{2n}(2).
Let M(n) = the Aurifeuillian M-part of 2^n+1, M(n) = 2^(n/2) + 2^((n+2)/4) + 1 for n congruent to 2 (mod 4).
Let M*(n) = GCD(M(n), J*(n)), this sequence lists all n such that M*(n) is prime.

Examples

			14 is in this sequence because the right Aurifeuillian primitive part of 2^14+1 is 29, which is prime.
26 is not in this sequence because the right Aurifeuillian primitive part of 2^26+1 is 8321, which equals 53 * 157 and is not prime.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[2000], Mod[#, 4] == 2 && PrimeQ[GCD[2^(#/2) + 2^((#+2)/4) + 1, Cyclotomic[2*#, 2]]] &]
  • PARI
    isok(n) = isprime(gcd(2^(n/2) + 2^((n+2)/4) + 1, polcyclo(2*n, 2))); \\ Michel Marcus, Jan 27 2015

A250208 Ratio of the primitive part of 2^n-1 to the product of primitive prime factors of 2^n-1.

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 5, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 1, 1, 11, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Eric Chen, Mar 02 2015

Keywords

Comments

As with A178764, it can be shown that all terms are either 1 or prime.
a(2*3^n) = 3 (n>=1).
a(4*5^n) = 5 (n>=1).
a(3*7^n) = 7 (n>=1).
a(10*11^n) = 11 (n>=1).
a(12*13^n) = 13 (n>=1).
a(8*17^n) = 17 (n>=1).
a(18*19^n) = 19 (n>=1).
...
a(A014664(k)*prime(k)^n) = prime(k).
For other n (while Phi_n(2) is squarefree), a(n) = 1.
a(n) != 1 for n = {6, 18, 20, 21, 54, 100, 110, 136, 147, 155, 156, 162, ...}.
At least, a(A049093(n)) = 1. (In fact, since Phi_n(2) is not completely factored for n = 991, 1207, 1213, 1217, 1219, 1229, 1231, 1237, 1243, 1249, ..., so it is unknown whether they are squarefree or not, but it is likely that Phi_n(2) is squarefree for all n except 364 and 1755 (because it is likely 1093 and 3511 are the only two Wieferich primes), so a(991), a(1207), a(1213), ..., are likely to be 1.)

Examples

			a(11) = 1 since Phi_11(2) = (2^11-1)/(2-1) = 2047, and the primitive prime factors of 2^11-1 are 23 and 89, so a(11) = 2047/(23*89) = 1.
a(18) = 3 since Phi_18(2) = 2^6 - 2^3 + 1 = 57, and the only primitive prime factor of 2^18-1 is 19, so a(18) = 57/19 = 3.
		

Crossrefs

Programs

  • Mathematica
    a250208[n_] = If[n == 364, 1093, If[n == 1755, 3511, GCD[Cyclotomic[n, 2], n]]]; Table[a250208[n], {n, 0, 200}]
  • PARI
    a(n) = if (n==364, 1093, if (n==1755, 3511, gcd(polcyclo(n, 2), n)));
    
  • PARI
    isprimitive(p, n) = {for (r=1, n-1, if (((2^r-1) % p) == 0, return (0)); ); return (1); }
    ppf(n) = {my(pf = factor(2^n-1)[,1]); prod(k=1,#pf, if (isprimitive(pf[k], n), pf[k], 1));}
    a(n) = if (issquarefree(m=polcyclo(n,2)), gcd(m, n), m/ppf(n)); \\ Michel Marcus, Mar 06 2015

Formula

a(n) = A019320(n) / A064078(n) while Phi_n(2) is squarefree.
a(n) = GCD(Phi_n(2), n) while Phi_n(2) is squarefree.
Notice: a(364) = 1093, a(1755) = 3511. (See A001220.)

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, ", "))
Showing 1-8 of 8 results.