Stupid cat & Doge （分形图）

【题目描述】

　　　　

【题目链接】

　　　　http://noi.openjudge.cn/ch0204/8463/

【算法】

　　　　为求等级N下的点的坐标可由几何关系找到其与等级N-1下对应点的关系，然后递归直至所有点的祖先（等级0）即可计算出坐标。

【代码】

 #include <bits/stdc++.h>
 using namespace std;
 typedef long long ll;
 int t;
 ll n,s,d;
 pair<ll,ll> calc(int state,ll num)
 {
     if(!state) return make_pair(,);
     ll side=<<(state-),sum=side*side;
     ll  rec=num/sum,pre_num=num%sum;
     pair<ll,ll> pr=calc(state-,pre_num);
     switch(rec) {
     case : swap(pr.first,pr.second); break;
     case : pr.first+=side; break;
     case : pr.first+=side,pr.second+=side; break;
     case : swap(pr.first,pr.second),pr.first=side-pr.first+,pr.second=side*-pr.second+;
     }
     return pr;
 }
 int main()
 {
     scanf("%d",&t);
     while(t--) {
         scanf("%d%lld%lld",&n,&s,&d);
         pair<ll,ll> p1=calc(n,s-);
         pair<ll,ll> p2=calc(n,d-);
         ll dx=p1.first-p2.first,dy=p1.second-p2.second;
         printf("%.0f\n",sqrt(dx*dx*1.0+dy*dy*1.0)*);
     }
     return ;
 }

【《算法竞赛进阶指南》大佬代码，侵删】

 /*
   Author: Yufei Du
   本程序仅供参考
 */
 #include <cmath>
 #include <cstdio>
 #include <cstring>
 #include <cstdlib>
 #include <algorithm>
 using namespace std;

 pair<__int64, __int64> recur (int stage, __int64 id)
 {
     if (stage == ) return make_pair(, );
     __int64 max =  << (stage - ), s = max * max;
     __int64 z = id / s, idsp = id % s;
     pair<__int64, __int64> psp = recur(stage - , idsp);
     if (z ==  || z == ) swap(psp.first, psp.second);
     if (z == )
     {
         psp.first = max - psp.first - ;
         psp.second = max - psp.second - ;
     }
     if (z ==  || z == ) psp.first += max;
     if (z ==  || z == ) psp.second += max;
     return psp;
 }
 int main ()
 {

     int kase; for (scanf("%d", &kase); kase; --kase)
     {
         int n; __int64 h, o;
         scanf("%d %I64d %I64d", &n, &h, &o);
         pair<__int64, __int64> hp = recur(n, h - );
         pair<__int64, __int64> op = recur(n, o - );
         __int64 dx = abs(hp.first - op.first), dy = abs(hp.second - op.second);
         printf("%.0f\n", (double)sqrt(dx * dx + dy * dy) * );
     }
     return ;
 }