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

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A064138 Sum of non-unitary divisors of n!. 1
0, 0, 0, 24, 144, 1398, 11184, 97200, 973296, 10950696, 131408352, 1593191808, 22304685312, 333297226080, 5103130001760, 81686161277280, 1470350902991040, 26490792085668288, 529815841713365760, 10635027891469974720 (list; graph; refs; listen; history; internal format)
OFFSET

1,4

LINKS

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

FORMULA

a(n) = sigma(n!)-usigma(n!) = A000203(n!)-A034448(A000142(n)) = A062569(n)-A034448(n!) = A048105(n!)

EXAMPLE

n=6,6!=720, sum of its 30 divisors is 2418, sum of 8 unitary divisors is 1020, so remaining 22 divisors gives a(6)=1398.

PROG

(PARI) usigma(n)= { local(f, s=1); f=factor(n); for(i=1, matsize(f)[1], s*=1 + f[i, 1]^f[i, 2]); return(s) } { n=0; f=1; for (n=1, 100, f*=n; write("b064138.txt", n, " ", sigma(f) - usigma(f)); ) } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Sep 08 2009]

CROSSREFS

Cf. A034448, A003444, A048105, A046656, A056657, A056171, A056172, A000203, A000142, A062569, A063955, A063960.

Sequence in context: A001342 A029572 A129923 * A158047 A042116 A042118

Adjacent sequences:  A064135 A064136 A064137 * A064139 A064140 A064141

KEYWORD

nonn

AUTHOR

Labos E. (labos(AT)ana.sote.hu), Sep 11 2001

EXTENSIONS

Term corrected and more terms added by Harry J. Smith (hjsmithh(AT)sbcglobal.net), Sep 08 2009

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
Recent Additions | More pages | Superseeker | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified February 15 07:58 EST 2012. Contains 205717 sequences.