#define _CRT_SECURE_NO_WARNINGS #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; bool printb(bool f) { if (f)printf("Yes\n"); else printf("No\n"); return f; } template void prt(T t, string sep = "\n") { cout << t << sep; return; } template void printl(vector a, string sep = " ") { for (int i = 0; i < a.size(); i++) { cout << a[i] << sep; } cout << "\n"; return; } #define all(a) a.begin(),a.end() #define rep(i, n) for (int i = 0; i < (int)(n); i++) using llong = long long; using pii = pair; using pll = pair; using pli = pair; using pil = pair; template using vec2 = vector>; template inline bool chmin(T& a, T b) { return (a > b) ? (a = b, true) : false; } template inline bool chmax(T& a, T b) { return (a < b) ? (a = b, true) : false; } bool bitIn(llong a, int b) { return ((a >> b) & 1); } int bitCnt(llong a) { int re = 0; while (a > 0) { if (a & 1)re++; a >>= 1; } return re; } llong powL(llong n, llong i) { llong re = 1; while (i >= 1) { if (i & 1) re *= n; n *= n; i >>= 1; } return re; } llong powL_M(llong n, llong i, llong m) { llong re = 1; while (i >= 1) { if (i & 1) { re *= n; re %= m; } n *= n; n %= m; i >>= 1; } return re; } struct point { llong x = 0, y = 0; }; //lからrまでの和を返す template T sum_num(T l, T r) { if (((l + r) & 1) == 0) { return (l + r) / 2 * (l - r + 1); } else { return (l - r + 1) / 2 * (l + r); } } //int dx[4] = { 0,1,0,-1 }, dy[4] = { 1,0,-1,0 }; struct edge { int to, co; }; static const double pi = 3.14159265358979; /* modintクラス。四則演算と累乗が定義されている。 割り算はmodが素数でない時にも使える。(逆元の存在条件注意) extGCD(),GCD()を含む */ /* template T extGCD(T a, T b, T& x, T& y) { if (b == 0) { x = 1; y = 0; return a; } T gcd = extGCD(b, a % b, y, x); y -= a / b * x; return gcd; } template T GCD(T a, T b) { T x, y; return extGCD(a, b, x, y); } static const int mod = 1e9 + 7; //問題文に合わせて変更すること class modint { public: long long x; modint(long long x = 0) :x((x% mod + mod) % mod) {} modint operator-() const { return (-x); } modint& operator+=(const modint& a) { if ((x += a.x) >= mod)x -= mod; return *this; } modint& operator-=(const modint& a) { if ((x += mod - a.x) >= mod) x -= mod; return *this; } modint& operator*=(const modint& a) { (x *= a.x) %= mod; return *this; } modint operator+(const modint& a) const { modint res(*this); return res += a; } modint operator-(const modint& a) const { modint res(*this); return res -= a; } modint operator*(const modint& a) const { modint res(*this); return res *= a; } modint inv() const { long long y, c; extGCD(x, (long long)mod, y, c); return y; } modint& operator/=(const modint& a) { return (*this) *= a.inv(); } modint operator/ (const modint& a) const { modint res(*this); return res /= a; } friend ostream& operator<<(ostream& os, const modint& a) { os << a.x; return os; } }; //pow(a , n) modint型aのn乗のmodを求める template modint powM(modint a, T n) { modint re(1); while (n > 0) { if (n & 1)re *= a; a *= a; n >>= 1; } return re; } //*/ int main() { int h, w; cin >> h >> w; vec2 co(h - 2, vector(w + 2)); rep(i, h - 2) { co[i][0] = 0; rep(j, w) { scanf("%lld", &co[i][j + 1]); if (co[i][j + 1] == -1)co[i][j + 1] = 1e18; } co[i][w + 1] = 0; } vec2 dis(h - 2, vector(w + 2, 1e18)); dis[0][0] = 0; using T = tuple; priority_queue, greater> q; q.emplace(0, 0, 0); int dx[8] = { -1,-1,0,1,1,1,0,-1 }, dy[8] = { 0,-1,-1,-1,0,1,1,1 }; while (!q.empty()) { auto i = q.top(); q.pop(); auto d = get<0>(i); auto x = get<1>(i), y = get<2>(i); if (dis[x][y] < d)continue; //cout << x << " " << y << " " << d << "\n"; rep(k, 8) { int nx = x + dx[k], ny = y + dy[k]; if (nx < 0 || ny < 0)continue; if (nx >= h - 2 || ny >= w + 2)continue; //cout << " " << nx << " " << ny << "\n"; if (chmin(dis[nx][ny], d + co[nx][ny])) { q.emplace(d + co[nx][ny], nx, ny); } } } //for (auto i : dis)printl(i); if (dis[h - 3][w + 1] == 1e18)prt(-1); else prt(dis[h - 3][w + 1]); }