UVA - 12183 ：Top Secret（N^2的循环矩阵乘法）

pro：N个数排成一圈。一次操作为，每个位置的数+=L*左+R*右，保留x为整数。 问S轮操作后每个位置的值。 N<=1000，S<=2^30，x<=9 。

sol：不难想到矩阵乘法，但是N为1000，显然普通的N^3矩阵乘法的复杂度不能AC。 不难发现这是经典的循环矩阵（第二行为第一行右移一格....依次），所以我们只需要求第一行的矩阵，其他每一行都可以在第一行找到对应的值。

（我的代码1500ms，VJ有神仙150ms，暂时不知道怎么搞的。

#include<bits/stdc++.h>
#define ll unsigned long long
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
int Mod;
int p[]={,,,,,,,
,,};
int tt[maxn][],res[maxn][];
void read(int &x){
    x=; char c=getchar();
    while(c>''||c<'') c=getchar();
    while(c>=''&&c<='') {
        x=x*+c-''; c=getchar();
    }
}
void MOD(int &x){
    if(x>=Mod) x-=Mod;
}
struct mat{
    int mp[][maxn],C;
    mat(){}
    mat(int c){ rep(i,,) rep(j,,c) mp[i][j]=; C=c; }
    mat friend operator *(mat a,mat b){
        mat res(a.C);
        rep(k,,a.C)
         rep(i,,) //循环矩阵，求第一行即可
          if(a.mp[i][k])
           rep(j,,a.C){
             int t=(j-k+a.C+)%a.C; if(t==) t=a.C;
             MOD(res.mp[i][j]+=1LL*a.mp[i][k]*b.mp[][t]%Mod);
        }
        return res;
    }
    mat friend operator ^(mat a,int x){
        mat res(a.C); rep(i,,) res.mp[i][i]=;
        while(x){
            if(x&) res=res*a;
            a=a*a; x>>=;
        }
        return res;
    }
};
int main()
{
    int T,N,L,R,S,X;
    scanf("%d",&T);
    while(T--){
        read(N); read(S); read(L);
        read(R); read(X);
        Mod=p[X];
        mat a(N),base(N);
        rep(i,,N) read(tt[i][]),res[i][]=;
        rep(i,,) {
            base.mp[i][i-==?N:i-]=L;
            base.mp[i][i+==N+?:i+]=R;
            base.mp[i][i]=;
        }
        base=(base^S);
        rep(k,,N)
         rep(i,,N)
           rep(j,,){
             int t=(k-i+N+)%N; if(t==) t=N;
             MOD(res[i][j]+=1LL*base.mp[][t]*tt[k][j]%Mod);
        }
        rep(i,,N-) printf("%d ",res[i][]);
        printf("%d\n",res[N][]);
    }
    return ;
}