pku 3150(修改后的矩阵乘法)

起初直接有矩阵乘法做，结果TLE。后来参考网上代码，发现他们巧妙利用转移矩阵的性质。从而降低复杂度。

转移矩阵的下一列是前一列的下移而来。因次计算出第一列，其余的就可推出。

code：

#include<stdio.h>
#include<string.h>

int n,m,k,d;
__int64 ans[510]={0};
__int64 E[501];

int multiply(__int64 (*A),__int64 (*B))
{  
 int i,j,t;
    __int64 C[502]={0};
 for(i=0;i<n;i++)
 {  
  j=(n-i)%n;
  for(t=0;t<n;t++)
  {  
      C[i]+=A[t]*B[j];
   j=(j+1)%n;
  }
 }
 for(i=0;i<n;i++)  
   A[i]=C[i]%m; 

 return 1;
}
int erfen()
{
   while(k)
   {
    if(k&1)multiply(ans,E);
    multiply(E,E);
    k>>=1;
   }
   return 0;
}
int main()
{
    int i,cnt=0;
 while(EOF!=scanf("%d%d%d%d",&n,&m,&d,&k))
 {     
        for(i=0;i<n;i++)
   scanf("%I64d",&ans[i]);
  memset(E,0,sizeof(E));
  for(E[0]=i=1;i<=d;i++)
  { 
   E[i]=E[n-i]=1;
  }
        erfen();
  for(i=0;i<n;i++)
   printf("%I64d ",ans[i]);
  printf("/n");
 }
 return 0;
}