OFFSET
1,1
LINKS
John Cerkan, Table of n, a(n) for n = 1..10000
PROG
(Python)
from sympy import factorint
def is_A064902(n):
f = factorint(n)
return (sum([f[i] for i in f]) == 2) and (max(f) % min(f) == 4)
def first_A064902(n):
x = 1
an = []
while len(an) < n:
if is_A064902(x):an.append(x)
n += 2
return an # John Cerkan, Apr 14 2018
(PARI) isok(n) = my(f = factor(n)); (#f~ == 2) && (vecmax(f[, 2]) < 2) && ((f[2, 1] % f[1, 1]) == 4); \\ Michel Marcus, Apr 16 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Patrick De Geest, Oct 13 2001
EXTENSIONS
Offset changed by John Cerkan, Apr 12 2018
STATUS
approved