#include #define PPque priority_queue, vector>, greater>> #define Pque priority_queue, vector>, greater>> #define pque priority_queue, greater> #define umap unordered_map #define uset unordered_set #define rep(i, s, f) for(ll i = s; i <= f; i++) #define per(i, s, f) for(ll i = s; i >= f; i--) #define all0(x) (x).begin() ,(x).end() #define all(x) (x).begin() + 1, (x).end() #define vvvvl vector #define vvvi vector>> #define vvvl vector>> #define vvvc vector>> #define vvvd vector>> #define vvvm vector>> #define vvi vector> #define vvl vector> #define vvs vector> #define vvc vector> #define vvp vector>> #define vvb vector> #define vvd vector> #define vvm vector> #define vp vector> #define vi vector #define vl vector #define vs vector #define vc vector #define vb vector #define vd vector #define vm vector #define P pair #define TU tuple #define rrr(l, r) mt()%(r-l+1)+l #define ENDL '\n' #define ull unsigned long long typedef long long ll; using namespace std; //////////////////////////////////////////////////////////////////////////////////////////////////////////// //これが本当の組み込み関数ってね(笑) template T or_less(vector &A, T x) { //x以下で最大要素の添字 前提: sort済み 存在しない: -1 return distance(A.begin(), upper_bound(A.begin(), A.end(), x)-1); } template T under(vector &A, T x) { //x未満の最大要素の添字 前提: sort済み 存在しない: -1 return distance(A.begin(), lower_bound(A.begin(), A.end(), x)-1); } template T or_more(vector &A, T x) { //x以上で最小要素の添字  前提: sort済み 存在しない: N . //distanceのA.beginは添字を出すために常にA.begin() NG: A.begin() + 1 return distance(A.begin(), lower_bound(A.begin(), A.end(), x)); } template T over(vector &A, T x) { //xより大きい最小要素の添字前提: sort済み 存在しない: N return distance(A.begin(), upper_bound(A.begin(), A.end(), x)); } template vector vec_shift_right(vector A, ll step, ll dir = 1, ll indexed = 1) {//dir = 1 : 右シフト dir = -1 : 左シフト ll N = A.size() - indexed; vector res(N+1); rep(i, indexed, N) { ll idx = i - step * dir; if(idx < indexed) idx += N; if(idx > N) idx -= N; res.at(i) = A.at(idx); } return res; } template void UNIQUE(vector &A) { sort(all0(A)); return A.erase(unique(A.begin(), A.end()), A.end()); } template void rev90(vector> &A, int indexed = 1) { reverse(A.begin() + indexed, A.end()); int n = A.size(); rep(i, indexed, n-1) { rep(j, i+1, n-1) { swap(A.at(i).at(j), A.at(j).at(i)); } } } void chmin(ll &a, ll b) { a = min(a, b); } void chmax(ll &a, ll b) { a = max(a, b); } ////////////////////////////////////////////////////////////////////// //数学系 /////////////////////////////////////////////////////////////////////// ll round(ll x, ll i) {return ll(x + 5 * pow(10, i-1))/ll(pow(10, i)) * ll(pow(10, i));} vp insu_bunkai(ll N) { vp res; for (ll i = 2; i * i <= N; i++) { ll cnt = 0; while(N % i == 0) { cnt++; N /= i; } if(cnt != 0) res.push_back(P(i, cnt)); } if(N != 1) res.push_back(P(N, 1)); return res; } ll extgcd (ll a, ll b, ll &x, ll &y) { if(b == 0) { x = 1;y = 0;return a;} ll d = extgcd(b, a%b, y, x); y -= a/b * x; return d; } template T ceil(T a, T b) { assert(b != 0); if(a % b == 0) return a / b; if((a <= 0 && b < 0) || (a >= 0 && b > 0)) return a/b + 1; else return a / b; } template T floor(T a, T b) { assert(b != 0); if(a % b == 0) return a / b; if((a <= 0 && b < 0) || (a >= 0 && b > 0)) return a/b; else return a/b - 1; } ll modpow(ll x, ll y, ll mod) { if(x > mod) x %= mod; if(y == 0) return 1; ll res = modpow(x, y >> 1, mod); res = res * res % mod; if(y & 1) res *= x, res %= mod; return res; } ll sqrt_(ll a) { ll l = 0; ll r = 3000000000LL; while(l < r) { ll mid = (l + r + 1) / 2; if(mid * mid <= a) l = mid; else r = mid - 1; } return l; } ////////////////////////////////////////////////////////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// //グローバル変数を置くところ(情報工学意識高め) #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") const ll int_max = 1001001001; const ll ll_max = 1001001001001001001LL; const double pi = 3.141592653589793; vl dx{0, 1, 0, -1, 0, 1, 1, -1, -1}; // (番兵) → ↓ ← ↑ ※ 右から時計回り  vl dy{0, 0, -1, 0, 1, 1, -1, -1, 1}; //const ll mod = 1000000007; //const ll mod = 998244353; ////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////// /* 外側のマスはせいぜい200選べばよう マスごとに9通り (200 * 10) * (200 * 10) = 4 * 10^6 処理に10^4かけるとNGなんだが、使うマスはせいぜい200マス程度なので気を使ってやる。 */ void solve() { ll H, W; cin >> H >> W; vvl A(H+1, vl(W+1, 0)); rep(i,1,H)rep(j,1,W) cin >> A.at(i).at(j); vvb used(H+1, vb(W+1, false)); ll sub = -ll_max; rep(i1, 0, H+1) { rep(j1, 0, W+1) { if(j1 >= 1 && (i1 != H+1 && i1 != 0)) break; rep(i2, i1, H+1) {//高さのみi1以下として良い、横には制約入れらんないはず。 rep(j2, 0, W+1) { if(j2 >= 1 && (i2 != H+1 && i2 != 0)) break; rep(dx1, -1, 1) rep(dy1, -1, 1) { if(dx1 == 0 && dy1 == 0) continue; rep(dx2, -1, 1) rep(dy2, -1, 1) { if(dx2 == 0 && dy2 == 0) continue; ll ans = 0; vp cnt; //1つ目のをまずたしあげる。 ll y1 = i1; ll x1 = j1; ll y2 = i2; ll x2 = j2; while(0 < y1 + dy1 && y1 + dy1 <= H && x1 + dx1 <= W &&0 < x1 + dx1) { y1 += dy1; x1 += dx1; if(!used[y1][x1]) ans += A.at(y1).at(x1), used.at(y1).at(x1) = true, cnt.push_back(P(y1, x1)); } while(0 < y2 + dy2 && y2 + dy2 <= H && x2 + dx2 <= W &&0 < x2 + dx2) { y2 += dy2; x2 += dx2; if(!used.at(y2).at(x2)) ans += A.at(y2).at(x2), used.at(y2).at(x2) = true, cnt.push_back(P(y2, x2)); } for(auto [y, x] : cnt) used.at(y).at(x) = false; chmax(sub, ans); } } } } } } cout << sub << endl; } //無断で0を平方数にカウントする人もいる //”部分文字列”と”連続部分文字列”は違うので確認すること //一般のグラフと、有向辺かつその貼り方に制約がある(多くの場合:番号がで解放に伸びる)はだいぶ違うので確認すること //座標を2で割った時の”切り捨て側(左側)”を求めるには 誤:(x / 2) マイナスの時!!! 正:floor(x, 2); //stringでの数字の下から1桁目は 正:S.at(N-1) 誤:S.at(0) //if(S.at(i) == 1) ← charなのに1...? // modは取りましたか...?(´・ω・`) //sortの比較関数は、 a == b ならば falseを返す必要がある(そうで無いとRE(発生しない場合もある)) int main() { ios::sync_with_stdio(false); std::cin.tie(nullptr); cout << fixed << setprecision(15); ll T = 1; //cin >> T; rep(i, 1, T) { solve(); } return 0; }