OFFSET
1,1
COMMENTS
A squarefree subsequence of heptagonal numbers.
EXAMPLE
32890 = 2*5*11*13*23 = 115(5*115-3)/2.
48790 = 2*5*7*17*41 = 140(5*140-3)/2.
102718 = 2*7*11*23*29 = 203(5*203-3)/2.
167314 = 2*7*17*19*37 = 259(5*259-3)/2.
MATHEMATICA
Select[Table[n*(5*n - 3)/2, {n, 1, 1000}], FactorInteger[#][[;; , 2]] == {1, 1, 1, 1, 1} &] (* Amiram Eldar, Apr 17 2022 *)
PROG
(Python)
from sympy import factorint
from itertools import count, islice
def agen():
for h in (n*(5*n-3)//2 for n in count(1)):
f = factorint(h, multiple=True)
if len(f) == len(set(f)) == 5: yield h
print(list(islice(agen(), 34))) # Michael S. Branicky, May 28 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Massimo Kofler, Apr 17 2022
STATUS
approved