########################################################################SVT.txt:将此文件另存为SVT.txt###要使用它,请留在###同一目录,进入Maple(输入:Maple)###然后键入:read`SVT.txt`###然后按照那里给出的说明进行操作### ###罗格斯大学Doron Zeilberger撰稿##zeilberg在数学dot rutgers dot edu########################################################################创建日期:2016年4月26日#本版本:2016年5月13日带有(linalg):打印(`创建日期:2016年4月26日`):打印(`本版本:2016年5月13日`):打印(`这是SVT.txt`):打印(“文章附带”):print(`On the Sequences Enumerating Singular Vector Tuples of d-dimension n by n…by n tensors)(关于用n张量枚举d维n的奇异向量元组的序列):印刷品(由Shalosh B.Ekhad和Doron Zeilberger创作):打印(“也可从Zeilberger的网站上获得”):打印(``):打印(`Please report bug to zeilberg at math dot rutgers dot edu`):打印(``):打印(“此软件包和纸张的最新版本”):打印(“可从以下位置获得”):打印(`http://www.math.rutgers.edu/~zeilberg/.`):打印(`--------------------------------------------`):print(`有关支持过程的列表,请键入ezra1();,有关“)的帮助:print(`特定过程,键入ezra(procedure_name);.`):打印(``):打印(`----------------------------------------`):打印(`----------------------------------------`):print(`有关MAIN过程的列表,请键入ezra();,有关“)的帮助:print(`特定过程,键入ezra(procedure_name);.`):打印(``):打印(`----------------------------------------`):使用(组合):数字:=20:ezra1:=处理()如果args=NULL,则print(`支持程序为:Asy、AsyC、A3f、B3f、A4f、B4f、CharRoots、Cnk、EmpirAsyC,findrec、findrec、MMMT,`):打印(`SVTseqSlow,SeqFromRec,Zinn`):打印(``):其他的ezra(参数):图1:结束时间:ezra:=进程()如果args=NULL,则print(`主要程序是:Adf、Bdf、Bdr、SVT、SVTseq、SVTsq3、SVTeq4、SVTseqPC、SVTpaper`):打印(``):elif nops([args])=1,op(1,[args')=A3f,则print(`A3f(a1,a2,a3):0≤b1≤a1,0≤b2≤a2,0≤B3≤a3`的B3(b1,b2,B3)之和:打印(`Try:`):打印(`A3f(3,3,3);`):elif nops([args])=1,op(1,[args')=A4f,则print(`A4f(a1,a2,a3,a4):0≤b1≤a1,0≤b2≤a2,0≤b3≤a3,0≤B4≤a4`)的B4(b1,b2,b3,B4)之和:打印(`Try:`):打印(`A4f(3,3,3,1);`):elif nops([args])=1,op(1,[args')=Adf,则print(`Adf(a):c(a),其中a是长度为k的非负整数列表,例如,输出x1^a[1]*的系数。。。xd^a[d]`):print(`在有理函数1/(1-e_2(x1,…,xd)-2*e_3(x1、…、xd)-(d-1)e_d(x1…,xt))*(1/((1-x1)**(1-xd))`):打印(`Try:`):打印(`Adf([3,3,3]);`):elif nops([args])=1,op(1,[args')=Asy,则print(`Asy(ope,n,n,mu,M):输入多项式系数的线性递归算子,`):print(`在变量n中,向前移位运算符用n表示,即Nf(n):=f(n+1)`):print(`及其特征根之一mu和正整数M找到渐近级数`):print(`annihating the requency,to order M`):打印(“不考虑初始条件”):打印(`Try:`):打印(`Asy(n^2*(n-1)*(n-2)*(n-3)+n*(n^2+5*n+1)+(n^2+7*n+1,n,n,3,4);`):elif nops([args])=1,op(1,[args')=AsyC,则print(`AsyC(ope,n,n,mu,M,INI,K):输入多项式系数的线性递归算子,`):print(`在变量n中,向前移位运算符用n表示,即Nf(n):=f(n+1)`):print(`及其特征根之一mu和正整数M以及列表`):print(`数字的INI,其长度是递归的顺序(即次数(ope,N)),`):print(`它还输入一个大整数K,告诉我们要估计/猜测前面的常数需要走多远`):print(`它输出渐近序列,包括前面的常数`):print(`annilihating the returnal,to order M,and saving the given initial conditions`):print(`如果mu错误,则返回FAIL `)。print(`除估计的常数外,其他一切都很严格。`):打印(`Try:`):打印(`AsyC(n ^2*(n-1)*(n-2)*(n-3)+n*(n ^2+5*n+1)+(n ^ 2+7*n+1),n,n,3,4,[1,1],1000);`):elif nops([args])=1,op(1,[args')=B3f,则打印(`B3f(a1,a2,a3):x1^a1*x2^a2*x3^a3在1/(1-(x1*x2+x1*x3+x2*x3)-2*x1*x2*x3`)中的系数:打印(`Try:`):打印(`B3f(3,4,6);`):elif nops([args])=1,op(1,[args')=B4f,则打印(`B4f(a1,a2,a3,a4):x1^a1*x2^a2*x3^a3*x4^a4在1/(1-e2(x1,x2,x3,x4)-2*e3(x1、x2、x3、x4)-3*e4打印(`Try:`):打印(`B4f(3,3,3,1);`):elif nops([args])=1,op(1,[args')=Bdf,则print(`Bdf(a):f(a),其中a是长度为d的非负整数列表,例如,输出x1^a[1]*的系数。。。xd^a[d]`):print(`在有理函数1/(1-e_2(x1,…,xd)-2*e_3(x1、…、xd)-(d-1)e_d打印(`Try:`):打印(`Bdf([3,3,3]);`):elif nops([args])=1,op(1,[args')=Bdr,则print(`Bdr(a):f(a),其中a是长度为d的非负整数列表,而正整数r输出系数x1^a[1]*。。。xd^a[d]`):print(`在有理函数1/(1-e_2(x1,…,xd)-2*e_3(x1,…,xd)-(d-1)e_d(x1,…,xd))^r`中:打印(`Try:`):打印(`Bdr([3,3,3],2);`):elif nops([args])=1,op(1,[args')=CharRoots,然后print(`(ope,n,n):输入具有多项式系数的线性递归运算符,ope`):print(`在变量n中,向前移位运算符用n表示,即Nf(n):=f(n+1)`):打印(`并输出有序列表`):print(其特征根的`。如果它有复杂根,则返回FAIL。尝试:a`):打印(`字符根((n+4)*n^2-n*n-1,n,n);`):打印(``):elif nargs=1且args[1]=Cnk,则print(`Cnk(n,k):长度为n且k为1的所有0-1向量`):elif nargs=1且args[1]=EmpirAsyC,则print(`EmpirAsyC(L1,n,mu,t,M):输入一个正数列表L1,该列表被认为给出了形式的渐近性`):print(`mu^n*n^t,使用ofer M的sympotics查找常量。尝试:`):打印(`EmpirAsyC([seq(5*3^i/i*(1+4/i),i=1..40)],n,3,-1,5);`):elif nargs=1且args[1]=findrec,则print(`findrec(f,DEGREE,ORDER,n,n):猜测递归运算符的湮灭`):print(`degree和order order.`的顺序f):print(`例如,try:findrec([seq(i,i=1..10)],0,2,n,n);`):elif nargs=1和args[1]=Findrec,然后print(`Findrec(f,n,n):给定一个列表f,尝试用`求线性递归方程:print(`poly-coffs.ope(n,n),其中n是离散变量,n是移位运算符`):print(`消除序列f。如果失败,则返回FAIL。这意味着要么没有重复出现`):print(`或其复杂性使您需要更多术语。它不处理常系数递归的解决方案`):打印(`例如,尝试Findrec([seq(i,i=1..20)],n,n);`):elif nargs=1且args[1]=MMMT,则print(`MMMT(A,x):变量x列表中MacMahon主定理给出的生成函数。`):打印(`Try:`):打印(`MMMT([[0,1,1],[1,0,1],[1,1,0]],[x1,x2,x3]);`):elif nargs=1且args[1]=SeqFromRec,则print(`SeqFromRec(ope,n,n,Ini,K):给定第一个L-1`):打印(`序列项Ini=[f(1),…,f(L-1)]`):print(`满足递归操作(n,n)f(n)=0`):print(`将其扩展到前K个值`):print(`例如,try:`):打印(`SeqFromRec(N-N-1,N,N,[1],10);`):elif nops([args])=1,op(1,[args')=SVT,则print(`SVT(N):输入一个非负的列表。整数N=[n1,…,nd],并输出奇异向量元组的数量`):打印(`对于一般的n1 x n2 x…x nd张量`):打印(`它使用文章中给出的常量公式`):print(Shmuel Friedland和Giorgio Ottaviani,奇异向量元组的数量和张量的最佳秩一近似的唯一性,):印刷品(《Found.Comput.Math.14》(2014),第6期,1209-1242.):打印(`尝试:SVT([3,3,3]);`):elif nops([args])=1,op(1,[args')=SVT纸张,然后print(`SVTpaper(d,n,Kama,M,K):输入整数d、符号n、正整数Kama、正整数M和大正整数K`):print(`输出一篇文章,首先告诉您序列枚举的第一个Kama项`):打印(`一个n乘以n乘以…乘以n(d倍)张量的奇异向量元组的数目,使用本文中的常数项公式:`):print(“Shmuel Friedland and Giorgio Ottaviani,奇异向量元组的数量和张量最佳秩一近似的唯一性”):印刷品(《Found.Comput.Math.14》(2014),第6期,1209-1242.):print(“然后它试图用那么多项找到一个由该序列满足的齐次线性递归方程。如果它失败了”):打印(`论文在那里完成了。如果成功了,它会陈述出来,并试图用它来找到一个简单的公式。`):打印(`如果成功,它会告诉您。请尝试:`):打印(`SVT纸张(3,n,70,52000):`):elif nops([args])=1,op(1,[args')=SVTseq,则print(`SVTseq(d,N):输入正整数d和正整数N,输出SVT([N$d])从N=1到N=N的前N项:打印(`警告,非常慢,直接从定义中完成。如果d=3,则执行SVTseq3(N),如果d=4,则执行SVTseq4(N)`):打印(`Try:`):打印(`SVTseq(3,30);`):elif nops([args])=1,op(1,[args')=SVTseqSlow thenprint(`SVTseqSlow(d,N):输入正整数d和正整数N,输出SVT([N$d])的前N项,从N=1到N=N。`):打印(`警告,非常慢,直接从定义中完成。如果d=3,则执行SVTseq3(N),如果d=4,则执行SVTseq4(N)`):打印(`Try:`):打印(`SVTseqSlow(3,30);`):elif nops([args])=1,op(1,[args')=SVTseq3,则print(`SVTseq3(N):输入一个正整数N,输出SVT([N,N,N])从N=1到N=N的前N项:打印(`Try:`):打印(`SVTseq3(30);`):elif nops([args])=1,op(1,[args')=SVTseq4,则print(`SVTseq4(N):输入一个正整数N,输出SVT([N,N,N])的前N项,从N=1到N=N.`):打印(`Try:`):打印(`SVTseq4(6);`):elif nops([args])=1,op(1,[args')=SVTseqPC,然后打印(`SVTseqPC(d,N):与SVTseq(d,N)类似,但对于d=3和N<=70预先计算以节省时间。`):打印(`如果不是d=3且N<=70,则使用SVTseq(d,N)`):打印(`Try:`):打印(`SVTseqPC(3,30);`):elif nops([args])=1,op(1,[args')=Zinn,则print(`Zinn(L):输入整数列表,使用Zinn-Justin的方法进行估算`):print(“μ和θ使得”):打印(`resh[i]约为常数*mu^i*i^theta`):print(`例如,try:`):打印(`Zinn([seq(5*i*2^i,i=1..30)]);`):其他的print(`没有ezra for `,args):图1:结束时间:######从FindRec开始sn:=进程(resh,n1):-1/log(op(n1+1,resh)*op(n1-1,resh结束时间:#Zinn(resh):Zinn-Justin的估算方法#C、μ和θ使得#resh[i]是近似值。常数*mu^i*i^theta#例如,尝试:#Zinn([seq(5*i*2^i,i=1..30)]);Zinn:=进程(重新)局部s1,s2,θ,mu,n1,i:如果nops({seq(符号(resh[i]),i=1..nops(resh))})<>1,则返回(失败):图1:n1:=nops(重新)-1:s1:=sn(resh,n1):s2:=sn(resh,n1-1):θ:=evalf(2*(s1+s2)/(s1-s2)^2):mu:=evalf(sqrt(op(n1+1,resh)/op(n1-1,resw))*exp(-(s1+s2)/((s1-s2)*s1)):[θ,μ]:结束时间:#findrec(f,DEGREE,ORDER,n,n):猜测一个递归算子的湮灭#度和序的顺序f#例如,尝试:findrec([seq(i,i=1..10)],0,2,n,n);findrecVerbose:=proc(f,度,顺序,n,n)本地操作,var,eq,i,j,n0,kv,var1,eq1,mu,a:如果(1+度)*(1+阶)+3+阶>nops(f),则错误(`数据不足,无法重复出现顺序`,order,`degree`,degree):图1:操作:=0:变量:={}:对于i从0到ORDER do对于0到度do的jope:=ope+a[i,j]*n^j*n^i:var:=var联合{a[i,j]}:日期:日期:等式:={}:对于n0,从1到(1+度)*(1+阶)+2 do等式1:=0:对于i从0到ORDER doeq1:=eq1+subs(n=n0,系数(ope,n,i))*op(n0+i,f):日期:eq:=eq联合{eq1}:日期:var1:=求解(eq,var):kv:={}:对于i从1到nops(var1)domu:=op(i,var1):如果op(1,mu)=op(2,mu),则kv:=kv联合{op(1,mu)}:图1:日期:ope:=子(var1,ope):如果ope=0,则返回(失败):图1:ope:={seq(系数(展开(ope),kv[i],1),i=1..nops(kv))}减去{0}:如果nops(ope)>1,则print(`There is some slack,There`,nops(ope)):打印(操作):返回(Yafe(操作[1],N)[2]):elif nops(ope)=1,则返回(Yafe(操作[1],N)[2]):其他的返回(失败):图1:结束时间:#findrec(f,DEGREE,ORDER,n,n):猜测一个递归算子的湮灭#度和序的顺序f#例如,尝试:findrec([seq(i,i=1..10)],0,2,n,n);findrec:=proc(f,度,阶,n,n)本地操作,var,eq,i,j,n0,kv,var1,eq1,mu,a:选项记住:如果不是findrecEx(f,DEGREE,ORDER,ithprime(20)),则返回(失败):图1:如果不是findrecEx(f,DEGREE,ORDER,ithprime(40)),则返回(失败):图1:如果不是findrecEx(f,DEGREE,ORDER,ithprime(80)),则返回(失败):图1:如果(1+度)*(1+阶)+5+阶>nops(f),则错误(`数据不足,无法重复出现顺序`,order,`degree`,degree):图1:操作:=0:变量:={}:对于i从0到ORDER do对于从0到DEGREE的j doope:=ope+a[i,j]*n^j*n^i:var:=var联合{a[i,j]}:日期:日期:等式:={}:对于从1到(1+度)*(1+阶)+4 do的n0等式1:=0:对于i从0到ORDER doeq1:=eq1+subs(n=n0,系数(ope,n,i))*op(n0+i,f):日期:eq:=eq并集{eq1}:日期:var1:=求解(eq,var):kv:={}:对于i从1到nops(var1)domu:=op(i,var1):如果op(1,mu)=op(2,mukv:=kv联合{op(1,mu)}:图1:日期:ope:=子(var1,ope):如果ope=0,则返回(失败):图1:ope:={seq(系数(展开(ope),kv[i],1),i=1..nops(kv))}减去{0}:如果nops(ope)>1,则返回(Yafe(操作[1],N)[2]):elif nops(ope)=1,则返回(Yafe(操作[1],N)[2]):其他的返回(失败):图1:结束时间:Yafe:=proc(ope,N)局部i,ope1,coe1,L:如果ope=0,则返回(1,0):图1:ope1:=展开(ope):五十: =度(ope1,N):coe1:=系数(ope1,N,L):ope1:=正常(ope1/coe1):ope1:=正常(ope1):操作1:=转换([seq(系数(系数(ope1,N,i))*N^i,i=度(ope2,N)..度(ope 1,N))],`+`):系数(coe1),ope1:结束时间:#FindrecOld(f,n,n,MaxC):给定一个列表f,尝试用#聚咖啡。#最大阶数+阶数<=最大C#例如,尝试Findrec([1,1,2,3,5,8,13,21,34,55,89],n,n,2);FindrecOld:=进程(f,n,n,MaxC)本地学位、订单、操作、L:对于从0到MaxC的L do对于从0到L的订单学位:=L级:如果(2+度)*(1+阶)+4>=nops(f),则打印(`degree,degree,`andorder,order`数据不足):返回(失败):图1:ope:=findrec([op(1..(2+度)*(1+阶)+4,f)],度,阶,n,n):如果操作失败,则返回(操作):图1:日期:日期:失败:结束时间:#Findrec(f,n,n):给定一个列表f,尝试用#多边形系数。Findrec:=proc(f,n,n)局部DEGREE,ORDER,ope:选项记住:对于来自1的ORDER,而2*(1+ORDER)+5+ORDER<=nops(f)do对于从1开始的度,而(1+DEGREE)*(1+ORDER)+5+ORDER<=nops(f)doope:=findrec(f,度,阶,n,n):如果操作失败,则返回(操作):图1:日期:日期:失败:结束时间:#SeqFromRec(ope,n,n,Ini,K):给定第一个L-1#序列Ini=[f(1),…,f(L-1)]的项#满足递归操作(n,n)f(n)=0#将其扩展到第一个K值SeqFromRec:=进程(操作,n,n,Ini,K)本地ope1、gu、L、n1、j1:ope1:=Yafe(ope,N)[2]:五十: =度(ope1,N):如果nops(Ini)<>L,则错误(`Ini应为长度`,L):图1:ope1:=展开(子项(n=n-L,ope1)/n^L):gu:=Ini:对于n1,从nops(Ini)+1到K dogu:=[op(gu),-add(gu[nops(gu)+1-j1]*subs(n=n1,coeff(ope1,n,-j1)),j1=1..L)]:日期:古:结束时间:#结束Findrec带有(linalg):#findrecEx(f,DEGREE,ORDER,m1):探究thre#很有可能再次出现学位#并使用素数m1对order进行排序#例如,尝试:findrecEx([seq(i,i=1..10)],0,2,n,n,1003);findrecEx:=proc(f,度,顺序,m1)本地操作,var,eq,i,j,n0,eq1,a,A1,D1、E1、等式、变量、f1、n、n:选项记住:f1:=f模m1:如果(1+度)*(1+阶)+5+阶>nops(f),则错误(`数据不足,无法重复出现顺序`,order,`degree`,degree):图1:操作:=0:变量:={}:对于i从0到ORDER do对于从0到DEGREE的j doope:=ope+a[i,j]*n^j*n^i:var:=var联合{a[i,j]}:日期:日期:等式:={}:对于从1到(1+度)*(1+阶)+4 do的n0等式1:=0:对于i从0到ORDER doeq1:=eq1+subs(n=n0,系数(ope,n,i))*op(n0+i,f1)mod m1:日期:eq:=eq联合{eq1}:日期:等式:=转换(等式,列表):变量:=转换(变量,列表):D1:=nops(变量):E1:=nops(等式):如果是E1然后是0返回(false):图1:如果E1-D1>=1,则对于j从1到nops(Var)doA1[D1,j]:=系数(等式[D1+1],Var[j]):日期:如果det(A1)mod m1<>0,则返回(假):图1:图1:如果E1-D1>=2,则对于j从1到nops(Var)doA1[D1,j]:=系数(等式[D1+2],Var[j]):日期:如果det(A1)mod m1<>0,则返回(假):图1:图1:正确:结束时间:###从FindRec结束#CharRoots(ope,n,n):输入具有多项式系数的线性递归算子,#在变量n中,并且其中前进档运算符由n表示,即Nf(n):=f(n+1)#并输出有序列表#它的特征根。尝试:#特征根((n+4)*n^2-n*n-1,n,n);CharRoots:=proc(ope,n,n)局部ope1,lu:ope1:=lcoeff(展开(ope),n):如果度(ope1,N)=0,则返回(失败):图1:lu:=[求解(ope1,N)]:lu:=排序(lu):结束时间:#Asy(ope,n,n,mu,M):输入多项式系数的线性递归算子,#变量n中,向前移位运算符用n表示,即Nf(n):=f(n+1)#它的一个特征根mu发现#一个渐近级数,其前导项为μ^n,满足递归性,与初始条件无关。尝试:#Asy((n+5)*n^2-(3*n+4)*n+2*n+2,n,n,2);Asy:=proc(ope,n,n,mu,M)局部ope,ope1,gu,t,i,x,d,eq,var,c,f:Ope:=展开(数字(Ope)):d: =度(Ope,n):ope1:=系数(展开(Ope),n,d):如果简化(subs(N=mu,ope1))<>0,则返回(失败):图1:gu:=加(x^d*subs(n=1/x,系数(Ope,n,i))*mu^i*(1+i*x)^t,i=ldegree(Ope,n)。。度(Ope,N):gu:=展开(gu):gu:=泰勒(gu,x=0.2):如果展开(系数(gu,x,1))=0,则返回(失败):图1:t: =求解(系数(gu,x,1),t):f: =1+加(c[i]/n^i,i=1..M):var:={seq(c[i],i=1..M)}:gu:=加(x^d*subs(n=1/x,系数(Ope,n,i))*mu^i*(1+i*x)^t*subs。。度(Ope,N):gu:=展开(gu):gu:=泰勒(gu,x=0,M+3):eq:={seq(展开(系数(gu,x,i)),i=2..M+1)}:var:=求解(eq,var):如果var=NULL,则返回(失败):图1:mu^n*n^t*subs(var,f):结束时间:#AsyC(ope,n,n,mu,M,INI,K):见ezraAsyC:=进程(ope,n,n,mu,M,INI,K)局部ope,gu,L1,C:Ope:=展开(数字(Ope)):gu:=Asy(ope,n,n,mu,M):如果gu=失败,则返回(失败):图1:如果nops(INI)<>度(Ope,N),则打印(INI,`应为长度`,度(Ope,N)):返回(失败):图1:L1:=序列源记录(Ope,n,n,INI,K):如果abs(evalf(L1[K]/subs(n=K,gu))<10^(-10),则返回(失败):图1:C: =评估(L1[K]/subs(n=K,gu)):C: =识别(evalf(C,10)):C*gu:结束时间:#SVT(N):输入非负列表。整数N=[n1,…,nd],并输出奇异向量元组的数量#对于一般n1 x n2 x。。。x nd张量。#警告,非常缓慢,直接从定义开始。如果N为3-dim,则执行SVT3SVT:=proc(N)局部z,d,i,gu,lu:选项记住:d: =无(N):lu:=添加(z[i],i=1..d):gu:=1:对于i从1到d dogu:=gu*正常(((lu-z[i])^N[i]-z[i]^N[i])/(lu-2*z[i)):日期:gu:=展开(gu):因为我从1到d做gu:=系数(gu,z[i],N[i]-1):日期:古:结束时间:#SVTseqSlow(d,N):输入正整数d和正整数N,输出B([N$d])从N=1到N=N的前N项。#警告,非常缓慢,直接从定义开始。如果N是3维,则执行SVTseq3(N);如果N是4维,则进行SVTseq 4(N),#尝试:#SVTseqSlow(3,30);SVTseqSlow:=进程(d,N)本地n1:[seq(SVT([n1$d]),n1=1..N)]:结束时间:#SVTseq(d,N):输入正整数d和正整数N,输出B([N$d])从N=1到N=N的前N项。##尝试:#SVTseq(3,30);SVTseq:=进程(d,N)局部n1:[seq(Adf([n1$d]),n1=0..N-1)]:结束时间:#SVTseqPC(d,N):与SVTseq(d,N)类似,但如果d=3且N<=70,则它会给出预先计算的值#尝试:#SVTseqPC(3,30);SVTseqPC:=进程(d,N)局部n1,L1,L2:第一层:=[1, 6, 37, 240, 1621, 11256, 79717, 572928, 4164841, 30553116, 225817021, 1679454816, 12556853401, 94313192616, 711189994357, 5381592930816, 40848410792017, 310909645663332, 2372280474687277, 18141232682656320, 139010366280363601, 1067160872528170536, 8206301850166625797, 63203453697218605440, 487480961825645988721, 3764885111909676266856, 29112650384540312124397, 225377607037294291831008, 1746648730442554074446041, 13549903631370307542264936, 105214484529950220234997141, 817709228332592281063815168, 6360391464220968384566019745, 49511928791684034124645730676, 385707501237368849206626328717, 3006846684379174364465426255616, 23455977069620091841962275705089, 183092319857177553953377636176456, 1430037723632571918813195121704997, 1117556244135989347042802862284928、873846304864847451749353666860367841,683632453962320044584305182576200288, 5350886598153793514481537897337176493, 41902005901674406571033924405839308768, 328277426887267454292261947436072998521, 2572978040988090019983161573900671390056, 2017年4938455027696505833506309037744673717,158256711967850325867904932175164613805568, 1241878121565109721305598151979300621187041, 9748914697534511954263668793857706936885416, 76557533012847845401553482298831558586858781, 601406624861928916400901634219448551509045216, 4725977743774598911951207466985442875802004041, 37149518272412304070221079042862933387873130216, 292110984170941385701186299397195372784867345221, 2297582777856431871486896289498367921392694659456,18076679242941480643112276737449814240450399614961, 142261109059091890572028251448021584831777493283976, 1119875738825722199969846898012482760095108120469421, 8817913547145642408782946965788762982054339519598816, 69449739062865201034991808469623104398486763332361881, 547117715028973597496117057639950244728453193275378536, 4311149329690424172211655089600627818487132234309850581, 33978514528762393389117490950553726381113420890825711616, 267862418462871260646089896898133632801592296425890350401, 2112090051496950649407979640746003519860385829377421719956, 16657262165394749047223410973405780093084150560578941606221, 131396189046352933535318043187955359702871390804208939532416, 1036686394280916864442473886582846818530956264890966401359201, 8180787959749319986573210152825554401372254551057541041233416]:L2级:=[1, 24, 997, 51264, 2940841, 180296088, 11559133741, 765337680384, 51921457661905, 3590122671128664, 252070718210663749, 17922684123178825536, 1287832671004683373753, 93368940577497932331288, 6821632357294515590873917, 501741975445243527381995520, 37121266623211130111114816929, 2760712710223967190110979892824, 206267049696409355312012281872181, 15475371710788303634337687336412224, 1165407938272751946314416465834086889, 88061985234092428001001009986265842520, 6674841972181830152646239523454436828797, 507365393175531718935697292222219542209024, 38665843085507634007413709865882990263469041,2953740030556494206297508806673140156327404248, 226139608318389776420643582181314917081622957061, 17348847768732971801911980996414434839471856232768, 1333495999756767180730069621847180744764824179044825,102679492043788645453302773152678289049925015355480088, 7919467782979625176350608079488125284094134882129213277, 611762406707957516785060447056716898028654120149205549056, 47326297794171750329395108432285966368588072502073467990849, 3666210684055121798038265960682272734448162887732387739437656, 284375989192008571736288450556047894658675111432554006369198101, 22085002503311452557061359446412536526728124547277566177787892288, 1717125381944152148262069150947173852066288121156424979958627271625, 133653371349899710144611275023346214165051471055586217662726722907736,10413736715423917925250315399202563233072956051319130146752533719981053, 812193134071724372030028233213580873621215660428760806777376397599783424, 63404056487163136402118230809098571181603513366964484046606874882278962769, 4954048702957761157113524005575838540485231482644154318297501579169017363736, 387410031814835435731535807396509767481143902823175117339431737579621537848213,3032019104953631935927245935785181752925383265226972671661881411766405320238624\0, 2374803834313774545027493708891468517565412197712400933613898364537188007266\415065, 18614165945764264383509075035019318709521380225496240889482096963404716\0142280221656, 1460043170995964708856156528528938389222377695391471169729078084\4366152779284725947757, 1145990988900840641498974352772435574812691486580998832\888840442789761192955239298207744, 90007356802866260340681675743343295423352567\419796418449737235391358299924670403492251105, 70736693805083529397480221136055\92581392110691450778592637430119392203264791502783256731864, 556250953151679190\326199484492492324591698905227339332845928022594121900029904848181878033749, 43\7669196230078776788737597365736295916093100444453781362348341558871850881876505\01694445489216, 344556869088661514142383463598133274533590051098179797634568099\3602708821591585237841811565673065, 2713979783990190110602974150069380468418799\26951597708635047601874185952373678877541675978258404568, 213882183493150066670\6483254859097343585216990790743086036561160356374638779426179264517198748656322\9, 1686385703886231464111507816879306616123474300541011700087381363447966207871\114366608478451366840144384, 13302888122049843048583783811778374771716769379195\1549233571332531946972464580592254702500946111712214577, 1049867536127299365636\5497227516437477855862285648134225237865420533609968529304268060706932598043796\535256, 82892607577994651456907533699515307250159154916797481720621027970596587\8763162178927441457954079168327435141、65476182751544172368076277493509617188339\0561330168239058377228182666146594166925063393822640826631362216652096, 5174054\6109669857407254864558010313843670271715836656387781610740275475557105777415128\17639228972097875063728153, 409027887749565473877643124644902175375544568079990\109075128262610563868460599234291869165223723204244596556580376,32347708141131\6703815112633244396816051373975539285755764374214440614410184428935385384474359\17662519366143191410589, 255916030760866495394527698951579075886076340724451964\0658122483760011731366878142308998366812224871870032223390662656, 2025395240895\1603246698248194238532846559543259009791531066210961323946905292332059644843530\9068417769429529956502685313, 1603524387663098386262074342517422638718342567678\2167996633336356019580148925904327992960651893088013216529811030500811352, 1269\9605587989209055001298402406398177669502575479970411995455982570766601035104530\61245176163841999507880978168128227321301, 100611904723621504989197065102687600\56534655553050036319143930926830585636462050949830520738795678147410970775185526\58278518336, 797349578873080205697751714231745229204003541172569345165375061133\1139698481136957917921419063990733623143125922769418501460361, 6320979719089295\5360382815502363297770861920831751714025080216909873677318083716624321039453174\6401261452355825000537762507812568, 5012477014633297839861201240291319860394600\7439132120509864627446835189634880231488902214692421032853360665076770619392321\171964157, 39760246007434002263122771664553011899184436308132273433144559413030\99745331143603268930380461032113850457743160958373952262009784832, 315479292291\3956577203778130688262135870108448404747784377827959627670582691460350012918017\35333187359438736270047927662791312330794513, 250388500221019531312713421084697\4108947429561673585282648831960094871813676160195925602880976709018133779466713\5887587705734069647224728, 1987816823728556596123625241672421825758958175776974\8598555938762939620589203386825268137188450332215916917816807530042965920535983\12894517, 157853302741847454273439691565911926574130298546204629589813585802776\508147757855151240838895127044522108017538593152066002261339870960194880, 12538\4340752346336853109321776495015359570646919793393138881734833630652099581700596\32146740282752565168203907489021983261477646662277381157849, 996190396276612251\1705647985861203267338110511349611130284466638206704819483647150478457066892315\24090521025556311814462365392336365072956287320, 791677018490697051619705534671\0226992773793074176044840611548011015860153453335364163605644990813188457321020\0459832117005673118892374654134729741, 6292999305076221635289521103045702396556\7376258264950085614014090402431439242581649692447928776432672043049984187261251\24412095339186826385979822080, 500344014007587203400001211826809049233422877132\7641155252146099007741056540797243004321781335218618839157977192712385019667895\41245938813176876065441, 397904295110177421949295675871828872974954018166295543\918640865248936056564428559978261186946670801542126458929974693277169748109675\4808754698102758680, 3165082579580981019251700832149775222812744145017531963735\1662364725334284986046134953645539434430463186131518819367838854558031474818175\64155533105479557, 251817353557598225180603802212268857541397484886337857197230\22795519864524807088276887689166290567591663002232193443179862722539742397\66085268872234560, 200391008125539684200079917324872930227709966655258060255824\6613861124071281344613410787763539498777520450453286440611775838083231201347796\8123653989200572169, 1594999756584197937091376922485656814469152689497999394432\8298668163227068670231581456699380322714150785801699035367785590624403429191411\25154923097178694101784, 126978664303800230048248365830511062083289620698244402\6029228239668793253643263609196501277414985310316356995046620271553280116065881\03981425934300726451012014445, 101108257767235959829205622808187626030087725574\0083953925780215620136332179178655408085422816846337203330708366098913467279944\8372382157581384048668873433076978176、8052417864994305539040572420568131583499\9730642478434734308748734687047297220429693220672311182622622531466083639080822\1475211272627518645358168693593115314639650865, 6414280001915575133736154854176\2126269147211399946165883923705134557821649378820742943693615808568491307650353\408642911901481644815713508372697983466547381624580088984, 51103382899294247295\3254172102565644737684958492382612372875509431010005371969091516202908232303526\7467941530833349312464046507487025779488297520521107193830389633201749, 4072206\2143744200374032801101071020066037586589924393678952402729912329661888558484672\2264183193238403129380380290080150371748126599234960329602439950125606394240547\435328, 32455369529124983521698348415787175273222918527213635107991789844112393\2936474026634592819859895767430163620142373287545847142914231553647287608110560\38156251482959017915577, 258713088005783581345986148647833214410935763755425067\8489673329018435731802025930947928694299728822248249096458582782759970039327392\397322391481123153017756601550869712717016, 20626412580361576753969781335412891\1746799490826642473542333196922308931176801844910842778549891743970867718151171\563738206902856248717386966184856153669522720001199393298593101, 16447538968034\7639937459460201287562707660861565560865034009374275547406662591987035247125543\1890673538137993003726605995976897147118933198139273849479582463692850041026480\0493568, 1311742475245668886089440516090460879217389738410928738340599513010903\3702245389167597855634031489176850524751528156463035826886634819842524312077986\12394112462467087561914459972545, 104632169460506745543037354025742426410639372\5611297068124324998083522514289486856007612900419270319119779715818648862820481\94024504765585587462025050901069089290189485679471464178904, 834736628061662143\0072352616934982796497241260060485638513673536445597251986895370586080624266250\2316702764017941400119600532512199664174644505379064037808273162077272192460332\39561461、6660393438770326033747326351168176537246337607561298874676981936\3305591339385803585921408064482807876740189356283370522312105693971947064238639\81366352404470368586638132792642025024, 531514773440676360168130393354967776826\7528050086635907987786001741901310843264343107577234932454090569312661525257204\3833627467420865066902434459184172282414583999076165755373790928695049、4242231\6052795413113172460739819368710758076902717876886789795582992846251243404121309\8386912180274694681833123444117101483385099407291468896319209536042224479872146\5346137039825333940805720, 3386380938240000283548844348656961586868383333620671\8053814194190770673931801651648103684148227782716806488695026791990480815721112\6302747981626037303254014002121929535185789724389130405980445, 2703575130703808\7514099317793796184145755757661734868175305795591792813036749399071326545031768\1600515636653074082234098163769402728927754852637057359871864909102466484033943\26786860528631069184, 215874379089302750363434099339956856165439763543601478254\82540812957994711597300985731935013799824966670402592422254595871368004176624119\623107904933101646550373409097151531725469605047871702334289, 17239419984865383\2578851474687227714626127688851198213618562023982474832941210578405158001998691\5507889321175240855981695826859982216933686899920181463978651507754984377988758\57094539817536031491224, 137689868153063928099873333533128517858367534535863772\7292604837460531139474545320078959898100562293194691352804447704120110227180638\5791978778226851094950539722271728058633084683603811083199796810933, 1099861575\2100832605395344771452142626486926763354545797087130977669518521307026996298176\5578888814716562889155614549732572477892130623640328748890689675678139417224122\8077994858349916559930497196726592, 8786779802522957580011353673943727787972540\6144613920377738046535236628937979712714677511935324724006177134550964707202769\0155528580474501859184149621594961451829759048495015594607481368258200517339691\77, 702063036987849353607209476891253472021486502536546242088832391600897387166\5558806008653685082823757772734359939779427221609358458166972237222212991391831\674267898210779250381540724993215430710809091637080, 56101715022433697024041069\7517664062038838213270295968405053561236911744923709706505124711432911080090546\8860981638271595028171871478531471802047757125291790458131645292406682330370277\99221994448419271992141, 448362074377990766195021286896356851920786439358786160\1237509388045184203863238979103610871462453373402268947216757552120565243118166\3516144458210435782115810038749325395856156029804972765105063247818630696960, 3583714208831485840578896294599155651628619947433246571053255848726194672796341\3797051012092580679764655712393855284036697524319773457805027103324323042817776\31855744946231369871312625188116653155965568635636321, 286476315985921240993141\4870223879677228862456446998356026407059298275383926229624749994926096428088666\7062656625857394777545994850016634523095135187169761380329534674073331629828580\8885978063233356671176929732056、229030882279005185475467799062701482655958335\2616908139295263144823788880951028738965830482073085086439042037053313923924317\2789644978769687267557733981383524338071233603547722299112132835801809387507985\80068155797, 183125352334330098712122561662980426921208778350383462145410463826\4680179881143231661424395130380893017405128221705229819522338336181968817515748\090199076848180385444685918581083991272970802462712881475178279087266880, 14643\7114518221480891687503635204462755058110770323827130354247029644748925830135536\3856305109230321708681548409807336076079089089009058628912165372672757248446099\39587950862100666386129499533089258783162245410036572009, 117111971969687057427\0210254881556032962353495981146853207702229157516642153370685912737642204566606\2622071936538669116702420910986977180382935134321807589362243345077789896456362\637854345505219281498852113714539128858584, 93669470859584928227586717292422075\4015683270581229363183561062069985279944377170279872066780284103646223763907752\953453865977461565322242701421070129816838698224377149984515004313945260364804197\62859144179496033914773546461, 749274068105905168876166859985374134357379312133\1831205101121921103303505234356840444257709229051392165161228233718375196837262\6331008832925675985575471372334720937227594844354248627281156257248212956975928\225249370721513984, 59941614933679463626504967326006653987145863699362348228589\6743612342900420147848389655133228300405590320581518634329459320050047004324130\1102778842845465216882867492527080317380210940911812087410683452703985854156802\265119409, 47957941878972464875437332732908289794308358487242825334257123056236\4961222307155241434265194530226984658295860898472846937787947315279893702655155\6400918466165360488699334143543149321203782032890089033510719518011351811193983\60, 383739293435178647672194769351789761649874413148601781731438951369594823344\8372925665904260246149082142480192679543186361849077910944584442810671478465015\4930274080324523686547229226676317207514738179545473053440563981693915851397, 3070824078000560388022474386115571855815928223911271993834520574732080220966078\7046194765660857372819020127474436174740212536376118680458802696278073129149363\29550638122220941587080990682733076498066391855408106775628122493604885824, 245\7626295136293652765756383256173311519115250867946831355236480001623985468918024\7541761193098509227952733181243998219759311118764070837071077317520055598504976\0363175987350680787500961016980930193166001019553647945143818583004710041, 1967\0632412423664537517753243083900445150435222365273043230045469432966150066154363\7065432946903757213919813514866071320953991101387069752544874939092570501758286\11133644152352531051232026965651420729669238667800882644205567654807233048, 157\4569061362861436030957981531576384763064351327180209417047625839369250751169633\9781484516929509891676918490128446446910000111440785998433443895571639861378908\13011820117846846999782962934533532469073189759728293070979551374424127301661,126050733451628912225951110363658495326752259575892831594755614780956578401661\6505972324038623764096423084659721998549557305618833211923474534729772635090193\8913017468751723915570829631330348441794960905063977418799793973100218202380566\528, 10091801196146249416397185495065889785550803661735153048983813791346720333\8380607669465873482085458855696141620594927297202455419978448727416955356621451\4341931018279416583887832077261173021621426832566164280973840827279782709915881\9768634625, 8080365962873125590289273560015133494621139501528831264195251674800\6888030607018726196294803444801424711814516634898530441502287861589234133381499\6730407840168867516239914473544772709590295699252067342620591514214888711947056\411007154106831448, 64704101621165712779466560043844027692907736782934228665968\7998933635985154050447370955969795533625097164457548305899477153353204089003209\3412798336632007585629459073800745306112029288863444260384461706308392819496918\1989962985248175802394362197, 5181678289483604607342187357659336100505698542310\980016938585959552079597119494109550848237661697125071413229910351846096470159945\7963785291252644667698397591227305293023591760499041844767822044043674634784240\8235077273434718194193183444541313910336, 4149983426538839925133935842999297763\8321569805787844674008655401509410967557984593713085153516812441295507274678731\3211194750275985284178215063522807185289251960128205437333660064686803923170520\765926183381303059270496138382638986953643802160067849, 33239846015519822683306\8874882951732051053161566176721634834861840558723304827047817831190746729377835\6573697283804832766085983266341857238633994660389821971658637682706085613816956\7901709342224678679131697482990930124402494491175659764349171116001496, 2662611\6206961049933674054716288921263260412313763748185026192913979934486383794171939\6907653397366631306602753694635742255777641921451745626927350232025498947145950\7349731896247377018272225061932350358881326025049628055049294153620360071973035\452359421, 21330070383484011501444379604071250784758660040858944722511781659462\0578628346743411220016706711041207616564731106269275314118505045096308106516098\0601768807052950841353100452877348965666483894408952911542197950547516070896596\31710194229211637340011145728, 170888122260835095803664806950007981332418643767\0130899871808284558483456341155972790761191787313800753165517615682325305499655\8560488602713501945634537752655290914067406207106327571417424687745498036674456\436026110175479092971218817558151261014852444227985, 13691974461984292604199859\2500440237429016722969734674682495200305948512125237053588655918399266558912759\8336003885301047641364280502532848875893985387816005904252155780050180296708526\016302056059801200934369849295608699506394112867790002072071901145603212440, 10\9712066997091643087962492545894433468584647525595467147826243328462876878091882\5409126630769399879068702181681647729352218057874287486921539330989309124073742\7443109107098338767259136560209011851133702916365999653851610392257962829105947\6501012566648113551733、8791770441383329722422199075589339407048153964657858407\3775793131444315340887615388658638101539598535665876552050544437365623809450881\7481046492625387329122489476265958925319991125055825626022122266552693941575094\2211546581403277872249755763269032479650216035648, 7045818148128715655179262090\1338146532299012569656685251683136139034245273043730899393841325305302890537626\0476364594782294008116351395074540185565936439590960504061975751806020623499179\362072794647713891290957677976864379507873513872925671761608835968975067155417,5647018355506196211433150770749725858723988418945725075385722208891431825215048\9752855200010786526875595270408856086659398979705714350497798116225995647807008\0838986531134345875770558448108272203316727690403693767351327613317059299970544\85321548080578040126389617752, 452625679941866052315328612569957965184259976613\1877776905462517883943306453127866097289465750974810223077276845361631712830147\8681838860550309485737742949554244089961965853414689344606781346106403917691208\476689989275400516505234394464957549927933137954761845092737357、362819778431962\0877695455275676808393485252814576790645350032975740618415506466947640710082246\9432673200842956422108576525938915517082032333501784600396942523683415565652038\8635790560831843097652659750221236352067763772160994002321949668588618666075833\200362408246759424, 29085333690455109056766890138687902816039160548377242020743\4388108747579284908367718561383750634783163628965271409549448420642955971093270\9552682935487208040464690460915410135527014321908620429267355555141907299749320\6864596474705930825520367253889353230242128077516636705, 2331782565963045777760\3800369868519761525181093510663144577018641685951004226867998683776135919890032\9773543582121822432143662209424211630823604475459521461424760096282119831077025\5917307559573105175340155404403425256546729207824197918042671122884712090440825\177006091934616, 18695304054438320801911884354338651452197392631824802633379680\5417290118378543551933439431406495648859209795166356250843178090409614705590611\2564858066201480115169784820626494457415687659766136161570406185020689717821199\13890382859134348205509249965939756878112474861232710757, 149901881352968983939\5522293397990426683721885460698662768129402438301423576720565319872014536955837\3066659006489047588654498373565283644775207970467963168570785059349283102910438\2762707794775849496738009350800508088399259623672920147234099489517128259954712\60838186647916890688, 120201890614396043292152088357771575561484315957331258199\8232001840450729584949447275001419398692935370254667854056646885050369077616181\1985913731122408814782312585636352682983951465375670988396362007945190747369011\31999482651052487460650274053445226538455297816341370792168019785, 963928412350\6855464449809920875436150950158978724735593078362062580367081608324773260416447\5537240915090867972917971623929500252975810303606489303752372534473043627885386\4436039209313509777765737707229688131987871984640935269854081481874480275117805\84740034918665999836819456968088、77304921561973471173747996636737093551836158\4887957327539893447470040317300038938247284501679175881750994473384343419957541\0241424107825147360729741931073643856081804626725297103435345215948991331623351\9359081081520431387822205346992144971806148707416350933510784091736410452709031\7、62009015580241457496300218636373206269814727106792288578063930961372304129\2210339914169979839281968071600352254683404531216408454522578254397460637422561\6103916174017158718088276984133755409680874526577812779758964226036196777184045\558546786851600776591224549613246799646881870762496, 49729831827414204843565862\0056484566753505119393736165263259760101209882504232142360486868719906919764798\2901463034901710678025294227349396119423337305744185047609606360154982540346481\4997710236263669257861037135299341744902323155194040188376598749728900583148262\455830225353552724863025, 39889972375891613617936323283599955783877853073347496\3316202836114740806316737052061221172237832297217362159812075964566983373912705\9734323381673162519969545306726055994300165533323128690406330285832859378045967\203810200116252772157122341467047827164471568791050174169914928119764117138456,3199910862166606912277039254677820822698798048582346913021255972216009631801265\7076644870788999105453345753166533590120923924011680346876235422526513632842936\3680518530708988284274394878347930726829330777235558713466009270450380601887413\641581818448251982158977339359134500219882489714959477, 25670780528212513771497\5786239322500941241903338947911715035803494535811745645970033249431923813601477\9397804399480031510560882509634808745840071780396315018739988775523286994877341\0144298625174323900333314829399166065848597608588698203562236466485304667688100\948299959508762765145378143752000, 20595245326726584758610229212329448286354824\5862087620972762811149921147073211253255414204396928994324544086125148245694571\4967366135155771605346006752890211231889808217867801046918983421527497699318215\6639628472974269408006773362412833486225942350295065998951745853259364139791080\47329992029945, 165242299274620869147709542011007038985041685512469910763517631\2666259162833465228639426153414841113494546131429806952390996776391367704742357\9959493293373707538351484171796994680284131697107609412683574659372894239132405\6340683245461602523496132126071400038458249536013876688274467954743726131288, 1325871760716780556386436615872058537411437268448479124741252070527340936586750\6067893265399350698138529271459899897810509808714415544124534247214517903067705\3919990127850311769321814382469491028176443145312166226059110810057814106466896\44956319603742937831916338068889473123005840642313569999898157, 106391643517342\7606703119661956209811986342604111601457596843271127467762084662827790711336845\0033727607117475889803378026638962040474731621640569943321407714941476540708075\6734955812396284518795976259267514901874237355504728421602907948342743988391929\0358715337467434742316043994769648852923099283456]:如果d=3并且N<=nops(L1),则返回([操作(1..N,L1)]):elif d=4且N<=nops(L2),则返回([op(1..N,L2)]):其他的SVTseq(d,N):图1:结束时间:#EmpirAsyC(L1,n,mu,t,M):输入一组正数L1,这些正数被认为给出了形式的非对称性#mu^n*n^t,使用ofer M的sympotics查找常量。尝试:#EmpirAsyC([seq(5*3^i/i*(1+4/i),i=1..40)],3,-1,5);EmpirAsyC:=proc(L1,n,mu,t,M)局部eq,var,c,lu,i:如果nops(L1)<20或M>nops(L2)/2,则打印(“列表太短”):返回(失败):图1:lu:=mu^n*n^t*add(c[i]/n^i,i=0..M):var:={seq(c[i],i=0..M)}:eq:={seq(L1[i]-subs(n=i,lu),i=nops(L1)-M.nops(L2))}:var:=求解(eq,var):如果var=NULL,则返回(失败):图1:evalf(子项(var,c[0]),10):结束时间:#SVTpaper(d,n,Kama,M,K):输入整数d、符号n、正整数Kama、正整数M和大正整数K`):#print(`输出一篇文章,首先告诉您序列枚举的第一个Kama项#n乘n乘n的奇异向量元组的个数。。。通过n(d倍)张量,使用本文中的常数项公式:#Shmuel Friedland和Giorgio Ottaviani,奇异向量元组的数目和张量最佳秩一近似的唯一性#已找到。计算。数学。14(2014),第6期,1209-1242。#然后,它试图用这么多的项找到由该序列满足的齐次线性递归方程。如果它失败了#论文在那里完成了。如果它成功了,它会陈述它,并尝试用它来找到一个非对称公式。#如果它成功了,它就会告诉你。SVTpaper:=proc(d,n,Kama,M,K)局部gu,a,ope,t0,ope,ru,mu,lu,n,i,ku,C:t0:=时间():打印(``):gu:=SVTseqPC(d,卡玛):print(`The Number of Singular Vector Tuples of a general`,[n$d],`张量`):打印(``):打印(`By Shalosh B.Ekhad`):打印(``):print(`定理:设a(n)为a`,d,'-维张量的奇异向量元组的个数,其维数均为`,n):打印(``):print(`The first`,Kama,`该序列的术语是`):打印(``):打印(gu):ope:=Findrec(gu,n,n):如果ope=FAIL,则print(`虽然我们知道序列是P-递归的,即满足一些具有常系数的齐次线性递归方程`):print(`input中的词条数`,Kama,`不够。如果您真的感兴趣,请尝试放大。`):打印(``):print(‘不过,让我们对渐近性进行一些非严格估计。’):ku:=Zinn(gu):print(`序列的行为似乎类似于`,C*ku[2]^n*n^ku[1],`对于某些常数C`):如果evalf(ku[2]/(d-1)^d-1)<0.1且evalf(ku[1]+(d-1)/2)<0.2,则print(`这些接近我们的推测渐近`):打印(C*((d-1)^d)^n*n^(-(d-1,/2)):打印(``):打印(`对于某些常量C`):#print(`估计得到的常数可能近似`):#打印(``):#C0:=EmpirAsyC(gu,n,(d-1)^d,-(d-1)/2,M):#打印(C0):#打印(``):#print(`因此渐近性的非严格估计为`):#打印(``):#打印(C0*((d-1)^d)^n*n^(-(d-1,/2)):#打印(``):图1:print(`This ends This paper that taked `,time()-t0,` seconds.`):打印(“祝你今天愉快”):返回():图1:print(`定理1:序列a(n)满足以下多项式系数线性递归方程`):打印(``):打印(添加(系数(ope,N,i)*a(N+i),i=0..度(ope、N))=0):打印(``):打印(“根据初始条件”):打印(``):lprint(seq(a(i)=gu[i],i=1..度(ope,N)):打印(``):打印(`and in Maple notification`):打印(``):lprint(加上(系数(ope,N,i)*a(N+i),i=0..度(ope、N))=0):打印(``):print(`只是为了好玩,这个序列的第`,K,`项,即`的奇异向量元组的数量):print(`a`,d,`-维张量,其维数均等于`,K,`等于`):打印(``):lprint(SeqFromRec(操作,n,n,[op(1..度(操作,n),gu)],K)[K]):打印(``):Ope:=展开(数字(Ope)):ru:=排序([求解(lcoeff(Ope,n),n)]):mu:=ru[nops(ru)]:lu:=AsyC(Ope,n,n,mu,M,[op(1..度(Ope、n),gu)],K):如果lu=失败且nops(ru)>1,则mu:=ru[nops(ru)-1]:lu:=AsyC(Ope,n,n,mu,M,[op(1..度(Ope、n),gu)],K):图1:如果lu=失败,则print(`我们无法找到渐近展开,但我们知道它存在。欢迎您亲自尝试。'):打印(``):print(`This ends This paper that taked `,time()-t0,` seconds.`):打印(“祝你今天愉快”):返回():图1:print(`定理2:a(n)的渐近公式,阶为`,M,`由`给出):打印(lu):打印(``):print(`和以Maple表示法`):打印(``):l打印(lu):打印(``):打印(``):print(`This ends This paper that taked `,time()-t0,` seconds.`):打印(“祝你今天愉快”):结束时间:#MMMT(A,x):麦克马洪主定理在变量x列表中给出的生成函数。#尝试:#MMMT([[0,1,1],[1,0,1],[1,1,0]],[x1,x2,x3]):MMMT:=程序(A,x)局部M,M1,i,j:M: =[]:对于i从1到nops(A)doM1:=[]:对于从1到nops(A[i])的j,do如果i=j,则M1:=[op(M1),1-A[i][j]*x[i]]:其他的M1:=[操作(M1),-A[i][j]*x[i]]:图1:日期:M: =[op(M),M1]:日期:1/det(M):结束时间:SVTseq3:=进程(N)本地i:[序列(A3f(i,i,i),i=0..N-1)]:结束时间:SVTseq4:=进程(N)本地i:[序列(A4f(i,i,i),i=0..N-1)]:结束时间:S1:=proc():op(排序([args])):结束:#B3f(a1,a2,a3):x1^a1*x2^a2*x3^a3在1/(1-2*(x1*x2+x1*x3+x2*x3)-2*x1*x2*x3*)中的系数B3f:=过程(a1,a2,a3)选项记住:如果a1<0、a2<0或a3<0,则0:elif a1=0、a2=0和a3=0,则1:其他的B3f(S1(a1-1,a2-1,a3))+B3f图1:结束时间:#A3f(a1,a2,a3):0的B3(b1,b2,B3)之和,A3f:=程序(a1,a2,a3)选项记住:如果a1<0、a2<0或a3<0,则0:elif a1=0、a2=0和a3=0,则1:其他的A3f(S1(a1-1,a2,a3))+A3f(S2(a1,a2-1,a3-A3f(S1(a1-1,a2-1,a3))-A3f+A3f(S1(a1-1,a2-1,a3-1))+B3f(S1(a1,a2,a3)):图1:结束时间:#B4f(a1,a2,a3,a4):x1^a1*x2^a2*x3^a3*x4^a4在1/(1-2*e_2(x1,x2,x3,x4)-3*e3(x1、x2、x3、x4)-4e4(x1和x2、x3和x4)中的系数B4f:=进程(a1、a2、a3、a4)选项记住:如果a1<0或a2<0或a3<0或r4<0,则0:elif a1=0,a2=0,以及a3=0和a4=0,然后1:其他的B4f(S1(a1-1,a2-1,a3,a4))+B4f+2*(B4f(S1(a1-1,a2-1,a3-1,a4))+B4f+3*B4f(S1(a1-1,a2-1,a3-1,a4-1)):图1:结束时间:#A4f(a1,a2,a3,a4):0的B4(b1,b2,b3,B4)之和A4f:=进程(a1,a2,a3,a4)选项记住:如果a1<0或a2<0或a3<0或r4<0,则0:elif a1=0且a2=0且a3=0且a4=0,则1:其他的A4f(S1(a1-1,a2,a3,a4))+A4f-A4f(S1(a1-1,a2-1,a3,a4))-A4f+A4f(S1(a1-1,a2-1,a3-1,a4))+A4f-A4f(S1(a1-1,a2-1,a3-1,a4-1))+B4f(S1(a1,a2,a3,a4)):图1:结束时间:#Cnk(n,k):长度为n且k为1的所有0-1向量Cnk:=proc(n,k)局部gu,gu1,gu11:选项记住:如果n=0,则如果k=0,那么返回({[]}):其他的返回({}):图1:图1:如果k>n或k<0,则返回({}):图1:gu1:=Cnk(n-1,k):gu:={seq([op(gu11),0],gu1中的gu11)}:gu1:=Cnk(n-1,k-1):gu:=gu联合{seq([op(gu11),1],gu11在gu1中)}:结束时间:#Bdf(a):f(a),其中a是非负整数列表,系数x1^k1*。。。xd ^kd(扩展^kd)#有理函数1/(1-e2(x1,…,xd)-2*e3(x1、…、xd)-(d-1)e_d(x1#尝试:#Bdf([3,3,3]);Bdf:=proc(a)局部d,gu,gu1,k,su:选项记住:d: =无(a):如果min(op(a))<0,则返回(0):图1:如果a=[0$d],则返回(1):图1:su:=0:对于k从2到d dogu:=Cnk(d,k):su:=su+(k-1)*add(Bdf(排序(a-gu1)),gu1 in gu):日期:苏:结束时间:#Adf(a):c(a),其中a是非负整数列表,系数x1^k1*。。。xd ^kd(扩展^kd)#有理函数1/(1-e_2(x1,…,xd)-2*e_3(x1、…、xd)-(d-1)e_d(x1))*(1/((1-x1)**(1倍)#尝试:#Adf([3,3,3]);Adf:=proc(a)局部d,gu,gu1,k,su:选项记住:d: =无(a):如果min(op(a))<0,则返回(0):图1:如果a=[0$d],则返回(1):图1:su:=Bdf(a):对于k从1到d dogu:=Cnk(d,k):su:=su+(-1)^(k-1)*添加(Adf(排序(a-gu1)),gu1 in gu):日期:苏:结束时间:#Bdr(a,r):f(a),其中a是长度d的非负整数列表,例如,输出x1^a[1]*的系数。。。xd^a[d]#有理函数1/(1-e2(x1,…,xd)-2*e3(x1、…、xd)-(d-1)e_d(x1#尝试:#Bdr([3,3,3]);Bdr:=proc(a,r)局部d,gu,gu1,k,su:选项记住:如果r=1,则返回(Bdf(a)):图1:d: =无(a):如果min(op(a))<0,则返回(0):图1:如果a=[0$d],则返回(1):图1:su:=Bdr(a,r-1):对于k从2到d dogu:=Cnk(d,k):su:=su+(k-1)*add(Bdr(排序(a-gu1),r),gu1 in gu):日期:苏:结束时间: