(PARI)an=矢量(200);对于(n=1,4,an[n]=[1,1,1,1][n]);对于(n=5,长度(an),an[n]=(an[n-1]*an[n-3]-an[n-2]^2)/an[n-4]);a(n)=符号(n)*an[abs(n)+(n==0)]
(PARI){a(n)=my(an);如果(n<0,-a(-n),如果(n==0,0,an=向量(max(4,n),i,1);an[4]=-1;对于(k=5,n,an[k]=(an[k-1]*an[k-3]-an[k-2]^2)/an[k-4]);an[n])}/*迈克尔·索莫斯2014年7月7日*/
(哈斯克尔)
a050512 n=a050512_list!!n个
a050512_list=0:1:1:(-1):zipWith div(zipWise(-)(zipWith(*))
(删除4 a050512_list)(删除2 a050512-list)
(地图(^2)(放置3 a050512_list))(尾部a050512-list)
|