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