usaco Section 3.1 Shaping Regions -- 矩形切割

      听说这道题用矩形切割做，于是学习了2004薛矛的oi论文。

     

---

 

      我的做法是：弄一个白纸集合，开始里面只有一张白纸(0,0) (A,B)。然后把输入的N个矩形倒序地，每一个都与白纸集合中的所有白纸依次切割。对于输入的矩形rec[i]与每张白纸，用rec[i]去覆盖它。如果有交集，则交集为rec[i]的颜色，把元素个数统计到rec[i]的颜色中，然后把白纸中没被覆盖的地方切成矩形加入到白纸集合，并且把原来被切的白纸从集合中删除。没有交集就什么都不用做。N个矩形都处理过后，白纸集合中矩形就都是颜色为1的矩形，把它们的元素个数统计到颜色1中去即可。

      其中矩形切割的函数是按着上面的伪代码写的。

      /* ID: lingxiu1 LANG: C++ TASK: rect1 */ #include<vector> #include<stdio.h> #include<string.h> #include<iostream> #include<algorithm> using namespace std; struct RECT { int x1, y1, x2, y2, col ; //x1,y1为左下角。x2,y2为右上角。col为颜色 bool cut;//标记是否被切割 } ; RECT rec[1001], temp; //rec存储输入的矩形 vector<RECT> set; //set存储被切出的白色矩形数集合 int cnt[2501] ; //cnt用来记录每个颜色有多少个元素 void count(int x1, int y1, int x2, int y2, int col) { cnt[col] += (y2 - y1) * (x2 - x1) ; //记录颜色为col的元素个数 } void cut(int j, RECT b, int dir) //set[j]为白纸，若b与set[j]有交集则切割 { if(set[j].y2 <= b.y1 || set[j].x1 >= b.x2 || set[j].x2 <= b.x1 || set[j].y1 >= b.y2 ) //没有交集 return ; //直接返回 set[j].cut = 1 ; //标记被切割 int k1, k2 ; if(dir == 1) //x方向（横） { k1 = max(set[j].x1, b.x1) ; k2 = min(set[j].x2, b.x2) ; //k1, k2为x方向交集 if(set[j].x1 < k1) { temp.x1 = set[j].x1 ; temp.y1 = set[j].y1 ; temp.x2 = k1 ; temp.y2 = set[j].y2 ; temp.cut = 0 ; set.push_back(temp) ; //切出白纸temp,加入集合 } if(k2 < set[j].x2) { temp.x1 = k2 ; temp.y1 = set[j].y1 ; temp.x2 = set[j].x2 ; temp.y2 = set[j].y2 ; temp.cut = 0 ; set.push_back(temp) ;//切出白纸temp,加入集合 } RECT c; c.x1 = k1 ; c.y1 = set[j].y1 ; c.x2 = k2 ; c.y2 = set[j].y2 ; set.push_back(c) ; //c为set[j]切掉x方向后的剩余 cut(set.size() - 1, b, dir + 1) ; //继续y方向的切割 } else //y方向 （竖） { k1 = max(set[j].y1, b.y1) ; k2 = min(set[j].y2, b.y2) ; //k1, k2为y方向交集 if(set[j].y1 < k1) { temp.x1 = set[j].x1 ; temp.y1 = set[j].y1 ; temp.x2 = set[j].x2 ; temp.y2 = k1 ; temp.cut = 0 ; set.push_back(temp) ; //切出白纸temp,加入集合 } if(k2 < set[j].y2) { temp.x1 = set[j].x1 ; temp.y1 = k2 ; temp.x2 = set[j].x2 ; temp.y2 = set[j].y2 ; temp.cut = 0 ; set.push_back(temp) ; //切出白纸temp,加入集合 } count(max(b.x1, set[j].x1), k1, min(b.x2, set[j].x2), k2, b.col) ; //计算交集的元素个数，颜色为b的颜色 } } int main(){ freopen("rect1.in","r",stdin); freopen("rect1.out","w",stdout); int a, b, n, i, j, k ; scanf("%d%d%d", &a, &b, &n); for(i = 0; i < n; i ++) { scanf("%d%d%d%d%d", &rec[i].x1, &rec[i].y1, &rec[i].x2, &rec[i].y2, &rec[i].col); } memset(cnt, 0, sizeof(cnt)) ; set.clear() ; temp.x1 = 0 ; temp.y1 = 0 ; temp.x2 = a ; temp.y2 = b ; temp.cut = 0; set.push_back(temp) ; //第一张白纸 (0,0) (a,b) for(i = n-1; i >= 0; i --)//倒序 { int tmp = set.size() ; for(j = 0; j < tmp; j ++) { cut(j, rec[i], 1); //每次覆盖切出的白纸 } for(vector<RECT>::iterator j = set.begin(); j != set.end(); ) if((*j).cut) { //被切过的白纸删除 j = set.erase(j) ; } else j ++ ; } for(i = 0; i < set.size(); i ++) { count(set[i].x1, set[i].y1, set[i].x2, set[i].y2, 1); //最后计算所有白纸的个数 } for(i = 1; i < 2501; i ++) { if(cnt[i]) printf("%d %d/n", i, cnt[i]); } return 0; }

USER: li xiu [lingxiu1]  TASK: rect1  LANG: C++

Compiling...   Compile: OK   Executing...

Test 1: TEST OK [0.000 secs, 3052 KB]

Test 2: TEST OK [0.000 secs, 3052 KB]

Test 3: TEST OK [0.000 secs, 3052 KB]

Test 4: TEST OK [0.000 secs, 3052 KB]

Test 5: TEST OK [0.000 secs, 3052 KB]

Test 6: TEST OK [0.000 secs, 3052 KB]

Test 7: TEST OK [0.000 secs, 3052 KB]

Test 8: TEST OK [0.000 secs, 3052 KB]

Test 9: TEST OK [0.000 secs, 3052 KB]

Test 10: TEST OK [0.000 secs, 3052 KB]

Test 11: TEST OK [0.011 secs, 3052 KB]

All tests OK.

Your program ('rect1') produced all correct answers! This is yoursubmission #12 for this problem. Congratulations!

速度还不错啊~