P1955 [NOI2015]程序自动分析(并查集+离散化

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;
    }
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值