OFFSET
1,1
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
filter:= n -> n mod numtheory:-phi(n) <> 0:
select(filter, [$1..100]); # Robert Israel, Feb 04 2025
MATHEMATICA
Select[Range[100], ! Divisible[#, EulerPhi[#]] &] (* Amiram Eldar, Jan 08 2025 *)
PROG
(Python)
from sympy import integer_log
def A379851(n):
def f(x): return n-(m:=integer_log(x, 3)[0])+sum((x//3**i).bit_length() for i in range(m+1))
m, k = n, f(n)
while m != k: m, k = k, f(k)
return m # Chai Wah Wu, Jun 08 2026
CROSSREFS
KEYWORD
nonn,changed
AUTHOR
Franz Vrabec, Jan 04 2025
STATUS
approved
