矩形的交集如果非空,那么仍然是一个矩形。考虑求出经过每个格子 (i,j) 的矩形数量 fi,j,f 可以通过二维前缀和在 O(n + m2) 的时间内求出,则 ans =∑C(fi,j,3)。
上述方法会重复统计,这些多出来的部分一定是交集覆盖了 (i,j −1) 或者 (i−1,j)。减去 交集同时覆盖 (i,j −1) 和 (i,j) 的答案,再减去交集同时覆盖 (i−1,j) 和 (i,j) 的答案,那么 交集同时覆盖 (i−1,j−1) 和 (i,j) 的答案会被多减一次,再加回来即可。按照该方法,每个方 案都将在交集的左上角格子处被恰好统计到一次。 对于计算交集同时覆盖 (i,j −1) 和 (i,j) 的答案,只需要把所有矩形的边长减小 1,其它 几项答案的计算同理。 时间复杂度 O(n + m2)。
#include<cstdio>
typedef long long ll;
const int N=100010,M=1010;
int Case,n,m,i,j,e[N][4],s[M][M];ll C3[N],ans;
ll solve(int dx,int dy)
{
for(i=1;i<=m;i++)
for(j=1;j<=m;j++)
s[i][j]=0;
for(i=1;i<=n;i++)
{
int xl=e[i][0]+dx,yl=e[i][1]+dy,xr=e[i][2],yr=e[i][3];
if(xl>xr||yl>yr)
continue;
s[xl][yl]++;
s[xr+1][yl]--;
s[xl][yr+1]--;
s[xr+1][