矩阵运算
矩阵运算:
/**=========================================
矩阵运算
复杂度:O(N^3)
=========================================*/
const int maxn = 105;
const int maxm = 105;
struct Matrix
{
int n, m;
int a[maxn][maxm];
void clear()
{
n = m = 0;
memset(a, 0, sizeof(a));
}
Matrix operator + (const Matrix &b) const
{
Matrix tmp;
tmp.n = n;
tmp.m = m;
for (int i = 0; i < n; i++)
for (int j = 0; j < m; j++)
tmp.a[i][j] = a[i][j] + b.a[i][j];
return tmp;
}
Matrix operator - (const Matrix &b) const
{
Matrix tmp;
tmp.n = n;
tmp.m = m;
for (int i = 0; i < n; i++)
for (int j = 0; j < m; j++)
tmp.a[i][j] = a[i][j] - b.a[i][j];
return tmp;
}
Matrix operator * (const Matrix &b) const
{
Matrix tmp;
tmp.clear();
tmp.n = n;
tmp.m = b.m;
for (int i = 0; i < n; i++)
for (int j = 0; j < b.m; j++)
for (int k = 0; k < m; k++)
tmp.a[i][j] += a[i][k] * b.a[k][j];
return tmp;
}
};
矩阵加法
对应位置相加。(等行等列)
代码如下:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <deque>
#include <list>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <numeric>
#include <iomanip>
#include <bitset>
#include <sstream>
#include <fstream>
#include <limits.h>
#define debug "output for debug\n"
#define pi (acos(-1.0))
#define eps (1e-6)
#define inf (1<<28)
#define sqr(x) (x) * (x)
#define mod 1000000007
using namespace std;
typedef long long ll;
typedef unsigned long long ULL;
const int MAX=105;
struct Matrix
{
int n,m;
int a[MAX][MAX];
void clear()
{
n=m=0;
memset(a,0,sizeof(a));
}
};
Matrix operator +(Matrix a,Matrix b)
{
Matrix c;
c.clear();
c.n=a.n;
c.m=a.m;
for(int i=0;i<a.n;i++)
{
for(int j=0;j<a.m;j++)
{
c.a[i][j]=a.a[i][j]+b.a[i][j];
}
}
return c;
}
int main()
{
int i,j,k,n1,m1;
Matrix a,b,c;
printf("输入矩阵的行(n)列(m):\n");
scanf("%d%d",&n1,&m1);
a.n=n1;a.m=m1;
b.n=n1;b.m=m1;
printf("输入矩阵a:\n");
for(i=0;i<a.n;i++)
{
for(j=0;j<a.m;j++)
scanf("%d",&a.a[i][j]);
}
printf("输入矩阵b:\n");
for(i=0;i<b.n;i++)
{
for(j=0;j<b.m;j++)
scanf("%d",&b.a[i][j]);
}
c=a+b;
printf("矩阵a+b的和:\n");
for(i=0;i<c.m;i++)
{
for(j=0;j<c.m;j++)
printf("%d ",c.a[i][j]);
printf("\n");
}
return 0;
}
矩阵减法
对应位置相减(等行等列)
代码如下:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <deque>
#include <list>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <numeric>
#include <iomanip>
#include <bitset>
#include <sstream>
#include <fstream>
#include <limits.h>
#define debug "output for debug\n"
#define pi (acos(-1.0))
#define eps (1e-6)
#define inf (1<<28)
#define sqr(x) (x) * (x)
#define mod 1000000007
using namespace std;
typedef long long ll;
typedef unsigned long long ULL;
const int MAX=105;
struct Matrix
{
int n,m;
int a[MAX][MAX];
void clear()
{
n=m=0;
memset(a,0,sizeof(a));
}
};
Matrix operator -(Matrix a,Matrix b)
{
Matrix c;
c.clear();
c.n=a.n;
c.m=a.m;
for(int i=0;i<a.n;i++)
{
for(int j=0;j<a.m;j++)
{
c.a[i][j]=a.a[i][j]-b.a[i][j];
}
}
return c;
}
int main()
{
int i,j,k,n1,m1;
Matrix a,b,c;
printf("输入矩阵的行(n)列(m):\n");
scanf("%d%d",&n1,&m1);
a.n=n1;a.m=m1;
b.n=n1;b.m=m1;
printf("输入矩阵a:\n");
for(i=0;i<a.n;i++)
{
for(j=0;j<a.m;j++)
scanf("%d",&a.a[i][j]);
}
printf("输入矩阵b:\n");
for(i=0;i<b.n;i++)
{
for(j=0;j<b.m;j++)
scanf("%d",&b.a[i][j]);
}
c=a-b;
printf("矩阵a-b的差:\n");
for(i=0;i<c.m;i++)
{
for(j=0;j<c.m;j++)
printf("%d ",c.a[i][j]);
printf("\n");
}
return 0;
}
/*
4 4
1 2 3 4
4 3 2 1
4 3 2 1
1 2 3 4
1 2 3 4
4 3 2 1
4 3 2 1
1 2 3 4
*/
矩阵乘法
若一矩阵的列数与另一矩阵的行数相等,则可定义这两个矩阵的乘积。如 A 是 m×n 矩阵和 B 是 n×p矩阵,它们的乘积 AB 是一个 m×p 矩阵。
代码如下:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <deque>
#include <list>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <numeric>
#include <iomanip>
#include <bitset>
#include <sstream>
#include <fstream>
#include <limits.h>
#define debug "output for debug\n"
#define pi (acos(-1.0))
#define eps (1e-6)
#define inf (1<<28)
#define sqr(x) (x) * (x)
#define mod 1000000007
using namespace std;
typedef long long ll;
typedef unsigned long long ULL;
const int MAX=105;
struct Matrix
{
int n,m;
int a[MAX][MAX];
void clear()
{
n=m=0;
memset(a,0,sizeof(a));
}
};
Matrix operator *(Matrix a,Matrix b)
{
Matrix c;
c.clear();
c.n=a.n;
c.m=b.m;
for(int i=0;i<a.n;i++)
{
for(int j=0;j<b.m;j++)
{
for(int k=0;k<a.m;k++)
c.a[i][j]+=a.a[i][k]*b.a[k][j];
}
}
return c;
}
int main()
{
int i,j,k,n1,m1;
Matrix a,b,c,d;
printf("输入矩阵A的行(n1)列(m1):\n");
scanf("%d%d",&a.n,&a.m);
printf("输入矩阵a:\n");
for(i=0;i<a.n;i++)
{
for(j=0;j<a.m;j++)
scanf("%d",&a.a[i][j]);
}
printf("输入矩阵B的行(n2)列(m2):\n");
scanf("%d%d",&b.n,&b.m);
printf("输入矩阵b:\n");
for(i=0;i<b.n;i++)
{
for(j=0;j<b.m;j++)
scanf("%d",&b.a[i][j]);
}
if(a.m==b.n)
{
c=a*b;
printf("矩阵a*b的积:\n");
for(i=0;i<c.n;i++)
{
for(j=0;j<c.m;j++)
printf("%d ",c.a[i][j]);
printf("\n");
}
}
else
printf("a*b is Unlawful!\n");
if(b.m==a.n)
{
d=b*a;
printf("矩阵b*a的积:\n");
for(i=0;i<d.n;i++)
{
for(j=0;j<d.m;j++)
printf("%d ",d.a[i][j]);
printf("\n");
}
}
else
printf("b*a is Unlawful!\n");
return 0;
}
1、给定n个点,m个操作,构造O(m+n)的算法输出m个操作后各点的位置。操作有平移、缩放、翻转和旋转
这里的操作是对所有点同时进行的。其中翻转是以坐标轴为对称轴进行翻转(两种情况),旋转则以原点为中心。如果对每个点分别进行模拟,那么m个操作总共耗时O(mn)。利用矩阵乘法可以在O(m)的时间里把所有操作合并为一个矩阵,然后每个点与该矩阵相乘即可直接得出最终该点的位置,总共耗时O(m+n)。假设初始时某个点的坐标为x和y,下面5个矩阵可以分别对其进行平移、旋转、翻转和旋转操作。预先把所有m个操作所对应的矩阵全部乘起来,再乘以(x,y,1),即可一步得出最终点的位置。
矩阵快速幂
1、A^k
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define MAX 105
int n;
struct Matrix
{
int a[MAX][MAX];
Matrix()
{
memset(a,0,sizeof(a));
}
};
Matrix operator *(Matrix a,Matrix b)
{
Matrix c;
int i,j,k;
memset(c.a,0,sizeof(c.a));
for(k=0;k<n;k++)
{
for(i=0;i<n;i++)
{
if(a.a[i][k]==0)
continue;
for(j=0;j<n;j++)
{
if(b.a[k][j]==0)
continue;
c.a[i][j]=c.a[i][j]+a.a[i][k]*b.a[k][j];
}
}
}
return c;
}
//a^k
Matrix operator ^(Matrix a,int k)
{
Matrix c;
int i;
//初始化为单位矩阵
for(i=0;i<n;i++)
c.a[i][i]=1;
//
while(k)
{
if(k&1)
c=c*a;
a=a*a;
k>>=1;
}
return c;
}
int main()
{
int i,j,k;
Matrix a,b,c;
//printf("输入矩阵的行列(n)和幂值(k):\n");
scanf("%d%d",&n,&k);
//printf("输入n*n的矩阵A:\n");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
scanf("%d",&a.a[i][j]);
}
b=a^k;
//printf("A^k的矩阵:\n");
for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
printf("%d ",b.a[i][j]);
printf("\n");
}
return 0;
}
2、给定矩阵A,请快速计算出A^n(n个A相乘)的结果,输出的每个数都mod p。
由于矩阵乘法具有结合律,因此A^4 = A * A * A * A = (A*A) * (A*A) = A^2 * A^2。我们可以得到这样的结论:当n为偶数时,A^n = A^(n/2) * A^(n/2);当n为奇数时,A^n = A^(n/2) * A^(n/2) * A (其中n/2取整)。这就告诉我们,计算A^n也可以使用二分快速求幂的方法。例如,为了算出A^25的值,我们只需要递归地计算出A^12、A^6、A^3的值即可。根据这里的一些结果,我们可以在计算过程中不断取模,避免高精度运算。