login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A038186 Numbers divisible by the sum and product of their digits. 15
1, 2, 3, 4, 5, 6, 7, 8, 9, 12, 24, 36, 111, 112, 132, 135, 144, 216, 224, 312, 315, 432, 612, 624, 735, 1116, 1212, 1296, 1332, 1344, 1416, 2112, 2232, 2916, 3132, 3168, 3276, 3312, 4112, 4224, 6624, 6912, 8112, 9612, 11112, 11115, 11133, 11172, 11232 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
The property "numbers divisible by the sum and product of their digits" leads to the Diophantine equation t*x1*x2*...*xr=s*(x1+x2+...+xr), where t and s are divisors of n; xi is from [1...9]. This corresponds to some arithmetic problems in geometry, see Sándor, 2002. - Ctibor O. Zizka, Mar 04 2008
LINKS
David A. Corneth, Table of n, a(n) for n = 1..10352 (first 1000 terms from T. D. Noe)
József Sándor, Geometric Theorems, Diophantine Equations and Arithmetic Functions, American Research Press, Rehoboth, 2002.
FORMULA
A188641(a(n)) * A188642(a(n)) = 1. - Reinhard Zumkeller, Apr 07 2011
MATHEMATICA
dspQ[n_]:=Module[{idn=IntegerDigits[n], t}, t=Times@@idn; t!=0 && Divisible[n, Total[idn]] && Divisible[n, t]]; Select[Range[11500], dspQ] (* Harvey P. Dale, Jul 11 2011 *)
PROG
(Haskell)
import Data.List (elemIndices)
a038186 n = a038186_list !! (n-1)
a038186_list = map succ $ elemIndices 1
$ zipWith (*) (map a188641 [1..]) (map a188642 [1..])
-- Reinhard Zumkeller, Apr 07 2011
(PARI) for(n=1, 10^4, d=digits(n); s=sumdigits(n); p=prod(i=1, #d, d[i]); if(p&&!(n%s+n%p), print1(n, ", "))) \\ Derek Orr, Apr 29 2015
(Python)
from math import prod
def sd(n): return sum(map(int, str(n)))
def pd(n): return prod(map(int, str(n)))
def ok(n): return n%sd(n) == 0 and pd(n) and n%pd(n) == 0
def aupto(limit): return [m for m in range(1, limit+1) if ok(m)]
print(aupto(11233)) # Michael S. Branicky, Jan 28 2021
CROSSREFS
Intersection of A005349 and A007602. - Reinhard Zumkeller, Apr 07 2011
Sequence in context: A001102 A051004 A032575 * A118575 A327453 A289791
KEYWORD
nonn,base,nice,look
AUTHOR
EXTENSIONS
More terms from Patrick De Geest, Jun 15 1999
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 01:35 EDT 2024. Contains 371964 sequences. (Running on oeis4.)