POJ3889 Fractal Streets

我对分治的理解：https://www.cnblogs.com/AKMer/p/9728574.html

题目传送门：http://poj.org/problem?id=3889

据说这种图叫分形图……一般找找规律就好了，像我这种菜鸡也只会做变化小一点的分形图了。

时间复杂度：\(O(能过)\)

空间复杂度：\(O(1)\)

代码如下：

#include <cmath>
#include <cstdio>
using namespace std;

int read() {
    int x=0,f=1;char ch=getchar();
    for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
    for(;ch>='0'&&ch<='9';ch=getchar())x=(x<<1)+(x<<3)+ch-'0';
    return x*f;
}

struct fuckpps {
    int x,y;
    fuckpps() {}
    fuckpps(int _x,int _y) {
        x=_x,y=_y;
    }
};

fuckpps find(int n,int angle,int id) {
    if(n==1) {
        if(angle==0) {
            if(id==1)return fuckpps(1,1);
            if(id==2)return fuckpps(1,2);
            if(id==3)return fuckpps(2,2);
            if(id==4)return fuckpps(2,1);
        }
        if(angle==1) {
            if(id==1)return fuckpps(1,1);
            if(id==2)return fuckpps(2,1);
            if(id==3)return fuckpps(2,2);
            if(id==4)return fuckpps(1,2);
        }
        if(angle==2) {
            if(id==1)return fuckpps(2,2);
            if(id==2)return fuckpps(2,1);
            if(id==3)return fuckpps(1,1);
            if(id==4)return fuckpps(1,2);
        }
        if(angle==3) {
            if(id==1)return fuckpps(2,2);
            if(id==2)return fuckpps(1,2);
            if(id==3)return fuckpps(1,1);
            if(id==4)return fuckpps(2,1);
        }
    }
    fuckpps res;
    int a=1<<(2*(n-1)),b=a*2,c=a*3;
    if(angle==0) {
        if(id<=a)res=find(n-1,1,id);
        if(id>a&&id<=b) {
            res=find(n-1,0,id-a);
            res.y+=1<<(n-1);
        }
        if(id>b&&id<=c) {
            res=find(n-1,0,id-b);
            res.x+=1<<(n-1);
            res.y+=1<<(n-1);
        }
        if(id>c) {
            res=find(n-1,3,id-c);
            res.x+=1<<(n-1);
        }
    }
    if(angle==1) {
        if(id<=a)res=find(n-1,0,id);
        if(id>a&&id<=b) {
            res=find(n-1,1,id-a);
            res.x+=1<<(n-1);
        }
        if(id>b&&id<=c) {
            res=find(n-1,1,id-b);
            res.x+=1<<(n-1);
            res.y+=1<<(n-1);
        }
        if(id>c) {
            res=find(n-1,2,id-c);
            res.y+=1<<(n-1);
        }
    }
    if(angle==2) {
        if(id<=a) {
            res=find(n-1,3,id);
            res.x+=1<<(n-1);
            res.y+=1<<(n-1);
        }
        if(id>a&&id<=b) {
            res=find(n-1,2,id-a);
            res.x+=1<<(n-1);
        }
        if(id>b&&id<=c)res=find(n-1,2,id-b);
        if(id>c) {
            res=find(n-1,1,id-c);
            res.y+=1<<(n-1);
        }
    }
    if(angle==3) {
        if(id<=a) {
            res=find(n-1,2,id);
            res.x+=1<<(n-1);
            res.y+=1<<(n-1);
        }
        if(id>a&&id<=b) {
            res=find(n-1,3,id-a);
            res.y+=1<<(n-1);
        }
        if(id>b&&id<=c)res=find(n-1,3,id-b);
        if(id>c) {
            res=find(n-1,0,id-c);
            res.x+=1<<(n-1);
        }
    }
    return res;
}//巨恶心的分情况讨论分治

int ABS(int a,int b) {
    if(a>b)return a-b;
    return b-a;
}

int sqr(int a) {
    return a*a;
}

int main() {
    int T=read();
    while(T--) {
        int n=read(),a=read(),b=read();
        fuckpps a1=find(n,0,a),a2=find(n,0,b);
        double gey=sqrt(sqr(ABS(a1.x,a2.x))+sqr(ABS(a1.y,a2.y)));
        int ans=trunc(10*gey+0.5);
        printf("%d\n",ans);
    }
    return 0;
}

原文地址：https://www.cnblogs.com/AKMer/p/9734953.html