HDU3926 Hand in Hand, 图形同构判断

图形比较特殊, 每个点的度最大是2, 故图形由链(list)和圆(circle), 其中度为0的点归为链.

由此, 我们可以对每个点进行DFS, 判断这个点所在图形是链还是圆.

然后就可以根据两个图的链的数目是否相等, 链的节点是否相等; 圆的数目是否相等, 圆的节点是否相等.

/******************************************************************************* # Author : Neo Fung # Email : [e-mail] # Last modified: 2012-02-05 21:14 # Filename: HDU3926 Hand in Hand.cpp # Description :  ******************************************************************************/#ifdef _MSC_VER#define DEBUG#endif#include <fstream>#include <stdio.h>#include <iostream>#include <string.h>#include <string>#include <limits.h>#include <algorithm>#include <math.h>#include <numeric>#include <functional>#include <ctype.h>#include <vector>#define MAX 10010using namespace std;bool inline cmp(const vector<int> &lhs,const vector<int> &rhs){  return lhs.size()<rhs.size();}struct GRAPH{  vector<int>  adj[MAX],list[MAX],cirle[MAX];  bool used[MAX];  bool visit[MAX];  int n,lcnt,ccnt;  void init(const int &t)  {    n=t;    lcnt=ccnt=0;    memset(visit,false,sizeof(visit));    memset(used,false,sizeof(used));    for (int i=0;i<MAX;++i)    {      list[i].clear();      cirle[i].clear();      adj[i].clear();    }  }  void add_adj(const int &x,const int &y)  {    used[x]=true;    bool  flag=true;    for(size_t i=0;i<adj[x].size();++i)      if(adj[x][i]==y)        flag=false;    if(flag)      adj[x].push_back(y);  }  bool dfs(const int &u)  {    visit[u]=true;    list[lcnt].push_back(u);    cirle[ccnt].push_back(u);    bool flag=true;    for(size_t j=0;j<adj[u].size();++j)      if(!visit[adj[u][j]])        flag=false;    if(flag)      if(adj[u].size()==1)        return false;      else        return true;    bool ans=false;    for(size_t j=0;j<adj[u].size();++j)      if(!visit[adj[u][j]])        ans = ans||dfs(adj[u][j]);    return ans;  }  void func()  {    for(int i=0;i<MAX;++i)    {      if(!visit[i] && used[i])      {        bool ans=dfs(i);        if(!ans)        {          cirle[ccnt].clear();          ++lcnt;        }        else        {          list[lcnt].clear();          ++ccnt;        }      }    }    sort(list,list+lcnt,cmp);    sort(cirle,cirle+ccnt,cmp);  }};bool check(const vector<int> *lhs,const int &lcnt,const vector<int> *rhs,const int &rcnt){  if(lcnt!=rcnt)    return false;  for(int i=0;i<lcnt;++i)    if(lhs[i].size()!=rhs[i].size())      return false;  return true;}int main(void){#ifdef DEBUG    freopen("../stdin.txt","r",stdin);  freopen("../stdout.txt","w",stdout); #endif    int ncases,n,m;  int a,b;  GRAPH lhs,rhs;  scanf("%d",&ncases);  for(int nc=1;nc<=ncases;++nc)  {        scanf("%d%d",&n,&m);    lhs.init(n);    while(m--)    {      scanf("%d%d",&a,&b);      lhs.add_adj(a,b);      lhs.add_adj(b,a);    }    lhs.func();    scanf("%d%d",&n,&m);    rhs.init(n);    while(m--)    {      scanf("%d%d",&a,&b);      rhs.add_adj(a,b);      rhs.add_adj(b,a);    }    rhs.func();    bool ans=check(lhs.list,lhs.lcnt,rhs.list,rhs.lcnt) &&       check(lhs.cirle,lhs.ccnt,rhs.cirle,rhs.ccnt);    printf("Case #%d: ",nc);    if(!ans)      printf("NO\n");    else      printf("YES\n");  }  return 0;}