q[x_]:=x+1;d=2;
p[n,x_]:=((x+d)^n+(x-d)^n)/2(*类似于在A161516号*)
表[展开[p[n,x]],{n,0,6}]
约简规则={x^y_?EvenQ->q[x]^(y/2),
x^y_?奇数q->xq[x]^((y-1)/2)};
t=表[Last[Most[FixedPointList[Expand[#1/.reductionRules]&,p[n,x]]],{n,0,30}]
表[系数[部分[t,n],x,1],{n,1,30}](*1923年*)