POJ_2749题解

14 篇文章 0 订阅
7 篇文章 0 订阅

解法:二分+2-SAT

二分最长路径长度,得出边,将边加上,再加上之前由于hate和friend的边,2-SAT判断下即可。

具体见代码


//time:454ms
//memory:3848KB
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<stack>
#include<queue>
#include<vector>
#include<iostream>
#define BUG printf("hehe\n");
using namespace std;

const int N=1010;
int n,a,b,D;
int pre[N],low[N],scc[N],cnt,num;
int x[N],y[N],sx1,sx2,sy1,sy2;
int d1[N],d2[N];
int fri1[N],fri2[N],hae1[N],hae2[N];
vector<int> G[N];
stack<int> s;

void DFS(int u)
{
	//BUG;
	pre[u]=low[u]=++num;
	s.push(u);
	for(int i=0;i<G[u].size();++i) {
		int v=G[u][i];
		if(!pre[v]) {
			DFS(v);
			low[u]=min(low[u],low[v]);
		}
		else if(!scc[v]) low[u]=min(low[u],pre[v]);
	}
	if(pre[u]==low[u]) {
		int tmp=-1;
		++cnt;
		while(tmp!=u) {
			tmp=s.top();
			s.pop();
			scc[tmp]=cnt;
		}
	}
}

bool Targin(int mid)
{
	cnt=0,num=0;
	int i,p,q,j;
	for(i=1;i<=(n<<1|1);++i) {
		G[i].clear();
	}
	memset(pre,0,sizeof pre);
	memset(scc,0,sizeof scc);
	memset(low,0,sizeof low);
	while(!s.empty()) s.pop();
	for(i=1;i<=a;++i) {               //hate产生的边
		p=hae1[i];
		q=hae2[i];
		G[p<<1].push_back(q<<1|1);
		G[q<<1].push_back(p<<1|1);
		G[q<<1|1].push_back(p<<1);
		G[p<<1|1].push_back(q<<1);
	}
	for(i=1;i<=b;++i) {                 //friend产生的边
		p=fri1[i];
		q=fri2[i];
		G[p<<1].push_back(q<<1);
		G[q<<1].push_back(p<<1);
		G[q<<1|1].push_back(p<<1|1);
		G[p<<1|1].push_back(q<<1|1);
	}

	for(i=1;i<=n;++i)                         //两个牛棚由于路径长度大于二分枚举长度而产生的边
		for(j=i+1;j<=n;++j) {
			if(d1[i]+d1[j]>mid) {
				G[i<<1].push_back(j<<1|1);
				G[j<<1].push_back(i<<1|1);
			}
			if(d2[i]+d2[j]>mid) {
				G[i<<1|1].push_back(j<<1);
				G[j<<1|1].push_back(i<<1);
			}
			if(d1[i]+d2[j]+D>mid) {
				G[i<<1].push_back(j<<1);
				G[j<<1|1].push_back(i<<1|1);
			}
			if(d2[i]+d1[j]+D>mid) {
				G[i<<1|1].push_back(j<<1|1);
				G[j<<1].push_back(i<<1);
			}
		}
	for(i=1;i<=(n<<1|1);++i)
		if(!pre[i]) DFS(i);
	for(i=1;i<=n;++i)
		if(scc[i<<1]==scc[i<<1|1])
			return false;
	return true;
	//BUG;
}

int main()
{
	int L,R,mid,ans;
	while(scanf("%d%d%d",&n,&a,&b)==3) {
		scanf("%d%d%d%d",&sx1,&sy1,&sx2,&sy2);
		D=abs(sx1-sx2)+abs(sy1-sy2);
		for(int i=1;i<=n;++i) {
			scanf("%d%d",x+i,y+i);
			d1[i]=abs(x[i]-sx1)+abs(y[i]-sy1);
			d2[i]=abs(x[i]-sx2)+abs(y[i]-sy2);
		}
		for(int i=1;i<=a;++i)
			scanf("%d%d",hae1+i,hae2+i);
		for(int i=1;i<=b;++i)
			scanf("%d%d",fri1+i,fri2+i);
		L=1,R=12000000;
		ans=12000000;
		while(L<=R) {
			mid=(R+L)>>1;
			if(Targin(mid)) R=mid-1,ans=mid;
			else L=mid+1;
		}
		if(12000000<=ans) printf("-1\n");
		else printf("%d\n",ans);
	}
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值