結果

問題 No.2509 Beam Shateki
ユーザー hibit_athibit_at
提出日時 2023-10-20 22:11:28
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
WA  
実行時間 -
コード長 5,877 bytes
コンパイル時間 1,575 ms
コンパイル使用メモリ 136,292 KB
実行使用メモリ 4,348 KB
最終ジャッジ日時 2023-10-20 22:11:39
合計ジャッジ時間 4,192 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge9
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 1 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 9 ms
4,348 KB
testcase_04 AC 9 ms
4,348 KB
testcase_05 AC 8 ms
4,348 KB
testcase_06 AC 9 ms
4,348 KB
testcase_07 AC 8 ms
4,348 KB
testcase_08 AC 8 ms
4,348 KB
testcase_09 AC 8 ms
4,348 KB
testcase_10 AC 9 ms
4,348 KB
testcase_11 AC 9 ms
4,348 KB
testcase_12 AC 9 ms
4,348 KB
testcase_13 AC 27 ms
4,348 KB
testcase_14 AC 28 ms
4,348 KB
testcase_15 AC 28 ms
4,348 KB
testcase_16 AC 28 ms
4,348 KB
testcase_17 AC 28 ms
4,348 KB
testcase_18 AC 28 ms
4,348 KB
testcase_19 AC 28 ms
4,348 KB
testcase_20 AC 28 ms
4,348 KB
testcase_21 AC 27 ms
4,348 KB
testcase_22 AC 28 ms
4,348 KB
testcase_23 AC 27 ms
4,348 KB
testcase_24 AC 28 ms
4,348 KB
testcase_25 AC 28 ms
4,348 KB
testcase_26 AC 28 ms
4,348 KB
testcase_27 AC 28 ms
4,348 KB
testcase_28 AC 29 ms
4,348 KB
testcase_29 AC 30 ms
4,348 KB
testcase_30 AC 28 ms
4,348 KB
testcase_31 AC 28 ms
4,348 KB
testcase_32 AC 27 ms
4,348 KB
testcase_33 AC 3 ms
4,348 KB
testcase_34 AC 12 ms
4,348 KB
testcase_35 AC 18 ms
4,348 KB
testcase_36 AC 8 ms
4,348 KB
testcase_37 AC 5 ms
4,348 KB
testcase_38 AC 8 ms
4,348 KB
testcase_39 AC 7 ms
4,348 KB
testcase_40 AC 12 ms
4,348 KB
testcase_41 AC 2 ms
4,348 KB
testcase_42 AC 8 ms
4,348 KB
testcase_43 AC 20 ms
4,348 KB
testcase_44 AC 22 ms
4,348 KB
testcase_45 AC 7 ms
4,348 KB
testcase_46 AC 10 ms
4,348 KB
testcase_47 AC 21 ms
4,348 KB
testcase_48 AC 10 ms
4,348 KB
testcase_49 AC 14 ms
4,348 KB
testcase_50 AC 15 ms
4,348 KB
testcase_51 AC 4 ms
4,348 KB
testcase_52 AC 12 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 2 ms
4,348 KB
testcase_57 AC 1 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 WA -
testcase_63 AC 2 ms
4,348 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <algorithm>
#include <iostream>
#include <iomanip>
#include <limits.h>
#include <map>
#include <math.h>
#include <numeric>
#include <queue>
#include <set>
#include <sstream>
#include <string>
#include <utility>
#include <vector>
#include <stack>
#include <complex>

using namespace std;
#define rep(i, n) for (ll i = 0; i < n; i++)
#define rep1(i, n) for (ll i = 1; i < n + 1; i++)
#define rev(i, n) for (ll i = n - 1; i >= 0; i--)
#define all(A) A.begin(), A.end()
#define itr(A, l, r) A.begin() + l, A.begin() + r
#define debug(var) cout << #var << " = " << var << endl;
typedef long long ll;

template <typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p)
{
    os << "(" << p.first << "," << p.second << ")";
    return os;
}

template <typename T1, typename T2>
istream &operator>>(istream &is, pair<T1, T2> &p)
{
    is >> p.first >> p.second;
    return is;
}

template <typename T>
ostream &operator<<(ostream &os, const vector<T> &v)
{
    for (ll i = 0; i < (ll)v.size(); i++)
    {
        os << v[i] << (i + 1 != (ll)v.size() ? " " : "");
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, const vector<vector<T>> &v)
{
    for (ll i = 0; i < (ll)v.size(); i++)
    {
        os << v[i] << endl;
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, const vector<vector<vector<T>>> &v)
{
    ll n = v.size();
    ll m = v[0].size();
    ll p = v[0][0].size();
    rep(k, p)
    {
        os << "k = " << k << endl;
        rep(i, n)
        {
            rep(j, m)
            {
                os << v[i][j][k];
                if (j < m - 1)
                {
                    os << " ";
                }
                else
                {
                    os << endl;
                }
            }
        }
    }
    return os;
}

template <typename T>
istream &operator>>(istream &is, vector<T> &v)
{
    for (T &in : v)
        is >> in;
    return is;
}

template <typename T, typename S>
ostream &operator<<(ostream &os, map<T, S> &mp)
{
    for (auto &[key, val] : mp)
    {
        os << key << ":" << val << " ";
    }
    cout << endl;
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, set<T> st)
{
    auto itr = st.begin();
    for (ll i = 0; i < (ll)st.size(); i++)
    {
        os << *itr << (i + 1 != (ll)st.size() ? " " : "");
        itr++;
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, multiset<T> st)
{
    auto itr = st.begin();
    for (ll i = 0; i < (ll)st.size(); i++)
    {
        os << *itr << (i + 1 != (ll)st.size() ? " " : "");
        itr++;
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, queue<T> q)
{
    while (q.size())
    {
        os << q.front() << " ";
        q.pop();
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, deque<T> q)
{
    while (q.size())
    {
        os << q.front() << " ";
        q.pop_front();
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, stack<T> st)
{
    while (st.size())
    {
        os << st.top() << " ";
        st.pop();
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, priority_queue<T> pq)
{
    while (pq.size())
    {
        os << pq.top() << " ";
        pq.pop();
    }
    return os;
}

template <typename T>
ostream &operator<<(ostream &os, priority_queue<T, vector<T>, greater<T>> mpq)
{
    while (mpq.size())
    {
        os << mpq.top() << " ";
        mpq.pop();
    }
    return os;
}
vector<vector<ll>> maze;
vector<ll> way;
vector<map<ll, ll>> sum(4);
map<ll, pair<ll, ll>> lim;
vector<ll> lefts;
vector<ll> rights;
ll ans = -2e9;
ll h, w;
void dfs()
{
    if (way.size() == 2)
    {
        ll p = way[0];
        ll q = way[1];
        // debug(vector<ll>({p, q}));
        auto [p_min, p_max] = lim[p];
        auto [q_min, q_max] = lim[q];
        for (ll i = p_min; i < p_max; i++)
        {
            for (ll j = q_min; j < q_max; j++)
            {
                if (p == q && i == j)
                {
                    continue;
                }
                ll tmp = sum[p][i] + sum[q][j];
                ll a, b, c, d;
                a = lefts[p];
                b = rights[p];
                c = lefts[q];
                d = rights[q];
                // debug(vector<ll>({i, j}));
                // debug(vector<ll>({a, b, c, d}));
                // debug(tmp);
                ll disc = a * d - b * c;
                // debug(disc);
                if (disc)
                {
                    ll u = (i * d - j * b) / disc;
                    ll v = (-i * c + j * a) / disc;
                    // debug(vector<ll>({u, v}));
                    if (0 <= u && u < h)
                    {
                        if (0 <= v && v < w)
                        {
                            tmp -= maze[u][v];
                        }
                    }
                }
                // debug(tmp);
                auto chmax = [](auto &a, auto b)
                { a = max(a, b); };
                chmax(ans,tmp);
            }
        }
        return;
    }
    rep(i, 4)
    {
        way.push_back(i);
        dfs();
        way.pop_back();
    }
}

int main()
{
    cin >> h >> w;
    lim[0] = {0, h};
    lim[1] = {0, w};
    lim[2] = {0, h + w - 1};
    lim[3] = {-w + 1, h};
    lefts = {1, 0, 1, 1};
    rights = {0, 1, 1, -1};
    maze.resize(h, vector<ll>(w));
    cin >> maze;
    // cout << maze << endl;
    rep(i, h)
    {
        rep(j, w)
        {
            sum[0][i] += maze[i][j];
            sum[1][j] += maze[i][j];
            sum[2][i + j] += maze[i][j];
            sum[3][i - j] += maze[i][j];
        }
    }
    // rep(i, 4)
    // {
    //     cout << sum[i] << endl;
    // }
    dfs();
    cout << ans << endl;
}
0