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.

A286499 Primes which divide a term of A073935.

This page as a plain text file.
%I A286499 #22 Feb 23 2022 23:07:06
%S A286499 2,3,5,7,17,19,43,101,163,257,487,1459,14407,26407,39367,62501,65537,
%T A286499 77659,1020101,1336337,86093443,242121643,258280327,3103616899,
%U A286499 4528177054183,15258789062501,411782264189299,21108889701347407,953735353027359375062501
%N A286499 Primes which divide a term of A073935.
%C A286499 A prime p is in this sequence if and only if p-1 = Product_{i} (p_i)^(a_i) with p_j - 1 = Product_{j<i} (p_j)^(a_j).
%C A286499 This sequence contains all Fermat primes (A019434).
%H A286499 John Machacek, <a href="https://arxiv.org/abs/1706.01008">Egyptian Fractions and Prime Power Divisors</a>, arXiv:1706.01008 [math.NT], 2017.
%e A286499 p = 43 is in the sequence because 43-1 = 42 = 2*3*7, 7-1 = 6 = 2*3, 3-1 = 2.
%t A286499 upTo[mx_] := Block[{ric}, ric[n_, p_] := If[n < mx, Block[{m = n p}, If[PrimeQ[n + 1], Sow[n+1]; ric[n (n + 1), n+1]]; If[IntegerExponent[n, p] == 1, While[m < mx, ric[m, p]; m *= p]]]]; Sort[Reap[ric[1, 2]][[2, 1]]]]; upTo[10^20] (* _Giovanni Resta_, May 27 2017 *)
%Y A286499 Cf. A073935.
%K A286499 nonn
%O A286499 1,1
%A A286499 _John Machacek_, May 27 2017
%E A286499 a(20)-a(29) from _Giovanni Resta_, May 27 2017