面积交(贪心)
for (int i = 1; i < (1 << n); ++i) {
int xl = 0, yl = 0, xr = W, yr = H;
for (int j = 0; j < n; ++j) {
if (i >> j & 1) {
xl = max(xl, x[j]);
xr = min(xr, xx[j]);
yl = max(yl, y[j]);
yr = min(yr, yy[j]);
}
}
area[i] = 1ll * max(0, (xr - xl)) * max(0, (yr - yl));
}
ps:area[i]代表状态为i的几个矩形的面积交,这里是状态压缩表示,若i的2进制表示中某位是1,则表示含有第i个矩形,具体求法画个图就很容易贪心求出,求相交矩形的四个边界即可(这里给出的是矩形的左下角的点和右上角的点)
面积并(容斥)
for (int i = 1; i < (1 << n); ++i) {
dp[i] = 0;
for (int s = i; s; s = (s - 1) & i) {
dp[i] += area[s] * (__builtin_parity(s) ? 1 : -1);
}
}
要求面积的并首先要求出任意矩形的面积交然后用容斥原理即可,dp[i]代表即使状态为i的多个矩形面积并,这里同样用的是状态压缩,__builtin_parity(s)函数的作用是返回 s这个状态中2进制表示1的个数模2的值,也就是表示状态中1的个数的奇偶,这里就可用来判断容斥当中是加还是减,第二重for循环目的是枚举求出i这个状态的所有子集,以便用于容斥求并,这里的做法很明显优于从0开始一直枚举到当前状态容斥,对于每一个状态的求并都要容斥一遍