有趣的分形图形-递归和数学方法解决-POJ 2083

来源：互联网发布：淘宝家居用品类目编辑：程序博客网时间：2024/05/21 09:17

Description

A fractal is an object or quantity that displays self-similarity, in a somewhat technical sense, on all scales. The object need not exhibit exactly the same structure at all scales, but the same "type" of structures must appear on all scales.
A box fractal is defined as below :

A box fractal of degree 1 is simply
X
A box fractal of degree 2 is
X X
X
X X
If using B(n - 1) to represent the box fractal of degree n - 1, then a box fractal of degree n is defined recursively as following
```
B(n - 1)        B(n - 1)        B(n - 1)B(n - 1)        B(n - 1)
```

Your task is to draw a box fractal of degree n.

Input

The input consists of several test cases. Each line of the input contains a positive integer n which is no greater than 7. The last line of input is a negative integer −1 indicating the end of input.

Output

For each test case, output the box fractal using the 'X' notation. Please notice that 'X' is an uppercase letter. Print a line with only a single dash after each test case.

Sample Input

1234-1

Sample Output

X-X X XX X-X X   X X X     XX X   X X   X X    X   X XX X   X X X     XX X   X X-X X   X X         X X   X X X     X           X     XX X   X X         X X   X X   X X               X X    X                 X   X X               X XX X   X X         X X   X X X     X           X     XX X   X X         X X   X X         X X   X X          X     X         X X   X X            X X             X            X X         X X   X X          X     X         X X   X XX X   X X         X X   X X X     X           X     XX X   X X         X X   X X   X X               X X    X                 X   X X               X XX X   X X         X X   X X X     X           X     XX X   X X         X X   X X

一递归：

先看下较为常规的递归解决：

由于图形是重复的，小的图形只是把大图形的左上角一部分输出，，只计算最大的图形，打表可加快速度。

#include <stdio.h>#include <string.h>int p[8] = {1,3,9,27,81,243,729};char map[730][730];//n当前的图形大小，x，y图形所在的坐标void print(int n,int x,int y){if(n == 0){map[x][y] = 'X'; //return;}print(n-1, x, y); //左上print(n-1, x+2*p[n-1], y); //右上print(n-1, x+p[n-1], y+p[n-1]); //中间print(n-1, x, y+2*p[n-1]);print(n-1, x+2*p[n-1],  y+2*p[n-1]);}int n;int main(){for(int i=0; i<p[6]; i++) memset(map[i], 32, p[6]);print(6, 0, 0); //打表while(scanf("%d", &n) && n-- >= 0){for(int i=0; i<p[n]; i++){map[i][p[n]]=0;puts(map[i]);map[i][p[n]]=' ';}puts("-");}return 0;}

二数学方法

在discuss里面看到这个短小精悍的程序。

#include"stdio.h"#include"math.h"main(){int i,j,n,ii,jj,k;while(scanf("%d",&n)&&n--!=-1){for(i=0;i<pow(3,n);i++,printf("\n"))for(j=0;j<pow(3,n);j++){for(ii=i,jj=j,k=0;k<n&&(ii%3+jj%3)%2==0;ii/=3,jj/=3,k++);printf("%c",32+56*(k==n));}printf("-\n");} }

以下摘自 http://www.matrix67.com 博客：

关于这个图形，还可以用来证明

寻找1/5 + 1/25 + 1/125 + .. = 1/4的图形证明

前段时间，网上涌现出一大批关于1/4 + 1/16 + 1/64 + ... = 1/3的图形证明(1) (2)。不过，有多少人想过，为什么这些图形都是证明底数为1/4的情况呢？同样是几何级数求和，能否构造一个图形来证明1/5 + 1/25 + 1/125 + .. = 1/4呢？

无妨让我们来尝试一下。绝大多数人的第一想法便是画一个正五边形，然后把它分成五等分。接下来，把其中一份涂上颜色，表示整个面积的1/5。然后呢？然后怎么办？我们需要有一种办法把剩下的四块中的其中一块再次分成五份，选取一份表示1/25，并且递归地做下去。我们似乎发现了问题：这个“递归”是做不下去的。为了能够递归地表示出1/5、1/25、1/125，你必须把一个正五边形分成五个小的正五边形，再把小正五边形分成更小的五个正五边形，而这似乎是做不到的。这也就是上面的那些图形都可以用来证明Σ(1/4)^n=1/3的根本原因：它们可以把自己分成更小的四个自己，从而能够递归地表示(1/4)^n。
现在，我们的问题就出来了：为了弄出一个Σ(1/5)^n=1/4的图形证明，我们必须寻找这样一种图形，它能把自身分割为五块相同的部分，每一部分都和自身相似。

我们能否找到这样的图形呢？不能！根本原因在于：一个平面图形是二维的。平面图形的维度决定了这样一个性质：两个相似图形的相似比和面积比成平方关系。如果一个正方形的边长是另一个正方形边长的两倍，那前者的面积就是后者的四倍；圆的面积必然与半径r成平方关系，最多再乘上一个常系数。类似地，把一条曲线放大一倍，其长度也会跟着变大一倍；把一个三维图形放大一倍，其体积将变大到原来的八倍。因此，把一个平面图形分成四个和自身相似的图形是相当和谐的——面积变为原来的1/4，边长就应该变为原来的1/2，而两个1/2边长正好就拼成一个原边长。但是，如果要想把一个平面图形分成五等份，每一份的面积就应该是1/5，则对应边应该变为原来的1/sqrt(5)，这样显然无法既无重复又无遗漏地填满原图形的边。看来，图形证明Σ(1/5)^n=1/4似乎就不可能了。

思考问题时要善于用想象来替代放弃。二维空间中可以证明Σ(1/4)^n=1/3。类似地，一维空间中还可以轻易证明Σ(1/2)^n=1。要是有一个什么东西是log(5)/log(2)维的就好了——这样的话，对应边扩大到原来的2倍，“面积”就会变成原来的2^(log(5)/log(2))倍，也就是5倍。这就是说，在log(5)/log(2)维的“空间”里，一个“图形”恰好可以包含5个与自身相似的“图形”，并且新的“边长”正好能整除原来的“边长”。看到这里有人会哈哈大笑起来——这种维度会有吗？
有。很多分形图形都有一些极其怪异的性质，相似比的变化和其所占空间的变化不成整次幂的关系。例如，大名鼎鼎的分形图形Sierpinski三角形就是这样——相似比为1/2，所占空间之比为1/3。一个图形有二维图形的样子，却没有二维图形的性质。Hausdorff维度就是专门用来处理这类问题的。我们常常用Hausdorff维度来描述一个分形图形，比如Sierpinski三角形的Hausdorff维度就是log(3)/log(2)——所占面积变为原来的三倍时，对应边变为原来的两倍。

嘿！那么，是不是Sierpinski三角形就可以用来证明Σ(1/3)^n=1/2呢？对！一个Sierpinski三角形由三个与自身相似的小Sierpinski三角形组成，因此你可以递归地表示出(1/3)^n。画出整个Sierpinski三角形的1/3，以及另外1/3的其中1/3，以及1/3的1/3的1/3，这样无限画下去，你会很快看出来，总区域正好占了原来的一半。

类似地，把Sierpinski地毯放大到原来的三倍，整个图形所占空间就变成了原来的八倍。它的Hausdorff维度就是log(8)/log(3)。于是，我们就能用它来证明Σ(1/8)^n=1/7。

回到本文最初的问题，为了证明Σ(1/5)^n=1/4，我们只需要找到这样一个分形图形，其Hausdorff维度的分子为log(5)。这样的图形不仅存在，而且还不止一个。大家可以在这个页面里找到各种不同Hausdorff维度的图形。我找到了一个有趣的分形图形叫做Vicsek雪花，它的Hausdorff维度是log(5)/log(3)。

不断选出更小意义上的那个1/5，你会一眼看出，选出部分的总和就是整个图形的1/4。