矩阵十题1（nyoj298）

继续深造矩阵类题目，看到大佬们都有矩阵十题博客，我也学习了，只有题目，没有解析系列1；
通过矩阵变换将所有操作都存到矩阵中，最后输出答案即可。
矩阵的构造方法
直接上代码

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<iostream>
#include<string>
#include <set>
#include<time.h>
//a&3==a%4
using namespace std ;
#define ll double
#define mem(a) memset(a,0,sizeof(a))
const double eps = 1e-8;
const int maxn = 110010;
const int inf = 0x3f3f3f3f;
const ll mod=10007ll;
const double PI=acos(-1.00);

struct matrix
{
    ll mat[10][10];
    matrix()
    {
        for(int i=0;i<10;i++)
        {
            for(int j=0;j<10;j++)
            {
                mat[i][j]=0;
            }
        }
    }
};

void out(matrix a,int n)
{
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            cout<<a.mat[i][j]<<"   ";
        }
        cout<<endl;
    }
}

matrix multiply(matrix a,matrix b,int n)
{
    matrix res;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            for(int k=0;k<n;k++)
            {
                res.mat[i][j]+=a.mat[i][k]*b.mat[k][j];
            }
        }
    }
    return res;
}

struct point
{
    double x,y;
}p[maxn];

int main()
{
    int N,m;
    int n=3;
    while(scanf("%d%d",&N,&m)!=EOF)
    {
        matrix res;
        for(int i=0;i<n;i++)
        {
            res.mat[i][i]=1;
        }

        for(int i=0;i<N;i++)
        {
            scanf("%lf%lf",&p[i].x,&p[i].y);
        }

        matrix X;X.mat[0][0]=1;X.mat[1][1]=-1;X.mat[2][2]=1;

        matrix Y;Y.mat[0][0]=-1;Y.mat[1][1]=1;Y.mat[2][2]=1;

        while(m--)
        {
            char c[10];
            scanf("%s",c);
            switch(c[0])
            {
                case 'M':{
                    double a,b;
                    scanf("%lf%lf",&a,&b);
                    matrix M;M.mat[0][0]=1;M.mat[1][1]=1;M.mat[2][2]=1;M.mat[0][2]=a;M.mat[1][2]=b;
                    res=multiply(M,res,n);}
                    break;
                case 'X':res=multiply(X,res,n);break;
                case 'Y':res=multiply(Y,res,n);break;
                case 'S':{
                    double P;
                    scanf("%lf",&P);
                    matrix S;S.mat[0][0]=P;S.mat[1][1]=P;S.mat[2][2]=1;
                    res=multiply(S,res,n);}
                    break;
                case 'R':{
                    double A;
                    scanf("%lf",&A);
                    A=A*PI/180.00;
                    matrix R;R.mat[0][0]=cos(A);R.mat[1][1]=cos(A);R.mat[0][1]=-sin(A);R.mat[1][0]=sin(A);R.mat[2][2]=1.0;
                    res=multiply(R,res,n);}
                    break;
            }
        }
        for(int i=0;i<N;i++)
        {
            double resx=res.mat[0][0]*p[i].x+res.mat[0][1]*p[i].y+res.mat[0][2];
            double resy=res.mat[1][0]*p[i].x+res.mat[1][1]*p[i].y+res.mat[1][2];
            printf("%.1f %.1f\n",resx,resy);
        }
    }
    return 0;
}