NYOJ298 点的变换 矩阵乘法

关于矩阵乘法的知识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;
}