AtCoder Problems in C++ #E - 魚の生息範囲 (Fish)

第１２回日本情報オリンピック 予選（過去問）から、その提出コードの解答例です。
AtCoderとは、コンテストを通じて、プログラミングやアルゴリズムを学習するサービスです。
atcoder.jp
プログラミングコンテストとは、「与えられた問題をいかに素早く、正確に」解くことができるかを競うものです。
「競技プログラミング」を略して、「競プロ」などと呼ばれています。

#E - 魚の生息範囲 (Fish)

atcoder.jp

僕が作成、提出したコードは、以下のとおりです。

/*
AtCoder Problems in C++
#E - 魚の生息範囲 (Fish)
https://atcoder.jp/contests/joi2013yo/tasks/joi2013yo_e
提出コードの解答例
https://neguse-atama.hatenablog.com
*/
#include<bits/stdc++.h>
using namespace std;
using ll=long long;
vector<ll> vec_op(vector<ll>& x){
  vector<ll> y=x;
  sort(y.begin(),y.end());
  y.erase(unique(y.begin(),y.end()),y.end());
  for(ll& i : x){
    i=lower_bound(y.begin(),y.end(),i)-y.begin();
  }
  return y;
}
int main(void){
  ll n,K;
  cin>>n>>K;
  vector<ll> x(2*n),y(2*n),z(2*n);
  for(ll i=0;i<n;i++){
    cin>>x[i]>>y[i]>>z[i]>>x[n+i]>>y[n+i]>>z[n+i];
  }
  vector<ll> s=vec_op(x),t=vec_op(y),u=vec_op(z);
  ll a=s.size(),b=t.size(),c=u.size();
  vector<vector<vector<ll>>> r(a-1,vector<vector<ll>>(b-1,vector<ll>(c-1)));
  for(ll i=0;i<n;i++){
    for(ll j=x[i];j<x[n+i];j++){
      for(ll k=y[i];k<y[n+i];k++){
        for(ll l=z[i];l<z[n+i];l++){
          r[j][k][l]++;
        }
      }
    }
  }
  ll ans=0;
  for(ll i=0;i<a-1;i++){
    for(ll j=0;j<b-1;j++){
      for(ll k=0;k<c-1;k++){
        if(r[i][j][k]>=K){
          ans+=(s[i+1]-s[i])*(t[j+1]-t[j])*(u[k+1]-u[k]);
        }
      }
    }
  }
  cout<<ans<<endl;
  return 0;
}