旅行者

旅行者

时间限制: 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;
}