P1955 [NOI2015] 程序自动分析

题目链接

luogu P1955

题面

题目描述

在实现程序自动分析的过程中,常常需要判定一些约束条件是否能被同时满足。

考虑一个约束满足问题的简化版本:假设 x 1 , x 2 , x 3 , ⋯ x_1,x_2,x_3,\cdots x1,x2,x3, 代表程序中出现的变量,给定 n n n 个形如 x i = x j x_i=x_j xi=xj x i ≠ x j x_i\neq x_j xi=xj 的变量相等/不等的约束条件,请判定是否可以分别为每一个变量赋予恰当的值,使得上述所有约束条件同时被满足。例如,一个问题中的约束条件为: x 1 = x 2 , x 2 = x 3 , x 3 = x 4 , x 4 ≠ x 1 x_1=x_2,x_2=x_3,x_3=x_4,x_4\neq x_1 x1=x2,x2=x3,x3=x4,x4=x1,这些约束条件显然是不可能同时被满足的,因此这个问题应判定为不可被满足。

现在给出一些约束满足问题,请分别对它们进行判定。

输入格式

输入的第一行包含一个正整数 t t t,表示需要判定的问题个数。注意这些问题之间是相互独立的。
对于每个问题,包含若干行:
第一行包含一个正整数 n n n,表示该问题中需要被满足的约束条件个数。接下来 n n n 行,每行包括三个整数 i , j , e i,j,e i,j,e,描述一个相等/不等的约束条件,相邻整数之间用单个空格隔开。若 e = 1 e=1 e=1,则该约束条件为 x i = x j x_i=x_j xi=xj。若 e = 0 e=0 e=0,则该约束条件为 x i ≠ x j x_i\neq x_j xi=xj

输出格式

输出包括 t t t 行。
输出文件的第 k k k 行输出一个字符串 YES 或者 NO(字母全部大写),YES 表示输入中的第 k k k 个问题判定为可以被满足,NO 表示不可被满足。

样例 #1

样例输入 #1

2
2
1 2 1
1 2 0
2
1 2 1
2 1 1

样例输出 #1

NO
YES

样例 #2

样例输入 #2

2
3
1 2 1
2 3 1
3 1 1
4
1 2 1
2 3 1
3 4 1
1 4 0

样例输出 #2

YES
NO

提示

【样例解释1】

在第一个问题中,约束条件为: x 1 = x 2 , x 1 ≠ x 2 x_1=x_2,x_1\neq x_2 x1=x2,x1=x2。这两个约束条件互相矛盾,因此不可被同时满足。
在第二个问题中,约束条件为: x 1 = x 2 , x 1 = x 2 x_1=x_2,x_1 = x_2 x1=x2,x1=x2。这两个约束条件是等价的,可以被同时满足。

【样例说明2】

在第一个问题中,约束条件有三个: x 1 = x 2 , x 2 = x 3 , x 3 = x 1 x_1=x_2,x_2= x_3,x_3=x_1 x1=x2,x2=x3,x3=x1。只需赋值使得 x 1 = x 2 = x 3 x_1=x_2=x_3 x1=x2=x3,即可同时满足所有的约束条件。
在第二个问题中,约束条件有四个: x 1 = x 2 , x 2 = x 3 , x 3 = x 4 , x 4 ≠ x 1 x_1=x_2,x_2= x_3,x_3=x_4,x_4\neq x_1 x1=x2,x2=x3,x3=x4,x4=x1。由前三个约束条件可以推出 x 1 = x 2 = x 3 = x 4 x_1=x_2=x_3=x_4 x1=x2=x3=x4,然而最后一个约束条件却要求 x 1 ≠ x 4 x_1\neq x_4 x1=x4,因此不可被满足。

【数据范围】

注:实际上 n ≤ 1 0 6 n\le 10^6 n106

题意

t t t个问题,每个问题给出 n n n条约束条件,其格式为 x x x y y y e e e,分两种情况

  • x x x y y y 1 1 1 → \to a x = a y a_x=a_y ax=ay
  • x x x y y y 0 0 0 → \to a x ≠ a y a_x\not=a_y ax=ay

求这 n n n条约束条件是否可以互相满足 ( Y E S , N O ) (YES,NO) (YES,NO)
example:
a 1 = a 2 , a 2 = a 3 , a 1 ≠ a 3 a_1=a_2,a_2=a_3,a_1\not=a_3 a1=a2,a2=a3,a1=a3是不可以互相满足的
a 1 = a 2 , a 1 = a 3 , a 2 = a 3 a_1=a_2,a_1=a_3,a_2=a_3 a1=a2,a1=a3,a2=a3是可以互相满足的

思路

对于 a x = a y a_x=a_y ax=ay,我们可以让 a y a_y ay 的父亲指向 a x a_x ax
当我们判断 a x a_x ax是否等于 a y a_y ay时,就可以 f i n d find find两个变量最老的“祖先”,判断二者是不是同一个家族的成员(同一个量)
好像已经完成了,但又好像没有……


是的,没有完成
我们回到题目看一下数据大小
在这里插入图片描述
W O W ! WOW! WOW!
1 0 9 10^9 109大小的编号,直接用并查集的父亲数组储存肯定会 M L E MLE MLE
可是 n n n只有 1 0 6 10^6 106,这就说明了父亲数组里有 1 0 9 − 1 0 6 10^9-10^6 109106个空内存,纯属浪费
所以我们只要将编号都压缩排序后再储存进父亲数组就可以了
那怎么压缩编号呢?很简单,离散化
比如说有这样一组问题

3
1 999999999 1
1 888888888 0
888888888 999999999 1
//压缩前^^^
3
1 3 1
1 2 0
2 3 1
//压缩后^^^

如果不用离散化就要开 1 0 9 10^9 109的父亲数组
可是用了离散化压缩编号后就成了非常小的编号了,就可以不用开大数组了
这就结束了吗?


不,还没有
还有一个重要的细节
我们在对压缩后的约束条件进行并查集处理前
需要将每一个约束条件的 e e e进行排序,将 e e e 1 1 1的排在前面, e e e 0 0 0的排在后面
为什么呢?
因为所有的约束条件没有先后之分
换句话来说,就是先将相等的条件 ( e = = 1 ) (e==1) (e==1)用并查集连接起来,然后用不相等的条件 ( e = = 0 ) (e==0) (e==0)进行判断当前的条件是否成立
其实就是将相等的条件看成已知,将不相等的条件看成要判断的条件


离散化细节见代码注释

#include<bits/stdc++.h>
using namespace std;
#define int long long
int t,n,cnt,len;
bool flag;
struct Node{
	int x,y,e;
}node[1000010];
int rak[2000020],dx[1000010],dy[1000010],fa[2000010];
int find(int x){//并查集 find"父亲" 
	if(fa[x]!=x)fa[x]=find(fa[x]);
	return fa[x];
}
void merge(int x,int y){//并查集连接 
	fa[find(x)]=find(y);
	return;
}
bool cmp(Node a,Node b){//对e进行排序 
	return a.e>b.e;
}
signed main(){
	scanf("%lld",&t);
	while(t--){
		scanf("%lld",&n);
		flag=true,cnt=0;//初始化 
		for(int i=1;i<=n;i++)scanf("%lld%lld%lld",&node[i].x,&node[i].y,&node[i].e),rak[++cnt]=node[i].x,rak[++cnt]=node[i].y;
		//初始化rak数组,rak数组的作用是排序去重,方便二分查找元素压缩后对应的编号 
		sort(rak+1,rak+1+cnt);//先排序,因为unique(去重)函数需要对有序数组才能起作用 
		len=unique(rak+1,rak+1+cnt)-(rak+1);//去重,得到去重后的数组长度 
		for(int i=1;i<=len;i++)fa[i]=i;//初始化父亲数组 
		for(int i=1;i<=n;i++)node[i].x=lower_bound(rak+1,rak+1+len,node[i].x)-rak,node[i].y=lower_bound(rak+1,rak+1+len,node[i].y)-rak;
		//用二分查找元素压缩的对应编号,放入结构体,方便对e排序 
		sort(node+1,node+1+n,cmp);//对e进行排序 
		for(int i=1;i<=n;i++){
			if(node[i].e)merge(node[i].x,node[i].y);//若当前约束条件为 x等于y,连接x和y 
			else if(find(node[i].x)==find(node[i].y)){
				//若当前约束条件为 x不等于y,但是在前面的约束条件中 直接 或 间接 表达了 x等于y
				//推翻了当前约束条件,则当前问题的约束条件不可以互相满足 
				flag=false;
				break;
			}
		}
		if(flag)printf("YES\n");
		else printf("NO\n");
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值