旅行者
时间限制: 2 Sec 内存限制: 512 MB
题目描述
小Y来到了一个新的城市旅行。她发现了这个城市的布局是网格状的,也就是有n条从东到西的道路和m条从南到北的道路,这些道路两两相交形成n×m个路口(i,j)(1≤i≤n,1≤j≤m)。
她发现不同的道路路况不同,所以通过不同的路口需要不同的时间。通过调查发现,从路口(i,j)到路口(i,j+1)需要时间ri,j,从路口(i,j)到路口(i+1,j)需要时间ci,j。注意这里的道路是双向的,也就是从路口(i,j+1)到路口(i,j)需要时间同样是ri,j。
小Y有q个询问,她想知道从路口(x1,y1)到路口(x2,y2)最少需要花多少时间。
输入
输入文件tourist.in第一行包含2个正整数n,m,表示城市的大小。
接下来n行,每行包含m-1个整数,第i行第j个正整数表示从一个路口到另一个路口的时间ri,j。
接下来n-1行,每行包含m个整数,第i行第j个正整数表示从一个路口到另一个路口的时间ci,j。
接下来一行,包含1个正整数q,表示小Y的询问个数。
接下来q行,每行包含4个正整数x1,y1,x2,y2,表示两个路口的位置。
输出
输出文件为tourist.out。
输出共q行,每行包含一个整数表示从一个路口到另一个路口最少需要花的时间。
样例输入
样例输入
2 2
2
3
6 4
2
1 1 2 2
1 2 2 1
样例输出
6
7
数据范围
n*m<=20000,q<=100000。
来源
zjoi2016
题解
考虑进行离线分治。
每次操作将当前矩形中的大边长度减半,枚举中轴线上的每个点跑最短路,并更新答案。
对于两点跨过中轴线的询问,最短路必然经过中轴线。
对于两点不在同一块的询问,最短路可能经过中轴线,也可能不经过。
故可以递归分治。
这题卡spfa,还卡常。。
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<algorithm>
#define _(d) while(d((ch=getchar()-48)>=0))
#define inf 1050000000
#define N 20010
#define M 100010
using namespace std;
int n,m,Q,ans[M],w[N],len[N][4],flag[N];
struct info{int x1,y1,x2,y2,pos;}t[M],s1[M],s2[M],s3[M];
struct node{
int x,y,val;
inline bool operator<(const node &p)
const{return val>p.val;}
};
priority_queue<node>q;
inline int next(int x){return ++x>20005?1:x;}
inline int pos(int x,int y){return m*(x-1)+y;}
const int t1[]={0,0,1,-1};
const int t2[]={1,-1,0,0};
inline int get(){
char ch;_(!);int x=ch;
_() x=(x<<3)+(x<<1)+ch;return x;
}//by wbs
inline void solve(int x1,int y1,int x2,int y2,int l,int r)
{
if(l>r)return;
int l1=0,l2=0,l3=0,tot=0,mid,cnt,pos1,pos2;
node st[150];
if((x2-x1)>(y2-y1))
{
mid=x1+x2>>1;
for(int i=y1;i<=y2;i++)st[++tot]=(node){mid,i};
for(int i=l;i<=r;i++)
{
if(t[i].x1<mid&&t[i].x2<mid)s1[++l1]=t[i];
else if(t[i].x1>mid&&t[i].x2>mid)s3[++l3]=t[i];
else s2[++l2]=t[i];
}
}
else
{
mid=y1+y2>>1;
for(int i=x1;i<=x2;i++)st[++tot]=(node){i,mid};
for(int i=l;i<=r;i++)
{
if(t[i].y1<mid&&t[i].y2<mid)s1[++l1]=t[i];
else if(t[i].y1>mid&&t[i].y2>mid)s3[++l3]=t[i];
else s2[++l2]=t[i];
}
}
cnt=l;
for(int j=1;j<=l1;j++)t[cnt++]=s1[j];
pos1=cnt-1;
for(int j=1;j<=l2;j++)t[cnt++]=s2[j];
pos2=cnt;
for(int j=1;j<=l3;j++)t[cnt++]=s3[j];
if((x2-x1)>(y2-y1))
{
solve(x1,y1,mid,y2,l,pos1);
solve(mid,y1,x2,y2,pos2,r);
}
else
{
solve(x1,y1,x2,mid,l,pos1);
solve(x1,mid,x2,y2,pos2,r);
}
for(int j=1;j<=tot;j++)
{
for(int a=x1;a<=x2;a++)
for(int b=y1;b<=y2;b++)w[pos(a,b)]=inf,flag[pos(a,b)]=0;
q.push((node){st[j].x,st[j].y,w[pos(st[j].x,st[j].y)]=0});
flag[pos(st[j].x,st[j].y)]=1;
while(!q.empty())
{
node p=q.top();q.pop();
int x=p.x,y=p.y,val=p.val,now=pos(x,y);
flag[now]=1;if(val!=w[now])continue;
for(int i=0;i<4;i++)
{
int X=x+t1[i],Y=y+t2[i],to=pos(X,Y);
if(X<x1||Y<y1||X>x2||Y>y2)continue;
if(w[to]>w[now]+len[now][i])
{
w[to]=w[now]+len[now][i];
if(!flag[to])q.push((node){X,Y,w[to]});
}
}
}
for(int i=l;i<=r;i++)
{
int id=t[i].pos,X1=t[i].x1,Y1=t[i].y1;
int X2=t[i].x2,Y2=t[i].y2;
ans[id]=min(ans[id],w[pos(X1,Y1)]+w[pos(X2,Y2)]);
}
}
}
int main()
{
int x,x1,y1,x2,y2;
n=get();m=get();
for(int i=1;i<=n;i++)
for(int j=1;j<m;j++)
{
x=get();
len[pos(i,j)][0]=x;
len[pos(i,j+1)][1]=x;
}
for(int i=1;i<n;i++)
for(int j=1;j<=m;j++)
{
x=get();
len[pos(i,j)][2]=x;
len[pos(i+1,j)][3]=x;
}
Q=get();
for(int i=1;i<=Q;i++)
{
x1=get();y1=get();x2=get();y2=get();
t[i]=(info){x1,y1,x2,y2,i};
}
memset(ans,63,sizeof(ans));
solve(1,1,n,m,1,Q);
for(int i=1;i<=Q;i++)printf("%d\n",ans[i]);
return 0;
}