图论当中的术语,假设G=(V,E)和G1=(V1,E1)是两个图,如果存在一个双射m:V→V1,使得对所有的x,y∈V均有xy∈E等价于m(x)m(y)∈E1,则称G和G1是同构的,这样的一个映射m称之为一个同构,如果G=G1,则称他为一个自同构。
简单来说,同构图的结点数必须相同,结构必须相同。
如图3.6,第一个图形和第二个图形的区别在于环的数量。第一个图形为一个环,第二个为两个环,所以不是同构图。
若删去z1和u1,删去v1和w1,连接z1和w1,成为一个v1u1的链和z1w1x1y1的环,依旧不是同构图,因为必须环数相同,链数相同。
但这还是缺少一个条件,比如图形A存在两个环a1和a2,a1有3个结点,a2有5个结点,图形B也有两个环,b1有4个结点,b2有4个结点,依旧不是同构图,这里的条件就是环上或链上的借点数相同,和结点顺序无关。
引入例题,HDU3926-Hand in Hand ,判断两次组成的图形是否是同构图。
思路之一:通过并查集确定环数/链数,和环内/链内的人数,再排序进行比较。
排序时按照人数排序,若人数相同要按照状态排序。注意这几点或许会比较容易过。
请先自己进行尝试,尝试后再参考代码。
1 #include<iostream> 2 #include<cstring> 3 #include<cstdio> 4 #include<math.h> 5 #include<vector> 6 #include<algorithm> 7 #include<queue> 8 #include<set> 9 using namespace std; 10 int pre[10100]; 11 struct e{ 12 int a,b; 13 }; 14 e s1[10010]; 15 e s2[10010]; 16 int find(int x) 17 { 18 while(x!=pre[x]) 19 x=pre[x]; 20 return x; 21 } 22 int cmp(e a,e b){ 23 if(a.a==b.a) return a.b>b.b; 24 else return a.a>b.a; 25 } 26 void init(int n) 27 { 28 for(int i=1;i<=n;i++) 29 pre[i]=i; 30 } 31 int main() 32 { 33 int t,cas=1;; 34 scanf("%d",&t); 35 while(t--) 36 { 37 for(int i=1;i<10010;i++) 38 { 39 s1[i].a=1;s1[i].b=0; 40 s2[i].a=1;s2[i].b=0;//最开始每个都是独立的,默认为链 41 } 42 bool flag=false; 43 int n1,m1,n2,m2; 44 45 scanf("%d%d",&n1,&m1); 46 init(n1); 47 for(int i=0;i<m1;i++) 48 { 49 int a,b; 50 scanf("%d%d",&a,&b); 51 int dx=find(a); 52 int dy=find(b); 53 if(dx!=dy) 54 { 55 pre[dx]=dy; 56 s1[dy].a+=s1[dx].a; 57 s1[dx].a=0;//把拉手的孩子数量加起来,下同 58 } 59 else s1[dy].b=1;//成环 60 } 61 62 scanf("%d%d",&n2,&m2); 63 init(n2); 64 for(int i=0;i<m2;i++) 65 { 66 int a,b; 67 scanf("%d%d",&a,&b); 68 int dx=find(a); 69 int dy=find(b); 70 if(dx!=dy) 71 { 72 pre[dx]=dy; 73 s2[dy].a+=s2[dx].a; 74 s2[dx].a=0; 75 } 76 else s2[dy].b=1; 77 } 78 if(n1==n2){ 79 80 sort(s1+1,s1+n1+1,cmp); 81 sort(s2+1,s2+n2+1,cmp);//排序,若孩子的数量相同则对是否是环进行排序,这里要注意 82 83 for(int i=0;i<n1;i++) 84 if(s1[i].a!=s2[i].a||s1[i].b!=s2[i].b) {//判断数量,状态 85 flag=true; 86 break; 87 } 88 } 89 if(n1!=n2) flag=true; 90 91 if(flag) printf("Case #%d: NO\n",cas++); 92 else printf("Case #%d: YES\n",cas++); 93 } 94 return 0; 95 }