Problem:
平面に長方形のシートが何枚か置かれている. これらのシートで覆われている部分の面積や周長を求めるプログラムを作成せよ.

ただし, 平面を座標平面とみなしたとき, シートの配置は次の条件 (1), (2) を満たすものとする.

(1) 各シートの長方形の 4 頂点の x, y 座標は全て 0 以上 10000 以下の整数であり, 長方形の各辺は x 軸, または y 軸に平行である.

(2) シートの枚数は高々 10000 枚以下である.

入力データ の 1 行目に長方形の個数 n と問題の種類を表す整数 r が空白で区切って書かれている. 2 行目以降の各行には, 各シートの左下頂点座標 (x1, y1) と右上頂点座標 (x2, y2) の座標値が x1, y1, x2, y2 の順で空白で区切って書いてある.

出力には r = 1 のとき 1 行目に面積を, r = 2 のとき 1 行目に面積, 2 行目に周囲の長さを出力する. いずれの場合も最後に改行を入れること.

また, テストデータの 40% では長方形の頂点の座標が 0 以上 100 以下で, そのうちの 1/2 では面積だけを求める問題である. さらに, 全体の 1/2 は面積だけを求める問題である.

入力例１

入力例２

入力例３

入力例４

5 1

5 2

2 2

3 2

0 0 3 2

0 0 3 2

0 0 8 9

2 2 8 8

1 1 2 5

1 1 2 5

0 0 9 8

3 0 4 9

0 4 6 5

0 4 6 5

5 0 7 9

3 3 5 6

3 3 5 6

5 0 7 6

5 0 7 6

　

出力例１

出力例２

出力例３

出力例４

29

29

80

45

38

36

36

入力

入力は複数のデータセットからなる．n, r がともに 0 のとき入力が終了する．データセットの数は 10 を超えない．

出力

データセットごとに, r = 1 のとき 1 行目に面積を, r = 2 のとき 1 行目に面積, 2 行目に周囲の長さを出力する.

入力例

5 1
0 0 3 2
1 1 2 5
0 4 6 5
3 3 5 6
5 0 7 6
5 2
0 0 3 2
1 1 2 5
0 4 6 5
3 3 5 6
5 0 7 6
2 2
0 0 8 9
0 0 9 8
3 2
2 2 8 8
3 0 4 9
5 0 7 9
0 0

出力例

29
29
38
80
36
45
36

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

Solution in C++11 (judge verdict: Memory Limit Exceeded)

#include <iostream>
#include <algorithm>
#include <utility>
#include <vector>
#include <list>
#include <map>
#include <set>
#include <queue>
#include <stack>
#include <string>
#include <sstream>
#include <cstring>
#include <cstdlib>

using namespace std;

#define ALL(c) c.begin(),c.end()
#define RALL(c) c.rbegin(),c.rend()
#define SORT(x) sort((x).begin(),(x).end())
#define REP(i,x,y) for(int i=(x);i<(y);++i)
#define MP(a,b) make_pair((a),(b))
#define F_ first
#define S_ second
#define debug(x) #x << "=" << (x)
 
#ifdef DEBUG
#define dump(x) std::cerr << debug(x) << " (L:" << __LINE__ << ")" << std::endl
#else
#define dump(x)
#endif
 
template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec){
	os << "[";
	for(const auto &v : vec){
		os << v << ",";
	}
	os << "]";
	return os;
}

typedef long long int lli;
typedef pair<short,short> P;

const int INF=1<<30;
const long long int INF_=1LL<<58;
const double EPS=1e-9;
const int dx[]={1,0,-1,0},dy[]={0,1,0,-1};

struct Rect{
	int x1,x2,y1,y2;
};

class Compare{
	public:
	bool operator()(const Rect &rect1,const Rect &rect2){
		return rect1.x1<rect2.x1;
	}
};

inline short max(short a,short b){
	if(a>b) return a;
	else return b;
}
	
int CalcSize(vector<vector<P> > &range,vector<Rect> &rects){
	int result=0;
	int index=0;
	REP(x,0,10010){
		while(index<rects.size()&&rects[index].x1<=x){
			REP(i,rects[index].x1,rects[index].x2) range[i].push_back(MP(rects[index].y1,rects[index].y2));
			++index;
		}
		sort(range[x].begin(),range[x].end());
		int y=0;
		REP(i,0,range[x].size()){
			if(range[x][i].S_>y) result+=range[x][i].S_-max(y,range[x][i].F_);
			y=max(y,range[x][i].S_);
		}
		//range[x].clear();
		vector<P>().swap(range[x]);
	}
	return result;
}

int Count(vector<vector<P> > &range,vector<vector<bool> > &flag,int x){
	int index=0,y=0,result=0;
	while(index<range[x].size()){
		for(y=max(y,range[x][index].F_); y<range[x][index].S_; ++y){
			if(!flag[(x-1)%3][y]) ++result;
			if(!flag[(x+1)%3][y]) ++result;
			if(!flag[x%3][y-1]) ++result;
			if(!flag[x%3][y+1]) ++result;
		}
		++index;
	}
	return result;
}

int CalcLength(vector<vector<P> > &range,vector<Rect> &rects){
	int result=0,index=0;
	vector<vector<bool> > flag(3,vector<bool>(10010));
	REP(x,1,10009){
		while(index<rects.size()&&rects[index].x1<=x+1){
			REP(i,rects[index].x1,rects[index].x2) range[i].push_back(MP(rects[index].y1,rects[index].y2));
			++index;
		}
		sort(range[x+1].begin(),range[x+1].end());
		if(x-2>=0) fill(flag[(x-2)%3].begin(),flag[(x-2)%3].end(),false);
		int y=0;
		REP(i,0,range[x+1].size()){
			if(range[x+1][i].S_>y) fill_n(flag[(x+1)%3].begin()+max(y,range[x+1][i].F_),range[x+1][i].S_-max(y,range[x+1][i].F_),true);
			y=max(range[x+1][i].S_,y);
		}
		result+=Count(range,flag,x);
		if(x-2>=0){
			//range[x-2].clear();
			vector<P>().swap(range[x-2]);
		}
	}
	return result;
}

int main(){
	std::ios::sync_with_stdio(false);
	std::cin.tie(0);
	int N,R;
	while(cin >> N >> R&&!(N==0&&R==0)){
		vector<vector<P> > range(10010);
		vector<Rect> rects(N);
		REP(i,0,N){
			int x1,y1,x2,y2;
			cin >> x1 >> y1 >> x2 >> y2;
			rects[i].x1=x1+4;
			rects[i].y1=y1+4;
			rects[i].x2=x2+4;
			rects[i].y2=y2+4;
		}
		sort(rects.begin(),rects.end(),Compare());
		cout << CalcSize(range,rects) << endl;
		if(R==2) cout << CalcLength(range,rects) << endl;
	}
	return 0;
}