GCJ Round 1C 2013 Problem B. Pogo

这题写了好久，大数据略可怕，其实是可以被优化的
甚至可以说，这题用最差复杂度为 1e9 级别的算法也无所谓
因为实际没那么多

算法：

显然步数的上界是 2(|dx| + |dy|) //来回折返，每次往目标方向前进1步
但是这个上界太差了
其实步数的上界可以减少到 3*sqrt(2(|dx|+|dy|)) //稍后再说
枚举最后一步的大小 : i
假设最后一步在y轴
将原本的 (dx, dy) 问题转化为 (dx, dy+i) 或 (dx, dy-i)
贪心的将问题转到 (dx, dy-i)
同理也可以转到 (dx-i, dy)
这两种情况如何选择呢？贪心地选择 |dx| 和 |dy| 中最大的那个
但是这题的贪心证明我不会，只是觉得应该如此

关于上界的优化：

我们假设最后一步的大小是 i，即用 i 步走到终点
可以从 1...k (k<=i) 不断的从原点不后退地往终点前进，直到不可前进
这时候，再多走两步一定会超过终点
所以就有：

$\frac{k(k+1)}{2} \leq |dx|+|dy|< \frac{(k+2)(k+3)}{2}$

即 $k<=\sqrt{2(|dx|+|dy|)}, k>=\sqrt{2(|dx|+|dy|)} -3$

然后在 k+1...i 的时候通过一次次往返，每次前进1个单位，最终逼近终点
可以得出

$\frac{k(k+1)}{2} + \frac{i-k}{2} \leq |dx|+|dy|$

$i <= 2(|dx|+|dy|) - k^2 <= 6\sqrt{2(|dx|+|dy|)} - 9$

其实从 1...k (k<=i) 前进是很浪费的
应该从 n, n-1, ... n-k 从大到小的前进，用小的数进行往返操作
这样能够减小一定的系数，不过没差
这题我用的上界是 $3\sqrt{2(|dx|+|dy|)}$

#include <iostream>
#include <map>
#include <set>
#include <algorithm>
#include <string>
#include <cstring>
#include <queue>
#include <cmath>
#include <cstdio>

using namespace std;

typedef long long ll;

const int maxn = 1000000+5;

int T;

int main() {
    freopen("a.in", "r", stdin);
    freopen("b.out", "w", stdout);

    cin>>T;
    for (int kase = 1; kase <= T; ++kase)
    {
        int tx, ty;
        cin>>tx>>ty;

        int dx, dy;
        int maxstep = 3*sqrt(2*(abs(tx)+abs(ty)));
        string path;

        int i;
        for (i = 1; i <= maxstep; ++i)
        {
            int k = i;
            path = "";
            dx = tx, dy = ty;
            while(k)
            {
                if (abs(dx) > abs(dy)){
                    if (dx > 0){
                        dx -= k;
                        path += 'E';
                    }
                    else
                    {
                        dx += k;
                        path += 'W';
                    }
                }
                else{
                    if (dy > 0){
                        dy -= k;
                        path += 'N';
                    }
                    else
                    {
                        dy += k;
                        path += 'S';
                    }   
                }
                k--;
            }
            if (dx==0 && dy==0)
                break;
        }
        printf("Case #%d: ", kase);
        for (int ii = i-1; ii >= 0; ii--)
            cout<<path[ii];
        puts("");
    }
}