2021湖南多校对抗赛第四场
排名
| 第一 | 第二 | 第三 |
|---|---|---|
| 中南1队 | 中南3队 | 长沙学院一队 |
团体成绩
| 学校 | 总题数 | 总罚时 |
|---|---|---|
| 中南大学 | 17 | 2526 |
题解
D
设
A=(a1a2…anana1…an−1an−1an…an−2…………a2a3…a1)A = \begin{pmatrix} a_{1} & a_{2} & \dots & a_{n} \\
a_{n} & a_{1} & \dots & a_{n-1} \\
a_{n-1} & a_{n} & \dots & a_{n-2}\\
\dots &\dots & \dots & \dots \\
a_{2} & a_{3} & \dots & a_{1}
\end{pmatrix}A=⎝⎜⎜⎜⎜⎛a1anan−1…a2a2a1an…a3……………anan−1an−2…a1⎠⎟⎟⎟⎟⎞
设多项式 f(x)=a1+a2⋅x+⋯+an⋅xn−1f(x) = a_1 + a_2\cdot x + \dots + a_{n} \cdot x^{n-1}f(x)=a1+a2⋅x+⋯+an⋅xn−1
设 wn0w_n^{0}wn0 , wn1w_n^1wn1 , …\dots… , wnn−1w_n^{n-1}wnn−1 是 111 的 nnn 次单位根 .
构造矩阵
w=(111…1(wn0)1(wn1)1(wn2)1…(wnn−1)1(wn0)2(wn1)2(wn2)2…(wnn−1)2……………(wn0)n−1(wn1)n−1(wn2)n−1…(wnn−1)n−1)w = \begin{pmatrix}
1 & 1 & 1 &\dots & 1 \\
(w_n^{0})^1 & (w_n^{1})^1 & (w_n ^ {2})^1 & \dots & (w_n^{n-1})^1 \\
(w_n^{0})^2 & (w_n^{1})^2 & (w_n ^ {2})^2 & \dots & (w_n^{n-1})^2 \\
\dots & \dots & \dots & \dots & \dots \\
(w_n^{0})^{n-1} & (w_n^{1})^{n-1} & (w_n ^ {2})^{n-1} & \dots & (w_n^{n-1})^{n-1} \\
\end{pmatrix}w=⎝⎜⎜⎜⎜⎛1(wn0)1(wn0)2…(wn0)n−11(wn1)1(wn1)2…(wn1)n−11(wn2)1(wn2)2…(wn2)n−1……………1(wnn−1)1(wnn−1)2…(wnn−1)n−1⎠⎟⎟⎟⎟⎞
AW=(f(wn0)f(wn1)f(wn2)…f(wnn−1)wn0f(wn0)wn1f(wn1)wn2f(wn3)…wnn−1f(wnn−1)(wn0)2f(wn0)(wn1)2f(wn1)(wn2)2f(wn3)…(wnn−1)2f(wnn−1)……………(wn0)n−1f(wn0)(wn1)n−1f(wn1)(wn2)n−1f(wn3)…(wnn−1)n−1f(wnn−1))AW = \begin{pmatrix}
f(w_n^0) & f(w_n^1) & f(w_n^2) & \dots & f(w_n^{n-1}) \\
w_n^0f(w_n^0) & w_n^{1}f(w_n^1) & w_n^{2}f(w_n^3) & \dots & w_n^{n-1}f(w_n^{n-1}) \\
(w_n^0)^2f(w_n^0) & (w_n^{1})^2f(w_n^1) & (w_n^{2})^2f(w_n^3) & \dots & (w_n^{n-1})^2f(w_n^{n-1}) \\
\dots & \dots & \dots & \dots & \dots \\
(w_n^0)^{n-1}f(w_n^0) & (w_n^{1})^{n-1}f(w_n^1) & (w_n^{2})^{n-1}f(w_n^3) & \dots & (w_n^{n-1})^{n-1}f(w_n^{n-1}) \\
\end{pmatrix}AW=⎝⎜⎜⎜⎜⎛f(wn0)wn0f(wn0)(wn0)2f(wn0)…(wn0)n−1f(wn0)f(wn1)wn1f(wn1)(wn1)2f(wn1)…(wn1)n−1f(wn1)f(wn2)wn2f(wn3)(wn2)2f(wn3)…(wn2)n−1f(wn3)……………f(wnn−1)wnn−1f(wnn−1)(wnn−1)2f(wnn−1)…(wnn−1)n−1f(wnn−1)⎠⎟⎟⎟⎟⎞
所以det(AW)=det(A)det(W)=∏i=0n−1f(wni)det(W)det(AW) = det(A)det(W) = \prod_{i =0}^{n-1}f(w_n^{i}) det(W)det(AW)=det(A)det(W)=i=0∏n−1f(wni)det(W)
det(A)=∏i=0n−1f(wni)det(A) = \prod_{i=0}^{n-1}f(w_n^i)det(A)=i=0∏n−1f(wni)
如果不考虑精度问题 , 我们可以使用 CZTCZTCZT 来实现任意长度的 DFTDFTDFT .
但是该算法无法通过本题 .
多项式的 结式 被定义为 :
μ1\mu_1μ1 , μ2\mu_2μ2 , …\dots… , μn\mu_nμn 是 nnn 次多项式 A(x)A(x)A(x) 的 nnn 个根 .
λ1\lambda_1λ1 , λ2\lambda_2λ2 , …\dots… , λm\lambda_mλm 是 mmm 次多项式 B(x)B(x)B(x) 的 mmm 个根 .
则 R(A,B)=bmn∏i=1nA(μi)=bmnamn∏j=1n∏i=1m(μi−λj)=(−1)nmanm∏i=1mB(λi)R(A,B) = b_m^n\prod_{i=1}^{n}A(\mu_i) = b_m^na_m^n\prod_{j=1}^{n}\prod_{i=1}^{m}(\mu_i - \lambda_j) = (-1)^{nm}a_n^m\prod_{i=1}^{m}B(\lambda_i)R(A,B)=bmn∏i=1nA(μi)=bmnamn∏j=1n∏i=1m(μi−λj)=(−1)nmanm∏i=1mB(λi)
我们根据定义可以得出如下性质 :
R(A,B)=(−1)nmR(B,A)R(A,B) = (-1)^{nm}R(B,A)R(A,B)=(−1)nmR(B,A).
$R(A,B) = a_nmb_mn \iff m = 0 $ 或者 n=0n = 0n=0 .
若bm=1b_m = 1bm=1 , 则 R(A,B)=R(A−CB,B)R(A,B) = R(A - CB,B)R(A,B)=R(A−CB,B) , 其中 C(x)C(x)C(x) 是一个任意多项式 .
取 C(x)=A(x)B(x)C(x) = \dfrac{A(x)}{B(x)}C(x)=B(x)A(x) , 上式变成 R(A,B)=R(A−AB⋅B,B)=R(Amod B,B)R(A,B) = R(A - \dfrac{A}{B}\cdot B,B) = R(A \mod B,B)R(A,B)=R(A−BA⋅B,B)=R(AmodB,B) .
我们可以用类似辗转相除法的方法实现这个算法 .
我们注意到 wn0w_n^0wn0 , wn1w_n^1wn1 , …\dots… , wnn−1w_n^{n-1}wnn−1 是 xn−1x^n - 1xn−1 的 nnn 个根 .
设 g(x)=−1+xng(x) = -1 + x^ng(x)=−1+xn , xnx^nxn 前面的系数是 111 .
所以答案是 R(f,g)R(f,g)R(f,g) .
这篇博客探讨了多项式和傅里叶变换在数学竞赛中的应用。文章通过2021湖南多校对抗赛第四场的问题,介绍了如何利用单位根和傅里叶矩阵计算多项式的卷积,以及卷积在多项式乘法和根的关系。讨论了快速傅里叶变换(FFT)的局限性,并提出了一种类似辗转相除法的算法来处理特定的多项式运算。最后,博客提到了多项式卷积的性质及其在数值计算中的重要性。
5642

被折叠的 条评论
为什么被折叠?



