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=10010;
const double eps=1e-8;
int n,m;
double a[maxn][2]={};
double b[5][5]={};//初始
double c[5][5]={};//添加
double e[5][5]={};//结束
void doit(int k){
    if(k==1){
        for(int i=1;i<=3;i++){
            for(int j=1;j<=3;j++){
                b[i][j]=c[i][j];
            }
        }
    }else{
        memset(e,0,sizeof(e));
        for(int i=1;i<=3;i++){
            for(int j=1;j<=3;j++){
                for(int w=1;w<=3;w++){
                    e[i][j]+=c[i][w]*b[w][j];
                }
            }
        }
        for(int i=1;i<=3;i++){
            for(int j=1;j<=3;j++){
                b[i][j]=e[i][j];
            }
        }
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++){
        scanf("%lf%lf",&a[i][0],&a[i][1]);
    }double x,y;
    for(int i=1;i<=m;++i){
        memset(c,0,sizeof(c));
        char ch;
        scanf("%s",&ch);
        if(ch=='M'){
            scanf("%lf%lf",&x,&y);
            c[1][1]=1.0;c[2][2]=1.0;
            c[3][3]=1.0;c[1][3]=x;c[2][3]=y;
        }else if(ch=='X'){
            c[1][1]=1.0;c[2][2]=-1.0;c[3][3]=1.0;
        }else if(ch=='Y'){
            c[1][1]=-1.0;c[2][2]=1.0;c[3][3]=1.0;
        }else if(ch=='S'){
            scanf("%lf",&x);
            c[1][1]=x;c[2][2]=x;c[3][3]=1.0;
        }else{
            scanf("%lf",&x);
            c[1][1]=cos(x/180*M_PI);c[1][2]=-sin(x/180*M_PI);
            c[2][1]=sin(x/180*M_PI);c[2][2]=cos(x/180*M_PI);
            c[3][3]=1.0;
        }
        doit(i);
    }
    for(int i=1;i<=n;i++){
        x=a[i][0]*b[1][1]+a[i][1]*b[1][2]+1*b[1][3];
        y=a[i][0]*b[2][1]+a[i][1]*b[2][2]+1*b[2][3];
        printf("%.1f %.1f
",x+eps,y+eps);
    }
    return 0;
}