Problem:
全探索お姉さんの休日

全探索お姉さんはとても優秀な女性である。
お姉さんは格子状の道の経路の数え上げを数千程度なら簡単に数え上げてしまう。

あなたと全探索お姉さんは今、六角形のタイルが敷き詰められた部屋にいる。
お姉さんは初めて見る六角形にとても興奮している様子である。
六角形の並びを座標に表すことに不慣れなお姉さんは図1のような座標系で部屋を表した。

図1

あなたはこの座標系上である地点からある地点まで移動したい。
しかし、お姉さんは 1 分ごとに動きたい方向をあなたに指示する。
いつものお姉さんなら同じ座標の場所を通らないようにあなたに移動の指示を出すだろう。
しかし、この座標系に不慣れなお姉さんは|x×y×t|(x：x 座標、y：y座標、t:最初の指示からの経過時間[分])を 6 で割った余り
に対応する方向（図で示す番号と対応）を指示するだけで、まったくのでたらめな方向にあなたを誘導する。

図2

お姉さんを傷つけたくないあなたは、お姉さんの指示をできるだけ守りつつ目的地までたどり着きたい。
あなたは許された行動は下の 7 つの行動である。

方向 0 へ 1 タイル移動

方向 1 へ 1 タイル移動

方向 2 へ 1 タイル移動

方向 3 へ 1 タイル移動

方向 4 へ 1 タイル移動

方向 5 へ 1 タイル移動

その場に留まる

お姉さんが指示を出した直後にあなたはこれら行動のうちの必ず1つを行う。

部屋には家具があり家具が配置されているタイルの中に移動することはできない。
また、y 座標の絶対値が ly を超えたり x 座標の絶対値が lx を超えるような移動は許されていない。
しかし、お姉さんはそのような移動を指示することがある。
指示を無視するとはお姉さんが示した方向と異なる方向へ移動するか、
もしくは、その場に留まることである。

目的地にたどり着くために最小で何度指示を無視すれば良いかを出力せよ。
目的地にたどり着くことが不可能な場合は -1 を出力せよ。

Input

入力は以下の形式で与えられる。

sx sy gx gy
n
x1 y1
...
xi yi
...
xn yn
lx ly

ここで、

sx, sy は出発地点の座標

gx, gy は目的地の座標

n は部屋に配置されている家具の数

xi, yi は家具があるタイルの座標

lx, ly は移動できる X, Y それぞれの座標の絶対値の上限
である。

Constraints

入力はすべて整数

-lx ≤ sx, gx ≤ lx

-ly ≤ sy, gy ≤ ly

(sx, sy) ≠ (gx, gy)

(xi, yi) ≠ (sx, sy)(1≤i≤n)

(xi, yi) ≠ (gx, gy)(1≤i≤n)

(xi, yi) ≠ (xj, yj)(i ≠ j)

0≤ n ≤ 1000

-lx ≤ xi ≤ lx (1≤i≤n)

-ly ≤ yi ≤ ly (1≤i≤n)

0 < lx, ly ≤ 100

Output

出力は1つの整数を含む1行で出力せよ。
目的地にたどり着ける場合は、最小の指示を無視する回数を出力せよ。
目的地にたどり着くことが不可能な場合は-1を出力せよ。

Sample Input 1

0 0 0 2
0
2 2

Output for the Sample Input 1

0

Sample Input 2

0 0 0 2
6
0 1
1 0
1 -1
0 -1
-1 -1
-1 0
2 2

Output for the Sample Input 2

-1

Sample Input 3

0 0 0 2
1
0 1
2 2

Output for the Sample Input 3

1

Solution in C++ (judge verdict: Time Limit Exceeded)

#include <iostream>
#include <algorithm>
#include <queue>
using namespace std;

const int INF = (1<<29);

const int dy[2][6] = {{1,0,-1,-1,-1,0},{1,1,0,-1,0,1}};
const int dx[6] = {0,1,-1,0,1,-1};
int wmin,hmin,wmax,hmax;

class state{
public:
  int y,x,t,cost;
  state(int y=0, int x=0, int t=0, int c=0):y(y),x(x),t(t),cost(c){}
  bool operator < (const state &s) const {
    return cost > s.cost;
  }
};

bool check(int y, int x){
  if(y < hmin || y > hmax) return false;
  if(x < wmin || x > wmax) return false;
  return true;
}

int main(){
  int sx,sy,gx,gy;
  while(cin >> sx >> sy >> gx >> gy){
    sx += 100;
    sy += 100;
    gx += 100;
    gy += 100;
    int n, h, w, a[6][201][201], ans=INF;
    state u(sy,sx,1,0), v;
    queue<state> Q;
    for(int i=0;i<6;i++) for(int j=0;j<201;j++) for(int k=0;k<201;k++) a[i][j][k] = INF;
    a[1][sy][sx] = 0;

    cin >> n;
    for(int i=0;i<n;i++){
      int x,y;
      cin >> x >> y;
      x += 100;
      y += 100;
      for(int j=0;j<6;j++) a[j][y][x] = -1;
    }
    cin >> w >> h;
    wmin = 100 - w;
    wmax = 100 + w;
    hmin = 100 - h;
    hmax = 100 + h;

    for(Q.push(u); !Q.empty(); Q.pop()){
      u = Q.front();
      if(u.y == gy && u.x == gx || ans <= u.cost){
	ans = min(ans,u.cost);
	continue;
      }
      
      int next = abs((u.x-100) * (u.y-100) * u.t) % 6;
      v = u;
      v.y += dy[u.x%2][next];
      v.x += dx[next];
      v.t++;
      if(check(v.y,v.x) && a[v.t%6][v.y][v.x] >= v.cost){
	a[v.t%6][v.y][v.x] = v.cost;
	Q.push(v);
      }

      v = u;
      v.t++;
      v.cost++;
      if(check(v.y,v.x) && a[v.t%6][v.y][v.x] > v.cost){
	a[v.t%6][v.y][v.x] = v.cost;
	Q.push(v);
      }

      for(int i=0;i<6;i++){
	v = u;
	v.y += dy[u.x%2][i];
	v.x += dx[i];
	v.t++;
	v.cost++;
	if(check(v.y,v.x) && a[v.t%6][v.y][v.x] > v.cost){
	  a[v.t%6][v.y][v.x] = v.cost;
	  Q.push(v);
	}
      }
    }

    if(ans == INF) cout << -1 << endl;
    else cout << ans << endl;
  }
}