圣诞树2

本文分享一个有趣的圣诞树小实例：

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define PI 3.14159265359
float sx,sy;
float sdCircle(float px,float py,float r)
{
    float dx=px-sx,dy=py-sy;
    return sqrtf(dx*dx+dy*dy)-r;
}
float opUnion(float d1,float d2)
{
    return d1<d2?d1:d2;
}
#define T px+scale*r*cosf(theta),py+scale*r*sin(theta)
float f(float px,float py,float theta,float scale,int n)
{
    float d=0.0f;
    for(float r=0.0f;r<0.8f;r+=0.02f)
    d=opUnion(d,sdCircle(T,0.05f*scale*(0.95f-r)));
    if(n>0)
    for(int t=-1;t<=1;t+=2)
    {
        float tt=theta+t*1.8f;
        float ss=scale*0.9f;
        for(float r=0.2df;r<0.8f;r+=0.1f)
        {
            d=opUnion(d,f(T,tt,ss*0.5f,n-1));
            ss*=0.8f;
        }
    }
    return d;
}
int ribbon()
{
    float x=(fmodf(sy,0.1f)/0.1f-1.5f)*0.5f;
    return sx>=x-0.05f&&sx<=x+0.05f;
}
int main(int argc,char* argv[])
{
    int n=argc>1?atoi(argv[1]):3;
    float zoom=argc>2?atof(argv[2]):1.0f;
    for(sy=0.8f;sy>0.0f;sy-=0.02f/zoom,putchar('\n'))
    for(sx=-0.35f;sx<0.35f;sx+=0.01f/zoom)
    {
        if(f(0,0,PI*0.5f,1.0f,n)<0.0f)
        {
            if(sy<0.1f)
            putchar('.');
            else{
                if(ribbon())
                putchar('=');
                else putchar("............................#j&o"[rand()%32]);            
                }
        }
        else putchar(' ');
    }
}
/*int main(int argc,char* argv[])
{
    int n=argc>1?atoi(argv[1]):3;
    for(sy=0.8f;sy>0.0f;sy-=0.02f,putchar('\n'))
    for(sx=-0.35f;sx<0.35f;sx+=0.01f)
    putchar(f(0,0,PI*0.5f,1.0f,n)<0?'*':' '); 
    getchar();
}*/
/*int main(int argc,char* argv[])
{
    int n=argc>1?atoi(argv[1]):4;
    for(int j=1;j<=n;j++)
    {
        int s=1<<j,k=(1<<n)-s,x;
        for(int y=s-j;y>=0;j--,putchar('\n'))
        {
            for(x=0;x<y+k;x++)
            printf(" ");
            for(x=0;x+y<s;x++)
            printf("%c ",'!'^y&x);
            for(x=1;x+y<s;x++)
            printf("%c ",'!'^y&(s-y-x-1));
        }
    }
}