結果
問題 | No.2889 Rusk |
ユーザー | yuusaan |
提出日時 | 2024-09-13 21:59:56 |
言語 | C++23 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 139 ms / 2,000 ms |
コード長 | 7,316 bytes |
コンパイル時間 | 5,726 ms |
コンパイル使用メモリ | 324,260 KB |
実行使用メモリ | 50,260 KB |
最終ジャッジ日時 | 2024-09-13 22:00:16 |
合計ジャッジ時間 | 10,969 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,248 KB |
testcase_01 | AC | 2 ms
5,248 KB |
testcase_02 | AC | 2 ms
5,376 KB |
testcase_03 | AC | 2 ms
5,376 KB |
testcase_04 | AC | 2 ms
5,376 KB |
testcase_05 | AC | 2 ms
5,376 KB |
testcase_06 | AC | 3 ms
5,376 KB |
testcase_07 | AC | 2 ms
5,376 KB |
testcase_08 | AC | 2 ms
5,376 KB |
testcase_09 | AC | 2 ms
5,376 KB |
testcase_10 | AC | 2 ms
5,376 KB |
testcase_11 | AC | 2 ms
5,376 KB |
testcase_12 | AC | 2 ms
5,376 KB |
testcase_13 | AC | 2 ms
5,376 KB |
testcase_14 | AC | 8 ms
5,376 KB |
testcase_15 | AC | 24 ms
10,752 KB |
testcase_16 | AC | 58 ms
22,784 KB |
testcase_17 | AC | 17 ms
8,704 KB |
testcase_18 | AC | 49 ms
19,584 KB |
testcase_19 | AC | 97 ms
35,040 KB |
testcase_20 | AC | 114 ms
42,344 KB |
testcase_21 | AC | 94 ms
35,200 KB |
testcase_22 | AC | 80 ms
30,464 KB |
testcase_23 | AC | 66 ms
24,576 KB |
testcase_24 | AC | 117 ms
41,984 KB |
testcase_25 | AC | 65 ms
24,320 KB |
testcase_26 | AC | 93 ms
34,176 KB |
testcase_27 | AC | 117 ms
43,648 KB |
testcase_28 | AC | 69 ms
26,752 KB |
testcase_29 | AC | 72 ms
27,776 KB |
testcase_30 | AC | 71 ms
27,520 KB |
testcase_31 | AC | 121 ms
43,008 KB |
testcase_32 | AC | 50 ms
20,224 KB |
testcase_33 | AC | 15 ms
7,808 KB |
testcase_34 | AC | 136 ms
50,168 KB |
testcase_35 | AC | 136 ms
50,228 KB |
testcase_36 | AC | 136 ms
50,192 KB |
testcase_37 | AC | 135 ms
50,116 KB |
testcase_38 | AC | 135 ms
50,140 KB |
testcase_39 | AC | 94 ms
45,800 KB |
testcase_40 | AC | 73 ms
36,352 KB |
testcase_41 | AC | 104 ms
47,744 KB |
testcase_42 | AC | 68 ms
33,928 KB |
testcase_43 | AC | 33 ms
17,408 KB |
testcase_44 | AC | 88 ms
43,120 KB |
testcase_45 | AC | 67 ms
33,536 KB |
testcase_46 | AC | 20 ms
11,520 KB |
testcase_47 | AC | 36 ms
18,816 KB |
testcase_48 | AC | 79 ms
39,024 KB |
testcase_49 | AC | 2 ms
5,376 KB |
testcase_50 | AC | 2 ms
5,376 KB |
testcase_51 | AC | 2 ms
5,376 KB |
testcase_52 | AC | 2 ms
5,376 KB |
testcase_53 | AC | 2 ms
5,376 KB |
testcase_54 | AC | 139 ms
50,260 KB |
ソースコード
#include <iostream> #include <string> #include <stdio.h> #include <algorithm> #include <vector> #include <tuple> #include <map> #include<set> #include<queue> #include<stack> #include <unordered_set> #include<thread> #include<bits/stdc++.h> #include <atcoder/all> #include <cstdio> // #pragma GCC target("avx") // #pragma GCC optimize("O3") // #pragma GCC optimize("unroll-loops") //if(a < 0 || h <= a || b < 0 || w <= b)return; using namespace std; using namespace atcoder; using ll = long long; using ld = long double; using ull = unsigned long long; using mint = modint998244353; //using mint = modint1000000007; //using VL = vector<ll>; template<typename T> using pq = priority_queue<T>;//降順?(最大取り出し) template<typename T> using pqg = priority_queue<T, vector<T>, greater<T>>;//昇順?(最小取り出し) template<typename T> using vector2 = vector<vector<T>>; template<typename T> using vector3 = vector<vector<vector<T>>>; template<typename T> using vector4 = vector<vector<vector<vector<T>>>>; template<typename T> using vector5 = vector<vector<vector<vector<vector<T>>>>>; template<typename T> using vector6 = vector<vector<vector<vector<vector<vector<T>>>>>>; template<typename T> using pairs = pair<T,T>; #define rep(i, n) for (ll i = 0; i < ll(n); i++) #define rep1(i,n) for(int i = 1;i <= int(n);i++) #define repm(i, m, n) for (int i = (m); (i) < int(n);(i)++) #define repmr(i, m, n) for (int i = (m) - 1; (i) >= int(n);(i)--) #define rep0(i,n) for(int i = n - 1;i >= 0;i--) #define rep01(i,n) for(int i = n;i >= 1;i--) // NのK乗根N < 2^64, K <= 64 uint64_t kth_root(uint64_t N, uint64_t K) { assert(K >= 1); if (N <= 1 || K == 1) return N; if (K >= 64) return 1; if (N == uint64_t(-1)) --N; auto mul = [&](uint64_t x, uint64_t y) -> uint64_t { if (x < UINT_MAX && y < UINT_MAX) return x * y; if (x == uint64_t(-1) || y == uint64_t(-1)) return uint64_t(-1); return (x <= uint64_t(-1) / y ? x * y : uint64_t(-1)); }; auto power = [&](uint64_t x, uint64_t k) -> uint64_t { if (k == 0) return 1ULL; uint64_t res = 1ULL; while (k) { if (k & 1) res = mul(res, x); x = mul(x, x); k >>= 1; } return res; }; uint64_t res; if (K == 2) res = sqrtl(N) - 1; else if (K == 3) res = cbrt(N) - 1; else res = pow(N, nextafter(1 / double(K), 0)); while (power(res + 1, K) <= N) ++res; return res; } // ユークリッドの互除法による最大公約数算出 ll GCD(ll a,ll b){ if(b == 0)return a; return GCD(b, a % b); } //拡張ユークリッドの互除法による(ax + by = GCD(a,b))を満たすx,yの算出 pair<long long, long long> extgcd(long long a, long long b) { if (b == 0) return make_pair(1, 0); long long x, y; tie(y, x) = extgcd(b, a % b); y -= a / b * x; return make_pair(x, y); } struct UnionFind { vector<int> par; // par[i]:iの親の番号 (例) par[3] = 2 : 3の親が2 UnionFind(int N) : par(N) { //最初は全てが根であるとして初期化 for(int i = 0; i < N; i++) par[i] = i; } int root(int x) { // データxが属する木の根を再帰で得る:root(x) = {xの木の根} if (par[x] == x) return x; return par[x] = root(par[x]); } void unite(int x, int y) { // xとyの木を併合 int rx = root(x); //xの根をrx int ry = root(y); //yの根をry if (rx == ry) return; //xとyの根が同じ(=同じ木にある)時はそのまま par[rx] = ry; //xとyの根が同じでない(=同じ木にない)時:xの根rxをyの根ryにつける } bool same(int x, int y) { // 2つのデータx, yが属する木が同じならtrueを返す int rx = root(x); int ry = root(y); return rx == ry; } }; ll n; struct UnionFind2 { vector2<int> par; // par[i]:iの親の番号 (例) par[3] = 2 : 3の親が2 UnionFind2(int N) : par(N,vector<int>(N)) { //最初は全てが根であるとして初期化 for(int i = 0; i < N; i++)for(int j = 0;j < N;j++) par[i][j] = i * N + j; } int root(int y,int x) { // データxが属する木の根を再帰で得る:root(x) = {xの木の根} if (par[y][x] == y * n + x) return y * n +x; return par[y][x] = root(par[y][x] / n,par[y][x] % n); } void unite(int y0, int x0,int y1,int x1) { // xとyの木を併合 int rx = root(y0,x0); //xの根をrx int ry = root(y1,x1); //yの根をry if (rx == ry) return; //xとyの根が同じ(=同じ木にある)時はそのまま par[rx / n][rx % n] = ry; //xとyの根が同じでない(=同じ木にない)時:xの根rxをyの根ryにつける } bool same(int y0, int x0,int y1,int x1) { // 2つのデータx, yが属する木が同じならtrueを返す int rx = root(y0,x0); int ry = root(y1,x1); return rx == ry; } }; //座標圧縮 vector<ll> Ccomp(vector<ll> a){ vector<ll> b = a; sort(b.begin(),b.end()); b.erase(unique(b.begin(),b.end()),b.end());//ダブり消去 vector<ll> rtn; rep(j,a.size()){ ll pb = lower_bound(b.begin(),b.end(),a[j]) - b.begin(); rtn.push_back(pb); } return rtn; } /// ここから//////////////////////////////////////////// using F = ll; using S = ll; string s; ll modPow(ll a, ll n, ll mod) { if(mod==1) return 0;ll ret = 1; ll p = a % mod; while (n) { if (n & 1) ret = ret * p % mod; p = p * p % mod; n >>= 1; } return ret; } void cincout(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); cout<< fixed << setprecision(15); } //seg,遅延segの設定-----ここから S op(S a,S b){return a + b;}//何を求めるか(最大値とか) S e(){return 0;}//モノイド(初期値) S bop(S a,S b){return a * b;}; S be(){return 1;} ll v; //bool f(ll a){return v > a;}//めんどいやつ //S mapping (F a,S b){return a + b;}//遅延処理 F composition (F a,F b){return a + b;}//遅延中の枝にさらに処理 F id(){return 0;}//遅延のモノイド vector<int> Op(vector<int> a,vector<int> b){a.insert(a.end(),b.begin(),b.end()); return a;} vector<int> E(){return vector<int> (0);} //segここまで string abc = "abcdefghijklmnopqrstuvwxyz"; //string abc = "ABCDEFGHIJKLMNOPQRSTUVWXYZ"; ll mod = 998244353; ll INF =ll(2e18); int main() { cincout(); cin >> n; vector<ll> a(n+1),b(n+1),c(n+1); rep1(j,n)cin >> a[j]; rep1(j,n)cin >> b[j]; rep1(j,n)cin >> c[j]; vector<ll> d(n+1),e(n+1); rep1(j,n)d[j] = b[j] - a[j],e[j] = c[j] - a[j]; vector3<ll> dp(n+1,vector2<ll>(3,vector<ll>(3,-INF))); rep(j,3)rep(i,3)dp[0][j][i] = 0; rep1(j,n){ dp[j][0][0] = dp[j-1][0][0] + a[j]; dp[j][1][0] = max(dp[j-1][1][0],dp[j-1][0][0]) + b[j]; dp[j][2][0] = max(dp[j-1][1][0],dp[j-1][2][0]) + a[j]; dp[j][1][1] = max({dp[j-1][0][0],dp[j-1][1][0],dp[j-1][1][1]}) + c[j]; dp[j][2][1] = max({dp[j-1][1][0],dp[j-1][1][1],dp[j-1][2][1],dp[j-1][2][0]}) + b[j]; dp[j][2][2] = max({dp[j-1][1][1],dp[j-1][2][1],dp[j-1][2][2]}) + a[j]; } ll ans = -INF; rep(j,3)rep(i,3)ans = max(ans,dp[n][j][i]); cout << ans << endl; return 0; }