[HNOI2019]白兔之舞 Bluestein's Algorithm FFT 生成函数 矩阵乘法

[HNOI2019]白兔之舞

题目传送门
名字不错 ^ o ^

分析

题意：较长，自己看！
听说要单位根反演？不会.jpg。
还好有非单位根反演的做法：-）
首先朴素Dp的话大概就是
$F(i,j,l)={\sum }_{a=0}^{i-1}{\sum }_{b=1}^{n}F(a,b,l-1)w[b][j]$
$l$那么大令人厌恶，给他生成函数掉！
$F(i,j)={\sum }_{a=0}^{i-1}{\sum }_{b=1}^{n}xF(a,b)w[b][j]$
$={\sum }_{b=1}^n({\sum }_{a=0}^{i-1}F(a,b))w[b][j]x$
设$S(i,b)={\sum }_{a=0}^{i}F(a,b)$
$F(i,j)={\sum }_{b=1}^{n}S(i-1,b)w[b][j]x$
$S(i,j)-S(i-1,j)={\sum }_{b=1}^{n}S(i-1,b)w[b][j]x$
式子化到这里如果不是$S$是多项式，我甚至想要矩乘！
那就试一试呗！
构造矩阵
[ w [ 1 ] [ 1 ] x + 1 w [ 1 ] [ 2 ] x w [ 1 ] [ 3 ] x w [ 2 ] [ 1 ] x w [ 2 ] [ 2 ] x + 1 w [ 2 ] [ 3 ] x w [ 3 ] [ 1 ] x w [ 3 ] [ 2 ] x w [ 3 ] [ 3 ] x + 1 ] \left [ \begin{matrix} w[1][1]x&