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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A064438 Numbers which are divisible by the sum of their quaternary digits. 4

%I

%S 1,2,3,4,6,8,9,12,16,18,20,21,24,28,30,32,33,35,36,40,42,48,50,52,54,

%T 60,63,64,66,68,69,72,76,78,80,81,84,88,90,91,96,100,102,108,112,114,

%U 120,126,128,129,132,136,138,140,144,148,150,154,156,160,162,168,171,180

%N Numbers which are divisible by the sum of their quaternary digits.

%C A good "puzzle" sequence -- guess the rule given the first twenty or so terms.

%H Harry J. Smith, <a href="/A064438/b064438.txt">Table of n, a(n) for n = 1..1000</a>

%H Paul Dalenberg, Tom Edgar, <a href="https://www.fq.math.ca/56-2.html">Consecutive factorial base Niven numbers</a>, Fibonacci Quart. (2018) Vol. 56, No. 2, 163-166.

%e Quaternary representation of 28 is 130, 1 + 3 + 0 = 4 divides 28.

%t Select[Range[200],Divisible[#,Total[IntegerDigits[#,4]]]&] (* _Harvey P. Dale_, Jun 09 2011 *)

%o (ARIBAS): maxarg := 190; for n := 1 to maxarg do if n mod sum(quaternarray(n)) = 0 then write(n," "); end; end; function quaternarray(n: integer): array; var k: integer; stk: stack; begin while n > 0 do k := n mod 4; stack_push(stk,k); n := (n - k) div 4; end; return stack2array(stk); end;.

%o (PARI) SumD(x)= { local(s); s=0; while (x>9, s+=x-10*(x\10); x\=10); return(s + x) } baseE(x, b)= { local(d,e,f); e=0; f=1; while (x>0, d=x-b*(x\b); x\=b; e+=d*f; f*=10); return(e) } { n=0; for (m=1, 10^9, if (m%(SumD(baseE(m, 4)))==0, write("b064438.txt", n++, " ", m); if (n==1000, break)) ) } \\ _Harry J. Smith_, Sep 14 2009

%o (PARI) isok(n) = !(n % sumdigits(n, 4)); \\ _Michel Marcus_, Jun 24 2018

%o (Python)

%o from sympy.ntheory.factor_ import digits

%o print [n for n in xrange(1, 201) if n%sum(digits(n, 4)[1:]) == 0] # _Indranil Ghosh_, Apr 24 2017

%Y Cf. A005349 (decimal), A049445 (binary), A064150 (ternary).

%K base,easy,nice,nonn

%O 1,2

%A _Len Smiley_, Oct 01 2001

%E More terms from _Matthew Conroy_, Oct 02 2001

%E Offset changed from 0 to 1 by _Harry J. Smith_, Sep 14 2009

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 20 01:17 EDT 2019. Contains 325168 sequences. (Running on oeis4.)