关于矩阵乘法的知识matrix67大牛的博客http://www.matrix67.com/blog/archives/276/已经讲得非常好了额,里面一共总结十个能利用矩阵乘法的解决的题目类型。这道题目是第一中类型的一个模板题。
之前连续做到多到矩阵乘法的题目,于是这几天开始学习矩阵乘法。
矩阵的乘法真的是很神奇的东西。
点的变换
时间限制:
2000 ms | 内存限制:
65535 KB
难度:
5
-
描述
-
平面上有不超过10000个点,坐标都是已知的,现在可能对所有的点做以下几种操作:
平移一定距离(M),相对X轴上下翻转(X),相对Y轴左右翻转(Y),坐标缩小或放大一定的倍数(S),所有点对坐标原点逆时针旋转一定角度(R)。
操作的次数不超过1000000次,求最终所有点的坐标。
提示:如果程序中用到PI的值,可以用acos(-1.0)获得。
-
输入
-
只有一组测试数据
测试数据的第一行是两个整数N,M,分别表示点的个数与操作的个数(N<=10000,M<=1000000)
随后的一行有N对数对,每个数对的第一个数表示一个点的x坐标,第二个数表示y坐标,这些点初始坐标大小绝对值不超过100。
随后的M行,每行代表一种操作,行首是一个字符:
首字符如果是M,则表示平移操作,该行后面将跟两个数x,y,表示把所有点按向量(x,y)平移;
首字符如果是X,则表示把所有点相对于X轴进行上下翻转;
首字符如果是Y,则表示把所有点相对于Y轴进行左右翻转;
首字符如果是S,则随后将跟一个数P,表示坐标放大P倍;
首字符如果是R,则随后将跟一个数A,表示所有点相对坐标原点逆时针旋转一定的角度A(单位是度)
输出
-
每行输出两个数,表示一个点的坐标(对结果四舍五入到小数点后1位,输出一位小数位)
点的输出顺序应与输入顺序保持一致
样例输入
-
2 5 1.0 2.0 2.0 3.0 X Y M 2.0 3.0 S 2.0 R 180
样例输出
-
-2.0 -2.0 0.0 0.0
-
只有一组测试数据
#include<iostream>
#include<cstdio>
#include<cmath>
#define MAXN 20
#define PI acos(-1.0)
using namespace std;
struct Matrix
{
int size;
double element[MAXN][MAXN];
void setSize(int);
Matrix operator* (Matrix);
};
void Matrix::setSize(int a)
{
for (int i=0; i<a; i++)
for (int j=0; j<a; j++)
element[i][j]=0;
size = a;
}
Matrix Matrix::operator* (Matrix param)
{
Matrix product;
product.setSize(size);
for (int i=0; i<size; i++)
for (int j=0; j<size; j++)
for (int k=0; k<size; k++)
product.element[i][j]+=element[i][k]*param.element[k][j];
return product;
}
Matrix m[5],beg;
double pos[11000][3];
int n,k;
void init()
{
for(int i=0;i<5;i++)
m[i].setSize(3);
m[0].element[0][0]=m[0].element[1][1]=m[0].element[2][2]=1;
m[1].element[2][2]=1;
m[2].element[0][0]=m[2].element[2][2]=1,m[2].element[1][1]=-1;
m[3].element[0][0]=-1,m[3].element[1][1]=m[3].element[2][2]=1;
m[4].element[2][2]=1;
}
int main()
{
char s[3];
init();
while(~scanf("%d%d",&n,&k))
{
beg.setSize(3);
beg.element[0][0]=beg.element[1][1]=beg.element[2][2]=1;
for(int i=0;i<n;i++)
scanf("%lf%lf",&pos[i][0],&pos[i][1]),pos[i][2]=1;
for(int i=0;i<k;i++)
{
scanf("%s",s);
if(s[0]=='M')
{
scanf("%lf%lf",&m[0].element[0][2],&m[0].element[1][2]);
beg=m[0]*beg;
}
if(s[0]=='S')
{
scanf("%lf",&m[1].element[0][0]);
m[1].element[1][1]=m[1].element[0][0];
beg=m[1]*beg;
}
if(s[0]=='X')
beg=m[2]*beg;
if(s[0]=='Y')
beg=m[3]*beg;
if(s[0]=='R')
{
double a;
scanf("%lf",&a);
a=(a/180)*PI;
m[4].element[0][0]=cos(a);
m[4].element[0][1]=-sin(a);
m[4].element[1][0]=sin(a);
m[4].element[1][1]=cos(a);
beg=m[4]*beg;
}
}
for(int i=0;i<n;i++)
{
for(int j=0;j<3;j++)
{
double x=0;
for(int k=0;k<3;k++)
x+=beg.element[j][k]*pos[i][k];
if(j<2)
printf("%.1lf%s",x,j==0?" ":"\n");
}
}
}
return 0;
}