【问题标题】:A fractal in c - Sierpinsky trianglec - 谢尔宾斯基三角形中的分形
【发布时间】:2012-10-02 13:54:07
【问题描述】:

我正在尝试用 c 编写一个项目,该项目显示一个名为 Sierpinski 分形的分形(其中节点由“#”表示)。所以一个 1-sierpinski 三角形看起来像:

##
#

2-谢尔宾斯基三角形

####
# #
##
# 

等等...这里有一个链接可以找到它的样子:http://fr.wikipedia.org/wiki/Triangle_de_Sierpiński

有人告诉我,它可以在没有任何循环的情况下完成,只需通过递归方法。所以我尝试了类似的东西:

//extracting the power of two's index
int puiss_2(int N){
    int i=0,j=1;
    for(i=0;i<N;i++){
        j=j*2;
        i++;
    }
    return j;
}

//the recursive method
void fractal(int N)
{
    int M;
    M= puiss_2(N);


    if(M==0){
        printf("##\n");
        printf("# ");
    }
    else{
        fractal(N-1);
        fractal(N-1);
        printf("\n");
        fractal(N-1);
        printf(" ");
    }
}

int main()
{
    int N;
    scanf("%d",&N);
    fractal(N);
}

当然它没有用,因为当我跳到一条线时,我无法反转它。所以当我调用它两次时:

分形(N-1); 分形(N-1);

两个连续的动机并没有聚集在一起......有没有人知道如何 做那个?或者我的算法设计完全错了?

【问题讨论】:

  • 谷歌谢尔宾斯基三角形了解它的实际外观。
  • 需要递归地执行此操作吗?如果您没有现成的代码可以迭代地执行此操作:rosettacode.org/wiki/Sierpinski_triangle
  • @Mike :是的,我想递归地这样做
  • @JimBalter,是的,我要编辑我的帖子
  • 这里(作为第一个答案)是 Java 中的递归解决方案:stackoverflow.com/questions/8448908/…

标签: c fractals


【解决方案1】:

这里有一些可能很复杂但递归的代码!

#include <stdio.h>
#include <stdlib.h>
#include <string.h>

void sierpinsky(int N, char c[1000]){
    int i=0,j,k,l,born;

    for(i=0;i<N;i++){printf("%c",c[i]);}
    printf("\n");

    if(N==1){}
    else{
        if((c[0]=='#')&&(c[1]=='#')&&(c[2]=='#')){
            for (j=0;2*j<N;j++){
                if(c[2*j]=='#'){
                    c[2*j]='#';c[2*j+1]=' ';
                }
                else{
                    c[2*j]=' ';c[2*j+1]=' ';
                }
            }
        }
        else if ((c[0]=='#')&&(c[1]!='#')&&(c[2]=='#')){
                for (j=0;4*j<N;j++){
                    if(c[4*j]=='#'){
                        c[4*j]='#';c[4*j+1]='#';c[4*j+2]=' ';c[4*j+3]=' ';
                    }
                    else{
                        c[4*j]=' ';c[4*j+1]=' ';c[4*j+2]=' ';c[4*j+3]=' ';
                    }
                }
            }
        else if ((c[0]=='#')&&(c[1]!='#')&&(c[2] !='#')){
                k=0;
            while(c[k+1] !='#'){k++;}
            born = k+1;
            j=0;

            while(j<N){
                if((c[j]=='#')&&(c[j+born]=='#')){
                for(l=0;l<born;l++){
                    c[j+l]='#';
                    }
                    j=j+born+1;
                }

                else if ((c[j]!='#')&&(c[j-1+born]=='#')&&(c[j-1+2*born] !='#'))
                {
                    c[j-1]='#';
                    for(l=0;l<born;l++){
                        c[j+l]='#';
                    }
                    j=j+born+1;
                }
                else{
                    c[j-1]= ' ';
                    c[j]=' ';
                    j++;
                }
            }
        }
        else if ((c[0] =='#')&&(c[1] =='#')&&(c[2] !='#')){
            for (j=0;4*j<N;j++){
                if(c[4*j]=='#'){
                    c[4*j]='#';c[4*j+1]=' ';c[4*j+2]=' ';c[4*j+3]=' ';
                }
                else{
                    c[4*j]=' ';c[4*j+1]=' ';c[4*j+2]=' ';c[4*j+3]=' ';
                }
            }

        }
            else{}

            sierpinsky(N-1, c);
        }
    }

int main()
{   int i,size;
    scanf("%d",&size);
    char c[1000];
    for(i=0;i<size;i++){c[i]='#';}
    for(i=size;i<1000;i++){c[i]='a';}
    sierpinsky(size, c);
}

【讨论】:

  • 事实上,sierpinsky 函数中的每个条件都与每行的三个第一个元素的开始方式有关,它规定了该行的其余部分,对吧?
  • 没错,但肯定应该存在更简单的东西:)
【解决方案2】:

我认为您不需要为此进行递归。让# 的三元组为一组。所以 n = No. of levels of the set 的值你应该在另一个下面打印一个。在第一行,打印集合 n 次。在下一行,n-1 次,以此类推。反复尝试。
编辑:如果您正在寻找递归解决方案,请忽略我的回答。

【讨论】:

  • 是的,我试图递归地找到一些东西
【解决方案3】:

您可能可以使用帕斯卡三角形对此进行编码。 1 1 1 1 2 1 1 3 3 1 1 4 6 4 1

如果你可以用循环将这个三角形作为一个整体打印出来,那么,也许你可以跳过偶数。

要打印普通三角形,只需计算与您想要的行数相关的空格数,并使用 for 循环(或其中几个)进行编码。检查哪个(帕斯卡)数字对应于哪个打印并跳过偶数。

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2015-10-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多