旅游胜地(二分+2-sat)

2695: 旅游胜地
题目描述
有一旅游胜地,可以看做是一个 n 个点、m 条边的连通无向图。每个点都有一个点权,编号为 i 的点的点权为 ai,表示这个点的景观值。
开发商为了能获取收益,准备选择一些点建设商业区,但是会让这些点的景观值大打折扣。具体地说,如果选择编号为 i 的点建设商业区,则会让该点的景观值变为 bi,其中 bi≤ai。
不过对于游客来讲,只要游览相邻的两个点时,景观值的差距不要太大就好。即确定好需要建设商业区的点后,记编号为 i 的点最终的点权为 wi,只要所有边 u,v 的 ∣wu−wv∣ 中最大值尽可能小即可。或者说,最小化有边相连的景观值差距的最大值。
请你帮助开发商计算所有可能的方案中,景观值最大差距的最小值。
输入
第一行,两个正整数 n,m(2≤n≤105, n−1≤m≤min(n(n−1)/2,105) ),表示图的点数和边数。
第二行,共 n 个正整数 ai (1≤ai≤109),表示编号为 i 的点不建设商业区的景观值。
第三行,共 n 个正整数 bi (1≤bi≤109, bi≤ai),表示编号为 i 的点建设商业区后的景观值。
接下来 m 行,每行两个正整数 u,v,表示 u 和 v 有一条边相连。
数据保证图是连通的,也保证没有重边或自环。
输出
一个整数,景观值最大差距的最小值。
样例输入 Copy

4 3
10 7 5 5
6 2 1 2
2 4
3 2
2 1

样例输出 Copy

2
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<stack>
using namespace std;
const int N=4e5+10;
const int inf=0x3f3f3f3f;
typedef long long ll;
int n,m; 
int a[N],b[N],d[N][2];
int dfn[N],low[N],col[N],stk[N],vis[N],tot,top,idx,cnt,co;
vector<int>ve[N];
void Tarjan(int x){
	dfn[x]=low[x]=++idx;
	vis[x]=1;
	stk[++top]=x;
	for(int i=0;i<ve[x].size();i++){
		int to=ve[x][i];
		if(!dfn[to]) {
			Tarjan(to);
			low[x]=min(low[x],low[to]);
		}
		else if(vis[to]){
			low[x]=min(low[x],dfn[to]);
		}
	}
	if(low[x]==dfn[x]){
		tot++;
		do{
			col[stk[top]]=tot;
			vis[stk[top]]=0;
		}while(stk[top--]!=x);
	}
}
int judge(int mid){
	idx=top=tot=0;
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
	memset(col,0,sizeof(col));
	memset(vis,0,sizeof(vis));
	memset(stk,0,sizeof(stk));
	for(int i=1;i<=n*2;i++){
		ve[i].clear();
	}
	
	for(int i=1;i<=m;i++){
		int x=d[i][0],y=d[i][1];
//		if(abs(a[x]-a[y])>=mid) ve[x].push_back(y),ve[y].push_back(x);
//		if(abs(a[x]-b[y])>=mid) ve[x].push_back(y+n),ve[y+n].push_back(x);
//		if(abs(b[x]-a[y])>=mid) ve[x+n].push_back(y),ve[y].push_back(x+n);
//		if(abs(b[x]-b[y])>=mid) ve[x+n].push_back(y+n),ve[y+n].push_back(x+n);
		if(abs(a[x]-a[y])>mid) ve[x].push_back(y+n),ve[y].push_back(x+n);
		if(abs(a[x]-b[y])>mid) ve[x].push_back(y),ve[y+n].push_back(x+n);
		if(abs(b[x]-a[y])>mid) ve[x+n].push_back(y+n),ve[y].push_back(x);
		if(abs(b[x]-b[y])>mid) ve[x+n].push_back(y),ve[y+n].push_back(x);
	}
	
	for(int i=1;i<=n*2;i++){
		if(!dfn[i]) Tarjan(i);
	}
	
	for(int i=1;i<=n;i++){
		if(col[i]==col[n+i]) return 0;
	}
	return 1;
}
int main()
{
	int i,j,k,x,y;
	scanf("%d %d",&n,&m);
	for(i=1;i<=n;i++) scanf("%d",&a[i]);
	for(i=1;i<=n;i++) scanf("%d",&b[i]);
	for(i=1;i<=m;i++){
		scanf("%d %d",&d[i][0],&d[i][1]);
	}
	int l=0,r=1e9,mid;
	while(l<=r){
		mid=l+r>>1;
		if(judge(mid)){
			r=mid-1;
		}
		else {
			l=mid+1;
		}
	}
	printf("%d",l);
	return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值