/*
    $$$$$$$\                      $$\           $$$$$$$\
    $$  __$$\                     \__|          $$  __$$\
    $$ |  $$ | $$$$$$\   $$$$$$\  $$\  $$$$$$$\ $$ |  $$ | $$$$$$\   $$$$$$\   $$$$$$$\ $$$$$$\
    $$$$$$$\ |$$  __$$\ $$  __$$\ $$ |$$  _____|$$$$$$$\ | \____$$\ $$  __$$\ $$  _____|\____$$\
    $$  __$$\ $$ /  $$ |$$ |  \__|$$ |\$$$$$$\  $$  __$$\  $$$$$$$ |$$ |  \__|$$ /      $$$$$$$ |
    $$ |  $$ |$$ |  $$ |$$ |      $$ | \____$$\ $$ |  $$ |$$  __$$ |$$ |      $$ |     $$  __$$ |
    $$$$$$$  |\$$$$$$  |$$ |      $$ |$$$$$$$  |$$$$$$$  |\$$$$$$$ |$$ |      \$$$$$$$\\$$$$$$$ |
    \_______/  \______/ \__|      \__|\_______/ \_______/  \_______|\__|       \_______|\_______|
*/
#include <bits/stdc++.h>

using namespace std;

#define PB push_back
#define MP make_pair
#define INS insert
#define LB lower_bound
#define UB upper_bound
#define pii pair <int,int>
#define pll pair <long long, long long>
#define si pair<string, int>
#define is pair<int, string>
#define X first
#define Y second
#define _ << " " <<
#define sz(x) (int)x.size()
#define all(a) (a).begin(),(a).end()
#define FOR(i, a, b) for (int i = (a); i < (b); ++i)
#define FORD(i, a, b) for (int i = (a); i > (b); --i)
#define FORR(i, l, r) for (int i = (l); i <= (r); ++i)
#define FORP(i, a, b) for ((i) = (a); (i) < (b); ++i)
#define FORA(i, x) for (auto &i : x)
#define REP(i, n) FOR(i, 0, n)
#define BITS(x) __builtin_popcount(x)
#define MSET memset
#define MCPY memcpy
#define SQ(a) (a) * (a)

typedef long long ll;
typedef long double ld;
typedef vector <int> vi;
typedef vector <pii> vpi;
typedef vector <ll> vll;
typedef vector <pll> vpl;
typedef vector <double> vd;
typedef vector <ld> vld;
typedef vector<si> vsi;
typedef vector<is> vis;
typedef vector<string> vs;
//((float) t)/CLOCKS_PER_SEC

const int MOD = 20183;
const int CX = 16807;
const int CY = 48271;
const int PI = acos(-1);
const int LOG = 32;
const int INF = 1e9 + 10;
const ll INFL = 1e18 + 10;
const int ABC = 30;
const int dx[] = {-1, 1, 0, 0};
const int dy[] = {0, 0, -1, 1};
const int dox[] = {-1, 1, 0, 0, -1, -1, 1, 1};
const int doy[] = {0, 0, -1, 1, -1, 1, -1, 1};

inline int sum(int a, int b){
  if (a + b >= MOD)
    return a + b - MOD;
  if (a + b < 0)
    return a + b + MOD;
  return a + b;
}

inline void add(ll &a, int b){
  a = sum(a, b);
}

inline ll mul(int a, int b){
  return ((ll)a * (ll)b) % MOD;
}

inline int sub(int a, int b){
    return (a - b + MOD) % MOD;
}

inline int pot(ll pot, int n){
    ll ret = 1;
    while (n){
        if (n & 1)
            ret = (ret * pot) % MOD;
        pot = (pot * pot) % MOD;
        n >>= 1;
    }
    return ret;
}

inline int divide(int a, int b){
    return mul(a, pot(b, MOD - 2));
}

ll lcm(ll a, ll b){
    return abs(a * b) / __gcd(a, b);
}

inline double ccw(pii A, pii B, pii C){
    return (A.X * B.Y) - (A.Y * B.X) + (B.X * C.Y) - (B.Y * C.X) + (C.X * A.Y) - (C.Y * A.X);
}

inline int CCW(pii A, pii B, pii C){
    double val = ccw(A, B, C);
    double eps = max(max(abs(A.X), abs(A.Y)), max(max(abs(B.X), abs(B.Y)), max(abs(C.X), abs(C.Y)))) / 1e9;
    if (val <= -eps)
        return -1;
    if (val >= eps)
        return 1;
    return 0;
}

void to_upper(string &x){
    REP(i, sz(x))
        x[i] = toupper(x[i]);
}

void to_lower(string &x){
    REP(i, sz(x))
        x[i] = tolower(x[i]);
}

string its(ll x){
    if (x == 0)
        return "0";
    string ret = "";
    while (x > 0){
        ret += (x % 10) + '0';
        x /= 10;
    }
    reverse(all(ret));
    return ret;
}

ll sti(string s){
    ll ret = 0;
    REP(i, sz(s)){
        ret *= 10;
        ret += (s[i] - '0');
    }
    return ret;
}

const int N = 1001;

ll g[N][N], sol[N][N][3], d = 5913, x = 701, y = 8;
char mat[N][N], e[3] = {'.', '=', '|'};
bool bio[N][N][3];

struct state{
    int x, y, g, d;
    state(int _x, int _y, int _g, int _d){
        x = _x, y = _y, g = _g, d = _d;
    }
};

class CMP{
    public:
        bool operator ()(state &a, state &b){
            return a.d > b.d;
        }
};

// 0 - climbing gear
// 1 - torch
// 2 - neither
// rocky - 0 or 1
// wet - 0 or 2
// narrow 1 or 2

void part1(){
    ll sol = 0;
    REP(i, x + 1){
        REP(j, y + 1){
            sol += g[i][j];
        }
    }
    cout << sol << '\n';
}

bool ok(int x, int y){
    return x > -1 && y > -1 && x < N && y < N;
}

ll dijkstra(){
    priority_queue <state, vector <state>, CMP> q;
    q.push(state(0, 0, 1, 0));
    while (!q.empty()){
        state s = q.top();
        q.pop();
        if (bio[s.x][s.y][s.g])
            continue;
        if (s.x == x && s.y == y && s.g == 1)
            return s.d;
        bio[s.x][s.y][s.g]++;
        REP(i, 4){
            int nx = s.x + dx[i],
                ny = s.y + dy[i];
            if (ok(nx, ny) && !bio[nx][ny][s.g] && g[nx][ny] != s.g)
                q.push(state(nx, ny, s.g, s.d + 1));
        }
        REP(i, 3){
            int nx = s.x,
                ny = s.y;
            if (ok(nx, ny) && !bio[nx][ny][i] && g[nx][ny] != i)
                q.push(state(nx, ny, i, s.d + 7));
        }
    }
}

void part2(){
    cout << dijkstra() << '\n';
}

int main () {
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);

    REP(i, N){
        REP(j, N){
            if (i == 0){
                if (j == 0)
                    g[i][j] = 0;
                else
                    g[i][j] = mul(j, CX);
            } else {
                if (j == 0)
                    g[i][j] = mul(i, CY);
                else
                    g[i][j] = mul(g[i - 1][j], g[i][j - 1]);
            }
            add(g[i][j], d);
        }
    }
    REP(i, N){
        REP(j, N){
            mat[i][j] = e[g[i][j] % 3];
            g[i][j] %= 3;
        }
    }
    g[x][y] = 0;
    part2();

    return 0;
}