%I#42 2021年7月27日15:49:41

%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，

%电话：60,63,64,66,68,69,72,76,78,80,81,84,88,90,91,96100102108112114，

%电话：120126128129132136138140144148150154156160162168171180

%N可被四位数之和整除的数字。

%C一个好的“拼图”序列——根据前二十个左右的术语猜测规则。

%H Harry J.Smith，n表，n=1..1000的a（n）</a>

%H Paul Dalenberg和Tom Edgar，<a href=“https://www.fq.math.ca/56-2.html“>连续阶乘基数Niven数，Fibonacci Quart.（2018）第56卷，第2期，163-166。

%e第四纪代表性28为130，1+3+0=4除以28。

%t选择[Range[200]，Divisible[#，Total[Integer Digits[#，4]]&]（*哈维·P·戴尔，2011年6月9日*）

%o（ARIBAS）：最大参数：=190；对于n:=1到maxarg do，如果n模和（四元数组（n））=0，则写入（n，“”）；结束；结束；函数quaternarray（n：整数）：数组；var k：整数；stk：堆栈；当n>0时开始k:=n mod 4；堆栈推送（stk，k）；n：=（n-k）div 4；结束；返回堆栈2阵列（stk）；结束；。

%o（PARI）

%o和D（x）={局部；s=0；而（x>9，s+=x-10*（x\10）；x\=10）；返回（s+x）}

%o基e（x，b）={局部（d，e，f）；e=0；f=1；while（x>0，d=x-b*（x\b）；x\=b；e+=d*f；f*=10）；返回（e）}

%o{n=0；对于（m=1，10^9，如果（m%（SumD（baseE（m，4）））==0，则写入（“b064438.txt”，n++，“”，m）；如果（n==1000，break））}\\_Harry J.Smith_，2009年9月14日

%o（PARI）isok（n）=！（n%和数字（n，4））；\\_米歇尔·马库斯，2018年6月24日

%o（Python）

%o来自sympy.theory.factor _导入数字

%o打印（如果n%总和（数字（n，4）[1:]）==0]，则[n代表范围（1201）中的n）#_Indranil Ghosh_，2017年4月24日

%Y参考A005349（十进制）、A049445（二进制）、A064150（三进制）。

%K基础，简单，漂亮，不

%O 1,2号机组

%A _Len Smiley，2001年10月1日

%E马修·康罗伊2001年10月2日的更多条款

%2009年9月14日，Harry J.Smith将E偏移从0更改为1