【基本算法题-2022.7.31】11. 分形之城

Just Coding

已于 2022-07-31 00:03:53 修改

阅读量291

点赞数 1

分类专栏： # 每日亿题文章标签：算法贪心算法

于 2022-07-31 00:00:30 首次发布

本文链接：https://blog.csdn.net/weixin_46020266/article/details/126079591

版权

每日亿题专栏收录该内容

12 篇文章 1 订阅

订阅专栏

每日一题包含七大板块，现在从最基本的算法题开始，此类题型包括位运算、递推、递归、二分、排序、贪心等，从简单到复杂，跟我一起从点滴积累，到最终一举成名，打遍天下！✨

💻基本算法题11. 分形之城

城市的规划在城市建设中是个大问题。

不幸的是，很多城市在开始建设的时候并没有很好的规划，城市规模扩大之后规划不合理的问题就开始显现。

而这座名为 Fractal 的城市设想了这样的一个规划方案，如下图所示：

在这里插入图片描述

当城区规模扩大之后，Fractal 的解决方案是把和原来城区结构一样的区域按照图中的方式建设在城市周围，提升城市的等级。

对于任意等级的城市，我们把正方形街区从左上角开始按照道路标号。

虽然这个方案很烂，Fractal 规划部门的人员还是想知道，如果城市发展到了等级 N，编号为 A 和 B 的两个街区的直线距离是多少。

街区的距离指的是街区的中心点之间的距离，每个街区都是边长为 10 米的正方形。

输入格式
第一行输入正整数n，表示测试数据的数目。

以下n行，输入n组测试数据，每组一行。

每组数据包括三个整数 N,A,B, 表示城市等级以及两个街区的编号，整数之间用空格隔开。

输出格式
一共输出n行数据，每行对应一组测试数据的输出结果，结果四舍五入到整数。

数据范围
1≤N≤31,
1≤A,B≤22N,
1≤n≤1000

输入样例：

输出样例：

10 
30 
50

🍀题解 — 递归+分治+数学坐标系公式+找规律

在这里插入图片描述

// 该题解想法来自小xxx同学

📝代码展示

#include <iostream>
#include <cmath>
#define x first
#define y second
using namespace std;
typedef long long LL;
typedef pair<LL, LL> PLL;

PLL calc(LL n, LL m) {
    if(n == 0) return {0, 0};
    LL len = 1ll << (n - 1); // n - 1级图构成n级图平移时的单位长度
    LL cnt = 1ll << (n - 1) * 2; // n - 1级图中所含的元素个数
    LL cnk = m / cnt; // 在n级图中，编号为m的元素所属块的编号
    LL idx = m % cnt; // 在n级图中，编号为m的元素在所属块中的编号
    PLL pos = calc(n - 1, idx); // 在n级图中，编号为m的元素在所属块中的坐标
    LL x = pos.x, y = pos.y;
    // 根据n级图中某个块的编号和这个块中某个元素的坐标，确定n - 1级图的坐标变换
    // 注意：离散坐标系跟实数坐标系有些许差别，不考虑这些差别
    // 结果必然是抓耳挠腮。hh
    if(cnk == 0) return {y, x}; // 没有平移，虽是旋转，坐标关系容易确定
    if(cnk == 1) return {x, y + len}; // 单纯平移
    if(cnk == 2) return {x + len, y + len}; // 单纯平移
    if(cnk == 3) return {len * 2 - y - 1, len - x - 1}; // 旋转又平移
}

LL rounding(double a) {
    LL b;
    if(a > 0) b = (a * 2 + 1) / 2;
    else b = (a * 2 - 1) / 2;  
    return b;
}

int main() {
    int times;
    cin >> times;
    while(times--) {
        LL n, a, b;
        cin >> n >> a >> b;
        PLL pa = calc(n, a - 1);
        PLL pb = calc(n, b - 1);
        //double x = pa.x - pb.x, y = pa.y - pb.y;
        //printf("%.0lf\n", sqrt(x * x + y * y) * 10.00);
        LL x = pa.x - pb.x, y = pa.y - pb.y;
        double dist = sqrt(x * x + y * y);
        cout << rounding(dist) << endl;
    }
    return 0;
}

在这里插入图片描述

【🎯知识点：画图】

左上角：我们可以发现，左上角的N−1级矩阵其实就是等级为N−1，也就是上一个矩阵，顺时针旋转90°，(补充:旋转矩阵,属于大学的线性代数内容)，得出转移后的矩阵中的一点坐标从(x,y)变为(y,x)
左下角：同左上角，它则是逆时针旋转90°而且还要水平翻转，也即是沿着X轴对称，原本逆时针后为(y,−x)，然后要对称,x坐标不变，y坐标取反，所以坐标为(−y,−x) 也就是所谓的(2×len−1−y,len−1−x) 最难理解的坐标，具体可以画图理解
右上角和右下角：通过N=2级图发现，其实和N=1是一样的，并没有旋转，只是平移，则右上角坐标为(x,y+len)，右下角坐标为(x+len,y+len)
总的来说以上四种转移，都可以通过画图理解