OFFSET
1,4
COMMENTS
Fixed points of the map x -> gcd(x, A003415(x)), i.e., if we start iterating with A085731 from any x = n (>= 1), we will eventually reach a(n), after which the result does not change anymore. This was found by LODA miner (see C. Krause link), and is easily seen to be true by Eric M. Schmidt's multiplicative formula for A085731. Note also that this sequence is idempotent, meaning a(a(n)) = a(n) for all n. - Antti Karttunen, Apr 05 2021
The largest divisor of n that is a term of A072873. - Amiram Eldar, Sep 14 2023
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
FORMULA
Multiplicative with a(p^e) = p^(e-(e mod p)).
a(n) = n / A327938(n).
MATHEMATICA
f[p_, e_] := p^(e - Mod[e, p]); a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 14 2023 *)
PROG
(PARI) A327939(n) = { my(f = factor(n)); for(k=1, #f~, f[k, 2] = (f[k, 2]-(f[k, 2]%f[k, 1]))); factorback(f); };
CROSSREFS
KEYWORD
nonn,easy,mult
AUTHOR
Antti Karttunen, Oct 01 2019
STATUS
approved