蓝桥杯-皮亚诺曲线距离

观察皮亚诺曲线，我们可以得到一下结论。
将整个图形分成九个区域，编号如下：

我们规定一号区域的方向为标准方向，则所有区域的方向为
1：x不翻转，y不翻转
2：x翻转，y不翻转
3：x不翻转，y不翻转
4：x不翻转，y翻转
5：x翻转，y翻转
6：x不翻转，y翻转
7：x不翻转，y不翻转
8：x翻转，y不翻转
9：x不翻转，y不翻转
知道了这些性质如何求出答案呢？
我们可以先求出x，y所在区域，则答案就是该区域之前的所有点数加上x，y到该子区域起点的距离。递归求解即可，递归时利用上面的性质，在传参数时对x，y进行翻转。
如何求x，y所在区间：
如果k的范围比较小，我们可以用len=3^(k-1),cnt=len*len分别表示一个子区域的边长和点数。
则x/len，y/len就是所在子区域的列数(0/1/2)和行数(0/1/2)。以此就可求出子区域编号。
但是，本题的k的范围是100。3^100很明显会爆long long。所有上述方法不可行。又观察到x,y的范围小于等于1e18。则我们可以暴力或者二分分别找到最小的使3^k1大于等于x、使3^k2大于等于y的k1和k2。这样就不存在爆longlong的问题。接着利用k1和k2求得所在区域编号，具体看代码。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
ll mi[39];
void init(){
    mi[0]=1;
    for(int i=1;i<39;i++)mi[i]=mi[i-1]*3;
}
ll get_(ll x){
    int l=0,r=38;
    while(l<r){
        int mid=l+r>>1;
        if(mi[mid]>x)r=mid;
        else l=mid+1;
    }
    return l;
}
ll f(ll n,ll x,ll y){
    if(n==0)return 0;
    ll k1=get_(x),k2=get_(y);
    if(k1<n){//子区域在左侧
        //子区域在下侧
        if(k2<n)return f(n-1,x,y);
        //子区域在中间
        if(mi[k2-1]*2>y)return mi[n-1]*mi[n-1]+f(n-1,mi[n-1]-x-1,y-mi[n-1]);
        //子区域在上侧
        return mi[n-1]*mi[n-1]*2+f(n-1,x,y-2*mi[n-1]);
    }
    else if(mi[k1-1]*2>x){//子区域在中间
        x-=mi[n-1];
        //子区域在下侧
        if(k2<n)return mi[n-1]*mi[n-1]*5+f(n-1,x,mi[n-1]-y-1);
        //子区域在中间
        if(mi[k2-1]*2>y)return mi[n-1]*mi[n-1]*4+f(n-1,mi[n-1]-x-1,mi[n-1]-(y-mi[n-1])-1);
        //子区域在上侧
        return mi[n-1]*mi[n-1]*3+f(n-1,x,mi[n-1]-(y-2*mi[n-1])-1);
    }
    else{//子区域在右侧
        x-=2*mi[n-1];
        //子区域在下侧
        if(k2<n)return mi[n-1]*mi[n-1]*6+f(n-1,x,y);
        //子区域在中间
        if(mi[k2-1]*2>y)return mi[n-1]*mi[n-1]*7+f(n-1,mi[n-1]-x-1,y-mi[n-1]);
        //子区域在上侧
        return mi[n-1]*mi[n-1]*8+f(n-1,x,y-mi[n-1]*2);
    }
}
int main(){
    init();
    ll n,x1,y1,x2,y2;
    cin>>n>>x1>>y1>>x2>>y2;
    cout<<abs(f(n,x1,y1)-f(n,x2,y2));
    return 0;
}