%I#20 2018年6月16日19:09:18

%S 0,0,0,1,0,0,1，-3,0,1，-2,3,1，-11,0,0,1，-6,0,1，-5,12,1，-4,7，-9，-3,3，

%T-2，-2,0,1，-4,12,3,8，-15，-2,5，-7,7,3，-2,0，4,1，-2.4，-7，-1,2，-3,4,1，

%U-1,1，-1,0,0,1,0,1,1,0,1，-9,0,1，-8,21,1，-7,13，-18

%N四次剖分a（4n+k）给出了k=0..3时a的第k个差值，N<3时a（N）=0，a（3）=1。

%H Alois P.Heinz，n的表格，n=0..65536的a（n）</a>

%F a（4*n）=a（n），

%F a（4*n+1）=a（n+1）-a（n），

%Fa（4*n+2）=a（n+2）-2*a（n+1）+a（n），

%F a（4*n+3）=a（n+3”）-3*a（n+2）+3*a（n+1）-a（n）。

%p a:=proc（n）选项记忆；（m->`if`（n<4，`if`）（n=3,1,0），加(

%pa（q+m-j）*（-1）^j*二项式（m，j），j=0..m））（irem（n，4，'q'））

%p端：

%p序列（a（n），n=0..100）；

%ta[n_]：=a[n]=模[{m，q}，{q，m}=商余数[n，4]；如果[n<4，如果[n==3，1，0]，和[a[q+m-j]*（-1）^j*二项式[m，j]，{j，0，m}]]；

%t表[a[n]，{n，0，100}]（*Jean-François Alcover_，2018年6月9日，摘自Maple*）

%Y参见A005590、A229653、A22965、A22965%、A2296506、A2296.57、A229685、A229699、A229660。

%K符号，特征，看

%0、8

%A _Alois P.Heinz，2013年9月27日