>>> import heapq, itertools
>>> seq = heapq.merge((2**i for i in itertools.count(1, 2)),
(9*2**i for i in itertools.count(1, 2)))
>>> list(itertools.islice(seq, 0, 4))
[2, 8, 18, 32]
import itertools
def sigma(n):
return sum(i for i in range(1, n//2+1) if n % i == 0)
def seq():
for n in itertools.count(1):
if sigma(2*n) % 2 == 0: continue
if sigma(3*n) % 3 == 0: yield n
continue
>>> list(itertools.islice(seq(), 0, 4))
[2, 8, 18, 32]
The others are either degenerate for 4 points of the A001105 case, or are not so easy to implement in raw Python.
I believe you are referring to one of the many listed ways to define https://oeis.org/A001105 , which is even more succinctly expressed as a(n) = 2*n^2. This is the answer I thought you wanted, and which I was avoiding saying.
My underlying point is your sequence was under-specified, and with no clear reason to pick your preferred answer over other answers.
Even knowing it's chemistry, three finite sequences which match your specification are: