POJ 2438 Children's Dining （求哈密顿回路模板）

题目大意：

n个孩子要求围成一圈，有m对敌对关系，有敌对关系的孩子不能相邻，输出一种方案。

解题思路：

有神犇说是特殊图所以可以用哈密顿回路做，不懂，姑且当作写个模板。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
#include<set>
#include<ctime>
#include<bitset>
#define LL long long
#define db double
#define EPS 1e-15
#define inf 1e10
#define pa pair<int,int>

using namespace std;

const int maxn=420;
int n,m;
bool mp[maxn][maxn];
int S,T,cnt,ans[maxn];
bool vis[maxn];

void reverse(int l,int r){
    while (l<r){
        swap(ans[l],ans[r]);
        l++,r--;
    }
}
void expend(){
    for (;;){
        bool flag=0;
        for (int i=1;i<=n;i++){
            if (!vis[i] && mp[T][i]){
                ans[cnt++]=i;
                T=i;
                vis[i]=1;
                flag=1;
                break;
            }
        }
        if (!flag) break;
    }
}
void hamiltun(){
    memset(vis,0,sizeof(vis));

    S=1;
    for (T=2;T<=n;T++) if (mp[S][T]) break; 
    //任意找两个相邻的节点S和T
    
    cnt=2;
    ans[0]=S, ans[1]=T;
    vis[S]=1, vis[T]=1;
    while (1) {
        expend();  
        //在它们的基础上扩展出一条尽量长的没有重复节点的路径
        reverse(0,cnt-1);
        swap(S,T);
        expend();
        //在它们的基础上扩展出一条尽量长的没有重复节点的路径

        int mid=0;
        if (!mp[S][T]){
        //若S与T不相邻,可以构造出一个回路使新的S和T相邻

            for (int i=1;i<cnt-2;i++){
            // 设路径S到T上有k+2个节点，依次为S,v1,v2…… vk和T.
            // 可以证明存在节点vi,满足vi与T相邻,且vi+1与S相邻
                if (mp[ans[i]][T] && mp[ans[i+1]][S]){
                    mid=i+1;
                    break;
                }
            }
            reverse(mid,cnt-1);
            //把原路径变成S→vi→T→vi+1→S,即形成了一个回路 

            T=ans[cnt-1];
        }
        if (cnt==n) break;
        //现在我们有了一个没有重复节点的回路.如果它的长度为N,则汉密尔顿回路就找到了  
        //否则,由于整个图是连通的,所以在该回路上,一定存在一点与回路以外的点相邻  
        //那么从该点处把回路断开,就变回了一条路径,再按照步骤1的方法尽量扩展路径
        for (int i=1;i<=n;i++){
            if (!vis[i]){
                int j;
                for (j=1;j<cnt-1;j++)
                    if (mp[ans[j]][i]) break;
                if (mp[ans[j]][i]){
                    T=i; mid=j;
                    break;
                }
            }
        }
        S=ans[mid-1];
        reverse(0,mid-1);
        reverse(mid,cnt-1);
        ans[cnt++]=T;
        vis[T]=1;
    }
}
int main(){
    int u,v;
    while (scanf("%d%d",&n,&m)!=EOF){
        if (n==0 && m==0) break;
        n<<=1;
        memset(mp,1,sizeof(mp));
        for (int i=1;i<=n;i++) mp[i][i]=0;
        for (int i=1;i<=m;i++){
            scanf("%d%d",&u,&v);
            mp[u][v]=0;
            mp[v][u]=0;
        }
        hamiltun();
        printf("%d",ans[0]);
        for (int i=1;i<cnt;i++)
            printf(" %d",ans[i]);
        puts("");
    }   
    return 0;
}