usaco Shaping Regions
有太多不知道，又常常焦急地想去知道，因为这份心情，所以我愿意将我会的与大家分享
——这也是我创建这个博客的目的！
今天发现自己辛辛苦苦做出一个题目后，越来越懒得写说明了。
突然想起了自己当初创建这个博客时，是抱着份什么心情............
hehe！
这道真的很有意思！做了两天！今天终于做出来了，真的很高兴。。。。
INPUT FORMAT
The order of the input lines dictates the order of laying down the rectangles. The first input line is a rectangle "on the bottom".
| Line 1: | A, B, and N, space separated (1 <= A,B <= 10,000) |
| Lines 2-N+1: | Five integers: llx, lly, urx, ury, color: the lower left coordinates and upper right coordinates of the rectangle whose color is `color' (1 <= color <= 2500) to be placed on the white sheet. The color 1 is the same color of white as the sheet upon which the rectangles are placed. |
其实题目看起来不难，如果数据小点，开个二维数组， 直接模拟也是可以出答案的。
但这个想法是好的，出题者是无情的。
测试数据A、B都达到了10000，开个10000*10000的二维数组。肯定不成！
这也是这道题目的难点。如何实现数据压缩。
题目过程就像在纸上画各色的矩形一样
（矩形的边都平行于纸的轮廓，假设颜色不叠加，并且完全覆盖），
最后纸上图案会变成什么样呢？
有一点可以想像的的出——就是图案可以拆分出若干个矩形，每个矩形内着色完全相同。
好下图所示
我选择存储矩形
struct rectangle {
int x1, y1;
int x2, y2;
}rec[100000];
int x2, y2;
}rec[100000];
rec数组初始化存入（0， 0）， （A， B）；然后就是模拟过程；
我的方法是判断新画到纸上的矩阵有没有切割数组中的矩阵rec[i]，
分为2种情况：
1、当前画的矩阵完全包含rec[i]，表示新画的矩形完全覆盖rec[i],当删除rec[i].
2、当前画的矩阵的4条边分别是否与rec[i]相交，若相交则会将rec[i]分成两个矩形。
然后将当前画的矩阵加入到rec数组中。
细节自己体会了！！！hehe....
代码如下：
/*
ID: guo geer
PROG: rect1
LANG: C++
*/
#include<iostream>
#include<fstream>
using namespace std;
struct rectangle {
int x1, y1;
int x2, y2;
};
int c[100000];
rectangle rec[100000];
int line(int a, int b, int x, int y) // a < b, x < y, [a, b]&[x, y] != NULL
{
if(x > a && x < b)
return 1;
if(y > a && y < b)
return 1;
if(x <= a && y >= b)
return 1;
return 0;
}
int main()
{
ifstream fin("rect1.in");
ofstream fout("rect1.out");
int a, b, n;
int i,j,k;
while(fin>>a>>b>>n)
{
int count = 0;
rec[0].x1 = 0, rec[0].y1 = 0;
rec[0].x2 = a, rec[0].y2 = b;
c[0] = 1;
count ++;
int x1, y1, x2, y2, col;
for(k=0; k<n; k++)
{
fin>>x1>>y1>>x2>>y2>>col;
for(i=0; i<count; i++)
{
if(rec[i].x1 >= x1 && rec[i].x2 <= x2 && rec[i].y1 >= y1 && rec[i].y2 <= y2)
{
c[i] = 0;
}
else if(rec[i].y1 < y1 && rec[i].y2 > y1 && line(rec[i].x1, rec[i].x2, x1, x2))
{
rec[count].x1 = rec[i].x1;
rec[count].y1 = y1;
rec[count].y2 = rec[i].y2;
rec[count].x2 = rec[i].x2;
c[count] = c[i];
count ++;
rec[i].y2 = y1;
}
else if(rec[i].x1 < x2 && rec[i].x2 > x2 && line(rec[i].y1, rec[i].y2, y1, y2))
{
rec[count].x2 = x2;
rec[count].y2 = rec[i].y2;
rec[count].y1 = rec[i].y1;
rec[count].x1 = rec[i].x1;
c[count] = c[i];
count ++;
rec[i].x1 = x2;
}
else if(rec[i].y1 < y2 && rec[i].y2 > y2 && line(rec[i].x1, rec[i].x2, x1, x2))
{
rec[count].x1 = rec[i].x1;
rec[count].y1 = rec[i].y1;
rec[count].y2 = y2;
rec[count].x2 = rec[i].x2;
c[count] = c[i];
count ++;
rec[i].y1 = y2;
}
else if(rec[i].x1 < x1 && rec[i].x2 > x1 && line(rec[i].y1, rec[i].y2, y1, y2))
{
rec[count].x1 = x1;
rec[count].y1 = rec[i].y1;
rec[count].y2 = rec[i].y2;
rec[count].x2 = rec[i].x2;
c[count] = c[i];
count ++;
rec[i].x2 = x1;
}
}
int count_1 = 0;
for(j=0; j<count; j++)
if(c[j] != 0)
{
rec[count_1] = rec[j];
c[count_1] = c[j];
count_1 ++;
}
count = count_1;
rec[count].x1 = x1;
rec[count].y1 = y1;
rec[count].x2 = x2;
rec[count].y2 = y2;
c[count] = col;
count ++;
}
int res[2501];
for(i=0; i<= 2500; i++)
res[i] = 0;
for(i=0; i<count; i++)
res[c[i]] += (rec[i].y2 - rec[i].y1) * (rec[i].x2 - rec[i].x1);
for(i=1; i<= 2500; i++)
if(res[i] != 0)
fout<<i<<' '<<res[i]<<endl;
}
return 0;
}