Rectangles POJ 3695

来源：互联网发布：掌趣上游网络牟正文编辑：程序博客网时间：2024/05/29 04:10

这道题n比较小，可以用容斥原理来做

好像直接扫描线也可以

把所有查询的状态用二进制存起来，然后用dfs进行容斥，对于每个状态，找到其是哪个查询容斥项，把面积加上去

#include<stdio.h>#include<string.h>#include<algorithm>using namespace std;const int inf=1000001;struct Rectangle{    int x1,x2,y1,y2;};Rectangle rc[25];int n,m;int ask[100005];int area[1<<20];void dfs(int x1,int y1,int x2,int y2,int deep,int flag,int sta){//才发现原来矩形求交这么简单    if(x1>=x2 || y1>=y2) return;//没有公共面积    if(deep==n){        if(sta) {            for(int i=0;i<m;i++){//找到所有包含该状态的查询                if ((ask[i]|sta)==ask[i]){                    area[ask[i]]+=flag*(x2-x1)*(y2-y1);                }            }        }        return;    }    dfs(x1,y1,x2,y2,deep+1,flag,sta);    dfs(max(x1,rc[deep].x1),max(y1,rc[deep].y1),min(x2,rc[deep].x2),min(y2,rc[deep].y2),deep+1,-flag,sta|(1<<deep));    //n个矩形的公共面积只要找最上的下界，最右的左界，最下的上界，最左的右界}int main(){    #ifndef ONLINE_JUDGE    freopen("in.txt","r",stdin);    #endif // ONLINE_JUDGE    int cas=0;    while(scanf("%d%d",&n,&m)!=EOF,n+m){        memset(area,0,sizeof(area));        memset(ask,0,sizeof(ask));        for(int i=0;i<n;i++) scanf("%d%d%d%d",&rc[i].x1,&rc[i].y1,&rc[i].x2,&rc[i].y2);        for(int i=0;i<m;i++){            int t,x;            scanf("%d",&t);            for(int j=0;j<t;j++) {                scanf("%d",&x);                ask[i]|=1<<(x-1);            }        }        dfs(0,0,inf,inf,0,-1,0);        printf("Case %d:\n", ++cas);        for(int i = 0; i < m; i++)            printf("Query %d: %d\n", i+1, area[ask[i]]);        printf("\n");    }    return 0;}

0 0