OFFSET
1,1
LINKS
PROG
(PARI) is_A391333(n) = if(!(n%2) || 2!=bigomega(n), 0, my(f=factor(n), a = f[1, 1], b = f[#f~, 1], Pa=Pol(digits(a, 3))*Mod(1, 3), Pb=Pol(digits(b, 3))*Mod(1, 3), Pn=Pol(digits(n, 3))*Mod(1, 3)); ((0==lift(Pn % Pa)) || (0==lift(Pn % Pb))));
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Antti Karttunen, Dec 07 2025
STATUS
approved
