%I #22 Jul 10 2022 10:47:41
%S 7,10,13,19,22,25,28,34,37,40,49,52,67,70,73,76,82,85,88,97,100,112,
%T 130,133,136,145,148,160,193,196,208,259,262,265,268,274,277,280,289,
%U 292,304,322,325,328,337,340,352,385,388,400,448,514,517,520,529,532
%N Numbers whose base-4 sum of digits is 4.
%C Subsequence of A016777. - _Michel Marcus_, Sep 03 2013
%C In general, the set of numbers with sum of base-b digits equal to b is a subset of { (b-1)*k + 1; k = 2, 3, 4, ... }. - _M. F. Hasler_, Dec 23 2016
%H Michael S. Branicky, <a href="/A226969/b226969.txt">Table of n, a(n) for n = 1..10000</a>
%e The quaternary expansion of 13 is (3,1), which has sum of digits 4.
%e The quaternary expansion of 40 is (2,2,0), which has sum of digits 4.
%e 17 is not on the list since the quaternary expansion of 17 is (1,0,1), which has sum of digits 2 not 4.
%t Select[Range@ 540, Total@ IntegerDigits[#, 4] == 4 &] (* _Michael De Vlieger_, Dec 23 2016 *)
%o (Sage) [i for i in [0..1000] if sum(Integer(i).digits(base=4))==4]
%o (PARI) select( is(n)=sumdigits(n,4)==4, [1..999]) \\ _M. F. Hasler_, Dec 23 2016
%o (Python)
%o agen = A226636gen(sod=4, base=4) # generator of terms using code in A226636
%o print([next(agen) for n in range(1, 57)]) # _Michael S. Branicky_, Jul 10 2022
%Y Cf. A018900, A187813.
%Y Cf. A226636 (b = 3), A227062 (b = 5), A227080 (b = 6), A227092 (b = 7), A227095 (b = 8), A227238 (b = 9), A052224 (b = 10).
%K nonn,base
%O 1,1
%A _Tom Edgar_, Sep 01 2013