NYOJ 298

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
struct Matrax {
	double m[3][3];
};
double PI;
struct Point {
	double x,y;
}p[10100];


void init(Matrax &p){
	for(int i=0;i<3;i++)
	for(int j=0;j<3;j++)
	p.m[i][j]=0;
}

Matrax enM(double L,double P){
	Matrax c;
	init(c);
	c.m[0][0]=c.m[1][1]=c.m[2][2]=1;
	c.m[0][2]=L;c.m[1][2]=P;
	return c;
}

Matrax enX(){
	Matrax c;
	init(c);
	c.m[1][1]=-1; 
	c.m[0][0]=c.m[2][2]=1;
	return c;
}

Matrax enY(){
	Matrax c;
	init(c);
	c.m[0][0]=-1; 
	c.m[1][1]=c.m[2][2]=1;
	return c;
}

Matrax enS(double P){
	Matrax c;
	init(c);
	c.m[1][1]=P; 
	c.m[0][0]=P;c.m[2][2]=1;
	return c;
}

Matrax enR(double R){
	Matrax c;
	init(c);
	c.m[0][0]=cos(R);c.m[0][1]=-1*sin(R);
	c.m[1][0]=sin(R);c.m[1][1]=cos(R);
	c.m[2][2]=1; 
	return c;
}

Matrax multi(Matrax a,Matrax b){
	Matrax c;
	for(int i=0;i<3;i++){
		for(int j=0;j<3;j++){
			c.m[i][j]=0;
			for(int k=0;k<3;k++)
			c.m[i][j]+=a.m[i][k]*b.m[k][j];
		}
	}
	return c;
}

int main(){
	PI=acos(-1.0);
	int N,M;
	double ts[3];
	char t; double x,y,L,P,R;
	while(scanf("%d %d",&N,&M)!=EOF){
		Matrax ans;
		init(ans);
		ans.m[0][0]=ans.m[1][1]=ans.m[2][2]=1;
		for(int i=0;i<N;i++)
		scanf("%lf%lf",&p[i].x,&p[i].y);
		for(int i=0;i<M;i++){
			cin>>t;
			Matrax tt;
			if(t==‘X‘){
				tt=enX();
				ans=multi(tt,ans);
			}
			else if(t==‘Y‘){
				tt=enY();
				ans=multi(tt,ans);
			}
			else if(t==‘M‘){
				scanf("%lf%lf",&L,&P);
				tt=enM(L,P);
				ans=multi(tt,ans);
			}
			else if(t==‘S‘){
				scanf("%lf",&P);
				tt=enS(P);
				ans=multi(tt,ans);
			}
			else if(t==‘R‘){
				scanf("%lf",&R);
				R=PI*R/180.0;
				tt=enR(R);
				ans=multi(tt,ans);
			}
		}
		for(int i=0;i<N;i++){
			ts[0]=p[i].x; ts[1]=p[i].y; ts[2]=1;
			double ansx=ts[0]*ans.m[0][0]+ts[1]*ans.m[0][1]+ts[2]*ans.m[0][2];
			if(fabs(ansx-0)<1e-8) ansx=0;
			double ansy=ts[0]*ans.m[1][0]+ts[1]*ans.m[1][1]+ts[2]*ans.m[1][2];
			if(fabs(ansy-0)<1e-8) ansy=0;
			printf("%.1lf %.1lf\n",ansx,ansy);
		}
	}
	return 0;
}