The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A081805 a(n) = n minus (largest prime power in n factorization); a(1) = 0. 3
 0, 0, 0, 0, 0, 3, 0, 0, 0, 5, 0, 8, 0, 7, 10, 0, 0, 9, 0, 15, 14, 11, 0, 16, 0, 13, 0, 21, 0, 25, 0, 0, 22, 17, 28, 27, 0, 19, 26, 32, 0, 35, 0, 33, 36, 23, 0, 32, 0, 25, 34, 39, 0, 27, 44, 48, 38, 29, 0, 55, 0, 31, 54, 0, 52, 55, 0, 51, 46, 63, 0, 63, 0, 37, 50, 57, 66, 65, 0, 64, 0, 41, 0 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,6 COMMENTS a(n) = 0 when n is a prime power (A000961). - Michel Marcus, Dec 03 2013 LINKS Antti Karttunen, Table of n, a(n) for n = 1..65537 FORMULA a(n) = n - A034699(n). - Michel Marcus, Jul 24 2017 MATHEMATICA f[n_] := (n - (#[[1]]^#[[2]]) & /@ {FactorInteger[n][[-1]] })[[1]]; Array[f, 80] (* Robert G. Wilson v, Aug 07 2018 *) PROG (PARI) A081805(n) = if(1==n, 0, my(f = factor(n)); n - vecmax(vector(#f~, k, f[k, 1]^f[k, 2]))); \\ Michel Marcus Jul 24 2017 & Antti Karttunen, Aug 06 2018 CROSSREFS Cf. A034699, A052128, A284600. Sequence in context: A195938 A184762 A330734 * A181189 A327889 A221702 Adjacent sequences:  A081802 A081803 A081804 * A081806 A081807 A081808 KEYWORD nonn AUTHOR Benoit Cloitre, Apr 10 2003 EXTENSIONS Term a(1) = 0 prepended by Antti Karttunen, Aug 06 2018 STATUS approved

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

Last modified January 24 19:19 EST 2020. Contains 331211 sequences. (Running on oeis4.)