小截斷多元多項式行列式求算

2021-10-06 18:31:05 字數 951 閱讀 6082

對於域 a

\mathbb a

a 上的 n

nn 階二元多項式矩陣 a(x

,y

)\mathbf a(x, y)

a(x,y)

,其中 x

xx 的次數都小於 a

aa,y

yy 的次數都小於 b

bb。我們想要求 det⁡a

modx

amod

yb

\det \mathbf a \bmod x^a \bmod y^b

detamo

dxam

odyb

。其中主要問題在於消元的過程中,可能出現常數項為 0

00 的情況,如果是一元多項式我們還可以用低次消去高次的方法,但多元多項式就顯得要複雜了。

因此我們考慮構造佔位多項式 xiy

j→xi

yjti

+j

x^iy^j \rightarrow x^iy^jt^

xiyj→x

iyjt

i+j,在這樣的情況下我們做插值,問題就變為了 abab

ab個 a+b

a+ba+

b 次一元多項式矩陣的行列式問題,因此我們通過低次消去高次的方法就可以完成消元。樸素實現是 θ(n

3(a+

b)2a

b)

\theta(n^3(a+b)^2ab)

θ(n3(a

+b)2

ab),如果使用 fft,那麼就變成了 θ(n

3(a+

b)ab

log⁡(a

+b))

\theta(n^3(a+b)ab\log (a+b))

θ(n3(a

+b)a

blog(a

+b))

。這個方法也可以用於更多元的情況,其實如果上來就讓你解決子集卷積的形式,可能反而容易想到。

MATLAB多項式及多項式擬合

多項式均表示為陣列形式,陣列元素為多項式降冪係數 1.polyval函式 求多項式在某一點或某幾個點的值.p 1,1,1 x 2 x 1 x 1,0,1 y polyval p,x 另外求函式在某一點或某幾個點的值可以用函式feval.x 1,0,1 y feval x exp x x 注意用的乘法...

數學 多項式 多項式求逆

多項式求逆 前置知識 ntt 給定乙個多項式 f left x right 求乙個多項式 g left x right 使得 f left x right g left x right equiv 1 left bmod 998244353 right 考慮遞迴求解。假定現在已經求出了 g 0 le...

matlab多項式與非多項式擬合

擬合標準 1 原始資料向量與擬合向量之間的距離最小,該距離的度量一般使用誤差平方和表示,即均方誤差 r q y 22 2 當均方誤差最小時,說明構造的擬合向量與原始向量最為接近,這種曲線擬合的方法稱為最小二乘法 3 計算均方誤差最小時的擬合係數,可以通過微積分中求解極值的方法實現 多項式擬合 1 多...