离线分治,最短路
这一篇说得很清楚了,时间复杂度也有分析:
http://blog.csdn.net/neither_nor/article/details/51733997
UOJ上的点实在是恶心。普通最短路不能过,发现每一次分治下去要连续跑最短路,前一次最短路的结果可以帮助下一次。然后套上堆优化dijkstra可卡过。听说用SPFA+SLF,LLL也可过?
BZOJ网站昨天就挂了,有空再上去交一下4456
#include<cstdio>
#include<cstring>
#include<algorithm>
#define Q 100005
#define MAX 20005
#define reg register
#define cmin(_i,_j) _i>_j?_i=_j:0
#define pack(_i,_j) ((_i-1)*m+_j)
#define add2(_i,_j,_k) {add(_i,_j,_k);add(_j,_i,_k);}
#define getc() (S == T && (T = (S = B) + fread(B, 1, 1 << 15, stdin), S == T) ? EOF : *S++)
char B[1 << 15], *S = B, *T = B;
using namespace std;
int in()
{
reg int r=0;
reg char c=getc();
while(c<'0'||c>'9')c=getc();
while(c>='0'&&c<='9')r=r*10+c-'0',c=getc();
return r;
}
struct edge{int next,to,val;}e[MAX*5];
struct node
{
int x, dis;
node(){}
node(int a, int b):x(a),dis(b){}
friend bool operator < (node a, node b)
{
return a.dis < b.dis;
}
}heap[MAX*5];
const int INF = 1<<29;
int ecnt, n, m, q, rx[MAX], ry[MAX], last[MAX], dis[MAX], tot, ans[Q], vis[MAX], timer, qq[MAX];
bool inq[MAX];
struct query
{
int x1, y1, x2, y2, ans, id;
}que[Q], temp[Q];
void heap_add(node x)
{
heap[++tot]=x;
for(int i = tot, j = i>>1; j; i=j, j>>=1)
{
if(heap[i]<heap[j])
swap(heap[i],heap[j]);
else break;
}
}
void heap_pop()
{
heap[1]=heap[tot--];
for(int i = 1, j = i<<1; j<=tot; i=j, j<<=1)
{
if(j<tot && heap[j|1]<heap[j])j|=1;
if(heap[j]<heap[i])
swap(heap[i],heap[j]);
else break;
}
}
void add(int a, int b, int c)
{
e[++ecnt]=(edge){last[a],b,c};
last[a]=ecnt;
}
void init()
{
for(int i = 1; i <= n; i++)
for(int j = 1; j <= m; j++)
{
int k = pack(i,j);
rx[k]=i;
ry[k]=j;
}
}
void dijk(int s, int x1, int y1, int x2, int y2, bool first)
{
timer++;
qq[0]=s;
heap_add(node(s,dis[s]=0));
memset(dis,63,sizeof(dis));
dis[s]=0;
int cnt=((x2-x1+1)*(y2-y1+1));
while(tot)
{
int x=heap[1].x;
heap_pop();
if(vis[x]==timer)continue;
cnt--;
vis[x]=timer;
for(reg int i = last[x]; i; i=e[i].next)
{
int y=e[i].to;
if(vis[y]==timer||rx[y]<x1||rx[y]>x2||ry[y]<y1||ry[y]>y2)continue;
if(dis[x]+e[i].val<dis[y])
{
dis[y]=dis[x]+e[i].val;
heap_add(node(y,dis[y]));
}
}
if(!cnt)tot=cnt;
}
}
void solve(int x1, int y1, int x2, int y2, int l, int r)
{
if(l>r)return;
if(x2-x1 <= y2-y1)
{
int mid=(y1+y2)>>1;
for(reg int i = x1; i <= x2; i++)
{
dijk(pack(i,mid),x1,y1,x2,y2,i==x1);
for(reg int j = l; j <= r; j++)
cmin(que[j].ans, dis[pack(que[j].x1,que[j].y1)]+dis[pack(que[j].x2,que[j].y2)]);
}
int nl=l-1, nr=r+1;
for(reg int i = l; i <= r; i++)
{
if(que[i].y1<mid && que[i].y2<mid)temp[++nl]=que[i];
else if(que[i].y1>mid && que[i].y2>mid)temp[--nr]=que[i];
cmin(ans[que[i].id],que[i].ans);
}
for(reg int i = l; i <= r; i++)que[i]=temp[i];
solve(x1,y1,x2,mid-1,l,nl);
solve(x1,mid+1,x2,y2,nr,r);
}
else
{
int mid=(x1+x2)>>1;
for(reg int i = y1; i <= y2; i++)
{
dijk(pack(mid,i),x1,y1,x2,y2,i==y1);
for(int j = l; j <= r; j++)
cmin(que[j].ans, dis[pack(que[j].x1,que[j].y1)]+dis[pack(que[j].x2,que[j].y2)]);
}
int nl=l-1, nr=r+1;
for(reg int i = l; i <= r; i++)
{
if(que[i].x1<mid && que[i].x2<mid)temp[++nl]=que[i];
else if(que[i].x1>mid && que[i].x2>mid)temp[--nr]=que[i];
cmin(ans[que[i].id],que[i].ans);
}
for(reg int i = l; i <= r; i++)que[i]=temp[i];
solve(x1,y1,mid-1,y2,l,nl);
solve(mid+1,y1,x2,y2,nr,r);
}
}
int main()
{
n=in(), m=in();
init();
for(reg int i = 1, v; i <= n; i++)
for(reg int j = 1; j < m; j++)
{
v=in();
add2(pack(i,j),pack(i,j+1),v);
}
for(reg int i = 1, v; i < n; i++)
for(reg int j = 1; j <= m; j++)
{
v=in();
add2(pack(i,j),pack(i+1,j),v);
}
q=in();
for(reg int i = 1; i <= q; i++)
{
que[i].x1=in(); que[i].y1=in();
que[i].x2=in(); que[i].y2=in();
que[i].ans=INF;
que[i].id=i;
}
memset(ans,63,sizeof(ans));
solve(1,1,n,m,1,q);
for(reg int i = 1; i <= q; i++)
printf("%d\n",ans[i]);
}