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.

A361215 Intersection of A361073 and 2 * A361611.

Original entry on oeis.org

8, 20, 50, 1406, 1516, 1558, 1868, 1898, 1948, 1978, 1986, 5862, 5972, 6014, 7122, 7966, 7996, 8270, 8348, 8366, 8548, 8618, 21092, 31804, 31822, 32158, 33092, 33162, 33316, 33414, 37124, 37190, 37292, 37394, 39164, 39214, 39316, 39346, 39484, 39562, 39604, 39622, 39692, 39794, 45044, 45244
Offset: 1

Views

Author

Zak Seidov and Robert Israel, Apr 09 2023

Keywords

Comments

If A361073(j) = 2*A361611(k) then x = 2*A361611(k+1) has the property that x, x - A361073(j) and x + A361073(j) are triprimes, so x >= A361073(j+1), with equality if and only if A361073(j+1) is even.

Examples

			a(4) = 1406 is a term because 1406 = A361073(20) = 2*A361611(17).
		

Crossrefs

Programs

  • Maple
    A:= {8}: lasta:= 8:
    for i from 2 to 1000 do
      for x from lasta+8 do
        if numtheory:-bigomega(x) = 3 and numtheory:-bigomega(x-lasta) = 3 and numtheory:-bigomega(x+lasta) = 3 then
           A:= A union {x}; lasta:= x; break
        fi
    od od:
    R:= {8}: lastb:= 4:
    while 2*lastb < lasta do
    for x from lastb+4 do
      if numtheory:-bigomega(x) = 2 and numtheory:-bigomega(x-lastb) = 2 and numtheory:-bigomega(x+lastb) = 2 then
         if member(2*x,A) then R:= R union {2*x} fi;
         lastb:= x; break
      fi
    od od:
    sort(convert(R,list));