Problem:
魚の生息範囲 (Fish)

問題

オーストラリア大陸の西には，広いインド洋が広がっている．海洋研究者である JOI 氏は，インド洋に生息しているある N 種類の魚の性質について研究している．

それぞれの魚の種類に対して，海の中に直方体状の生息範囲が定まっている．魚は境界も含めて生息範囲の中のどの場所にも移動できるが，生息範囲の外に出ることは決してない．海の中の点は，3 つの実数 (x, y, d) によって表される： (x, y, d) は，上空から見たときにある地点を基準にして東に x，北に y 進んだ位置であり，海面からの深さが d の点を表す．ただし，海面は平面であるとする．

JOI 氏は，K 種類以上の魚の生息範囲が重なる場所がどのくらいあるかを知りたい．そのような場所全体の体積を求めるプログラムを作成せよ．

入力

入力は 1 + N 行からなる．

1 行目には，2 つの整数 N, K (1 ≦ K ≦ N ≦ 50) が空白を区切りとして書かれている．これは，魚が N 種類であり，K 種類以上の魚の生息範囲が重なる場所の体積を求めたいことを表す．

続く N 行のうちの i 行目 (1 ≦ i ≦ N) には，6 つの整数 Xi,1, Yi,1, Di,1, Xi,2, Yi,2, Di,2 (0 ≦ Xi,1 ＜ Xi,2 ≦ 1000000 (= 106)，0 ≦ Yi,1 ＜ Yi,2 ≦ 1000000 (= 106)，0 ≦ Di,1 ＜ Di,2 ≦ 1000000 (= 106)) が書かれている．これは，i 種類目の魚の生息範囲が

8 点 (Xi,1, Yi,1, Di,1), (Xi,2, Yi,1, Di,1), (Xi,2, Yi,2, Di,1), (Xi,1, Yi,2, Di,1), (Xi,1, Yi,1, Di,2), (Xi,2, Yi,1, Di,2), (Xi,2, Yi,2, Di,2), (Xi,1, Yi,2, Di,2) を頂点とする直方体であることを表す．

出力

K 種類以上の魚の生息範囲が重なる場所全体の体積を 1 行で出力せよ．

入出力例

入力例 1

3 2
30 50 0 50 70 100
10 20 20 70 90 60
40 60 20 90 90 70

出力例 1

49000

入出力例 1 において，例えば，点 (45, 65, 65) は 1 種類目の魚と 3 種類目の魚の生息範囲であるので，条件を満たす場所である．一方，点 (25, 35, 45) は 2 種類目の魚のみの生息範囲であるので，条件を満たす場所ではない．また，魚の生息範囲は下の図のようになっている．点 O は海面上の基準の地点を表す．

入力例 2

1 1
0 0 0 1000000 1000000 1000000

出力例 2

1000000000000000000

問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++11 (judge verdict: Accepted)

#include <bits/stdc++.h>
using namespace std;

typedef long long ll;
#define rep(i,n) for(int (i)=0;(i)<(int)(n);++(i))
#define each(itr,c) for(__typeof(c.begin()) itr=c.begin(); itr!=c.end(); ++itr)
#define all(x) (x).begin(),(x).end()
#define pb push_back
#define fi first
#define se second

struct area
{
    ll x1,y1,z1,x2,y2,z2;

    bool inside(ll a1, ll a2, ll b1, ll b2, ll c1, ll c2)
    {
        return (x1<=a1 && a2<=x2 && y1<=b1 && b2<=y2 && z1<=c1 && c2<=z2);
    }
};

int main()
{
    int n,m;
    scanf(" %d %d", &n, &m);

    vector<area> f(n);
    vector<ll> X,Y,Z;
    rep(i,n)
    {
        scanf(" %lld %lld %lld %lld %lld %lld", &f[i].x1, &f[i].y1, &f[i].z1, &f[i].x2, &f[i].y2, &f[i].z2);

        X.pb(f[i].x1);
        X.pb(f[i].x2);
        Y.pb(f[i].y1);
        Y.pb(f[i].y2);
        Z.pb(f[i].z1);
        Z.pb(f[i].z2);
    }

    sort(all(X));
    X.erase(unique(all(X)),X.end());
    sort(all(Y));
    Y.erase(unique(all(Y)),Y.end());
    sort(all(Z));
    Z.erase(unique(all(Z)),Z.end());

    ll ans=0;

    rep(i,X.size()-1)rep(j,Y.size()-1)rep(k,Z.size()-1)
    {
        int ct=0;
        rep(l,n) if(f[l].inside(X[i],X[i+1],Y[j],Y[j+1],Z[k],Z[k+1])) ++ct;
        if(ct>=m) ans+=(X[i+1]-X[i])*(Y[j+1]-Y[j])*(Z[k+1]-Z[k]);
    }

    printf("%lld\n", ans);
    return 0;
}