[BZOJ3161]孤舟蓑笠翁

bzoj

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;
}
posted @ 2018-07-13 20:49  租酥雨  阅读(176)  评论(0编辑  收藏  举报