The OEIS is supported by the many generous donors to the OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A351168 If n = p_1^e_1 * ... * p_k^e_k, where p_1 < ... < p_k are primes, then a(n) is obtained by replacing the last factor p_k^e_k by (p_k - 1)^e_k; a(1) = 1. 4
 1, 1, 2, 1, 4, 4, 6, 1, 4, 8, 10, 8, 12, 12, 12, 1, 16, 8, 18, 16, 18, 20, 22, 16, 16, 24, 8, 24, 28, 24, 30, 1, 30, 32, 30, 16, 36, 36, 36, 32, 40, 36, 42, 40, 36, 44, 46, 32, 36, 32, 48, 48, 52, 16, 50, 48, 54, 56, 58, 48, 60, 60, 54, 1, 60, 60, 66, 64, 66, 60, 70, 32, 72, 72, 48 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,3 COMMENTS First time a term appears four or more times in a row is when n = 1684. LINKS Amiram Eldar, Table of n, a(n) for n = 1..10000 Jan Misali, What should we call the other bases?, Web page, no date. Inspiration for this sequence. FORMULA a(n) = n*(1 - 1/p_k)^e_k where prime factorization n = p_1^e_1 * ... * p_k^e_k with ascending p_1 < ... < p_k. a(n) = n*(1 - 1/A006530(n))^A071178(n). EXAMPLE The prime factorization of 44 is 2^2 * 11^1, so a(44) = 2^2 * 10^1 = 40. The prime factorization of 50 is 2^1 * 5^2, so a(50) = 2^1 * 4^2 = 32. MATHEMATICA a[n_] := Module[{f = FactorInteger[n]}, n*(1 - 1/f[[-1, 1]])^f[[-1, 2]]]; a = 1; Array[a, 100] (* Amiram Eldar, Feb 04 2022 *) PROG (PARI) a(n) = my(f=factor(n), r=matsize(f)); if(r, f[r, 1]--); factorback(f); \\ Kevin Ryde, Feb 03 2022 CROSSREFS Cf. A006530 (greatest prime), A071178 (its exponent). Cf. A171462 (one instance of the decrement), A003958 (all primes decremented), A351419, A351425. Sequence in context: A193631 A190993 A326147 * A196082 A273724 A108755 Adjacent sequences: A351165 A351166 A351167 * A351169 A351170 A351171 KEYWORD nonn AUTHOR Ben Polson, Feb 03 2022 EXTENSIONS a(1) = 1 prepended by Michel Marcus, Feb 04 2022 Edited by N. J. A. Sloane, Feb 11 2022 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified March 31 21:40 EDT 2023. Contains 361673 sequences. (Running on oeis4.)