[BZOJ3161]孤舟蓑笠翁
sol
先预处理每个状态是否合法。记\(dis[i][j]\)表示从任意一个绝招点到达左手在\(i\)右手在\(j\)这个状态的最小步数,以及\(vis[i][j]\)表示是从哪一个绝招点转移过来的。这一步可以\(bfs\)实现,复杂度\(O((n+m)^2)\)。
然后再枚举两个相邻状态,如果两个状态的\(vis\)不同,那么就可以更新这两个绝招点的答案。
可以理解为从每个绝招点开始扩张,当遇到“接壤”时即意味着两个绝招点在最短路上相遇。
code
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
int gi(){
int x=0,w=1;char ch=getchar();
while ((ch<'0'||ch>'9')&&ch!='-') ch=getchar();
if (ch=='-') w=0,ch=getchar();
while (ch>='0'&&ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return w?x:-x;
}
#define cmin(a,b) (a>b?a=b:a)
const int N = 1005;
struct Graph{
int to[N*10],nxt[N*10],head[N],cnt;
void link(int u,int v){
to[++cnt]=v,nxt[cnt]=head[u],head[u]=cnt;
}
}G1,G2;
int n,m,k,mn,mx,x[N],y[N],ok[N][N],q[2][N*N],hd,tl,vis[N][N],dis[N][N],ans[N*100];
int main(){
n=gi();m=gi();mn=gi();mx=gi();
for (int i=1;i<=n;++i) x[i]=gi(),y[i]=gi();
for (int i=1;i<=n;++i)
for (int j=1;j<=n;++j){
int len=abs(x[i]-x[j])+abs(y[i]-y[j]);
if (len>=mn&&len<=mx) ok[i][j]=1;
}
k=gi();
for (int i=1;i<=k;++i){
q[0][i]=gi();q[1][i]=gi();
vis[q[0][i]][q[1][i]]=i;
}
for (int i=1;i<=m;++i){
int u=gi(),v=gi(),t=gi();
if (!t) G1.link(u,v),G1.link(v,u);
else G2.link(u,v),G2.link(v,u);
}
for (int i=1;i<=n;++i) G1.link(i,i),G2.link(i,i);
hd=1,tl=k;
while (hd<=tl){
int u=q[0][hd],v=q[1][hd++];
for (int i=G1.head[u];i;i=G1.nxt[i])
for (int j=G2.head[v];j;j=G2.nxt[j]){
int uu=G1.to[i],vv=G2.to[j];
if (!ok[uu][vv]||vis[uu][vv]) continue;
vis[uu][vv]=vis[u][v];dis[uu][vv]=dis[u][v]+1;
q[0][++tl]=uu;q[1][tl]=vv;
}
}
memset(ans,63,sizeof(ans));
for (int u=1;u<=n;++u)
for (int v=1;v<=n;++v)
for (int i=G1.head[u];i;i=G1.nxt[i])
for (int j=G2.head[v];j;j=G2.nxt[j]){
int uu=G1.to[i],vv=G2.to[j];
if (!vis[u][v]||!vis[uu][vv]) continue;
if (vis[u][v]!=vis[uu][vv]){
int len=dis[u][v]+dis[uu][vv]+1;
cmin(ans[vis[u][v]],len);
cmin(ans[vis[uu][vv]],len);
}
}
for (int i=1;i<=k;++i) printf("%d\n",ans[i]==ans[0]?-1:ans[i]);
return 0;
}