NYOJ 298 点的变换 矩阵乘法

http://acm.nyist.net/JudgeOnline/problem.php?pid=298

最好还是自己手推一下矩阵式子..不算太难..但是有一些小知识....

首先当然是矩阵的细节..矩阵是不支持交换率的..所以如图的式子乘进去时要放在左边...

还有的比如说:

cmath里的sin函数用的是弧度制..需要把度数/180*M_PI ( M_PI是cmath里定义的常数π );

double在取固定小数位的时候小负数四舍五入会出现-0.0之类的情况,可以自己const一个小数eps加上去.

代码如下

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<cmath>
 #include<queue>
 using namespace std;
 const int maxn=;
 const double eps=1e-;
 int n,m;
 double a[maxn][]={};
 double b[][]={};//初始
 double c[][]={};//添加
 double e[][]={};//结束
 void doit(int k){
     if(k==){
         for(int i=;i<=;i++){
             for(int j=;j<=;j++){
                 b[i][j]=c[i][j];
             }
         }
     }else{
         memset(e,,sizeof(e));
         for(int i=;i<=;i++){
             for(int j=;j<=;j++){
                 for(int w=;w<=;w++){
                     e[i][j]+=c[i][w]*b[w][j];
                 }
             }
         }
         for(int i=;i<=;i++){
             for(int j=;j<=;j++){
                 b[i][j]=e[i][j];
             }
         }
     }
 }
 int main(){
     scanf("%d%d",&n,&m);
     for(int i=;i<=n;i++){
         scanf("%lf%lf",&a[i][],&a[i][]);
     }double x,y;
     for(int i=;i<=m;++i){
         memset(c,,sizeof(c));
         char ch;
         scanf("%s",&ch);
         if(ch=='M'){
             scanf("%lf%lf",&x,&y);
             c[][]=1.0;c[][]=1.0;
             c[][]=1.0;c[][]=x;c[][]=y;
         }else if(ch=='X'){
             c[][]=1.0;c[][]=-1.0;c[][]=1.0;
         }else if(ch=='Y'){
             c[][]=-1.0;c[][]=1.0;c[][]=1.0;
         }else if(ch=='S'){
             scanf("%lf",&x);
             c[][]=x;c[][]=x;c[][]=1.0;
         }else{
             scanf("%lf",&x);
             c[][]=cos(x/*M_PI);c[][]=-sin(x/*M_PI);
             c[][]=sin(x/*M_PI);c[][]=cos(x/*M_PI);
             c[][]=1.0;
         }
         doit(i);
     }
     for(int i=;i<=n;i++){
         x=a[i][]*b[][]+a[i][]*b[][]+*b[][];
         y=a[i][]*b[][]+a[i][]*b[][]+*b[][];
         printf("%.1f %.1f\n",x+eps,y+eps);
     }
     return ;
 }