康托尔三分集

1883年,德国数学家康托(G.Cantor)提出了如今广为人知的三分康托集,或称康托尔集。三分康托集是很容易构造的,然而,它却显示出许多最典型的分形特征。它是从单位区间出发,再由这个区间不断地去掉部分子区间的过程。

  三分康托集的构造过程是:
  第一步,把闭区间[0,1]平均分为三段,去掉中间的 1/3 部分段,则只剩下两个闭区间[0,1/3]和[2/3,1]。
  第二步,再将剩下的两个闭区间各自平均分为三段,同样去掉中间的区间段,这时剩下四段闭区间:[0,1/9],[2/9,1/3],[2/3,7/9]和[8/9,1]。
  第三步,重复删除每个小区间中间的 1/3 段。如此不断的分割下去, 最后剩下的各个小区间段就构成了三分康托集。

后一个间隔的开端是前一个间隔的后端加一个单位

#include <stdio.h>


int Pow(int a,int n){
    long long sum=1;
    for(int i=1;i<=n;i++){
        sum *= a;
    }
    return sum;
}

void  Contract(int x,int y){
    int m = x,n = y;
    if(x%y==0) printf("%d",x/y);
    else {
        while (n != 0) {
            int temp = m % n;
            m = n;
            n = temp;
        }
        printf("%d/%d",x/m,y/m);
    }
}

int  Contract1(int x,int y){
    int m = x,n = y;
    if(x%y==0) return x/y;
    else {
        while (n != 0) {
            int temp = m % n;
            m = n;
            n = temp;
        }
        return y/m;
    }
}

void Cantor(int begin,int end,int n){
    long long sum=1;
    int temp = begin,cnt=0,num=0,N=4,judge=0,t=0;
    int length = end -temp;
    sum  = Pow(2,n);
    for(int i=0;i<sum;i++){
        if(i==0) {
            printf("[%d,", temp);
            begin = temp*Pow(3,n) + length;
            Contract(begin,Pow(3,n));
            printf("]\n");
            cnt++;
        }
        else  {
            printf("[");
            Contract(begin,Pow(3,n));
            printf(",");
            Contract(begin+length,Pow(3,n));
            if(cnt!=1 ) begin += length ;
            printf("]\n");
            cnt++;
        }
        t = Contract1(begin+length,Pow(3,n));
        if(cnt==2) {
            num+=2;
            if(num%N==0){
                begin += length*Pow(3,n)/t ;
            }
            else begin += length*3;
            cnt = 0;
        }
        if(cnt!=2) begin += length ;

    }

}
int main()
{
    int n,begin,end;
    scanf("%d %d %d",&begin,&end,&n);
    Cantor(begin,end,n);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值