OFFSET
1,1
COMMENTS
The least terms not divisible by 3 are 127250, 494750, 578830, ... .
LINKS
Karl-Heinz Hofmann, Table of n, a(n) for n = 1..10000
EXAMPLE
PROG
(PARI) a396672(upto=7500) = my(p1=15, p2); forstep(k=21, upto, 2, my(f=factor(k), m); if(omega(f)==2 && bigomega(f)==2, p2=k; m=p1+p2; if(gcd(p1, m)!=1 || gcd(p2, m)!=1, print1(m/2, ", ")); p1=p2))
(Python)
from sympy import factorint
def A046388_isok_and_factorization(n):
return [n % 2 == 1 and (len(f := factorint(n, multiple=True)) == 2 and f[0] != f[1]), f]
while len(A396672) < (aupto := 47):
k += 2
if (r:=A046388_isok_and_factorization(k))[0]:
A046388.append(k)
if (f := factorint(A046388[-2], multiple=True))[0] in r[1] or f[1] in r[1]:
print(A396672) # Karl-Heinz Hofmann, Jun 02 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Hugo Pfoertner, Jun 02 2026
STATUS
approved
