結果

問題 No.2114 01 Matching
ユーザー planesplanes
提出日時 2022-11-03 19:39:08
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 452 ms / 5,000 ms
コード長 2,388 bytes
コンパイル時間 2,106 ms
コンパイル使用メモリ 196,836 KB
実行使用メモリ 64,128 KB
最終ジャッジ日時 2024-11-24 12:26:58
合計ジャッジ時間 16,165 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
5,248 KB
testcase_01 AC 2 ms
5,248 KB
testcase_02 AC 236 ms
45,692 KB
testcase_03 AC 78 ms
15,872 KB
testcase_04 AC 114 ms
21,120 KB
testcase_05 AC 67 ms
18,640 KB
testcase_06 AC 230 ms
39,820 KB
testcase_07 AC 196 ms
29,812 KB
testcase_08 AC 210 ms
38,780 KB
testcase_09 AC 163 ms
25,096 KB
testcase_10 AC 203 ms
28,920 KB
testcase_11 AC 408 ms
62,720 KB
testcase_12 AC 425 ms
62,704 KB
testcase_13 AC 165 ms
41,820 KB
testcase_14 AC 213 ms
64,128 KB
testcase_15 AC 210 ms
64,124 KB
testcase_16 AC 291 ms
63,104 KB
testcase_17 AC 200 ms
29,780 KB
testcase_18 AC 128 ms
32,820 KB
testcase_19 AC 165 ms
42,416 KB
testcase_20 AC 418 ms
62,580 KB
testcase_21 AC 217 ms
51,928 KB
testcase_22 AC 140 ms
21,316 KB
testcase_23 AC 165 ms
41,028 KB
testcase_24 AC 43 ms
10,036 KB
testcase_25 AC 45 ms
12,480 KB
testcase_26 AC 452 ms
62,652 KB
testcase_27 AC 426 ms
62,660 KB
testcase_28 AC 229 ms
35,968 KB
testcase_29 AC 413 ms
62,592 KB
testcase_30 AC 252 ms
40,368 KB
testcase_31 AC 420 ms
62,720 KB
testcase_32 AC 205 ms
31,096 KB
testcase_33 AC 240 ms
45,704 KB
testcase_34 AC 189 ms
37,908 KB
testcase_35 AC 172 ms
27,784 KB
testcase_36 AC 421 ms
62,596 KB
testcase_37 AC 159 ms
28,156 KB
testcase_38 AC 214 ms
35,968 KB
testcase_39 AC 86 ms
15,012 KB
testcase_40 AC 155 ms
27,464 KB
testcase_41 AC 184 ms
28,432 KB
testcase_42 AC 152 ms
29,416 KB
testcase_43 AC 226 ms
39,956 KB
testcase_44 AC 264 ms
60,764 KB
testcase_45 AC 374 ms
62,520 KB
testcase_46 AC 60 ms
10,972 KB
testcase_47 AC 429 ms
62,464 KB
testcase_48 AC 420 ms
62,576 KB
testcase_49 AC 179 ms
42,168 KB
testcase_50 AC 74 ms
13,032 KB
testcase_51 AC 193 ms
30,424 KB
testcase_52 AC 216 ms
62,596 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

  #include <bits/stdc++.h> 
using namespace std;
using ll =long long;
#define all(v) v.begin(),v.end()
 #define rep(i,a,b) for(int i=a;i<b;i++)
#define rrep(i,a,b) for(int i=a;i>=b;i--)

ll INF=2e18;




int main() {

   

  ll N,M,K;cin>>N>>M>>K;
  vector<ll> B(N);
  vector<ll> A(M);
  for(ll i=0;i<N;i++) cin>>B[i];
  for(ll i=0;i<M;i++) cin>>A[i];

sort(all(B));
sort(all(A));

if(N<M) {
  swap(N,M);
vector<ll> C=B;
B=A;
A=C;
}


 vector<vector<ll>> b(N+M,vector<ll> (0));
 vector<vector<ll>> r(N+M,vector<ll> (0));

 ll ind=0;
 map<ll,ll> S;
 for(ll i=0;i<N;i++) {
  if(!S.count(B[i]%K)) {
    S[B[i]%K]=ind;
    ind++;
  }
    b[S[B[i]%K]].push_back(B[i]/K);
  }

ll ans=0;


for(ll i=0;i<M;i++) {
  if(!S.count(A[i]%K)) {
    S[A[i]%K]=ind;
    ind++;
  }
  r[S[A[i]%K]].push_back(A[i]/K);
}





for(ll w=0;w<N+M;w++) {
  if(b[w].size()<r[w].size()) {
    cout<<-1<<endl;
    return 0;
  }
  if(r[w].size()==0) continue;

  ll x=b[w].size();
  ll y=r[w].size();



ll now=y-x;
vector<pair<ll,ll>> vec(0);
for(auto p:b[w]) vec.push_back(make_pair(p,0));
for(auto p:r[w]) vec.push_back(make_pair(p,1));

vector<ll> R(x+y+1,0);
vector<ll> Rind(0);

sort(all(vec));




for(ll i=0;i<x+y;i++) {
  R[i]=now;
  if(vec[i].second==1){
    Rind.push_back(i);
     now--;
  }
  else now++;
}



vector<ll> d(x+y+1,0);

for(ll i=x+y-1;i>=0;i--) {
  d[i]=d[i+1];
  if(vec[i].second==1) d[i]-=vec[i].first;
  else d[i]+=vec[i].first;
}


ll g=max(x,y);
vector<vector<ll>> dp(y,vector<ll> (0));
vector<pair<ll,ll>> note(x+y+g+1,{INF,INF});

ll ri=-1;
now=y-x;




for(ll j=0;j<x+y;j++) {


note[now+g]=make_pair(j,ri);


  if(vec[j].second==1) {
    ll fb=0;
    for(ll h=j+1;h<x+y;h++) {
      if(vec[h].second==1) break;
      fb++;
    }

dp[ri+1]=vector<ll> (fb+1,INF);


ll k=R[j+1];


if(note[k+g].second!=INF) {
  dp[ri+1][0]=d[j+1]-d[note[k+g].first];
  if(note[k+g].second!=-1) {
    dp[ri+1][0]+=dp[note[k+g].second][note[k+g].first-Rind[note[k+g].second]-1];
  }
}



for(ll h=1;h<=fb;h++) {
   k=R[j+h+1];

   if(note[k+g].second!=INF) {
    dp[ri+1][h]=d[note[k+g].first]-d[j+h+1];


    if(note[k+g].second!=-1) {
      dp[ri+1][h]+=dp[note[k+g].second].back();
    }
   }


   dp[ri+1][h]=min(dp[ri+1][h],dp[ri+1][h-1]);
}


}




if(vec[j].second==1) {
  now--;
  ri++;

  }
else now++;
}



ans+=dp[y-1].back();

}







cout<<ans<<endl;
}
0