最强大脑
最强大脑
难度级别：E； 运行时间限制：3000ms； 运行空间限制：256000KB； 代码长度限制：2000000B
试题描述
zhb国是一个传说中的国度，国家的居民叫做最强(chang)大脑(diao)。Zhb国是一个N×M的矩形方阵，每个格子代表一个街区。
然而zhb国是没有交通工具的。居民大脑(diao)们完全靠地面的弹射装置来移动。
每个街区都装有弹射装置。使用弹射装置是需要支付一定费用的。而且每个弹射装置都有自己的弹射能力。
我们设第i行第j列的弹射装置有Aij的费用和Bij的弹射能力。并规定有相邻边的格子间距离是1。那么，任何居民大脑(diao)都只需要在(i,j)支付Aij的费用就可以任意选择弹到距离不超过Bij的位置了。如下图
（从红色街区交费以后可以跳到周围的任意蓝色街区。）
现在的问题很简单。有三个居民，她们分别是zhb的maze，分别叫做X，Y，Z。现在她们决定聚在一起玩找zhb玩(….)，于是想往其中一人的位置集合。但由于zhb太抠门了，不给报销路费，所以告诉你3个maze的坐标，求往哪里集合大家需要花的费用总和最低。
Zhb:如果你完美解决了这个问题，就授予你“最强大脑”称号。
输入
输入的第一行包含两个整数N和M，分别表示行数和列数。
接下来是2个N×M的自然数矩阵，为Aij和Bij
最后一行六个数，分别代表X，Y，Z所在地的行号和列号。
输出
第一行输出一个字符X、Y或者Z。表示最优集合地点。
第二行输出一个整数，表示最小费用。
如果无法集合，只输出一行NO
输入示例
4 4
0 0 0 0
1 2 2 0
0 2 2 1
0 0 0 0
5 5 5 5
5 5 5 5
5 5 5 5
5 5 5 5
2 1 3 4 2 2
输出示例
Z
15
其他说明
20% N, M ≤ 10; Bij ≤ 20
40% N, M ≤ 100; Bij ≤ 20
100% 1 ≤ N, M ≤ 150; 0 ≤ Bij ≤ 109; 0 ≤ Aij ≤ 1000
嗯，直接上代码
//////////////////// //版权所有：张煜暄// // 禁止抄袭 // //如有雷同纯属巧合// // ！！！！！！！ // //////////////////// #include<iostream> #include<cstdio> #include<cmath> #include<algorithm> #include<queue> #include<cstring> #define INF 2147483647 using namespace std; inline int read() { int x=0,t=1,c; while(!isdigit(c=getchar()))if(c==\'-\')t=-1; while(isdigit(c))x=x*10+c-\'0\',c=getchar(); return x*t; } bool is_in(int left,int down,int right,int up,int x,int y,int r) { if(left>right||down<up)return false; int X,Y,distx,disty; if(x<left)X=0; else if(x>=left&&x<=right)X=1; else if(x>right)X=2; if(y<up)Y=0; else if(y>=up&&y<=down)Y=1; else if(y>down)Y=2; if(X==0) { if(Y==0) { distx=left-x; disty=up-y; if(distx+disty<=r)return true; else return false; } if(Y==1) { if(left-x<=r)return true; else return false; } if(Y==2) { distx=left-x; disty=y-down; if(distx+disty<=r)return true; else return false; } } if(X==1) { if(Y==0) { if(up-y<=r)return true; else return false; } if(Y==1) { return true; } if(Y==2) { if(y-down<=r)return true; else return false; } } if(X==2) { if(Y==0) { distx=x-right; disty=up-y; if(distx+disty<=r)return true; else return false; } if(Y==1) { if(x-right<=r)return true; else return false; } if(Y==2) { distx=x-right; disty=y-down; if(distx+disty<=r)return true; else return false; } } } bool is_all_in(int left,int down,int right,int up,int x,int y,int r) { if(left>right||down<up)return false; int distx,disty; distx=abs(left-x); disty=abs(down-y); if(distx+disty>r)return false; disty=abs(y-up); if(distx+disty>r)return false; distx=abs(x-right); disty=abs(down-y); if(distx+disty>r)return false; disty=abs(y-up); if(distx+disty>r)return false; return true; } int first[140000],target[10000000],next[10000000],value[10000000],countn=1; void AddEdge(int u,int v,int length) { target[countn]=v; value[countn]=length; next[countn]=first[u]; first[u]=countn++; } int from,posx,posy,r,v,A[160][160],B[160][160]; void build(int o,int left,int down,int right,int up,bool cut) { if(left==right&&up==down)return; int lc=o<<1,rc=lc+1,m; AddEdge(o,lc,0); AddEdge(o,rc,0); if((cut||left==right)&&up!=down) { m=up+down>>1; build(lc,left,m,right,up,!cut); build(rc,left,down,right,m+1,!cut); } else { m=left+right>>1; build(lc,left,down,m,up,!cut); build(rc,m+1,down,right,up,!cut); } } void update(int o,int left,int down,int right,int up,bool cut) { if(is_all_in(left,down,right,up,posx,posy,r)) { AddEdge(from,o,v); return; } int lc=o<<1,rc=lc+1,m; if((cut||left==right)&&up!=down) { m=up+down>>1; if(is_in(left,m,right,up,posx,posy,r))update(lc,left,m,right,up,!cut); if(is_in(left,down,right,m+1,posx,posy,r))update(rc,left,down,right,m+1,!cut); } else { m=left+right>>1; if(is_in(left,down,m,up,posx,posy,r))update(lc,left,down,m,up,!cut); if(is_in(m+1,down,right,up,posx,posy,r))update(rc,m+1,down,right,up,!cut); } } int query(int o,int left,int down,int right,int up,bool cut) { if(left==right&&up==down)return o; int lc=o<<1,rc=lc+1,m; if((cut||left==right)&&up!=down) { m=up+down>>1; if(posy<=m)return query(lc,left,m,right,up,!cut); else return query(rc,left,down,right,m+1,!cut); } else { m=left+right>>1; if(posx<=m)return query(lc,left,down,m,up,!cut); else return query(rc,m+1,down,right,up,!cut); } } struct State { int pos,cost; bool operator < (const State &b)const { return cost>b.cost; } }take,put; priority_queue<State> PQ; int fastest[140000],N,M; int Dijkstra(int x0,int y0,int x1,int y1) { memset(fastest,127,sizeof(fastest)); while(!PQ.empty())PQ.pop(); int s,t; posx=x0; posy=y0; s=query(1,1,N,M,1,0); posx=x1; posy=y1; t=query(1,1,N,M,1,0); put.pos=s; put.cost=0; PQ.push(put); while(!PQ.empty()) { take=PQ.top(); PQ.pop(); if(fastest[take.pos]<take.cost)continue; if(take.pos==t) { return take.cost; } for(int i=first[take.pos];i;i=next[i]) { if(fastest[target[i]]>take.cost+value[i]) { fastest[target[i]]=take.cost+value[i]; put.pos=target[i]; put.cost=fastest[put.pos]; PQ.push(put); } } } return -1; } int main() { N=read(),M=read(); build(1,1,N,M,1,0); for(int i=1;i<=N;i++) { for(int j=1;j<=M;j++) { B[i][j]=read(); } } for(int i=1;i<=N;i++) { for(int j=1;j<=M;j++) { A[i][j]=read(); if(B[i][j]) { posx=j; posy=i; r=B[i][j]; v=A[i][j]; from=query(1,1,N,M,1,0); update(1,1,N,M,1,0); } } } int y0=read(),x0=read(),y1=read(),x1=read(),y2=read(),x2=read(); int XtoY,YtoX,XtoZ,ZtoX,YtoZ,ZtoY,ans=INF; char place=\'N\'; XtoY=Dijkstra(x0,y0,x1,y1); YtoX=Dijkstra(x1,y1,x0,y0); XtoZ=Dijkstra(x0,y0,x2,y2); ZtoX=Dijkstra(x2,y2,x0,y0); YtoZ=Dijkstra(x1,y1,x2,y2); ZtoY=Dijkstra(x2,y2,x1,y1); if(~YtoX&&~ZtoX) { if(YtoX+ZtoX<ans) { place=\'X\'; ans=YtoX+ZtoX; } } if(~ZtoY&&~XtoY) { if(ZtoY+XtoY<ans) { place=\'Y\'; ans=XtoY+ZtoY; } } if(~XtoZ&&~YtoZ) { if(XtoZ+YtoZ<ans) { place=\'Z\'; ans=XtoZ+YtoZ; } } if(place==\'N\')puts("NO"); else printf("%c\n%d",place,ans); }
P.S.话说这个跨度有些大吧！！？