c - Sierpinski三角形中的分形

时间:2012-10-02 13:54:08

标签: c fractals

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

##
#

2-sierpinski三角形

####
# #
##
# 

依此类推......这是一个找到它的样子的链接: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);

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

3 个答案:

答案 0 :(得分:2)

这是一些可能复杂但递归的代码!

#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);
}

答案 1 :(得分:0)

我认为你不需要递归。让#的三元组为1。所以n = No. of levels of the set的值应该打印一个低于另一个。在第一行中,打印设置n次。在下一行,n-1次,依此类推。迭代试试。
编辑:如果您正在寻找递归解决方案,请忽略我的答案。

答案 2 :(得分:0)

您可以使用pascal三角形编码。        1       1 1      1 2 1     1 3 3 1    1 4 6 4 1

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

打印普通三角形只计算与所需行数相关的空格数,并使用for循环(或其中几个)编码。检查哪个(pascal)编号对应哪个打印并跳过偶数。