OFFSET
1,1
LINKS
John Cerkan, Table of n, a(n) for n = 1..10000
MATHEMATICA
pmp11Q[n_]:=Module[{fi=FactorInteger[n][[All, 1]]}, Mod[fi[[2]], fi[[1]]] == 11]; Select[ Range[12000], PrimeNu[#]==PrimeOmega[#]==2&&pmp11Q[#]&] (* Harvey P. Dale, Jun 25 2018 *)
PROG
(Python)
from sympy import factorint
def is_A064909(n):
f = factorint(n)
return (sum([f[i] for i in f]) == 2) and (max(f) % min(f) == 11)
def first_A064909(n):
x = 1
an = []
while len(an) < n:
if is_A064909(x): an.append(x)
x += 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]) == 11);
CROSSREFS
KEYWORD
nonn
AUTHOR
Patrick De Geest, Oct 13 2001
EXTENSIONS
Offset and name fixed by John Cerkan, Apr 12 2018
STATUS
approved