Форум умных людей

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

typedef struct{double x,y;}cplx;

typedef struct{double size;cplx center,orient;}square;

cplx newCplx(double x,double y){
    cplx z;
    z.x=x;
    z.y=y;
    return z;
}

square newSquare(double size,cplx center,cplx orient){
    square sq;
    sq.size=size;
    sq.center=center;
    sq.orient=orient;
    return sq;
}

cplx add(cplx z1,cplx z2){return newCplx(z1.x+z2.x,z1.y+z2.y);}

cplx subt(cplx z1,cplx z2){return newCplx(z1.x-z2.x,z1.y-z2.y);}

cplx mult(cplx z1,cplx z2){
    return newCplx(z1.x*z2.x-z1.y*z2.y,z1.x*z2.y+z1.y*z2.x);
}

cplx mult1(cplx z,double k){
    return newCplx(z.x*k,z.y*k);
}

cplx w1,w2,z1,z2;
double xMin,xMax,yMin,yMax,dX,dY;
double k;

square fSq1(square sq){
    double size=sq.size*0.8;
    cplx orient=mult(sq.orient,w1);
    cplx center=add(sq.center,mult(mult1(sq.orient,sq.size),z1));
    return newSquare(size,center,orient);
}

square fSq2(square sq){
    double size=sq.size*0.6;
    cplx orient=mult(sq.orient,w2);
    cplx center=add(sq.center,mult(mult1(sq.orient,sq.size),z2));
    return newSquare(size,center,orient);
}

double fXMin(square sq,int l,int l0){
    if(l==l0||sq.center.x-k*sq.size>xMin)return sq.center.x;
    double x1=fXMin(fSq1(sq),l+1,l0);
    double x2=fXMin(fSq2(sq),l+1,l0);
    double x=fmin(sq.center.x,fmin(x1,x2));
    if(x<xMin)xMin=x;
    return x;
}

double fXMax(square sq,int l,int l0){
    if(l==l0||sq.center.x+k*sq.size<xMax)return sq.center.x;
    double x1=fXMax(fSq1(sq),l+1,l0);
    double x2=fXMax(fSq2(sq),l+1,l0);
    double x=fmax(sq.center.x,fmax(x1,x2));
    if(x>xMax)xMax=x;
    return x;
}

double fYMin(square sq,int l,int l0){
    if(l==l0||sq.center.y-k*sq.size>yMin)return sq.center.y;
    double y1=fYMin(fSq1(sq),l+1,l0);
    double y2=fYMin(fSq2(sq),l+1,l0);
    double y=fmin(sq.center.y,fmin(y1,y2));
    if(y<yMin)yMin=y;
    return y;
}

double fYMax(square sq,int l,int l0){
    if(l==l0||sq.center.y+k*sq.size<yMax)return sq.center.y;
    double y1=fYMax(fSq1(sq),l+1,l0);
    double y2=fYMax(fSq2(sq),l+1,l0);
    double y=fmax(sq.center.y,fmax(y1,y2));
    if(y>yMax)yMax=y;
    return y;
}

int main()
{
    w1=newCplx(0.8,0.6);
    w2=newCplx(0.6,-0.8);
    z1=mult(subt(newCplx(-0.5,0.5),mult(newCplx(-0.4,-0.4),w1)),newCplx(0,-1));
    z2=mult(subt(newCplx(0.5,0.5),mult(newCplx(0.3,-0.3),w2)),newCplx(0,-1));
    xMin=xMax=yMin=yMax=0;
    k=4.5*sqrt(2);
    square sq=newSquare(1,newCplx(0,0),newCplx(0,1));
    int l=0;
    double dX=0,dY=0,dX0,dY0;
    double eps=1e-14;
    do{
        l+=10;
        dX0=dX,dY0=dY;
        fXMin(sq,1,l),fXMax(sq,1,l),fYMin(sq,1,l),fYMax(sq,1,l);
        dX=xMax-xMin,dY=yMax-yMin;
        printf("L = %d\n",l);
        printf("xMin = %18.15f\n",xMin);
        printf("xMax = %18.15f\n",xMax);
        printf("yMin = %18.15f\n",yMin);
        printf("yMax = %18.15f\n",yMax);
        printf("  dX = %18.15f\n",dX);
        printf("  dY = %18.15f\n\n",dY);
    }while(fabs(dX-dX0)>eps||fabs(dY-dY0)>eps);
    return 0;
}