poj3233

令B={ A A}

{O E}

那么B^k={A A+A^2+A^3...A^K}

{O E}

代码：

#include<stdio.h> const int NUM=60; int n,k,m; struct Matrix { int data[NUM][NUM]; Matrix operator =(Matrix a) { for(int i=0;i<n;i++) { for(int j=0;j<2*n;j++) { data[i][j]=a.data[i][j]; } } copy(); } Matrix operator *(Matrix b) { Matrix c; for(int i=0;i<2*n;i++) { for(int j=0;j<2*n;j++) { int sum=0; for(int k=0;k<2*n;k++) { sum+=data[i][k]*b.data[k][j]; sum=sum%m; } c.data[i][j]=sum; } } c.copy(); return c; } void copy() { for(int i=n;i<2*n;i++) { for(int j=0;j<n;j++) { data[i][j]=0; } for(int j=n;j<2*n;j++) { if(i==j) data[i][j]=1; else data[i][j]=0; } } } void setE() { for(int i=0;i<2*n;i++) { for(int j=0;j<2*n;j++) { if(i==j) data[i][j]=1; else data[i][j]=0; } } } void print() { for(int i=0;i<n;i++) { for(int j=n;j<2*n;j++) { printf("%d ",data[i][j]); } printf("\n"); } //printf("\n"); } }input,output; int main() { freopen("input.txt","r",stdin); while(scanf("%d %d %d",&n,&k,&m)!=EOF) { for(int i=0;i<n;i++) { for(int j=0;j<n;j++) { scanf("%d",&input.data[i][j]); input.data[i][j+n]=input.data[i][j]; input.copy(); } } Matrix temp=input; output.setE(); //input.print(); //output.print(); while(k>0) { if(k&1) { output=(output*temp); //(output*temp).print(); } temp=temp*temp; //temp.print(); k=k>>1; //output.print(); } output.print(); } return 0; }