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.

A129473 Primes p of Erdos-Selfridge class 5+ with largest prime factor of p+1 not of class 4+.

This page as a plain text file.
%I A129473 #3 Mar 31 2012 13:48:24
%S A129473 15913,18541,22921,36353,47741,49201,52267,55333,60589,64969,66137,
%T A129473 66721,69203,72707,73291,74167,75773,78401,79861,80737,82051,84533,
%U A129473 90227,90373,95191,95483,95629,97673,99133,101323,103951,104681,104827
%N A129473 Primes p of Erdos-Selfridge class 5+ with largest prime factor of p+1 not of class 4+.
%C A129473 See A129470
%e A129473 a(1) = 15913 = -1+2*73*109 is a prime of class 5+ since 73 is of class 4+, but the largest divisor of 15913+1 is 109 which is only of class 2+.
%o A129473 (PARI) class(n,s=1)={n=factor(n+s)[,1];if(n[ #n]<=3,1,for(i=2,#n,n[1]=max(class(n[i],s)+1,n[1]));n[1])}; A129473(n=100,p=1,a=[])={ local(f); while( #a<n, until( f[ #f] > 3 & 4 > class(f[ #f]), f=factor(1+p=nextprime(p+1))[,1]); forstep( i=#f-1,2,-1, if( 5 < f[1] = max( f[1],1+class( f[i] )), next(2))); if( f[1] == 5, a=concat(a,p); /*print(#a," ",p)*/)); a}
%Y A129473 Subsequence of A129470; see also A129471-A129472, A129477-A129478, A129469, A005113, A005105-A005108, A081633.
%K A129473 easy,nonn
%O A129473 1,1
%A A129473 _M. F. Hasler_, Apr 17 2007