Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A081705
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A081705 k-tuple abundance of abundant numbers. +0
4
1, 1, 1, 2, 7, 1, 1, 6, 1, 5, 1, 2, 4, 1, 1, 3, 1, 3, 1, 2, 2, 1, 7, 1, 1, 1, 6, 8, 5, 3, 31, 2, 1, 30, 1, 1, 1, 28, 5, 2, 14, 4, 1, 1, 3, 1, 2, 1, 14, 4, 1, 29, 4, 1, 28, 7, 4, 5, 3, 1, 11, 2, 1, 6, 3, 12, 1, 11, 1, 1, 6, 5, 27, 18, 1, 1, 17, 1, 2, 3, 3, 1, 1, 14, 4, 4, 13, 1, 1, 12, 2, 3, 10, 1, 5, 1, 4 (list; graph; listen)
OFFSET

1,4

COMMENT

Note that only increasing steps at the beginning of an aliquot chain count toward k-tuple abundance. I wonder if there are an infinite number of k-tuply abundant numbers for all k? Another interesting question - are there any numbers that are completely abundant (that is, numbers whose aliquot chain increases forever). Though there are several numbers whose aliquot chains aren't yet fully determined, all the ones I've checked have had a finite k-tuple abundance.

Lenstra shows that there are in fact infinitely many k-tuply abundant numbers for every k > 0.

REFERENCES

H. W. Lenstra, Jr., "Advanced Problems and Solutions, 6064", The American Mathematical Monthly, Vol. 84, No. 7. (Aug. - Sep., 1977), p. 580.

FORMULA

a(n) = 0 if n is not abundant, otherwise 1 + (a(sigma(n)-n)) Note, however, that non-abundant numbers are excluded from this sequence.

a(n) = number of increasing steps at the start of the aliquot chain of A005101(n).

EXAMPLE

a(4)=2 because the 4th abundant number is 24 which has aliquot sequence 24->36->55->17->1, which has two increasing steps at the beginning.

MAPLE

A005101 := proc(n) option remember ; local a ; if n =1 then RETURN(12) ; else a := A005101(n-1)+1 ; while numtheory[sigma](a) <= 2*a do a := a+1 ; od ; RETURN(a) ; fi ; end: aliqRis := proc(n) local r, a, an ; r := 0 ; a := n; while true do an := numtheory[sigma](a)-a ; if an > a then r := r+1 ; a := an ; else RETURN(r) ; fi ; od ; end: A081705 := proc(n) aliqRis(A005101(n)) ; end: for n from 1 to 100 do printf("%d, ", A081705(n)) ; od ; - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Mar 07 2007

CROSSREFS

Cf. A081699, A081700.

Sequence in context: A084794 A116691 A011340 this_sequence A145057 A065254 A010592

Adjacent sequences: A081702 A081703 A081704 this_sequence A081706 A081707 A081708

KEYWORD

nonn

AUTHOR

Gabriel Cunningham (gcasey(AT)mit.edu), Apr 02 2003, Dec 15 2006

EXTENSIONS

More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Mar 07 2007

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified December 4 15:11 EST 2009. Contains 170347 sequences.


AT&T Labs Research