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

 

Logo


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.

License Agreements, Terms of Use, Privacy Policy. .

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