PKUSC 2021

求和变换

给定一大小为 $n\times n$ 的矩阵 $M$，每次变换将 $M_{i,j}$ 变为当前的 ${\sum }_{k=0}^{n-1}{M}_{i,k}+M_{k,j}$。求原矩阵变换 $t$ 次的结果，每个位置的答案对 $P$ 取模。
$1⩽n⩽1000,0⩽t⩽10^9,2⩽P⩽10^9$
时间限制1s，空间限制512MB。

设 $n\times n$ 的矩阵 $S$ 满足 $S_{i,j}=1$，则 $MS$ 即可得到将行求和后的矩阵。同理 $SM$ 即为将列求和后的矩阵。设 $F(M,t)$ 表示 $M$ 经过 $t$ 次变换后的结果，则有 $F(M,0)=M,F(M,t)=SF(M,t-1)+F(M,t-1)S$。发现规律：
$F(M,0)=M$
$F(M,1)=SM+MS$
$F(M,2)=S^{2}M+2SMS+M{S}^2$
$F(M,3)=S^{3}M+3S^{2}MS+3SM{S}^2+M{S}^3$
$F(M,4)=S^{4}M+4S^{3}MS+6S^{2}M{S}^2+4SM{S}^3+M{S}^4$
归纳可得：
$F(M,t)={\sum }_{i=0}^t\left(\genfrac{}{}{0ex}{}{t}{i}\right)S^{i}M{S}^{t-i}$
容易发现 $S^0=I,S^t=n^{t-1}S$。设 $B=S^{t_1}M{S}^{t_2}$，根据矩阵乘法的定义：
$B_{i,l}={\sum }_{j=0}^{n-1}{\sum }_{k=0}^{n-1}{n}^{t_1-1}{S}_{i,j}{M}_{j,k}{n}^{t_2-1}{S}_{k,l}=n^{t1+t2-2}{\sum }_{j=0}^{n-1}{\sum }_{k=0}^{n-1}{M}_{j,k}$
设矩阵 $M^{\prime }$ 中的任意一位均为 $M$ 中所有数的和，得到：
F ( M , t ) = S t M + M S t + ∑ i = 1 t − 1 ( t i ) n t − 2 M ′ = n t − 1 ( S M + M S ) + n t − 2 M ′ ∑ i = 1 t − 1 ( t i ) F(M,t)=S^tM+MS^t+\sum_{i=1}^{t-1} \binom{t}{i}n^{t-2} M'=n^{t-1}(SM+MS)+n^{t-2} M' \sum_{i=1}^{t-1} \binom{t}{i} F(M,t)=StM+MSt+∑i=1t−1​(it​)nt−2M′=nt−1(SM+MS)+nt−2M′∑