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.
%I A129472 #3 Mar 31 2012 13:48:24 %S A129472 3181,4513,4957,6067,7177,8731,9397,10433,13171,14947,15761,17389, %T A129472 19387,19609,22051,22273,22453,22717,23531,23753,24197,26161,27823, %U A129472 28711,37369,37591,38183,38923,39293,40993,41143,42697,43067,44621,44843 %N A129472 Primes p of Erdos-Selfridge class 4+ with largest prime factor of p+1 not of class 3+. %C A129472 See A129470 %e A129472 a(1) = 3181 = -1+2*37*43 is a prime of class 4+ since 37 is of class 3+, but the largest divisor of 3181+1 is 43 which is only of class 2+. %o A129472 (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])}; A129472(n=100,p=1,a=[])={ local(f); while( #a<n, until( f[ #f] > 3 & 3 > class(f[ #f]), f=factor(1+p=nextprime(p+1))[,1]); forstep( i=#f-1,2,-1, if( 4 < f[1] = max( f[1],1+class( f[i] )), next(2))); if( f[1] == 4, a=concat(a,p); /*print(#a," ",p)*/)); a} %Y A129472 Subsequence of A129470; see also A129471, A129473, A129477-A129478, A129469, A005113, A005105-A005107. %K A129472 easy,nonn %O A129472 1,1 %A A129472 _M. F. Hasler_, Apr 17 2007