100851K King’s Inspection

传送门

题目大意

给你一张图，求这张图的汉密尔顿回路。

分析

因为m≤n+20，所以如果存在回路一定是在一个环中加入了至多20条边。我们先考虑dfs，但我们发现如果出现图1这种情况就会是复杂度爆炸

图1图2

我们发现如果有很多这样的三角形程序就会爆炸。所以我们考虑优化。我们发现对于所有出度为1的点我们完全不必枚举它的出边，所以根据这些就写出了代码1，然后它居然能水过，但是对于像图2这种数据它就只有被卡掉的份啦，所以我们考虑继续优化。我们考虑将所有只有一条路径的点从图中真实的去掉，也就是说我们将所有点的相连关系提前记录下来，以便后续输出答案，而在搜索过程中只考虑出度不为1的点，听起来可能有些麻烦，参照代码2看看写写就好理解了。

代码1

#include<bits/stdc++.h>
using namespace std;
vector<int>v[];
int n,m,x[],y[],vis[],id[],od[],nxt[];
int ans[];
inline void bad(){
       puts("There is no route, Karl!");
       exit();
}
inline void dfs(int x){
       if(vis[x])return;
       vis[x]=;
       for(int i=;i<v[x].size();i++)dfs(v[x][i]);
}
inline void pr(){
      for(int i=;i<=n;i++)
        printf(" %d",ans[i]);
      puts("");
      exit();
}
inline void search(int x,int cnt){
      if(cnt==n){
          ans[cnt]=x;
          if(x==)pr();
          return;
      }
      if(vis[x])return;
      ans[cnt]=x;
      vis[x]=;
      if(nxt[x]){
        search(nxt[x],cnt+);
      }else {
          for(int i=;i<v[x].size();i++){
            search(v[x][i],cnt+);
          }
      }
      vis[x]=;
      return;
}
int main(){
      freopen("king.in","r",stdin);
      freopen("king.out","w",stdout);
      scanf("%d%d",&n,&m);
      if(m<n)bad();
      for(int i=;i<=m;++i){
          scanf("%d%d",&x[i],&y[i]);
          v[x[i]].push_back(y[i]);
      }
      memset(vis,,sizeof(vis));
      dfs();
      for(int i=;i<=n;i++)
        if(!vis[i])
          bad();
      for(int i=;i<=n;i++)v[i].clear();
      for(int i=;i<=m;i++){
          v[y[i]].push_back(x[i]);
      }
      memset(vis,,sizeof(vis));
      dfs();
      for(int i=;i<=n;i++)
        if(!vis[i])
          bad();
      for(int i=;i<=n;i++)v[i].clear();
      for(int i=;i<=m;i++){
          v[x[i]].push_back(y[i]);
        if(od[x[i]]==)od[x[i]]=y[i];
          else if(od[x[i]]!=y[i])od[x[i]]=-;
        if(id[y[i]]==)id[y[i]]=x[i];
          else if(id[y[i]]!=x[i])id[y[i]]=-;
      }
      for(int i=;i<=n;i++)
        if(od[i]==||id[i]==)
          bad();
      for(int i=;i<=n;i++)
        if(od[i]>&&id[i]>){
          nxt[i]=od[i];
          nxt[id[i]]=i;
        }
      memset(vis,,sizeof(vis));
      search(,);
      bad();
      return ;
}

代码2

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<cctype>
#include<cmath>
#include<cstdlib>
#include<queue>
#include<ctime>
#include<vector>
#include<set>
#include<map>
#include<stack>
using namespace std;
#define sp cout<<"---------------------------------------------------"<<endl
vector<int>v[];
vector<int>cpx;
int nxt[],vis[],arr[],n,m;
inline void bad(){
       puts("There is no route, Karl!");
       exit();
}
inline void check(){
      int i,x;
      memset(arr,,sizeof(arr));
      for(i=,x=;i<=n;i++,x=nxt[x]){
          if(arr[x])return;
          arr[x]=;
      }
      for(i=,x=;i<=n;i++,x=nxt[x]){
          printf("%d ",x);
      }
      cout<<<<endl;
      exit();
}
inline void work(int wh){
      if(wh==cpx.size()){
        check();
        return;
      }
      int x=cpx[wh];
      for(int i=;i<v[x].size();i++)
        if(!vis[v[x][i]]){
          nxt[x]=v[x][i];
          vis[v[x][i]]=;
          work(wh+);
          vis[v[x][i]]=;
        }
      return;
}
int main(){
      freopen("king.in","r",stdin);
      freopen("king.out","w",stdout);
      int i,j,k;
      scanf("%d%d",&n,&m);
      for(i=;i<=m;i++){
          int x,y;
          scanf("%d%d",&x,&y);
          v[x].push_back(y);
      }
      for(i=;i<=n;i++)
        if(v[i].empty())bad();
      for(i=;i<=n;i++)
        if(v[i].size()>){
          cpx.push_back(i);
        }else {
          if(vis[v[i][]])bad();
          vis[v[i][]]=;
        }
      for(i=;i<=n;i++)
        nxt[i]=v[i][];
      work();
      bad();
      return ;
}