[HDU2461] Rectangles [2008 Asia Hefei Regional Contest Online F]

题意

给定不多于20个矩形，每次询问这些矩形的子集的矩形面积并。

题解

容斥原理，直接2^20枚举显然会有很多无用的状态，那么用dfs来枚举状态，当并的面积为0时，之后的状态就没有用了，这样就可以提高效率。

代码

/****************************************\
* Author : ztx
* Title  : F-Rectangles
* ALG    : 容斥原理
* CMT    : g++  I64d
* Time   :
\****************************************/

#include <cstdio>
#define Rep(i,l,r) for(i=(l);i<=(r);i++)
#define rep(i,l,r) for(i=(l);i< (r);i++)
#define Rev(i,r,l) for(i=(r);i>=(l);i--)
#define rev(i,r,l) for(i=(r);i> (l);i--)
typedef long long ll ;
typedef double lf ;
int CH , NEG ;
template <typename TP>inline void read(TP& ret) {
    ret = NEG = 0 ; while (CH=getchar() , CH<'!') ;
    if (CH == '-') NEG = true , CH = getchar() ;
    while (ret = ret*10+CH-'0' , CH=getchar() , CH>'!') ;
    if (NEG) ret = -ret ;
}
template <typename TP>inline void readc(TP& ret) {
    while (ret=getchar() , ret<'!') ;
    while (CH=getchar() , CH>'!') ;
}
template <typename TP>inline void reads(TP *ret) {
    ret[0]=0;while (CH=getchar() , CH<'!') ;
    while (ret[++ret[0]]=CH,CH=getchar(),CH>'!') ;
    ret[ret[0]+1]=0;
}

#define  max(a,b) (a)>(b)?(a):(b)
#define  min(a,b) (a)<(b)?(a):(b)

struct Rec {
    int x1 , y1 , x2 , y2 ;
    inline Rec intersect(Rec &B) {
        Rec ret ;
        ret.x1 = max(x1,B.x1) ;
        ret.y1 = max(y1,B.y1) ;
        ret.x2 = min(x2,B.x2) ;
        ret.y2 = min(y2,B.y2) ;
        return ret ;
    }
    inline int area() {
        if (x1 >= x2 || y1 >= y2) return 0 ;
        return (x2-x1)*(y2-y1) ;
    }
    inline void in() {
        read(x1) , read(y1) , read(x2) , read(y2) ;
    }
}rec[21] ;

int n , m , kiss = 0 ;
int tot , id[21] , ans ;

inline void dfs(int dep , Rec now ,int idx) {
    // 直接用dfs找到所有的状态集合，可以同时剪去不必要的枚举
    Rec nxt ;
    int i ;
    Rep (i,idx,tot)
        if (nxt=now.intersect(rec[id[i]]),nxt.area()) {
            if (dep&1) ans -= nxt.area() ;
            else ans += nxt.area() ;
            dfs(dep+1,nxt,i+1) ;
        }
}

int main() {
int i , j ;
//    #define READ
    #ifdef  READ
        freopen(".in" ,"r",stdin ) ;
        freopen(".out","w",stdout) ;
    #endif
    while (++kiss) {
        read(n) , read(m) ;
        if (n+m == 0) break ;
        Rep (i,1,n) rec[i].in() ;
        printf("Case %d:\n",kiss) ;
        Rep (i,1,m) {
            ans = 0 ;
            read(tot) ;
            Rep (j,1,tot) read(id[j]) ;
            Rep (j,1,tot) {
                ans += rec[id[j]].area() ;
                dfs(1,rec[id[j]],j+1);
            }
            printf("Query %d: %d\n",i,ans) ;
        }
        puts("") ;
    }
    #ifdef  READ
        fclose(stdin) ; fclose(stdout) ;
    #else
        getchar() ; getchar() ;
    #endif
    return 0 ;
}