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

 

Logo

"Email this user" was broken Aug 14 to 9am Aug 16. If you sent someone a message in this period, please send it again.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A098771 Numbers n such that sigma(n)=sigma(d_1)*sigma(d_2)*...*sigma(d_k) where d_1 d_2 ... d_k is the decimal expansion of n. 2
1, 2, 3, 4, 5, 6, 7, 8, 9, 38, 58, 66, 87, 568, 766, 2799, 6478, 6785, 26486, 36685, 36893, 64788, 66758, 67586, 69664, 77686, 276895, 467683, 474668, 488973, 548678, 587667, 658648, 666369, 668954, 839666, 939669, 958968, 968765, 2386698 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Number of terms of this sequence up to 15000000 is 72. Has this sequence at least one multi-digit prime term?

LINKS

Table of n, a(n) for n=1..40.

EXAMPLE

97888745 is in the sequence because sigma(97888745) = 117936000 = sigma(9)*sigma(7)*sigma(8)*sigma(8)*sigma(8)*sigma(7)*sigma(4)*sigma(5).

MATHEMATICA

Do[If[h=IntegerDigits[n]; l=Length[h]; !MemberQ[h, 0]&&DivisorSigma[ 1, n]==Product[DivisorSigma[1, h[[k]]], {k, l}], Print[n]] {n, 100000000}]

CROSSREFS

Cf. A098684, A098685.

Sequence in context: A257554 A270393 A257787 * A276810 A024659 A129513

Adjacent sequences:  A098768 A098769 A098770 * A098772 A098773 A098774

KEYWORD

base,nonn

AUTHOR

Farideh Firoozbakht, Oct 03 2004

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 | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy .

Last modified August 19 01:41 EDT 2017. Contains 290789 sequences.