BZOJ4456 [Zjoi2016]旅行者

好神的分治做法

我们把矩形顺着比较短的一边从中间切成两半，从这条中线上每个点做一次到这个矩形内每个点的dijkstra，把两个端点都在这个矩形里的询问用中线上每个点到两个端点的最短路更新，然后分治两半边

复杂度分析：瞎YY一下我们知道这个算法每次更新答案要做比较小的一边的边长次dijkstra，所以当是一个正方形的时候分治的复杂度最高，而询问一定是传到底的时候复杂度最高

因为每次切一半两边都是一样的，所以每个询问只要传到底，不管往哪边传都可以，为了方便分析不妨设每个询问往两边各传一半

T(n,m,Q)表示短边长为n，长边长为m，矩形里有Q个询问的复杂度，T(n,Q)表示边长为n，里边有Q个询问的正方形的复杂度

把最后的式子化简一下，面积为S的话，复杂度是O(S*sqrt(S)*log S)的

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 40010
#define MAXM 100010
#define ll long long
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
char xB[1<<15],*xS=xB,*xT=xB;
#define getc() (xS==xT&&(xT=(xS=xB)+fread(xB,1,1<<15,stdin),xS==xT)?0:*xS++)
inline int read(){
    char ch=getc();
    int f=1,x=0;
    while(!(ch>='0'&&ch<='9')){if(ch=='-')f=-1;ch=getc();}
    while(ch>='0'&&ch<='9'){x=x*10+(ch-'0');ch=getc();}
    return x*f;
}
struct vec{
	int to;
	int fro;
	int v;
};
struct pt{
	int x;
	int y;
	int num;
	pt(){
		
	}
	pt(int _x,int _y);
};
struct que{
	pt x;
	pt y;
	int num;
};
struct sp{
	int to;
	int v;
	sp(){
		
	}
	sp(int _to,int _v){
		to=_to;
		v=_v;
	}
	friend bool operator <(sp x,sp y){
		return x.v>y.v;
	}
};
vec mp[MAXN*2];
int tai[MAXN],cnt;
int ans[MAXM];
que a[MAXM],ta[MAXM];
int dis[MAXN],vis[MAXN];
priority_queue<sp>q;
int T;
int srd[MAXN];
int n,m,Q;
int hx[MAXN],hy[MAXN];
pt::pt(int _x,int _y){
	x=_x;
	y=_y;
	num=(x-1)*m+y;
}
inline void be(int x,int y,int z){
	mp[++cnt].to=y;
	mp[cnt].fro=tai[x];
	tai[x]=cnt;
	mp[cnt].v=z;
}
inline void bde(int x,int y,int z){
	be(x,y,z);
	be(y,x,z);
}
void dijkstra(int X,int Y,int XX,int YY,int x){
	int i,y;
	T++;
	vis[x]=T;
	dis[x]=0;
	q.push(sp(x,0));
	while(!q.empty()){
		x=q.top().to;
		if(srd[x]==T){
			q.pop();
			continue ;
		}
		srd[x]=T;
		q.pop();
		for(i=tai[x];i;i=mp[i].fro){
			y=mp[i].to;
			if(hx[y]<X||hx[y]>XX||hy[y]<Y||hy[y]>YY){
				continue ;
			}
			if(vis[y]!=T){
				vis[y]=T;
				dis[y]=INF;
			}
			if(dis[y]>dis[x]+mp[i].v){
				dis[y]=dis[x]+mp[i].v;
				q.push(sp(y,dis[x]+mp[i].v));
			}
		}
	}
}
void cal(int x,int y,int xx,int yy,int l,int r){
	if(x>xx||y>yy||l>r){
		return ;
	}
	int i,j;
	if(xx-x<=yy-y){
		int mid=y+yy>>1;
		for(i=x;i<=xx;i++){
			dijkstra(x,y,xx,yy,pt(i,mid).num);
			for(j=l;j<=r;j++){
				ans[a[j].num]=min(ans[a[j].num],dis[a[j].x.num]+dis[a[j].y.num]);
			}
		}
		int tl=l-1,tr=r+1;
		for(i=l;i<=r;i++){
			if(a[i].x.y<mid&&a[i].y.y<mid){
				ta[++tl]=a[i];
			}
			if(a[i].x.y>mid&&a[i].y.y>mid){
				ta[--tr]=a[i];
			}
		}
		memcpy(a+l,ta+l,sizeof(que)*(r-l+1));
		cal(x,y,xx,mid-1,l,tl);
		cal(x,mid+1,xx,yy,tr,r);
	}else{
		int mid=x+xx>>1;
		for(i=y;i<=yy;i++){
			dijkstra(x,y,xx,yy,pt(mid,i).num);
			for(j=l;j<=r;j++){
				ans[a[j].num]=min(ans[a[j].num],dis[a[j].x.num]+dis[a[j].y.num]);
			}
		}
		int tl=l-1,tr=r+1;
		for(i=l;i<=r;i++){
			if(a[i].x.x<mid&&a[i].y.x<mid){
				ta[++tl]=a[i];
			}
			if(a[i].x.x>mid&&a[i].y.x>mid){
				ta[--tr]=a[i];
			}
		}
		memcpy(a+l,ta+l,sizeof(que)*(r-l+1));
		cal(x,y,mid-1,yy,l,tl);
		cal(mid+1,y,xx,yy,tr,r);
	}
}
int main(){
	int i,j,x,y,xx,yy;
	n=read();
	m=read();
	for(i=1;i<=n;i++){
		for(j=1;j<=m;j++){
			hx[pt(i,j).num]=i;
			hy[pt(i,j).num]=j;
		}
	}
	for(i=1;i<=n;i++){
		for(j=1;j<m;j++){
			x=read();
			bde(pt(i,j).num,pt(i,j+1).num,x);
		}
	}
	for(i=1;i<n;i++){
		for(j=1;j<=m;j++){
			x=read();
			bde(pt(i,j).num,pt(i+1,j).num,x);
		}
	}
	Q=read();
	memset(ans,0x3f,sizeof(ans));
	for(i=1;i<=Q;i++){
		x=read();
		y=read();
		xx=read();
		yy=read();
		a[i].x=pt(x,y);
		a[i].y=pt(xx,yy);
		if(x==xx&&y==yy){
			ans[i]=0;
		}
		a[i].num=i;
	}
	cal(1,1,n,m,1,Q);
	for(i=1;i<=Q;i++){
		printf("%d\n",ans[i]);
	}
	return 0;
}

/*

*/