POJ 2438 哈密顿回路模板题

题意：

n个人，m行描述n个人的关系，n个人两两都有边相连，m行表示这两个点没有边相连。输出一条哈密顿回路（每个点必须经过一次）

分析：

哈密顿回路模板：

#include<stdio.h>
#include<string.h>
#include<algorithm>

using namespace std;

const int maxn=405;

int ans[maxn],g[maxn][maxn];
bool vis[maxn];
int n,m;

void init()
{
    for(int i=0;i<=n;i++)
        for(int j=0;j<=n;j++)
            if(i==j)
                g[i][j]=0;
            else
                g[i][j]=1;
    memset(vis,0,sizeof(vis));
    memset(ans,0,sizeof(ans));
}

void revese(int ans[],int s,int t)//将ans数组中s到t的部分倒置
{
    while(s<t){
        swap(ans[s],ans[t]);
        s++;
        t--;
    }
}

void Hamilton()
{
    int s=1,t;//初始化s取1号点
    int ansi=2;
    int i,j,w;
    for(i=1;i<=n;i++)
        if(g[s][i]) break;
    t=i;//取任意连接s的点为t
    vis[s]=vis[i]=true;
    ans[0]=s;ans[1]=t;
    while(true){
        while(true){//从t向外扩展
            for(i=1;i<=n;i++)
                if(g[t][i]&&!vis[i]){
                    ans[ansi++]=i;
                    vis[i]=true;
                    t=i;
                    break;
                }
            if(i>n)break;
        }
        //将当前得到的序列倒置，s和t互换，从t继续扩展，相当于在原来的序列上从s扩展
        w=ansi-1;
        i=0;
        revese(ans,i,w);
        swap(s,t);
        //从新的t向外扩展，相当于在原来的序列上从s向外扩展
        while(true){
            for(i=1;i<=n;i++)
                if(g[t][i]&&!vis[i]){
                    ans[ansi++]=i;
                    vis[i]=true;
                    t=i;
                    break;
                }
            if(i>n)break;
        }
        //如果s和t不相邻，进行调整
        if(!g[s][t]){
            for(i=1;i<ansi-2;i++)
                if(g[ans[i]][t]&&g[s][ans[i+1]])//取序列中一点i，使得ans[i]与t相连接且ans[i+1]与s相连
                    break;
            //将从ans[i+1]到t部分的ans[]倒置
            w=ansi-1;
            i++;
            t=ans[i];
            revese(ans,i,w);
        }
        //如果当前s和t相连
        if(ansi==n)return;//如果当前序列中包含n个元素，算法结束
        //当前序列中的元素个数小于n，寻找点ans[i]，使得ans[i]与ans[]外一点相连
        for(j=1;j<=n;j++){
            if(vis[j])continue;
            for(i=1;i<ansi-2;i++)
                if(g[ans[i]][j])
                    break;
            if(g[ans[i]][j])
                break;
        }
        s=ans[i-1];
        t=j;
        revese(ans,0,i-1);//将ans[]中s到ans[i-1]部分的ans[]倒置
        revese(ans,i,ansi-1);//将ans[]中ans[i]到t的部分倒置
        ans[ansi++]=j;//将点j加入到ans[]的尾部
        vis[j]=true;
    }
}

int main()
{
    while(scanf("%d%d",&n,&m)&&(n||m))
    {
        n<<=1;
        init();
        int a,b;
        while(m--){
            scanf("%d%d",&a,&b);
            g[a][b]=g[b][a]=0;//巧妙的建立反图
        }
        Hamilton();
        printf("%d",ans[0]);
        for(int i=1;i<n;i++)
            printf(" %d",ans[i]);
        printf("\n");
    }
    return 0;
}