q=x^3;s=x^2+x+1;z=40;
p[n,x_]:=(2x+1)^n;
表[展开[p[n,x]],{n,0,7}]
减少[{p1,q,s,x}]:=
固定点[多项式商@@#1+
多项式余数@@#1&)[{#1,q,x}]&,p1]
t=表[reduce[{p[n,x],q,s,x}],{n,0,z}];
u1=表[系数[部分[t,n],x,0],{n,1,z}](*A192814号*)
u2=表[系数[部分[t,n],x,1],{n,1,z}](*A192815号*)
线性递归[{5,-3,7},{1,1,1},30](*G.C.格鲁贝尔,2019年1月3日*)