login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A065968 a(n) = n * Sum_{primes p dividing n} (1 - 1/p). 1
0, 1, 2, 2, 4, 7, 6, 4, 6, 13, 10, 14, 12, 19, 22, 8, 16, 21, 18, 26, 32, 31, 22, 28, 20, 37, 18, 38, 28, 59, 30, 16, 52, 49, 58, 42, 36, 55, 62, 52, 40, 85, 42, 62, 66, 67, 46, 56, 42, 65, 82, 74, 52, 63, 94, 76, 92, 85, 58, 118, 60, 91, 96, 32, 112 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

COMMENTS

For a prime p, a(p) = phi(p), where phi is Euler's totient function. - Alonso del Arte, Jun 16 2012

LINKS

Harry J. Smith, Table of n, a(n) for n = 1..1000

EXAMPLE

a(10) = 13 because 10 is divisible by the primes 2 and 5, and 10*(1/2 + 4/5) = 13.

PROG

(PARI) { for (n=1, 1000, s=0; p=2; while (p<=n, if(n%p == 0, s+=1 - 1/p); p=nextprime(p + 1)); write("b065968.txt", n, " ", n*s) ) } \\ Harry J. Smith, Nov 05 2009

CROSSREFS

Cf. A000010 (same definition except for the crucial difference of product instead of sum for the 1 - 1/p).

Sequence in context: A270227 A162251 A244011 * A105669 A256963 A019657

Adjacent sequences:  A065965 A065966 A065967 * A065969 A065970 A065971

KEYWORD

nonn

AUTHOR

N. J. A. Sloane, Dec 09 2001

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 April 20 16:17 EDT 2019. Contains 322310 sequences. (Running on oeis4.)