cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A114440 Numbers which divided by the sum of their digits (Harshad or Niven numbers) give integers which are also divisible by the sum of their digits (until a single-digit Harshad remains).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 12, 18, 21, 24, 27, 36, 42, 45, 48, 54, 63, 72, 81, 84, 108, 162, 216, 243, 324, 378, 405, 432, 486, 648, 756, 864, 972, 1296, 1458, 1944, 2916, 3402, 4374, 5832, 6804, 7290, 8748, 11664, 13122, 13608, 15552, 17496, 23328, 26244
Offset: 1

Views

Author

Piotr K. Olszewski (piotrkornelolszewski(AT)poczta.onet.pl), Feb 14 2006

Keywords

Comments

The sequence is finite with a(15095), a 1434-digit number, being the final term. - Hans Havermann and Ray Chandler, Jan 21 2014

Examples

			The number 216 is a term of the sequence because it is divisible by the sum of its digits: 2+1+6=9; 216/9=24. Also, the successive quotients are divisible by the sum of their digits, until a single-digit Harshad remains: 24: 2+4=6; 24/6=4 and 4: 4/4=1.
		

Crossrefs

Programs

  • Mathematica
    s=w={1}; Do[t={}; Do[v=s[[k]]; u={}; Do[If[Total[IntegerDigits[c*v]]==c, AppendTo[u,c*v]], {c,2,7000}]; t=Join[t,u], {k,Length[s]}]; s=Sort[t]; w=Join[w,s], {440}]; Union[w] (* Hans Havermann, Jan 21 2014 *)
  • PARI
    v=vector(118); for(n=1, 9, v[n]=n; print1(n ", ")); c=9; for(n=10, 10^9, d=length(Str(n)); m=n; s=0; for(j=1, d, s=s+m%10; m=m\10); if(s==1, next); if(n%s==0, m=n/s, next); forstep(j=c, 1, -1, if(v[j]<=m, if(v[j]==m, c++; v[c]=n; print1(n ", ")); next(2)))) /* Donovan Johnson, Apr 09 2013 */

Extensions

Offset corrected by Donovan Johnson, Apr 09 2013
a(54)-a(235) from Donovan Johnson, Apr 09 2013
a(236)-a(15095) from Hans Havermann and Ray Chandler, Jan 21 2014

A235601 Smallest number m such that repeated application of A235600 takes n steps to reach 1, where A235600(k) = k/A007953(k) if the digital sum A007953(k) divides k, A235600(k) = k otherwise.

Original entry on oeis.org

1, 2, 12, 108, 1944, 52488, 1102248, 44641044, 2008846980, 108477736920, 6508664215200, 421761441144960, 22142475660110400, 1793540528468942400, 160701231350817239040, 15909421903730906664960, 1874419162475932276162560
Offset: 0

Views

Author

N. J. A. Sloane and David W. Wilson, Jan 18 2014

Keywords

Comments

Numbers m > 1 which never reach 1 are not candidates for a(n).
There is no analog in base 2 (cf. A235602).
Comment from David W. Wilson, Jan 20 2013: let S(0) = {1}; for each n >= 1, compute the set S(n) of possible predecessors of elements of S(n-1). Then a(n) is the smallest element of S(n). Using this approach, I was able to compute up to a(100).
The sequence is finite with a(440), a 1434-digit number being the final term. - Hans Havermann and Ray Chandler, Jan 21 2014
Sequence A236338 gives the count of iterations of A235600 required to reach 1 when starting from any n. Otherwise said: This sequence is the RECORDS transform of A236338. - M. F. Hasler, Jan 22 2014
The terms are a proper subset of A114440. - Robert G. Wilson v, Jan 22 2014

Examples

			a(4) = 1944: 1944 ->1944/18 = 108 -> 108/9 = 12 -> 12/3 = 4 -> 4/4 = 1 in 4 steps.
		

Crossrefs

Programs

  • Mathematica
    s={1}; Print[s[[1]]]; Do[t={}; Do[v=s[[k]]; u={}; Do[If[Total[IntegerDigits[c*v]]==c, AppendTo[u,c*v]], {c,2,7000}]; t=Join[t,u], {k,Length[s]}]; s=Sort[t]; Print[s[[1]]], {440}] (* Hans Havermann, Jan 21 2014 *)

Extensions

a(8) from Hans Havermann, Jan 19 2014
a(9)-a(100) from David W. Wilson, Jan 21 2014
a(101)-a(440) from Hans Havermann and Ray Chandler, Jan 21 2014
Showing 1-2 of 2 results.