当前位置: 代码迷 >> 综合 >> [U]3.1.4 Shaping Regions 递归,计算几何
  详细解决方案

[U]3.1.4 Shaping Regions 递归,计算几何

热度:82   发布时间:2024-01-20 20:43:28.0

一道很有意思的题目。只是今天的课程比较紧张,中午敲了下题,下午上完课后接着敲了下题。改了些bug,终于还是过了。鼓励一下~

解题思路:

题目大意是矩形覆盖问题,求最终见到的各种颜色矩形的面积。由于矩形比较多1000个而且坐标范围用点阵表示会爆空间。题中给的Hints也是含混不清。在队友提点下,终于到了解决的办法。用递归,因为题中有很明显的子问题性质,当前色块的矩形被上面的矩形覆盖之后可以分割成小矩形,再通过小矩形被更上面的矩形覆盖后进而递推。

由于没有想出用好的方法来表示矩形的碰撞情况,若用枚举的话,考虑情况太多而且代码冗余量太大。

矩形覆盖有3种情况:

1.完全不相交:

2.部分相交:

1.>1个角覆盖;

2.>2个角覆盖;

3.完全被包含(4个角覆盖):

显然枚举太麻烦!!!

考虑一般情况:


这样蓝色矩形就在黄色矩形的切割下分为4个小矩形。

把4个小矩形的坐标定义好,再进行递归。

那么其他情况怎么办呢??


同样的在这种情况下,也分为4个小矩形,明显的S1和S4是不符合条件的。

通过我们的坐标定义可以看出S1.ury<S1.lly 这样明显的翻折过来了,这样的矩形不能用。

而S4.urx<S4.llx,同样这样的矩形也不能用。

这种表达方式把所有的情况都包含在内了?还缺了一种:完全相交的情况!

通过坐标定义可以切出符合题意的矩形。将这种情况滤除掉就好了。

另外细节方面要注意的就是细心,还有就是坐标的问题。urx和ury那个点时没有实心点的。这点要好好注意。

总算是切掉了!

Code:

/*
ID:13ysen
LANG:C++
PROG:rect1
*/
#include<stdio.h>
#include<algorithm>
#define MAXR 1001
using namespace std;struct node
{int llx,lly,urx,ury;
}rect[MAXR];struct Color
{int color;int cnt;
}rec[MAXR];bool cmp( Color a,Color b ){ return a.color<b.color; }int N;int min( int a,int b ){ return a<b?a:b; }
int max( int a,int b ){ return a>b?a:b; }bool judge( node a )
{if( a.llx>a.urx || a.lly>a.ury )return false;return true;
}
void recursion( node a,int layer,int index )
{if( judge(a)==false )return ;bool recur=false;for( int i=layer+1;i<=N;i++ ){if( a.llx>rect[i].urx || a.lly>rect[i].ury || a.urx<rect[i].llx || a.ury<rect[i].lly )continue;recur=true;node s1,s2,s3,s4;s1=s2=s3=s4=a;s1.lly=rect[i].ury+1;s3.ury=rect[i].lly-1;s2.urx=rect[i].llx-1;s4.llx=rect[i].urx+1;s4.lly=s2.lly=max( rect[i].lly,a.lly );s4.ury=s2.ury=min( rect[i].ury,a.ury );recursion( s1,i,index );recursion( s2,i,index );recursion( s3,i,index );recursion( s4,i,index );break;}if( recur==false )rec[index].cnt+=(a.urx-a.llx+1)*(a.ury-a.lly+1);return ;
}
int main()
{freopen( "rect1.in","r",stdin );freopen( "rect1.out","w",stdout );int A,B;scanf( "%d %d %d",&A,&B,&N );for( int i=0;i<=N;i++ )rec[i].cnt=rec[i].color=0;rect[0].llx=rect[0].lly=0;rect[0].urx=A-1;rect[0].ury=B-1;rec[0].color=1;for( int i=1;i<=N;i++ ){scanf( "%d %d %d %d %d",&rect[i].llx,&rect[i].lly,&rect[i].urx,&rect[i].ury,&rec[i].color );rect[i].urx-=1;rect[i].ury-=1;}for( int i=0;i<=N;i++ )recursion( rect[i],i,i );sort( rec,rec+N+1,cmp );int ans=0;for( int i=0;i<=N;i++ ){if( rec[i].color==rec[i+1].color )ans+=rec[i].cnt;else{if( rec[i].cnt+ans>0 )printf( "%d %d\n",rec[i].color,rec[i].cnt+ans );ans=0;}}return 0;
}