简单矩阵运算类_纯属练手,取mod的时候不怎么好用的说-_-

1. #include <iostream>
2. using namespace std;
3. #define Error {cout<<"Error!"<<endl;return (*this);}
4. class Mat
5. {
6. public:
7. const static int MAX_MAT_SIZE=10;
8. inline void init(int l,int c);
9. inline Mat operator*(const Mat& b);
10. inline void operator*=(const Mat& b);
11. inline Mat operator+(const Mat& b);
12. inline void operator+=(const Mat& b);
13. inline Mat operator-(const Mat& b);
14. inline void operator-=(const Mat& b);
15. inline Mat operator^(int k);
16. inline void show();
17. void in(int l,int c);
18. private:
19.     int mat[MAX_MAT_SIZE][MAX_MAT_SIZE];
20.     long long moudle;
21.     int size_l,size_c;
22. };
23. inline void Mat::init(int l,int c)
24. {
25.     int i,j;
26.     for(i=0;i<l;i++)for(j=0;j<c;j++)mat[i][j]=0;
27.     size_l=l;size_c=c;
28.     moudle=10000;
29. }
30. inline Mat Mat::operator *(const Mat& b)
31. {
32.     if(size_c!=b.size_l) Error
33.     int i,j,k;Mat tmp;
34.     tmp.init(size_l,b.size_c);
35.     for(i=0;i<size_l;i++)
36.         for(k=0;k<b.size_c;k++)
37.             for(j=0;j<size_c;j++)
38.             {tmp.mat[i][k]+=mat[i][j]*b.mat[j][k];if(tmp.mat[i][k]>=tmp.moudle)tmp.mat[i][k]%=tmp.moudle;}
39.     return tmp;
40. }
41. inline void Mat::operator *=(const Mat &b)
42. {
43.     (*this)=(*this)*b;
44. }
46. inline Mat Mat::operator +(const Mat& b)
47. {
48.     if(size_l!=b.size_l||size_c!=b.size_c) Error
49.     int i,j;Mat tmp;
50.     tmp.init(size_l,b.size_c);
51.     for(i=0;i<size_l;i++)
52.         for(j=0;j<size_c;j++)
53.         {tmp.mat[i][j]=mat[i][j]+b.mat[i][j];if(tmp.mat[i][j]>=tmp.moudle)tmp.mat[i][j]%=tmp.moudle;}
54.     return tmp;
55. }
56. inline void Mat::operator +=(const Mat& b)
57. {
58.     (*this)=(*this)+b;
59. }
61. inline Mat Mat::operator -(const Mat& b)
62. {
63.     if(size_l!=b.size_l||size_c!=b.size_c) Error
64.     int i,j;Mat tmp;
65.     tmp.init(size_l,b.size_c);
66.     for(i=0;i<size_l;i++)
67.         for(j=0;j<size_c;j++)
68.             tmp.mat[i][j]=mat[i][j]-b.mat[i][j];
69.     return tmp;
70. }
71. inline void Mat::operator -=(const Mat& b)
72. {
73.     (*this)=(*this)-b;
74. }
76. inline Mat Mat::operator ^(int k)
77. {
78.     if(k==1) return (*this);
79.     Mat tmp=(*this)*(*this);
80.     if(k&1) return (tmp^(k/2))*(*this);
81.     else return tmp^(k/2);
82. }
83. inline void Mat::show()
84. {
85.  int i,j;
86.   for(i=0;i<size_l;cout<<mat[i++][j]<<endl)
87.       for(j=0;j<size_c-1;j++) cout<<mat[i][j]<<' ';
88. }
89. void Mat::in(int l,int c)
90. {
91.     for(int i=0;i<l;i++)for(int j=0;j<c;j++)cin>>mat[i][j];
92. }
93. int main()
94. {
95.     int l,c,k;
96.     Mat a,b;
97.     while(cin>>l>>c)
98.         {
99.             a.init(l,c);
100.             a.in(l,c);
101.             cin>>k;
102.             b=a^k;
103.             b.show();
105.         }
106.     return 0;
107. }