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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A098189 Sum of unitary divisors minus Euler phi: a(n) = A034448(n) - A000010(n). 7
0, 2, 2, 3, 2, 10, 2, 5, 4, 14, 2, 16, 2, 18, 16, 9, 2, 24, 2, 22, 20, 26, 2, 28, 6, 30, 10, 28, 2, 64, 2, 17, 28, 38, 24, 38, 2, 42, 32, 38, 2, 84, 2, 40, 36, 50, 2, 52, 8, 58, 40, 46, 2, 66, 32, 48, 44, 62, 2, 104, 2, 66, 44, 33, 36, 124, 2, 58, 52, 120, 2, 66, 2, 78, 64, 64, 36, 144, 2 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

LINKS

Michael De Vlieger, Table of n, a(n) for n = 1..10000

FORMULA

a(n)>A063919(n) if n>1.

a(A000040(k))=2.

EXAMPLE

a(1)=1-1=0.

MATHEMATICA

Table[DivisorSum[n, # &, CoprimeQ[#, n/#] &] - EulerPhi@ n, {n, 120}] (* Michael De Vlieger, Mar 01 2017 *)

PROG

(PARI) a(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d)) - eulerphi(n); \\ Michel Marcus, Feb 25 2014

(PARI) a(n)=my(f=factor(n)); prod(k=1, #f[, 2], f[k, 1]^f[k, 2]+1) - eulerphi(f) \\ Charles R Greathouse IV, Mar 01 2017

CROSSREFS

Cf. A034448, A000010, A063919, A098190-A098195.

Sequence in context: A179592 A058671 A016002 * A162733 A016003 A087139

Adjacent sequences:  A098186 A098187 A098188 * A098190 A098191 A098192

KEYWORD

nonn

AUTHOR

Labos Elemer, Sep 03 2004

EXTENSIONS

Edited by R. J. Mathar, Mar 02 2009

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 February 19 18:22 EST 2019. Contains 320327 sequences. (Running on oeis4.)