【CSP】202309-2 坐标变换（其二）

2023年 第31次CCF计算机软件能力认证  202309-2 坐标变换（其二）

原题链接：坐标变换（其二）

时间限制： 2.0 秒

空间限制： 512 MiB

题目描述

对于平面直角坐标系上的坐标 (𝑥,𝑦)，小 P 定义了如下两种操作：

1. 拉伸 𝑘 倍：横坐标 𝑥 变为 𝑘𝑥，纵坐标 𝑦 变为 𝑘𝑦；

2. 旋转 𝜃：将坐标 (𝑥,𝑦) 绕坐标原点 (0,0) 逆时针旋转 𝜃 弧度（0≤𝜃<2𝜋）。易知旋转后的横坐标为 𝑥cos⁡𝜃−𝑦sin⁡𝜃，纵坐标为 𝑥sin⁡𝜃+𝑦cos⁡𝜃。

设定好了包含 𝑛 个操作的序列 (𝑡1,𝑡2,⋯,𝑡𝑛) 后，小 P 又定义了如下查询：

• i j x y：坐标 (𝑥,𝑦) 经过操作 𝑡𝑖,⋯,𝑡𝑗（1≤𝑖≤𝑗≤𝑛）后的新坐标。

对于给定的操作序列，试计算 𝑚 个查询的结果。

输入格式

从标准输入读入数据。

输入共 𝑛+𝑚+1 行。

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

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

接下来 𝑚 行依次输入 𝑚 个查询，每行包含空格分隔的四个整数 𝑖、𝑗、𝑥 和 𝑦，含义如前文所述。

输出格式

输出到标准输出。

输出共 𝑚 行，每行包含空格分隔的两个实数，表示对应查询的结果。

样例输入

10 5
2 0.59
2 4.956
1 0.997
1 1.364
1 1.242
1 0.82
2 2.824
1 0.716
2 0.178
2 4.094
1 6 -953188 -946637
1 9 969538 848081
4 7 -114758 522223
1 9 -535079 601597
8 8 159430 -511187

样例输出

-1858706.758 -83259.993
-1261428.46 201113.678
-75099.123 -738950.159
-119179.897 -789457.532
114151.88 -366009.892

样例解释

第五个查询仅对输入坐标使用了操作八：拉伸 0.716 倍。

横坐标：159430×0.716=114151.88

纵坐标：−511187×0.716=−366009.892

由于具体计算方式不同，程序输出结果可能与真实值有微小差异，样例输出仅保留了三位小数。

子任务

80% 的测试数据满足：𝑛,𝑚≤1000；

全部的测试数据满足：

• 𝑛,𝑚≤10^5；

• 输入的坐标均为整数且绝对值不超过 10^6；

• 单个拉伸操作的系数 𝑘∈[0.5,2]；

• 任意操作区间 𝑡𝑖,⋯,𝑡𝑗（1≤𝑖≤𝑗≤𝑛）内拉伸系数 𝑘k 的乘积在 [0.001,1000] 范围内。

评分方式

如果你输出的浮点数与参考结果相比，满足绝对误差不大于 0.1，则该测试点满分，否则不得分。

提示

• C/C++：建议使用 double 类型存储浮点数，并使用 scanf("%lf", &x); 进行输入，printf("%f", x); 输出，也可以使用 cin 和 cout 输入输出浮点数；#include <math.h> 后可使用三角函数 cos() 和 sin()。

• Python：直接使用 print(x) 即可输出浮点数 x；from math import cos, sin 后可使用相应三角函数。

• Java：建议使用 double 类型存储浮点数，可以使用 System.out.print(x); 进行输出；可使用 Math.cos() 和 Math.sin() 调用三角函数。    

解题思路

1.如果你的思路是一个个挨个处理：打住！会超时！

2.直接处理行不通，那说明我们肯定是要把这些数据预处理一下再下手坐标变换的。如何预处理：

• 首先把复杂的操作简单化：输入的操作太多了，一个个处理得变换多久啊！但我们可以发现，这些操作无非就是对坐标进行拉伸和旋转两种变换，那么把这些操作整合一下，一共拉伸多少倍、一共旋转多少度，再对坐标进行总的拉伸和旋转，这行不行呢？疑问提出，开始验证，（这里就不赘述计算过程了），发现把坐标按照操作序列一步一步变换和把坐标进行整合拉伸变换所得结果一样，计划通！
• 由于题目给出了i，j限定了序列范围，由此，可以很快想到前缀和，开两个对象，分别存储序列累计的度数和倍数变换，这样后面就可以快速对坐标进行处理了。

3.注意细节:

• 凡是涉及到小数点的数据处理，最好使用double来做，精度更好。
• 题目给出的提示一定别忽略，比如此题，你用cin、cout输入输出的话，会比scanf、printf慢两倍多。

AC代码

//直接处理会超时 用前缀和做
#include<bits/stdc++.h>

using namespace std;

const int N = 1e5 + 10;

int n, m;

struct hhh{
    double type;    //操作类型
    double ko;      //k或θ值
    double bei;     //倍数前缀和
    double dushu;   //度数度数前缀和
}t[N];

int main()
{
    scanf("%d%d", &n, &m);
    
    t[0].bei = 1, t[0].dushu = 0; //倍数初始值为1 度数初始值为0
    for(int i = 1; i <= n; i ++)  //注意 这里从1开始 后面i、j就不用再减一处理了 同时这样也好算前缀和
    {
        scanf("%lf%lf", &t[i].type, &t[i].ko); 
        if(t[i].type == 1)   t[i].bei = t[i-1].bei * t[i].ko, t[i].dushu = t[i-1].dushu;
        else    t[i].dushu = t[i-1].dushu + t[i].ko, t[i].bei = t[i-1].bei;
    }
  
    while(m --)
    {
        int i, j;
        double x, y;
        scanf("%d%d%lf%lf", &i, &j, &x, &y); 
       
        x *= t[j].bei / t[i-1].bei;
        y *= t[j].bei / t[i-1].bei;
        
        double _cosx = x*cos(t[j].dushu - t[i-1].dushu) , _sinx = x * sin(t[j].dushu - t[i-1].dushu);
        double _cosy = y*cos(t[j].dushu - t[i-1].dushu) , _siny = y * sin(t[j].dushu - t[i-1].dushu);
        x = _cosx - _siny, y = _sinx + _cosy;
      
        printf("%.3lf %.3lf\n", x, y);
    }
    
    return 0;
}