【CSP】坐标变换2(问题拆解,快速输入,知识补充)

最新推荐文章于 2024-09-08 14:55:21 发布

Orange分享

最新推荐文章于 2024-09-08 14:55:21 发布

阅读量3.6k

点赞数 41

分类专栏： C++实用技术分享 CSP-J真题解决文章标签：算法数学建模

本文链接：https://blog.csdn.net/weixin_72391681/article/details/141719951

版权

C++实用技术分享同时被 2 个专栏收录

16 篇文章 0 订阅

订阅专栏

CSP-J真题解决

6 篇文章 0 订阅

订阅专栏

1. 题目背景与任务分析

题目背景
本题要求对二维平面上的点进行指定角度的旋转，并输出旋转后的坐标，要求精确到小数点后六位。这类题目广泛用于考察选手对数学计算、坐标变换以及编程语言中浮点数处理的能力。

任务明确

输入：多个坐标点及旋转角度。
输出：旋转后的新坐标，精确到小数点后六位。

分析与难点

几何计算：利用旋转矩阵进行坐标变换。
精度控制：确保输出结果满足精度要求，避免浮点数误差。
高效输入输出：优化大量数据的输入输出效率。

2. 问题拆解与功能模块设计

总问题：如何准确、高效地实现坐标旋转并满足精度要求？

子问题划分：

坐标旋转计算公式推导
- 任务：如何通过旋转矩阵公式计算旋转后的新坐标？
- 方案：使用三角函数计算旋转矩阵。
浮点数精度输出
- 任务：如何确保输出精确到小数点后六位？
- 方案：使用printf的格式化输出控制输出精度。
高效输入输出优化
- 任务：如何处理大量数据的输入输出？
- 方案：通过自定义快速输入函数，并设定合理的缓冲区大小。

3. 子问题的详细解决方案

子问题1：坐标旋转公式的推导与实现

公式推导：旋转矩阵的推导基于三角函数的和差公式。假设旋转前的点坐标为 (x,y)(x, y)(x,y)，旋转角度为 θ\thetaθ，则旋转后的坐标可以通过以下公式计算：
```
double x_new = x * cos(theta) - y * sin(theta);
double y_new = x * sin(theta) + y * cos(theta);
```

子问题2：浮点数精度的输出

精度控制：C++中可以通过printf("%.6f", value)来保证输出精确到小数点后六位。
格式说明：%.6f表示输出一个浮点数，精确到小数点后六位。
代码实现：
```
printf("%.6f %.6f\n", x_new, y_new);
```

子问题3：高效输入输出优化

手动设定缓冲区：通过设定合理大小的缓冲区，并使用fread批量读取数据，提高了I/O效率。

代码实现：

char buf[1 << 24], *p1 = buf, *p2 = buf;
#define getchar() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1 << 24, stdin), p1 == p2) ? EOF : *p1++)

inline int read() {
    int x = 0;
    char ch = getchar();
    while (ch < '0' || ch > '9') ch = getchar();
    while (ch >= '0' && ch <= '9') {
        x = x * 10 + ch - '0';
        ch = getchar();
    }
    return x;
}

4. 弧度制与角度制的转换

弧度与角度的关系
在C++编程中，sin、cos等函数使用弧度制。因此，若角度以度数给出，需要转换为弧度：

double theta_radian = theta_degree * M_PI / 180.0;

5. 常用数学函数与头文件

数学函数总结

函数名	功能	用法	头文件
`sin`	计算正弦值	`sin(弧度)`	`<cmath>`
`cos`	计算余弦值	`cos(弧度)`	`<cmath>`
`atan2`	计算两点之间的角度	`atan2(y, x)`	`<cmath>`
`sqrt`	计算平方根	`sqrt(x)`	`<cmath>`
`pow`	计算幂次	`pow(base, exp)`	`<cmath>`

6. `double`与`float`的比较与选择

主要区别：

精度：double为64位，float为32位，因此double能提供更高的精度。
内存占用：double占用8字节，float占用4字节。
使用场景：double适用于需要高精度计算的场景，如科学计算、金融分析等。

选择建议：在竞赛题目中，为了减少因精度不足导致的错误，推荐使用double。

7. 浮点数快速输入的实现与分析

在竞赛中，处理大规模输入输出是一个关键的优化点。为了实现快速的浮点数输入，我们可以手动设置缓冲区大小，并通过自定义的输入函数来提高效率。这一部分将详细解释这种方法的实现和原理，并对代码中的每个变量进行详细的说明。

浮点数快速输入的代码实现

#include <cstdio> // 包含标准I/O函数，如fread, printf

char buf[1 << 24], *p1 = buf, *p2 = buf;
#define getchar() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1 << 24, stdin), p1 == p2) ? EOF : *p1++)

// 自定义快速输入函数
inline double read_double() {
    double x = 0, factor = 1;
    char ch = getchar();
    while (ch < '0' || ch > '9') {
        if (ch == '-') factor = -1;
        ch = getchar();
    }
    while (ch >= '0' && ch <= '9') {
        x = x * 10 + ch - '0';
        ch = getchar();
    }
    if (ch == '.') {
        double frac = 1;
        ch = getchar();
        while (ch >= '0' && ch <= '9') {
            frac /= 10;
            x += (ch - '0') * frac;
            ch = getchar();
        }
    }
    return x * factor;
}

代码详解与变量说明

char buf[1 << 24], *p1 = buf, *p2 = buf;
- buf: 定义了一个字符数组，用作输入缓冲区。1 << 24表示将整数1左移24位，这相当于2^24，即16777216字节（大约16MB）。这是一个非常大的缓冲区，可以一次性存储大量的输入数据。
- p1: 指向缓冲区的当前读取位置。
- p2: 指向缓冲区的末尾，或者说是数据的读取结束位置。
- 作用: 通过定义这样一个大的缓冲区，可以减少fread函数的调用次数，从而提高数据读取的效率。
#define getchar()
- 定义: 自定义一个getchar函数，用来从缓冲区读取下一个字符。
- 实现原理:
  - 当p1和p2相等时，表示缓冲区的数据已经读取完毕，此时需要调用fread函数从标准输入（通常是键盘输入）再读取一批数据到缓冲区。
  - 如果fread返回的数据量为0，表示输入结束，getchar函数返回EOF（End of File，表示文件结束）。
  - 否则，p1自增，从缓冲区读取下一个字符并返回。
- 作用: 高效地读取输入，避免频繁地进行I/O操作。
fread(buf, 1, 1 << 24, stdin)
- 功能: 从标准输入（stdin）中读取数据到缓冲区buf中。
- 参数:
  - buf: 数据存储的目标缓冲区。
  - 1: 每次读取的单位大小为1字节。
  - 1 << 24: 总共读取的字节数，即16777216字节。
  - stdin: 指定读取的来源为标准输入。
- 原理: 一次性读取大量数据到缓冲区，从而减少I/O操作的开销。
- 应用场景: 适用于处理大规模数据输入的场景，如竞赛题目或需要高效读取大文件的情况。
inline double read_double()
- 功能: 快速读取一个浮点数。
- 过程:
  - 首先通过getchar跳过非数字字符（如空格、换行符）。
  - 检测到负号时，将factor设为-1以处理负数输入。
  - 读取整数部分，通过累加的方式将字符转换为对应的数值。
  - 如果遇到小数点，则进入小数部分的处理，逐位将小数部分转化并累加到x中。
- 变量作用:
  - x: 用来存储读取到的浮点数的值。
  - factor: 用来处理负数情况。
  - frac: 用于计算小数部分的分数值。
- 优势: 与标准的输入函数相比，这种方法减少了I/O操作次数，并且直接处理了小数部分的输入，效率更高。

理解1 << 24 表达式的意义

1 << 24 是一个位操作符，它将1的二进制位向左移动24位，结果是16777216。这个值表示缓冲区的大小，即我们一次性可以读取的最大字节数。

fread()函数的工作原理

fread 是C标准库中的一个函数，用来从文件或标准输入中读取数据。其工作方式是直接从I/O流中读取指定数量的字节并存储到缓冲区中，从而减少频繁的I/O操作。

EOF的意义

EOF (End of File) 是一个常量，用来标识输入流的结束。它通常在读取到输入的末尾或者在fread函数返回0字节时被返回。

8.手动设定缓冲区大小的方案 vs 使用系统规定大小的方案

1. 手动设定缓冲区大小的方案

优点:
- 控制更精确：可以根据具体情况设定缓冲区的大小，最大化内存使用效率。
- 大数据处理效率高：对于非常大的数据量，手动设定较大的缓冲区可以减少I/O操作次数，提高输入输出效率。
缺点:
- 内存占用大：手动设置的缓冲区大小如果过大，可能会导致内存占用过高，不适合在内存资源有限的环境中使用。
适用场景:
- 竞赛场景或者大数据处理场景，输入输出数据量巨大，且需要极高的I/O效率。

2. 使用系统规定大小的方案

优点:
- 内存使用灵活：缓冲区大小由系统管理，通常不会占用过多内存。
- 适合通用场景：对于一般的程序，这种方式已经足够满足需求，无需额外调优。
缺点:
- 效率稍低：由于系统缓冲区的大小不一定最优，因此在处理大规模数据时，效率可能不如手动设置缓冲区。
适用场景:
- 通用的编程场景，不需要极致优化输入输出的场合。

9. 最终代码与详细注释

#include <iostream>
#include <cmath>    // 包含数学函数头文件，如sin, cos, atan2, sqrt
#include <cstdio>   // 用于printf, scanf等C风格的I/O操作
using namespace std;

int main() {
    double x, y, theta;
    
    // 输入坐标(x, y)和旋转角度theta
    scanf("%lf %lf %lf", &x, &y, &theta);

    // 将角度转化为弧度制
    theta = theta * M_PI / 180.0;  // M_PI为π的值，直接转换为弧度
    
    // 计算cos(theta)和sin(theta)
    double cos_theta = cos(theta);
    double sin_theta = sin(theta);

    // 使用旋转矩阵计算新的坐标(x', y')
    double x_new = x * cos_theta - y * sin_theta;
    double y_new = x * sin_theta + y * cos_theta;

    // 按照题目要求，精确到小数点后6位输出结果
    printf("%.6f %.6f\n", x_new, y_new);

    return 0;
}

代码分析：