Rectangles POJ 3695

这道题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;
}