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

• 問題: No.2509 Beam Shateki
• ユーザー: fastmath
• 提出日時: 2023-10-20 21:29:50
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 170 ms / 2,000 ms
• コード長: 5,102 bytes
• コンパイル時間: 1,567 ms
• コンパイル使用メモリ: 154,652 KB
• 実行使用メモリ: 4,348 KB
• 最終ジャッジ日時: 2023-10-20 21:30:00
• 合計ジャッジ時間: 7,101 ms
• ジャッジサーバーID （参考情報）: judge13 / judge15

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 4,348 KB
testcase_01	AC	2 ms; 4,348 KB
testcase_02	AC	2 ms; 4,348 KB
testcase_03	AC	22 ms; 4,348 KB
testcase_04	AC	24 ms; 4,348 KB
testcase_05	AC	23 ms; 4,348 KB
testcase_06	AC	22 ms; 4,348 KB
testcase_07	AC	21 ms; 4,348 KB
testcase_08	AC	22 ms; 4,348 KB
testcase_09	AC	22 ms; 4,348 KB
testcase_10	AC	23 ms; 4,348 KB
testcase_11	AC	23 ms; 4,348 KB
testcase_12	AC	23 ms; 4,348 KB
testcase_13	AC	149 ms; 4,348 KB
testcase_14	AC	167 ms; 4,348 KB
testcase_15	AC	158 ms; 4,348 KB
testcase_16	AC	160 ms; 4,348 KB
testcase_17	AC	152 ms; 4,348 KB
testcase_18	AC	152 ms; 4,348 KB
testcase_19	AC	150 ms; 4,348 KB
testcase_20	AC	152 ms; 4,348 KB
testcase_21	AC	154 ms; 4,348 KB
testcase_22	AC	170 ms; 4,348 KB
testcase_23	AC	152 ms; 4,348 KB
testcase_24	AC	148 ms; 4,348 KB
testcase_25	AC	150 ms; 4,348 KB
testcase_26	AC	148 ms; 4,348 KB
testcase_27	AC	155 ms; 4,348 KB
testcase_28	AC	151 ms; 4,348 KB
testcase_29	AC	167 ms; 4,348 KB
testcase_30	AC	148 ms; 4,348 KB
testcase_31	AC	153 ms; 4,348 KB
testcase_32	AC	158 ms; 4,348 KB
testcase_33	AC	5 ms; 4,348 KB
testcase_34	AC	41 ms; 4,348 KB
testcase_35	AC	80 ms; 4,348 KB
testcase_36	AC	17 ms; 4,348 KB
testcase_37	AC	3 ms; 4,348 KB
testcase_38	AC	16 ms; 4,348 KB
testcase_39	AC	17 ms; 4,348 KB
testcase_40	AC	45 ms; 4,348 KB
testcase_41	AC	2 ms; 4,348 KB
testcase_42	AC	22 ms; 4,348 KB
testcase_43	AC	93 ms; 4,348 KB
testcase_44	AC	129 ms; 4,348 KB
testcase_45	AC	15 ms; 4,348 KB
testcase_46	AC	26 ms; 4,348 KB
testcase_47	AC	111 ms; 4,348 KB
testcase_48	AC	33 ms; 4,348 KB
testcase_49	AC	54 ms; 4,348 KB
testcase_50	AC	64 ms; 4,348 KB
testcase_51	AC	5 ms; 4,348 KB
testcase_52	AC	46 ms; 4,348 KB
testcase_53	AC	1 ms; 4,348 KB
testcase_54	AC	2 ms; 4,348 KB
testcase_55	AC	2 ms; 4,348 KB
testcase_56	AC	1 ms; 4,348 KB
testcase_57	AC	2 ms; 4,348 KB
testcase_58	AC	2 ms; 4,348 KB
testcase_59	AC	1 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 <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <vector>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <queue>
#include <ctime>
#include <cassert>
#include <complex>
#include <string>
#include <cstring>
#include <chrono>
#include <random>
#include <bitset>
#include <fstream>
#include <array>
#include <functional>
#include <stack>
#include <memory>
using namespace std;
#define int long long
#define ii pair <int, int>
#define app push_back
#define all(a) a.begin(), a.end()
#define bp __builtin_popcountll
#define ll long long
#define mp make_pair
#define x first
#define y second
#define Time (double)clock()/CLOCKS_PER_SEC
#define munq(a) sort(all(a)); a.resize(unique(all(a))-a.begin())
#define sz(a) ((int)a.size())

#ifdef LOCAL
#define debug(x) do { cout << #x << ": " << x << endl; } while(0)
#define debug2(x, y) do { std::cerr << #x << ", " << #y << ": " << x << ", " << y << endl; } while (0)
#define debug3(x, y, z) do {std::cerr << #x << ", " << #y << ", " << #z << ": " << x << ", " << y << ", " << z << endl; } while(0)
#else
#define debug(x)
#define debug2(x, y) 
#define debug3(x, y, z) 
#endif

#define FORI(i,a,b) for (int i = (a); i < (b); ++i)
#define FOR(i,a) FORI(i,0,a)
#define ROFI(i,a,b) for (int i = (b)-1; i >= (a); --i)
#define ROF(i,a) ROFI(i,0,a)
#define rep(a) FOR(_,a)
#define each(a,x) for (auto& a: x)
#define FORN(i, n) FORI(i, 1, n + 1)

using vi = vector<int>;
template <typename T>
std::istream& operator >>(std::istream& input, std::pair <T, T> & data)
{
    input >> data.x >> data.y;
    return input;
}
template <typename T>
std::istream& operator >>(std::istream& input, std::vector<T>& data)
{
    for (T& x : data)
        input >> x;
    return input;
}
template <typename T>
std::ostream& operator <<(std::ostream& output, const pair <T, T> & data)
{
    output << "(" << data.x << "," << data.y << ")";
    return output;
}
template <typename T>
std::ostream& operator <<(std::ostream& output, const std::vector<T>& data)
{
    for (const T& x : data)
        output << x << " ";
    return output;
}
ll div_up(ll a, ll b) { return a/b+((a^b)>0&&a%b); } // divide a by b rounded up
ll div_down(ll a, ll b) { return a/b-((a^b)<0&&a%b); } // divide a by b rounded down 
ll math_mod(ll a, ll b) { return a - b * div_down(a, b); }
#define tcT template<class T
#define tcTU tcT, class U
tcT> using V = vector<T>; 
tcT> bool ckmin(T& a, const T& b) {
    return b < a ? a = b, 1 : 0; 
} // set a = min(a,b)
tcT> bool ckmax(T& a, const T& b) {
    return a < b ? a = b, 1 : 0; 
}
tcT> vector <T> presum(vector <T> &a) {
    vector <T> p(a.size() + 1);
    FOR (i, a.size()) {
        p[i + 1] = p[i] + a[i];
    }   
    return p;
}
tcT> vector <T> sufsum(vector <T> &a) {
    vector <T> p(a.size() + 1);
    for (int i = (int)a.size() - 1; i >= 0; --i) {
        p[i] = p[i + 1] + a[i];
    }
    return p;
}
ll gcd(ll a, ll b) {
    while (b) {
        tie(a, b) = mp(b, a % b);
    }
    return a;
}
int Bit(int mask, int bit) { return (mask >> bit) & 1; }
signed main() {
    #ifdef LOCAL
    #else
    #define endl '\n'
    ios_base::sync_with_stdio(0); cin.tie(0);
    #endif
    int n, m;
    cin >> n >> m;
    vector <vector <int> > a(n, vector <int> (m));
    for (int i = 0; i < n; ++i) {
        for (int j = 0; j < m; ++j) {
            cin >> a[i][j];
        }
    }
    vector <vector <pair <int, int > > > l;
    for (int i = 0; i < n; ++i) {
        vector <pair <int, int > > nl;
        for (int j = 0; j < m; ++j) {
            nl.push_back({i,j});
        }
        l.push_back(nl);
    }

    for (int j = 0; j < m; ++j) {
        vector <pair <int, int > > nl;
        for (int i = 0; i < n; ++i) {
            nl.push_back({i,j});
        }
        l.push_back(nl);
    }

    for (int sum = 0; sum <= n + m; ++sum) {
        vector <pair <int, int > > nl;
        for (int i = 0; i < n; ++i) {
            int j = sum - i;
            if (0 <= j && j < m) {
                nl.push_back({i,j});
            }
        }
        l.push_back(nl);
    }

    for (int dif = -m; dif <= n; ++dif) {
        vector <pair <int, int > > nl;
        for (int i = 0; i < n; ++i) {
            int j = i - dif;
            if (0 <= j && j < m) {
                nl.push_back({i,j});
            }
        }
        l.push_back(nl);
    }

    for (auto &e : l) {
        sort(all(e));
    }
    int ans = 0;

    //cout << "ok" << endl;
    for (auto l1 : l) {
        for (auto l2 : l) {
            int ptr = 0;
            int sz = (int)l2.size();

            int sum = 0;
            for (auto p : l1) {
                sum += a[p.x][p.y];
            }
            for (auto p : l2) {
                sum += a[p.x][p.y];
            }

            for (auto p : l1) {
                while(ptr < sz && l2[ptr]<p) {
                    ptr++;
                }
                if (ptr<sz&&l2[ptr]==p) {
                    sum -= a[p.x][p.y];
                }
            }
            ckmax(ans,sum);
        }
    }
    cout<<ans<<endl;
}