OFFSET
1,1
PROG
(PARI) fcnt(n, m) = {local(s); s=0; if(n == 1, s=1, fordiv(n, d, if(d > 1 & d <= m, s=s+fcnt(n/d, d)))); s};
f(n) = n - fcnt(n, n);
is(n) = j=1; x=2^(floor(log(n)/log(2))+2); while(j<x&&f(j)!=n, j++); f(j)!=n; \\ Jinyuan Wang, Feb 14 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Aidan Clarke, Jan 21 2019
EXTENSIONS
More terms from Jinyuan Wang, Feb 14 2019
STATUS
approved