poj 3233 Matrix Power Series

求矩阵和的问题。这题我参考罗牛的思路，把前K次和分解为2的幂和，比如25 ＝ 1+8+16。比较郁闷的是，我的代码和他的比，可以说除了变量名其它都一样，结果人家跑16ms，我的141ms。

#include <cstdio>
int m,n,K;
int pow[30][31][31],sum[30][31][31];
void add(int a[31][31],int b[31][31],int c[31][31])
{
    for(int i = 0; i < n; i++)
        for(int j = 0; j < n; j++)
            a[i][j] = (c[i][j]+b[i][j])%m;
}
void mul(int a[31][31],int b[31][31],int c[31][31])
{
    int i,j,k;
    for(i = 0; i < n; i++)
        for(j = 0; j < n; j++)
        {
            a[i][j] = 0;
            for(k = 0; k < n; k++)
                a[i][j] = (a[i][j]+b[i][k]*c[k][j])%m;
        }

}
int main()
{
    int i,j,t;
    int tmp[31][31],ans[31][31]={};
    scanf("%d%d%d",&n,&K,&m);
    for(i = 0; i < n; i++)
        for(j = 0; j < n; j++)
        {
            scanf("%d",&t);
            pow[0][i][j] = sum[0][i][j] = t%m;
        }
    for(i = 1; (1<<i) <= K; i++)
    {
        mul(pow[i],pow[i-1],pow[i-1]);
        mul(tmp,sum[i-1],pow[i-1]);
        add(sum[i],tmp,sum[i-1]);
    }
    for(i = 0; K; K>>=1,i++)
    {
        if(K&1)
        {
            mul(tmp,ans,pow[i]);
            add(ans,tmp,sum[i]);
        }
    }
    for(i = 0; i < n; i++)
    {
        printf("%d",ans[i][0]);
        for(j = 1; j < n; j++)
            printf(" %d",ans[i][j]);
        printf("\n");
    }
    return 0;
}