使用(linalg):P:=proc(q)局部a,b,c,d,j,k,n,P,t,x,y:
对于n从1到q dox:=n:对于p从1到2 do
d: =ilog10(x)+1:a:=转换(x,基数,10):c:=[]:
对于从1到nops(a)的k,做c:=[op(c),a[-k]]:od:t:=[op([]),c]:
对于从2到d的k,做b:=[op([]),c[nops(c)]]:
对于从1到nops(c)-1的j,dob:=[op(b),c[j]]:od:
c: =b:t:=[op(t),c]:od;x: =det(t):如果x=0,则中断:
否则,如果p=1,则y:=x:fi:fi:od:
如果n=x且y<>x,则打印(n);fi:od:end:P(10^8):
|