login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A144776 Define f(n) = 1 if n is prime, 2 * rad(n) if four divides n and rad(n) otherwise: then a(n) = 0 for composite n where f(n) is not less than n and otherwise equals the number of positive integers k less than n for which f(k) < f(n). 1
0, 0, 0, 0, 0, 0, 0, 5, 5, 0, 0, 0, 0, 0, 0, 8, 0, 12, 0, 0, 0, 0, 0, 17, 14, 0, 10, 0, 0, 0, 0, 14, 0, 0, 0, 22, 0, 0, 0, 28, 0, 0, 0, 0, 29, 0, 0, 26, 25, 26, 0, 0, 0, 24, 0, 42, 0, 0, 0, 0, 0, 0, 41, 21, 0, 0, 0, 0, 0, 0, 0, 35, 0, 0, 42, 0, 0, 0, 0, 46, 23, 0, 0, 0, 0, 0, 0, 64, 0, 58, 0, 0, 0, 0, 0 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,8
COMMENTS
For the given terms, nearly all n for which a(n) obtains a new maximum are multiples of eight. Only 18, 36 and 45 are not.
LINKS
EXAMPLE
f(8) = 4 and f(9) = 3. For 1, 2, 3, 5 and 7, f(k) = 1, so a(8) = a(9) = 5.
PROG
(PARI)
A007947(n) = factorback(factorint(n)[, 1]);
f(n) = if(isprime(n), 1, if(n%4, A007947(n), 2*A007947(n)));
A144776(n) = if(n<2, 0, my(x=f(n)); if(!isprime(n)&&(x>=n), 0, sum(k=1, n-1, (f(k)<x)))); \\ Antti Karttunen, Jul 03 2018
CROSSREFS
Sequence in context: A186639 A266668 A043299 * A371264 A065937 A197738
KEYWORD
easy,nonn
AUTHOR
Reikku Kulon, Sep 21 2008
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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 02:23 EDT 2024. Contains 371264 sequences. (Running on oeis4.)