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 A056809 #55 Mar 03 2024 16:02:44 %S A056809 33,85,93,121,141,201,213,217,301,393,445,633,697,841,921,1041,1137, %T A056809 1261,1345,1401,1641,1761,1837,1893,1941,1981,2101,2181,2217,2305, %U A056809 2361,2433,2461,2517,2641,2721,2733,3097,3385,3601,3693,3865,3901,3957,4285 %N A056809 Numbers k such that k, k+1 and k+2 are products of two primes. %C A056809 Each term is the beginning of a run of three 2-almost primes (semiprimes). No runs exist of length greater than three. For the same reason, each term must be odd: If k were even, then so would be k+2. In fact, one of k or k+2 would be divisible by 4, so must indeed be 4 to have only two prime factors. However, neither 2,3,4 nor 4,5,6 is such a run. - _Rick L. Shepherd_, May 27 2002 %C A056809 k+1, which is twice a prime, is in A086005. The primes are in A086006. - _T. D. Noe_, May 31 2006 %C A056809 The squarefree terms are listed in A039833. - _Jianing Song_, Nov 30 2021 %H A056809 D. W. Wilson, <a href="/A056809/b056809.txt">Table of n, a(n) for n = 1..10000</a> %F A056809 a(n) = A086005(n) - 1 = 2*A086006(n) - 1 = 4*A123255(n) + 1. - _Jianing Song_, Nov 30 2021 %e A056809 121 is in the sequence because 121 = 11^2, 122 = 2*61 and 123 = 3*41, each of which is the product of two primes. %t A056809 f[n_] := Plus @@ Transpose[ FactorInteger[n]] [[2]]; Select[Range[10^4], f[ # ] == f[ # + 1] == f[ # + 2] == 2 & ] %t A056809 Flatten[Position[Partition[PrimeOmega[Range[5000]],3,1],{2,2,2}]] (* _Harvey P. Dale_, Feb 15 2015 *) %t A056809 SequencePosition[PrimeOmega[Range[5000]],{2,2,2}][[;;,1]] (* _Harvey P. Dale_, Mar 03 2024 *) %o A056809 (PARI) forstep(n=1,5000,2, if(bigomega(n)==2 && bigomega(n+1)==2 && bigomega(n+2)==2, print1(n,","))) %o A056809 (PARI) is(n)=n%4==1 && isprime((n+1)/2) && bigomega(n)==2 && bigomega(n+2)==2 \\ _Charles R Greathouse IV_, Sep 08 2015 %o A056809 (PARI) list(lim)=my(v=List(),t); forprime(p=2,(lim+1)\2, if(bigomega(t=2*p-1)==2 && bigomega(t+2)==2, listput(v,t))); Vec(v) \\ _Charles R Greathouse IV_, Sep 08 2015 %Y A056809 Intersection of A070552 and A092207. %Y A056809 Cf. A045939, A039833, A086005, A086006, A123255. %K A056809 nonn %O A056809 1,1 %A A056809 Sharon Sela (sharonsela(AT)hotmail.com), May 04 2002 %E A056809 Edited and extended by _Robert G. Wilson v_, May 04 2002