結果

問題 No.3605 Grand Cross
コンテスト
ユーザー yaaya
提出日時 2026-07-31 23:11:39
言語 C++23
(gcc 15.2.0 + boost 1.90.0)
コンパイル:
g++-15 -O2 -lm -std=c++23 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 754 ms / 2,000 ms
+ 944µs
コード長 2,610 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 2,548 ms
コンパイル使用メモリ 346,600 KB
実行使用メモリ 30,656 KB
最終ジャッジ日時 2026-07-31 23:11:56
合計ジャッジ時間 14,583 ms
ジャッジサーバーID
(参考情報)
judge2_0 / judge3_0
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 1
other AC * 49
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code

#include <bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(ll i=a;i<b;i++)
#define rrep(i,a,b) for(ll i=a-1;i>=b;i--)
#define ll long long
#define ull unsigned ll
#define ld long double
#define bl __int128_t
#define fi first
#define se second
#define vel vector<ll>
#define vvel vector<vel>
#define pll pair<ll,ll>
#define vepll vector<pll>
#define vvepll vector<vepll>
#define ves vector<string>
#define vem vector<mint>
#define vvem vector<vem>
#define pmm pair<mint,mint>
#define cleout(i) cout<<fixed<<setprecision(i)
template<class T>using PQ=priority_queue<T,vector<T>,greater<T>>;
//               上  右 下 左
vector<int> di={-1, 0, 1, 0};
vector<int> dj={ 0, 1, 0,-1};

vector<int> dx={ 0, 1, 0,-1};
vector<int> dy={ 1, 0,-1, 0};


vector<int> ddx={ 1, 1, 1, 0, -1, -1, -1, 0 };
vector<int> ddy={ 1, 0, -1, -1, -1, 0, 1, 1 };

ll inf=1000000000000000000;//1e18
// LLONG_MAX

mt19937_64 rng((ull)chrono::steady_clock::now().time_since_epoch().count());

//[x^M]1/(1-x)^N=comb(N-1+M,M)

ll get(vel &sum,ll l,ll r){
    ll ret=0;
    ret+=sum[r];
    if(l>0)ret-=sum[l-1];
    return ret;
}

void _solve(){
    ll N,M;
    cin>>N>>M;
    vel a(N),b(M),c(N),d(M);
    rep(i,0,N)cin>>a[i];
    rep(i,0,M)cin>>b[i];
    rep(i,0,N)cin>>c[i];
    rep(i,0,M)cin>>d[i];
    ll ans=-1;
    auto check=[&](){//左端a固定でbの真ん中を全探索すると同じ色のmaxlenが一意に決まる
        ll N=a.size();
        ll M=b.size();
        map<ll,vel> col;
        rep(i,0,N)col[c[i]].push_back(i);
        vel sa=a;
        vel sb=b;
        rep(i,1,N)sa[i]+=sa[i-1];
        rep(i,1,M)sb[i]+=sb[i-1];
        rep(i,0,M){
            if(col.count(d[i])){
                auto it=upper_bound(col[d[i]].begin(),col[d[i]].end(),min((N-(N%2==0))/2,min(i,M-i-1)));
                if(it==col[d[i]].begin()){
                    continue;
                }
                it--;
                //a[0:*it*2+1)とb[i:i+2*it+1)
                ans=max(ans,get(sa,0,*it*2)+get(sb,i-*it,i+*it));
            }
        }
    };
    //左端a固定
    check();
    //cout<<ans<<endl;

    reverse(a.begin(),a.end());reverse(c.begin(),c.end());
    //右端a固定
    check();
    //cout<<ans<<endl;

    swap(a,b);
    swap(c,d);
    //右端b固定
    check();
    //cout<<ans<<endl;

    reverse(a.begin(),a.end());reverse(c.begin(),c.end());
    //左端b固定
    check();

    cout<<ans<<"\n";
}


int main(){
    cin.tie(nullptr);
 	ios_base::sync_with_stdio(false);

    
    ll _;
    bool multitest=1;
    if(multitest)cin>>_;
    else _=1;
    rep(__,0,_){
        _solve();
    }
}
0