3a三次方h c语言表达式,希尔伯特曲线——第八届蓝桥杯C语言B组（国赛）第三题...

原创

标题：希尔伯特曲线

希尔伯特曲线是以下一系列分形曲线 Hn 的极限。我们可以把 Hn 看作一条覆盖 2^n × 2^n 方格矩阵的曲线，曲线上一共有 2^n × 2^n 个顶点(包括左下角起点和右下角终点)，恰好覆盖每个方格一次。

[p1.png]

Hn(n > 1)可以通过如下方法构造：

1. 将 Hn-1 顺时针旋转90度放在左下角

2. 将 Hn-1 逆时针旋转90度放在右下角

3. 将2个 Hn-1 分别放在左上角和右上角

4. 用3条单位线段把4部分连接起来

对于 Hn 上每一个顶点 p ，我们定义 p 的坐标是它覆盖的小方格在矩阵中的坐标(左下角是(1, 1)，右上角是(2^n, 2^n)，从左到右是X轴正方向，从下到上是Y轴正方向)，

定义 p 的序号是它在曲线上从起点开始数第几个顶点(从1开始计数)。

以下程序对于给定的n(n <= 30)和p点坐标(x, y)，输出p点的序号。请仔细阅读分析源码，填写划线部分缺失的内容。

#include

long long f(int n, int x, int y) {

if (n == 0) return 1;

int m = 1 << (n - 1);

if (x <= m && y <= m) {

return f(n - 1, y, x);

}

if (x > m && y <= m) {

return 3LL * m * m + f(n - 1, ________________ , m * 2 - x + 1); // 填空

}

if (x <= m && y > m) {

return 1LL * m * m + f(n - 1, x, y - m);

}

if (x > m && y > m) {

return 2LL * m * m + f(n - 1, x - m, y - m);

}

}

int main() {

int n, x, y;

scanf("%d %d %d", &n, &x, &y);

printf("%lld", f(n, x, y));

return 0;

}

注意：只填写划线处缺少的内容，不要填写已有的代码或符号，也不要填写任何解释说明文字等。

(分别称下图三个网格称为A、B、C网格)

#include

long long f(int n, int x, int y) { //(x，y)只是用来判断在哪部分

if (n == 0) return 1;int m = 1 << (n - 1);if (x <= m && y <= m) { //左下角

returnf(n - 1, y, x); //返回给上一层的是此点在此层的序号

}if (x > m && y <= m) { //右下角

return 3LL * m * m + f(n - 1,m+1-y, 2 * m - x + 1); //填空

}if (x <= m && y > m) { //左上角

return 1LL * m * m + f(n - 1, x, y -m);

}if (x > m && y > m) { //右上角

return 2LL * m * m + f(n - 1, x - m, y -m);

}

}intmain() {intn, x, y;

scanf("%d %d %d", &n, &x, &y); //n为2的指数,(x,y)为顶点坐标

printf("%lld", f(n, x, y));return 0;

}/*得找出当前层的点对下一层相应点的映射关系*/

答案：m+1-y

要明确的是：每个网格中的1*1的方格都对应曲线上的一个顶点——即使有些1*1的方格内不存在转折点，但也是曲线上的一个顶点。

int m = 1 << (n - 1);/*左移1位相当于原数*2，左移2位相当于原数*pow(2,2);

所以 m==1*pow(2,n-1);

其实 m存储的就是边长为：pow(2,n) x pow(2,n) 的二分之一的长

比如：输入的 n = 3，网格规模为8x8，则 m==4*/

if (x <= m && y <= m) { //左下角

return f(n - 1, y, x); //返回给上一层的是此点在此层的序号

}if (x > m && y <= m) { //右下角

return 3LL * m * m + f(n - 1,m+1-y, 2 * m - x + 1); //填空

}if (x <= m && y > m) { //左上角

return 1LL * m * m + f(n - 1, x, y -m);

}if (x > m && y > m) { //右上角

return 2LL * m * m + f(n - 1, x - m, y -m);

}/*四个 if()中的 return 返回给上一层都是当前层顶点(x，y)的序号

根据曲线的走势，左下角是序号的第一部分，其实是左上角、右上角、右下角

比如我们要求左上角的顶点序号，所以我们应该要加上左下角的顶点数：1LL*m*m;

要求右下角的顶点序号，要加上前面三部分(左下角、左上角、右上角)的顶点数：3LL * m * m

XLL * m * m + f() (X=0/1/2/3,LL 代表长长整型) f()代表的是下一层传过来的点序号

所以 '+' 左边的部分加上右边的部分==(x，y)的序号*/

XLL * m * m + f()中 f()的两个参数的确定：

n-1传递给下一层，方便计算m，这里不解释；

另外两个参数分别对应下一层的(x，y)，即用这层的 x’，y‘表示下一层的 x，y ；

左上角、右上角、左下角中的参数容易确定，我们用待定系数法来确定右下角的两个参数；

if (x > m && y <= m) { //右下角

return 3LL * m * m + f(n - 1,m+1-y, 2 * m - x + 1); //填空

}

由 2 * m - x + 1 可知下一层的 y 和当前层的 x 是相关联的；

所以推测下一层的 x 和当前层的 y 是相关联的；

设当前层的坐标为 (x'，y')，下一层的坐标为(x，y)；

所以 y = 2 * m - x' + 1；即 x' = 2 * m - y + 1；(用下一层的坐标表示当前层的坐标)

可设 y’ = m*a + b*x + c；由于 m已知，所以可化简为 y' = b*x + F

我们从B、C两图中找出可以带入上式求出未知参数的坐标；

(1，1)——> (8，4)

(2，1)——> (8，3)

(3，1)——> (8，2)

左边的坐标从B图中任意取，右边的坐标要从C图的右下角部分取(我们的目标是右下角)；

比如B图中的序号为7的顶点的坐标为(2，4)；在C图右下角也找到序号为7的顶点，坐标为(5，3)；

将左边坐标的横坐标、右边坐标的纵坐标代入 y' = b*x + F得：

4=b+F;

3=2b+F;

解得 b = -1，F=5；

所以 m*a+c=5；从 x' 的表达式很容易就可以推测出  F = m + 1；

当然，我这种是不严格的解题方法！

22:31:53

2018-05-08