2019牛客暑期多校训练营（第八场）D-Distance 定期重构

题目传送门

题意：

　　在一个三维空间中，给出q次操作，每次操作可以在空间中加上一个固定点，或者询问一个点，对于一个询问操作，输出距离这个点最近的固定点的曼哈顿距离。

思路：

　　官方题解：先假设所有询问都在加标记之后，那么我们可以同过一次bfs求出网格中每个点离最近标记点的距离，询问就可以O（1）回答。

　　现在考虑定期重构处理增量标记，记一个新增标记队列，每次拿出bfs预处理之后的结果，在暴力枚举队列中每一个新标记，这些结果取个min即可。

　　当队列元素超过一个阈值E时，我们把队列中的标记也进行bfs，更新每个位置的答案，清空新增标记队列。

　　克制复杂度是O（qnmh/E+qE），当E=根号nmh时取最小。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define dep(i,b,a) for(int i=b;i>=a;i--)
#define clr(a,b) memset(a,b,sizeof(a))
#define pb push_back
#define pii pair<int,int >
using namespace std;
const int maxn=;
const int inf=0x3f3f3f3f;
int dis[maxn];
int n,m,h;
int getid(int x,int y,int z){
    return (x-)*m*h+(y-)*h+z;
}

int dir[][]={{,,},{-,,},{,,},{,-,},{,,},{,,-}};
struct node{
    int x,y,z;
};
int getdis(node a,node b){
    return abs(a.x-b.x)+abs(a.y-b.y)+abs(a.z-b.z);
}
vector<node >ve;
queue<node >q;
int main(){
    int qi;
    cin>>n>>m>>h>>qi;
    int sq=sqrt(1ll*n*m*h);
    int si=;
    clr(dis,inf);
    rep(i,,qi){
        int op,x,y,z;
        scanf("%d%d%d%d",&op,&x,&y,&z);
        if(op==){
            ve.pb({x,y,z});
            si++;
            if(si==sq){
                rep(i,,si-){
                    q.push(ve[i]);
                    dis[getid(ve[i].x,ve[i].y,ve[i].z)]=;
                }
                while(!q.empty()){
                    node st=q.front();
                    q.pop();
                    rep(i,,){
                        int xx=st.x+dir[i][];
                        int yy=st.y+dir[i][];
                        int zz=st.z+dir[i][];
                        if(xx<||xx>n||yy<||yy>m||zz<||zz>h)continue;
                        if(dis[getid(xx,yy,zz)]>dis[getid(st.x,st.y,st.z)]+){
                            dis[getid(xx,yy,zz)]=dis[getid(st.x,st.y,st.z)]+;
                            q.push({xx,yy,zz});
                        }
                    }
                }
                ve.clear();
                si=;
            }
        }else{
            int ans=dis[getid(x,y,z)];
            rep(i,,si-){
                ans=min(ans,getdis({x,y,z},ve[i]));
            }
            printf("%d\n",ans);
        }
    }
}