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!)
A063919 Sum of proper unitary divisors (or unitary aliquot parts) of n, including 1. 31
1, 1, 1, 1, 1, 6, 1, 1, 1, 8, 1, 8, 1, 10, 9, 1, 1, 12, 1, 10, 11, 14, 1, 12, 1, 16, 1, 12, 1, 42, 1, 1, 15, 20, 13, 14, 1, 22, 17, 14, 1, 54, 1, 16, 15, 26, 1, 20, 1, 28, 21, 18, 1, 30, 17, 16, 23, 32, 1, 60, 1, 34, 17, 1, 19, 78, 1, 22, 27, 74, 1, 18, 1, 40, 29, 24, 19, 90, 1, 22, 1, 44 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,6

COMMENTS

For definition of unitary divisor see A034448.

LINKS

Antti Karttunen, Table of n, a(n) for n = 1..65537 (first 1000 terms from Harry J. Smith)

FORMULA

a(n) = A034460(n), n>1. - R. J. Mathar, Oct 02 2008

For n > 1: a(n) = sum (A077610(n,k): k = 1 .. A034444(n) - 1). - Reinhard Zumkeller, Mar 12 2012

EXAMPLE

a(10) = 8 because the unitary divisors of 10 are 1, 2, 5 and 10, with sum 18 and 18-10 = 8.

MAPLE

A063919 := proc(n)

    if n = 1 then

        1;

    else

        A034448(n)-n ;

    end if;

end proc: # R. J. Mathar, May 14 2013

MATHEMATICA

a[n_] := Total[Select[Divisors[n], GCD[#, n/#] == 1&]]-n; a[1] = 1; Table[a[n], {n, 82}] (* Jean-Fran├žois Alcover, Aug 31 2011 *)

PROG

(PARI) usigma(n) = sumdiv(n, d, if(gcd(d, n/d)==1, d))

{ for (n=1, 1000, if (n>1, a=usigma(n) - n, a=1); write("b063919.txt", n, " ", a) ) } \\ Harry J. Smith, Sep 02 2009

(PARI)

A034460(n) = (sumdivmult(n, d, if(gcd(d, n/d)==1, d))-n); \\ From A034460

A063919(n) = if(1==n, n, A034460(n)); \\ Antti Karttunen, Jun 12 2018

(Haskell)

a063919 1 = 1

a063919 n = sum $ init $ a077610_row n

-- Reinhard Zumkeller, Mar 12 2012

CROSSREFS

The values of sequence are A034448(n)-n (for n > 1).

Cf. A001065, A034448, A034460.

Sequence in context: A275109 A320832 A034460 * A308135 A072815 A080304

Adjacent sequences:  A063916 A063917 A063918 * A063920 A063921 A063922

KEYWORD

easy,nonn,nice

AUTHOR

Felice Russo, Aug 31 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 January 29 04:46 EST 2020. Contains 331335 sequences. (Running on oeis4.)