CSP-坐标变换（其二）

问题描述
对于平面直角坐标系上的坐标 (x,y)，小 P 定义了如下两种操作：

拉伸 k 倍：横坐标 x 变为 kx，纵坐标 y 变为 ky；

旋转 θ：将坐标 (x,y) 绕坐标原点 (0,0) 逆时针旋转 θ 弧度（0≤θ<2π）。易知旋转后的横坐标为 xcos⁡θ−ysin⁡θ，纵坐标为 xsin⁡θ+ycos⁡θ。

设定好了包含 n 个操作的序列 (t1,t2,⋯,tn) 后，小 P 又定义了如下查询：

i j x y：坐标 (x,y) 经过操作 ti,⋯,tj（1≤i≤j≤n）后的新坐标。
对于给定的操作序列，试计算 m 个查询的结果。

输入格式
从标准输入读入数据。

输入共 n+m+1 行。

输入的第一行包含空格分隔的两个正整数 n 和 m，分别表示操作和查询个数。

接下来 n 行依次输入 n 个操作，每行包含空格分隔的一个整数（操作类型）和一个实数（k 或 θ），形如 1 k（表示拉伸 k 倍）或 2 θ（表示旋转 θ）。

接下来 m 行依次输入 m 个查询，每行包含空格分隔的四个整数 i、j、x 和 y，含义如前文所述。
————————————————
注意，计算旋转结果时要避免纵坐标受到横坐标已经计算从而变化的影响，此处加入代码

double temp=c;

用于规避该问题。

#include <bits/stdc++.h>
using namespace std;
int main() {
    int n,m;
    cin>>n>>m;
    struct op {
        int type;
        double rat;
    }ops[n];
    for (int i = 0; i < n; i++) {
        cin>>ops[i].type;
        scanf("%lf", &(ops[i].rat));
    }
    for (int i = 0; i < m; i++) {
        int a,b;
        double c,d;
        cin>>a>>b;
        scanf("%lf", &c);
        scanf("%lf", &d);

        for (int j=a-1;j<=b-1;j++) {
            if (ops[j].type==1) {
                c=c*ops[j].rat;
                d=d*ops[j].rat;

            }else {
                double temp=c;
                c=c*cos(ops[j].rat)-d*sin(ops[j].rat);
                d=temp*sin(ops[j].rat)+d*cos(ops[j].rat);

            }
        }
        cout<<fixed<<setprecision(3)<<c<<" "<<fixed<<setprecision(3)<<d<<endl;
    }

    return 0;
}