結果

問題 No.2114 01 Matching
ユーザー planes
提出日時 2022-11-03 19:39:08
言語 C++14
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 51
権限があれば一括ダウンロードができます

ソースコード

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