P1955 [NOI2015]程序自动分析(并查集+离散化
题目描述
传送门
在实现程序自动分析的过程中,常常需要判定一些约束条件是否能被同时满足。
考虑一个约束满足问题的简化版本:假设 x1,x2,x3,⋯ 代表程序中出现的变量,给定 n 个形如 xi=xj或 xi≠xj 的变量相等/不等的约束条件,请判定是否可以分别为每一个变量赋予恰当的值,使得上述所有约束条件同时被满足。例如,一个问题中的约束条件为:x1=x2,x2=x3,x3=x4,x4≠x1,这些约束条件显然是不可能同时被满足的,因此这个问题应判定为不可被满足。
现在给出一些约束满足问题,请分别对它们进行判定。
输入格式
输入的第一行包含一个正整数 t,表示需要判定的问题个数。注意这些问题之间是相互独立的。
对于每个问题,包含若干行:
第一行包含一个正整数 n,表示该问题中需要被满足的约束条件个数。接下来 n 行,每行包括三个整数 i,j,e,描述一个相等/不等的约束条件,相邻整数之间用单个空格隔开。若 e=1,则该约束条件为 xi=xj。若e=0,则该约束条件为 xi≠xj。
输出格式
输出包括 ttt 行。
输出文件的第 k 行输出一个字符串 YES 或者 NO(字母全部大写),YES 表示输入中的第 k 个问题判定为可以被满足,NO 表示不可被满足。
输入输出样例
输入 #1
2
2
1 2 1
1 2 0
2
1 2 1
2 1 1
输出 #1
NO
YES
输入 #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】
在第一个问题中,约束条件为:x1=x2,x1≠x2。这两个约束条件互相矛盾,因此不可被同时满足。
在第二个问题中,约束条件为:x1=x2,x1=x2。这两个约束条件是等价的,可以被同时满足。
【样例说明2】
在第一个问题中,约束条件有三个:x1=x2,x2=x3,x3=x1。只需赋值使得 x1=x2=x3,即可同时满足所有的约束条件。
在第二个问题中,约束条件有四个:x1=x2,x2=x3,x3=x4,x4≠x1。由前三个约束条件可以推出 x1=x2=x3=x4,然而最后一个约束条件却要求 x1≠x4,因此不可被满足。
【数据范围】
注:实际上 n≤10^6 。
这题i,j的范围是10的九次方,如果用并查集开数组,肯定会爆!
那么很容易想到:
离!散!化!
总得来说离散化有三步走战略:
1.去重(可以用到unique去重函数)
2.排序
3.二分索引(可以用到lower_bound函数)
这样就把数据范围缩小到10的五次方。
先排序,把所有e=1的操作放在前面,然后再进行e=0的操作,在进行e=1的操作的时候,我们只要把它约束的两个变量放在同一个集合里面即可。在e=0,即存在一条不相等的约束条件,对于它约束的两个变量,如果在一个集合里面,那就不可能满足!如不相等的约束条件都满足,那就YES。
之后就很简单了!
最后注意并查集的查询函数,要查完之后将每一个数的祖先再次优化。这里TLE了几次。枯了wuwu
int find(int n){
if(f[n]==n)
return n;
return f[n]=find(f[n]);
}
喜闻乐见的 代码:
#include <iostream>
#include <algorithm>
#include <cstring>
#include <math.h>
#include <string>
#include <vector>
#include <map>
#include <stack>
#include <cstdio>
#include <list>
#include <queue>
#include <set>
using namespace std;
#define ll long long
int f[1000300],num[2000030];//flag
int find(int n){
if(f[n]==n)
return n;
return f[n]=find(f[n]);
}
typedef struct{
int x,y,a;
}node;
node k[1003000];
bool cmp(node x,node y){
return x.a>y.a;
}
int main(){
int t,n,i,length,flag,xx,yy;
cin>>t;
while(t--){
cin>>n;
length=0;
memset(f,0, sizeof(f));
for(i=1;i<=n;i++){
cin>>k[i].x>>k[i].y>>k[i].a;
num[length++]=k[i].x;
num[length++]=k[i].y;
}
sort(num,num+length);
length=unique(num,num+length)-num;
for(i=1;i<=n;i++){
k[i].x=lower_bound(num,num+length,k[i].x)-num;
k[i].y=lower_bound(num,num+length,k[i].y)-num;
}
//离散化完成
flag=1;
for(i=1;i<=length;i++) f[i]=i;
sort(k+1,k+n+1,cmp);
for(i=1;i<=n;i++){
xx=find(k[i].x);
yy=find(k[i].y);
if(k[i].a==1) {
if(xx!=yy)
f[xx]=yy;
}
else {
if(xx==yy){
flag=0;
goto A;
}
}
}
A:
if(flag==1) cout<<"YES"<<endl;
else cout<<"NO"<<endl;
}
}