poj.org/problem?id=3248
这是一道比较复杂的题目,思路也不是很容易想出来,看了讨论,想了好几天才捋清思路,写了半天代码,调试了半天,终于给AC了。
思路:
1. 以某方格x,y为中心,Manhattan距离为d的所有的所有方格,整体旋转45度后呈现类似正方形的形状, 其中x+y为偶数的方格组成一个正方形, x+y为奇数的方格组成另一个正方形;如果通过旋转构造,使得方格个数=正方形的面积,那么问题就可以转化为分别由奇数和偶数方格变幻的得到的几个正方形在平面上所占的面积和,然后使用线段树求解就是一个可行方案
2. 但是,不是所有方格都是有效的,题目给出了矩形限制区域,需要从格子总和中减去矩形区域外部的格子;这时候思路需要恢复到旋转前,超出矩形某个边界的格子具有类似等腰三角形的形态,那么这个问题就归结为底座在同一水平线上的等腰三角形的“面积”和,这只需要排序,然后遍历一遍即可,当然需要在四个边界上都处理一次
菜鸡的实现
1. 旋转坐标,构造正方形坐标,这里需要分x+y为奇数和偶数两种情况处理
for (mm=0; mm<2; mm++) {
kk=0; for (i=0; i<k; i++) {
x=rs[i].x; y=rs[i].y; d=rs[i].d;
if (x<0||x>=n||y<0||y>=m) assert(0);
if (((x+y)&1)!=mm) {
c = d/2;
xx=x-y; if (xx<0) xx=-((-xx+1)/2); else xx=xx/2;
sx = xx-c;
xx=x+y; if (xx<0) xx=-((-xx+1)/2); else xx=xx/2;
sy = xx-c;
c = c+c+1;
} else {
if (d==0) continue;
c = (d+1)/2;
xx=x-y+1; if (xx<0) xx=-((-xx+1)/2); else xx=xx/2;
sx = xx-c;
xx=x+y-1; if (xx<0) xx=-((-xx+1)/2); else xx=xx/2;
sy = xx-c+1;
c = c+c;
}
qs[kk].