%I#19 2021年3月1日06:08:03
%S 1,1,0,-1,-1,0,1,0,-1,2,1,-3,-3,4,7,-4,-14,3,28,4,-53,-30,91,95,
%电话-142,-240191556,-165,-1192,-14023581215,-4322,-4194721111492,
%U-10300,-2783910194617703506,-127099,-58080242071218346,-419278
%N在逆Euler变换下向左移动。
%H Alois P.Heinz,n的表格,n=1..1000时的a(n)</a>
%H N.J.A.Sloane,转换</a>
%p with(numtheory):ietr:=proc(p)local a,c;c: =proc(n)选项记忆;局部j;n*p(n)-添加(c(j)*p(n-j),j=1..n-1)结束;a: =proc(n)选项记忆;局部d`如果`(n=0,1,add(mobius(n/d)*c(d),d=除数(n))/n)end-end:aa:=ietr('a'):a:=n->`if`(n<=1,1,aa(n-1)):seq(a(n),n=1.48);#_Alois P.Heinz,2008年9月8日
%t ietr[p_]:=模块[{a,c},
%tc[n]:=c[n]=n*p[n]-和[c[j]*p[n-j],{j,1,n-1}];
%t a[n_]:=a[n]=如果[n==0,1,总和[MoebiusMu[n/d]*c[d],
%t{d,除数[n]}]/n];a] ;
%t aa=ietr[a];
%t a[n]:=如果[n<=1,1,aa[n-1]];
%t数组[a,48](*_Jean-François Alcover_,2021年3月1日,在_Alois P.Heinz_*之后)
%K符号,特征
%O 1,11号
%A _克里斯蒂安·G·鲍尔,1999年1月4日
|