#921122 (C++17) No.2509 Beam Shateki

• 問題: No.2509 Beam Shateki
• ユーザー: shobonvip
• 提出日時: 2023-10-20 21:58:31
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 6 ms / 2,000 ms
• コード長: 3,033 bytes
• コンパイル時間: 4,859 ms
• コンパイル使用メモリ: 265,700 KB
• 実行使用メモリ: 4,348 KB
• 最終ジャッジ日時: 2023-10-20 21:59:01
• 合計ジャッジ時間: 6,773 ms
• ジャッジサーバーID （参考情報）: judge14 / judge9

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	3 ms; 4,348 KB
testcase_01	AC	2 ms; 4,348 KB
testcase_02	AC	2 ms; 4,348 KB
testcase_03	AC	2 ms; 4,348 KB
testcase_04	AC	3 ms; 4,348 KB
testcase_05	AC	2 ms; 4,348 KB
testcase_06	AC	3 ms; 4,348 KB
testcase_07	AC	3 ms; 4,348 KB
testcase_08	AC	3 ms; 4,348 KB
testcase_09	AC	3 ms; 4,348 KB
testcase_10	AC	3 ms; 4,348 KB
testcase_11	AC	3 ms; 4,348 KB
testcase_12	AC	2 ms; 4,348 KB
testcase_13	AC	5 ms; 4,348 KB
testcase_14	AC	5 ms; 4,348 KB
testcase_15	AC	6 ms; 4,348 KB
testcase_16	AC	5 ms; 4,348 KB
testcase_17	AC	5 ms; 4,348 KB
testcase_18	AC	5 ms; 4,348 KB
testcase_19	AC	5 ms; 4,348 KB
testcase_20	AC	5 ms; 4,348 KB
testcase_21	AC	5 ms; 4,348 KB
testcase_22	AC	4 ms; 4,348 KB
testcase_23	AC	5 ms; 4,348 KB
testcase_24	AC	5 ms; 4,348 KB
testcase_25	AC	5 ms; 4,348 KB
testcase_26	AC	5 ms; 4,348 KB
testcase_27	AC	5 ms; 4,348 KB
testcase_28	AC	5 ms; 4,348 KB
testcase_29	AC	5 ms; 4,348 KB
testcase_30	AC	5 ms; 4,348 KB
testcase_31	AC	5 ms; 4,348 KB
testcase_32	AC	5 ms; 4,348 KB
testcase_33	AC	2 ms; 4,348 KB
testcase_34	AC	3 ms; 4,348 KB
testcase_35	AC	4 ms; 4,348 KB
testcase_36	AC	3 ms; 4,348 KB
testcase_37	AC	2 ms; 4,348 KB
testcase_38	AC	2 ms; 4,348 KB
testcase_39	AC	2 ms; 4,348 KB
testcase_40	AC	3 ms; 4,348 KB
testcase_41	AC	2 ms; 4,348 KB
testcase_42	AC	3 ms; 4,348 KB
testcase_43	AC	4 ms; 4,348 KB
testcase_44	AC	5 ms; 4,348 KB
testcase_45	AC	2 ms; 4,348 KB
testcase_46	AC	3 ms; 4,348 KB
testcase_47	AC	5 ms; 4,348 KB
testcase_48	AC	3 ms; 4,348 KB
testcase_49	AC	3 ms; 4,348 KB
testcase_50	AC	3 ms; 4,348 KB
testcase_51	AC	2 ms; 4,348 KB
testcase_52	AC	4 ms; 4,348 KB
testcase_53	AC	2 ms; 4,348 KB
testcase_54	AC	2 ms; 4,348 KB
testcase_55	AC	2 ms; 4,348 KB
testcase_56	AC	2 ms; 4,348 KB
testcase_57	AC	2 ms; 4,348 KB
testcase_58	AC	2 ms; 4,348 KB
testcase_59	AC	2 ms; 4,348 KB
testcase_60	AC	2 ms; 4,348 KB
testcase_61	AC	2 ms; 4,348 KB
testcase_62	AC	2 ms; 4,348 KB
testcase_63	AC	2 ms; 4,348 KB

ソースコード

#include<bits/stdc++.h>
using namespace std;

//* ATCODER
#include<atcoder/all>
using namespace atcoder;
typedef modint998244353 mint;
//*/

/* BOOST MULTIPRECISION
#include<boost/multiprecision/cpp_int.hpp>
using namespace boost::multiprecision;
//*/

typedef long long ll;

#define rep(i, s, n) for (int i = (int)(s); i < (int)(n); i++)
#define rrep(i, s, n) for (int i = (int)(n)-1; i >= (int)(s); i--)

template <typename T> bool chmin(T &a, const T &b) {
	if (a <= b) return false;
	a = b;
	return true;
}

template <typename T> bool chmax(T &a, const T &b) {
	if (a >= b) return false;
	a = b;
	return true;
}

template <typename T> T max(vector<T> &a){
	assert(!a.empty());
	T ret = a[0];
	for (int i=0; i<(int)a.size(); i++) chmax(ret, a[i]);
	return ret;
}

template <typename T> T min(vector<T> &a){
	assert(!a.empty());
	T ret = a[0];
	for (int i=0; i<(int)a.size(); i++) chmin(ret, a[i]);
	return ret;
}

template <typename T> T sum(vector<T> &a){
	T ret = 0;
	for (int i=0; i<(int)a.size(); i++) ret += a[i];
	return ret;
}

int main(){
	int h, w; cin >> h >> w;
	
	vector a(h, vector<int>(w));
	
	rep(i,0,h){
		rep(j,0,w){
			cin >> a[i][j];
		}
	}

	vector<int> tat(w);
	vector<int> yok(h);
	vector<int> na1(h+w-1);
	vector<int> na2(h+w-1);

	rep(i,0,w){
		rep(j,0,h){
			tat[i] += a[j][i];
		}
	}
	
	rep(i,0,h){
		rep(j,0,w){
			yok[i] += a[i][j];
		}
	}

	rep(i,0,h+w-1){
		rep(j,0,w){
			if (i-j >= 0 && i-j < h && j >= 0 && j < w){
				na1[i] += a[i-j][j];
			}
		}
	}

	rep(i,0,h+w-1){
		rep(j,0,w){
			if (h-1-i+j >= 0 && h-1-i+j < h && j >= 0 && j < w){
				na2[i] += a[h-1-i+j][j];
			}
		}
	}

	int ans = 0;

	rep(x,0,h){
		rep(y,x+1,h){
			chmax(ans, yok[x] + yok[y]);
		}
	}
	
	rep(x,0,w){
		rep(y,x+1,w){
			chmax(ans, tat[x] + tat[y]);
		}
	}

	rep(x,0,h+w-1){
		rep(y,x+1,h+w-1){
			chmax(ans, na1[x] + na1[y]);
		}
	}
	
	rep(x,0,h+w-1){
		rep(y,x+1,h+w-1){
			chmax(ans, na2[x] + na2[y]);
		}
	}

	rep(x,0,h){
		rep(y,0,w){
			chmax(ans, yok[x] + tat[y] - a[x][y]);
		}
	}

	rep(x,0,h){
		rep(y,0,h+w-1){
			int del = 0;
			// (y-j, j)
			// (x, y-x)
			if (0 <= y-x && y-x < w) del = a[x][y-x];
			chmax(ans, yok[x] + na1[y] - del);
		}
		rep(y,0,h+w-1){
			int del = 0;
			// (h-1-y+j, j)
			// (x, x+y+1-h)
			// h-1-y+j = x
			// j = x+y+1-h
			if (0 <= x+y+1-h && x+y+1-h < w) del = a[x][x+y+1-h];
			chmax(ans, yok[x] + na2[y] - del);
		}
	}

	
	rep(x,0,w){
		rep(y,0,h+w-1){
			int del = 0;
			// (y-j, j)
			// (y-x, x)
			if (0 <= y-x && y-x < h) del = a[y-x][x];
			chmax(ans, tat[x] + na1[y] - del);
		}
		rep(y,0,h+w-1){
			int del = 0;
			// (h-1-y+j, j)
			// (h-1-y+x, x)
			if (0 <= h-1-y+x && h-1-y+x < h) del = a[h-1-y+x][x];
			chmax(ans, tat[x] + na2[y] - del);
		}
	}

	rep(x,0,h+w-1){
		rep(y,0,h+w-1){
			int del = 0;
			// (x-j, j)
			// (h-1-y+k, k)
			// x-k = h-1-y+k
			// 2k = x+y+1-h
			int k = (x+y+1-h)/2;
			if ((x+y+1-h)%2 == 0 && 0 <= x-k && x-k < h && 0 <= k && k < w){
				del = a[x-k][k];
			}
			chmax(ans, na1[x] + na2[y] - del);
		}
	}

	cout << ans << '\n';
}