OFFSET
1,4
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
PROG
(PARI)
A032742(n) = if(1==n, n, n/vecmin(factor(n)[, 1]));
(Python)
from sympy import factorint
def a(n): return 1 if n==1 else (n//min(factorint(n)))%4
print([a(n) for n in range(1, 106)]) # Michael S. Branicky, Apr 26 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Apr 22 2022
STATUS
approved