[caioj 1479，利用矩阵乘法解决的经典题目一]点的变换

给定n(n<=10000)个点，m次操作，所有点的坐标都是已知的，现在可能对n个点做以下几种操作：
平移一定距离(M)，相对X轴上下翻转(X)，相对Y轴左右翻转(Y),坐标缩小或放大一定的倍数(S),所有点对坐标原点逆时针旋转一定角度(R)。
操作的次数不超过1000000次，求最终所有点的坐标。

这道题最主要的算法是矩阵乘法，这些操作都可以转换成一个3*3的矩阵。转换后，首先先将总的操作相乘，组成一个最终矩阵，再用这个最终矩阵乘各个点的坐标组成的3*1的矩阵，如{x}{y}{1}（这是竖着的），就可以得出每个点的最终坐标了，这样时间复杂度仅为O（m+n），但暴力为O（mn），快了很多。

图片转自Matrix67大牛

#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
using namespace std;
struct node
{
    double a[5][5];
    node()
    {
        for(int i=1;i<=3;i++)
        {
            for(int j=1;j<=3;j++)a[i][j]=0.0;
        }
    }
}he[11000];
double dx[11000],dy[11000];
node pre1,pre2,pre3,pre4,pre5,ss;
node chengfa1(node a,node b)
{
    node c;
    for(int i=1;i<=3;i++)
    {
        for(int j=1;j<=3;j++)
        {
            for(int k=1;k<=3;k++)
            {
                c.a[i][j]+=a.a[i][k]*b.a[k][j];
            }
        }
    }
    return c;
}
node chengfa2(node a,node b)
{
    node c;
    for(int i=1;i<=3;i++)
    {
        for(int k=1;k<=3;k++)
        {
            c.a[i][1]+=a.a[i][k]*b.a[k][1];
        }
    }
    return c;
}
int main()
{
    char s[5];
    int n,m;
    double x,y;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)scanf("%lf%lf",&dx[i],&dy[i]);
    for(int i=1;i<=n;i++)
    {
        he[i].a[1][1]=dx[i];he[i].a[2][1]=dy[i];he[i].a[3][1]=1;
    }
    ss.a[1][1]=ss.a[2][2]=ss.a[3][3]=1.0;
    pre1.a[1][1]=pre1.a[2][2]=pre1.a[3][3]=1.0;
    pre2.a[3][3]=1.0;
    pre3.a[1][1]=pre3.a[3][3]=1;pre3.a[2][2]=-1.0;
    pre4.a[1][1]=-1.0;pre4.a[2][2]=pre4.a[3][3]=1.0;
    pre5.a[3][3]=1.0;
    for(int i=1;i<=m;i++)
    {
        scanf("%s",s+1);
        if(s[1]=='M')
        {
            scanf("%lf%lf",&x,&y);
            pre1.a[1][3]=x;pre1.a[2][3]=y;
            ss=chengfa1(pre1,ss);
        }
        if(s[1]=='S')
        {
            scanf("%lf",&x);
            pre2.a[1][1]=pre2.a[2][2]=x;
            ss=chengfa1(pre2,ss);
        }
        if(s[1]=='X')ss=chengfa1(pre3,ss);
        if(s[1]=='Y')ss=chengfa1(pre4,ss);
        if(s[1]=='R')
        {
            scanf("%lf",&x);
            double pp=acos(-1.0);
            double aa=pp*x/180.0;
            pre5.a[1][1]=pre5.a[2][2]=cos(aa);
            pre5.a[1][2]=-sin(aa);pre5.a[2][1]=sin(aa);
            ss=chengfa1(pre5,ss);
        }
    }
    for(int i=1;i<=n;i++)
    {
        node ans;
        ans=chengfa2(ss,he[i]);
        printf("%.1lf %.1lf\n",ans.a[1][1],ans.a[2][1]);
    }
    return 0;
}