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

A361899 a(n) = 3*(6858365065530*(2^45 - 1)*n + 153479820268467961)^2.

Original entry on oeis.org

70668165688923686196507258250492563, 174687593550891106640307045856561008882907291372256643, 698750373759134872171732581703201135992894186495330123, 1572188340624731296664944773228844067526467943619713003
Offset: 0

Views

Author

Arkadiusz Wesolowski, Mar 28 2023

Keywords

Comments

The "1/k" heuristic predicts that primes of the form k*2^m + 1 with k odd and m > 0 have almost a 1/k chance of being Fermat divisors (Dubner and Keller). This sequence yields a correction to the "1/k" heuristic, because it generates special values of k.
If:
1) k is of the form 3*a^2, where a is an odd positive integer not divisible by 3,
2) k is not a SierpiƄski number,
3) for all odd positive integers m the numbers k*2^m + 1 are composite,
then the probability that a Fermat number is divisible by a prime of the form k*2^m + 1 equal to 0.
Every term meets the first and third condition. For any n, at least one of the primes from A361898 (except 3) divides every integer in the sequence a(n)*2^m + 1 with m odd.
What is the smallest odd integer k such that every prime of the form k*2^m + 1 (m > 0) does not divide any Fermat number?

References

  • H. Suyama, A note on the factors of Fermat numbers II, Abstracts of Papers Presented to the Amer. Math. Soc., Vol. 5 (1984), p. 132.

Crossrefs

Programs

  • Magma
    [3*(6858365065530*(2^45-1)*n+153479820268467961)^2: n in [0..3]];
    
  • Mathematica
    Table[3 (6858365065530 (2^45 - 1) n + 153479820268467961)^2, {n, 0, 3}]
  • PARI
    a(n)=3*(6858365065530*(2^45-1)*n+153479820268467961)^2

Formula

G.f.: (70668165688923686196507258250492563 + 174687593550891106428302548789789950293385516620778954*x + 174687593106461552462815941200289167933694087130037883*x^2)/(1 - x)^3.
a(n) = 3*(2*(Product_{i=1..13} A361898(i))*n + 153479820268467961)^2.
a(n) = 3*((29062/1192737)*(2^48 - 1)*(2^45 - 1)*n + 153479820268467961)^2.

A360652 Primes of the form x^2 + 432*y^2.

Original entry on oeis.org

433, 457, 601, 1657, 1753, 1777, 1801, 2017, 2089, 2113, 2281, 2689, 2833, 2953, 3457, 3889, 4057, 4129, 4153, 4177, 4513, 4657, 4729, 5113, 5209, 5449, 5569, 5737, 5953, 6217, 6361, 6673, 6961, 7057, 7321, 7369, 7537, 7753, 7873, 8353, 8377, 8713, 8761, 8929
Offset: 1

Views

Author

Arkadiusz Wesolowski, Feb 15 2023

Keywords

Comments

Supersequence of A351332. Thus every prime congruent to 1 mod 3 that divides a Fermat number is in this sequence.
Every Fermat number that is a semiprime has a prime of this form as a factor.

Crossrefs

Programs

  • Magma
    [p: p in PrimesUpTo(8929) | NormEquation(432, p) eq true];
    
  • PARI
    select(p->my(m=Mod(2, p)^(p\12)); p>11 && (m==1||m==p-1), primes(1110))

A361900 Numbers k such that 3*153479820268467961^2*2^k + 1 is prime.

Original entry on oeis.org

600, 810, 1074, 7974, 22290, 43086
Offset: 1

Views

Author

Arkadiusz Wesolowski, Mar 28 2023

Keywords

Comments

Let p be a prime number of the form 3*153479820268467961^2*2^k + 1 with k > 0, then the multiplicative order of 2 modulo p is not of the form 2^(m+1), m >= 0. Hence, p does not divide any Fermat number F(m) = 2^(2^m) + 1.

Crossrefs

Programs

  • Mathematica
    Select[Range[2, 10^4, 2], PrimeQ[3*153479820268467961^2*2^# + 1] &]
Showing 1-3 of 3 results.