AcWing 刷题计划：0x00基础算法——AcWing 98. 分形之城

1.题目描述：


2.题意：
给一个n级城市，按照上述城市的建法，然后给你两个编号，问他们之间的直线距离（欧几里得距离）。

3.思路：
递归进行坐标变换。 我们以一级城市为例，我们以中心交界的地方作为原点建系（我们先以单位长度建系，最后再把系数乘回去），那么1、2、3、4号城市的坐标分别为：（-1,1）、（1,1）、（1，-1）、（-1，-1）。那么我们如果知道一级城市的坐标，对于二级城市，我们如法炮制。在二级城市中间建系，这时两个系的坐标是不一样的，我们如果分成4个象限来看，我们能发现，这四个象限和我们一级城市是相似的。（具体性质可以百度一下：分形）我们可以利用一定的变换使得我们的一级城市变成二级城市的一个象限。对于左上角的象限，我们可以看成是一级城市以原坐标系顺时针旋转再沿y轴翻转的结果。 那么变化公式就是：(-y,-x)。但是我们的坐标单位长度不一样怎么办呢？我们可以发现，坐标相差刚好是一号城市的二个城市的长度的差。那么我们就在原坐标下记录一下边长然后相加就好了。 其他的象限可以依次推导出来。所以我们只需要记录下当前城市的级数，以及对应城市的房屋数，房屋数方便我们把编号对应到n-1级城市，以免编号出错。这里注意的就是数据范围，可能爆int，所以常数要用long long类型。

4.代码：

//AcWing 98. 分形之城
//#include<bits/stdc++.h>
//#pragma GCC optimize(3,"Ofast","inline")
#include<iostream>
#include<string>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
//#include<random>
#include<cstdlib>
#include<ctime>
#include<map>
#include<stack>
#include<queue>
#define FAST ios::sync_with_stdio(false)
#define DEV_RND ((int)rand()*RAND_MAX+rand())
#define RND(L,R) (DEV_RND%((R)-(L)+1)+(L))
#define abs(a) ((a)>=0?(a):-(a))
#define sz(x) ((int)(x).size())
#define all(x) (x).begin(),(x).end()
#define mem(a,b) memset(a,b,sizeof(a))
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define rep(i,a,n) for(int i=a;i<n;++i)
#define per(i,n,a) for(int i=n-1;i>=a;--i)
#define pb push_back
#define mp make_pair
#define x first
#define y second
#define li inline
#define re register
using namespace std;
//typedef uniform_int_distribution<int> RNDI;
typedef pair<int,int> PII;
typedef vector<int> VI;
typedef double db;
typedef long long ll;
typedef long double ld;
const int maxn = 1e5+5;
const int maxm = 100000+5;
const int mod = 1e9+7;
const int inf=0x3f3f3f3f;
const double eps = 1e-7;
const double pi=acos(-1);
//int dir[4][2]={{-1,0},{1,0},{0,-1},{0,1}};
//li int f(int x){return x==par[x]?par[x]:par[x]=f(par[x]);}
//mt19937 eng(time(0));
//li int RND(int L,int R){RNDI rnd(L,R);return rnd(eng);}
li ll lowbit(ll x){return x&(-x);}
ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
void ex_gcd(ll a,ll b,ll &d,ll &x,ll &y){if(!b){d=a,x=1,y=0;}else{ex_gcd(b,a%b,d,y,x);y-=x*(a/b);}}//x=(x%(b/d)+(b/d))%(b/d);
li ll qpow(ll a,ll b,ll MOD=mod){ll res=1;a%=MOD;while(b>0){if(b&1)res=res*a%MOD;a=a*a%MOD;b>>=1;}return res%MOD;}
li ll qmul(ll a,ll b,ll MOD=mod){return (a*b-(ll)((long double)a/MOD*b)*MOD+MOD)%MOD;}
li ll inv(ll x,ll p){return qpow(x,p-2,p);}
li ll jos(ll n,ll k,ll s=1){ll res=0;rep(i,1,n+1) res=(res+k)%i;return (res+s)%n;}
db f(db x){return x;}
li db sim(db l,db r){return (f(l)+4.*f((l+r)/2.)+f(r))*(r-l)/6.;}
db asr(db l,db r,db ans,db eps){db m=l+(r-l)/2.,L=sim(l,m),R=sim(m,r);return fabs(L+R-ans)<=15.*eps?L+R+(L+R-ans)/15.:asr(l,m,L,eps/2)+asr(m,r,R,eps/2);}
db asr(db l,db r,db eps){return asr(l,r,sim(l,r),eps);}
namespace IO
{
	li ll read()
	{
		ll x=0,sign=1;char c=getchar();
		while(c>'9'||c<'0') {if(c=='-') sign=-1;c=getchar();}
		while('0'<=c&&c<='9') x=x*10+c-'0',c=getchar();
		return x*sign;
	}
	template<typename T>
	li void write(T x,char t='\n')
	{
		if(x<0){x=-x;putchar('-');};
		static int sta[25];int top=0;
		do{sta[top++]=x%10;}while(x/=10);
		while(top) putchar(sta[--top]+'0');
		putchar(t);
	}
}
using namespace IO;
/*-------------head-------------*/
//
typedef pair<ll,ll> PLL;
ll n,A,B;
PLL calc(int n,ll m)
{
	if(n==0) return mp(0,0);
	ll len=1ll<<(n-1),siz=1ll<<(2*n-2);
	//len:n-1级城市的边长;siz：n-1级城市的房屋数
	auto d=calc(n-1,m%siz);//当前编号在n-1级城市里面的坐标
	int x=d.x,y=d.y,u=m/siz;
	//x:n-1级坐标下的x坐标；y:n-1级坐标下的y坐标；u:第几象限
	if(u==0) return mp(-y-len,-x+len);//左上角
	if(u==1) return mp(x+len,y+len);//右上角
	if(u==2) return mp(x+len,y-len);//右下角
	if(u==3) return mp(y-len,x-len);//左下角
}
li void solve()
{
	n=read(),A=read(),B=read();
	auto a=calc(n,A-1);//-1：方便编号的运算
	auto b=calc(n,B-1);
	db x=a.x-b.x,y=a.y-b.y;
	db ans=sqrt(x*x+y*y)*5;
	printf("%.0f",ans);
	//puts("");
}
int main()
{
	//srand(time(0));
	//freopen("C:\\Users\\Administrator\\Desktop\\in.txt","r",stdin);
	for(int QwQ=read();QwQ;QwQ--) solve();
	//while(~scanf("%d",&n)) solve();
	return 0;
}