POJ 1182 食物链

Description

动物王国中有三类动物A,B,C,这三类动物的食物链构成了有趣的环形。A吃B, B吃C,C吃A。 
现有N个动物,以1-N编号。每个动物都是A,B,C中的一种,但是我们并不知道它到底是哪一种。 
有人用两种说法对这N个动物所构成的食物链关系进行描述: 
第一种说法是"1 X Y",表示X和Y是同类。 
第二种说法是"2 X Y",表示X吃Y。 
此人对N个动物,用上述两种说法,一句接一句地说出K句话,这K句话有的是真的,有的是假的。当一句话满足下列三条之一时,这句话就是假话,否则就是真话。 
1) 当前的话与前面的某些真的话冲突,就是假话; 
2) 当前的话中X或Y比N大,就是假话; 
3) 当前的话表示X吃X,就是假话。 
你的任务是根据给定的N(1 <= N <= 50,000)和K句话(0 <= K <= 100,000),输出假话的总数。 

Input

第一行是两个整数N和K,以一个空格分隔。 
以下K行每行是三个正整数 D,X,Y,两数之间用一个空格隔开,其中D表示说法的种类。 
若D=1,则表示X和Y是同类。 
若D=2,则表示X吃Y。

Output

只有一个整数,表示假话的数目。

Sample Input

100 7
1 101 1 
2 1 2
2 2 3 
2 3 3 
1 1 3 
2 3 1 
1 5 5

Sample Output

3
 
 

基础回顾:

find()函数找根结点的递归写法如下:

  1. int find(int x)  
  2. {  
  3.     return x == pre[x] ? x : find(pre[x]);  
  4. }

下面我们通过这道题来讲解一下并查集的深层次应用。

输入:动物个数n以及k句话,接着输入k行,每一行形式为:d x y,在输入时可以先判断题目所说的条件2和3,即:        1> 若(x>n||y>n):即当前的话中,x或y比n大,则假话数目sum加1.        2> 若(d==2&&x==y):即当前的话表示x吃x,则假话数目sum加1.

而不属于这两种情况外的话语要利用并查集进行判断当前的话是否与此前已经说过的话相冲突.

  1. struct node  
  2. {  
  3.     int parent;                     //p[i].parent表示节点i的父节点  
  4.     int relation;                   //p[i].relation表示节点i与其父节点(即p[i].parent)的关系  
  5. }p[50010];  

   此处relation有三种取值(假设节点x的父节点为rootx,即p[x].parent=rootx):    p[x].relation=0   ……表示节点x与其父节点rootx的关系是:同类    p[x].relation=1   ……表示节点x与其父节点rootx的关系是:x被根结点rootx吃    p[x].relation=2   ……表示节点x与其父节点rootx的关系是:x吃根结点rootx

   

初始化函数为:

  1. void init(int n)  
  2. {  
  3.     for(int i = 1;i <= n; ++i)  
  4.     {  
  5.         p[i].parent = i;          //初始时集合编号就设置为自身  
  6.         p[i].relation = 0;        //因为p[i].parent=i,即节点i的父亲节点就是自身,所以此时节点i与其父亲节点的关系为同类(即p[i].relation=0)  
  7.     }  
  8. }  
下面详细讲解并查集的两个重要操作:查找和合并.
查找操作:
   在查找时节点不仅有父亲节点域parent[i], 而且还有表示节点与其父亲节点的关系域relation[i].
查找过程中对父亲节点域的处理和简单的并查集处理一样,即在查找过程中同时实现路径压缩,但正是由于路径压缩,使得表示节点与其父亲节点的关系域发生了变化 .(因为如果路径压缩之前节点x的父亲节点为rootx的话,那么在路径压缩之后节点x的父亲节点就变为了节点rootx的父亲节点rootxx,所以此时p[x].relation存储的应该是节点x与现在父亲节点rootxx的关系),此处可以画图理解一下:
POJ 1182 食物链_第1张图片POJ 1182 食物链_第2张图片
很明显查找之前节点x的父亲节点为rootx,假设此时p[x].relation=1(即表示x的父亲节点rootx吃x)且p[rootx].relation=0(即表示rootx和其父亲节点rootxx是同类),由这两个关系可以推出rootxx吃x,而合并以后节点x的父亲节点为rootxx(实现了路径压缩),且节点x的父亲节点rootxx吃x,即查找之后p[x].relation=1。

合并操作:

     在将元素x与y所在的集合合并时,假设元素x所在的集合编号为rootx,元素y所在的集合编号为rooty,合并时直接将集合rooty挂到集合rootx上,即p[rooty].parent=rootx,此时原来集合rooty中的根节点rooty的关系域也应随之发生变化,因为合并之前rooty的父亲节点就是其自身,故此时p[rooty].relation=0,而合并之后rooty的父亲节点为rootx,所以此时需判断rootx与rooty的关系,即更新p[rooty]的值,同理画图理解:      POJ 1182 食物链_第3张图片POJ 1182 食物链_第4张图片    此时假设p[x].relation=0(即x与rootx的关系是同类),p[y].relation=1(即rooty吃y),则有:         1> 输入d=1时,即输入的x和y是同类,则有上述关系可以推出rooty吃rootx,即p[rooty].relation=2;         2> 输入d=2时,即输入的x吃y,则有上述关系可以推出rooty与rootx是同类(因为rooty吃y,x吃y,则rooty与x是同类,又rootx与x是同类),即p[rooty].relation=0;    当然,这只是一种可能,其它的可能情况和上面一样分析。

当元素x与元素y在同一集合时,则不需要合并,因为此时x与y的父亲节点相同,可以分情况讨论:
        1> d=1时,即x与y是同类时,此时要满足这个要求,则必须满足p[x].relation=p[y].relation,这很容易推出来.
        2> d=2时,即表示x吃y,此时要满足这个要求,则也必须满足一定的条件,如x和rootx是同类(即p[x].relation=0),此时要满足x吃y,则必须满足rootx吃y,即p[y].relation=1,可以像上面一样画图来帮助理解.

关系域更新:

仔细再想想,rootx-x x-yy-rooty,是不是很像向量形式?于是我们可以大胆的从向量入手:

rootx     rooty

|         |

x    ~     y

对于集合里的任意两个元素xy而言,它们之间必定存在着某种联系,因为并查集中的元素均是有联系的(这点是并查集的实质,要深刻理解),否则也不会被合并到当前集合中。那么我们就把这2个元素之间的关系量转化为一个偏移量

由上面可知: x->y 偏移量0时 xy同类

x->y 偏移量1时 xy

x->y 偏移量2时 xy

有了这个假设,我们就可以在并查集中完成任意两个元素之间的关系转换了。

不妨继续假设,x的当前集合根节点rootxy的当前集合根节点rootyx->y偏移值为d-1(题中给出的询问已知条件)。为什么是d-1呢?因为题中说d=1时代表同类而我们自己设的relation=0表示同类。

(1)如果rootxrooty不相同,那么我们把rooty合并到rootx上,并且更新relation关系域的值(注意:p[i].relation表示i的根结点到i的偏移量!!!!(向量方向性一定不能搞错)

    此时 rootx->rooty = rootx->x + x->y + y->rooty,这一步就是大牛独创的向量思维模式

    上式进一步转化为:rootx->rooty = (relation[x]+d-1+3-relation[y])%3 = relation[rooty],(模3是保证偏移量取值始终在[0,2]间)

(2)如果rootxrooty相同(xy已经在一个集合中,不需要合并操作了,根结点相同),那么我们就验证x->y之间的偏移量是否与题中给出的d-1一致

    此时 x->y = x->rootx + rootx->y

    上式进一步转化为:x->y = (3-relation[x]+relation[y])%3     若一致则为真,否则为假。

#include <iostream> #include <stdio.h> using namespace std; struct node {     int parent;     int relation; }P[50001]; int set_find(int x) {     if(x!=P[x].parent)     {         int temp=P[x].parent;         P[x].parent=set_find(P[x].parent);         P[x].relation=(P[x].relation+P[temp].relation)%3;         return P[x].parent;     }     return x; } int main() {     int N,K,d,x,y,rootx,rooty;     int sum=0;     cin>>N>>K;     for(int i=1;i<=N;i++)     {         P[i].parent=i;         P[i].relation=0;     }     for(int i=0;i<K;i++)     {         scanf("%d%d%d",&d,&x,&y);         if(x>N||y>N||(d==2&&x==y))         {             sum++;             continue;         }         rootx=set_find(x);         rooty=set_find(y);         if(rootx!=rooty)         {             P[rooty].parent=rootx;             P[rooty].relation=(P[x].relation+d-1+3-P[y].relation)%3;         }         else         {             if(d==1&&P[x].relation!=P[y].relation)             {                 sum++;                 continue;             }             if(d==2&&((3-P[x].relation+P[y].relation)%3!=d-1))             {                 sum++;                 continue;             }         }     }     cout<<sum<<endl;     return 0; }

你可能感兴趣的:(POJ 1182 食物链)