OFFSET
1,1
COMMENTS
Integers k such that k and k+1 are in A046305.
EXAMPLE
1215 = 3^5*5 is a product of 6 primes. The next integer, 1216 = 2^6*19, is a product of 7 primes. Thus, 1215 is in this sequence.
MAPLE
q:= n-> andmap(x-> numtheory[bigomega](x)>5, [n, n+1]):
select(q, [$1..30000])[]; # Alois P. Heinz, Jul 10 2021
MATHEMATICA
Select[Range[100000], Total[Transpose[FactorInteger[#]][[2]]] >= 6 && Total[Transpose[FactorInteger[# + 1]][[2]]] >= 6 &]
PROG
(Python)
from sympy import factorint
def prod6(n): return sum(factorint(n).values()) >= 6
def aupto(lim): return [k for k in range(lim+1) if prod6(k) and prod6(k+1)]
print(aupto(27135)) # Michael S. Branicky, Jul 10 2021
(PARI) isA346207(k) = (bigomega(k) >= 6) && (bigomega(k+1) >= 6) \\ Jianing Song, Jul 10 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Tanya Khovanova, Jul 10 2021
STATUS
approved