hust1017

/*
分析:
    LZ在经过两个星期的挣扎后,终于鼓起勇气开始看DancingLinks,于昨
日上午两个小时看懂思想,又用俩个小时加入自己的理解,敲出自己模板,
然后用了整整一个下午来de自己的那一处小小的bug,结果愣是整整一下午到
太阳都消失了都牟牛揪揣那与哦bug,头疼类跟要爆炸一样。于是于晚饭前打
开游戏玩两局,然后!于今天上午起床吃过饭再次狂奔而来debug!这次经过
一个小时的努力,终于de出了它,为以后写有双表头的dancinglinks打下了
雄厚的、不可磨灭的、坚实的、顽强的、……的基础。也终于可以在个人blog
中加入dancinglinks这一光听名字就觉得很唬B的一个分类,哦也~~~  -、-I
    有谁能比我还水么-、-III,专业课不开数据结构,十字链表还是学

dancinglinks学到的人,桑不起呀,囧~~~

    DancingLinks果题,俺就不啰嗦了,囧~~~

                                                             2012-10-23
*/








//题目链接http://acm.hust.edu.cn:8080/judge/problem/viewProblem.action?id=10702








#include"stdio.h"
#include"string.h"
#include"stdlib.h"
#include"math.h"
#define N 2222
#define V 1111111
int n,m;
int S[N];
int L[V],R[V],U[V],D[V],C[V],H[V];
int c_T[N],h_T[N];
int head,base;
int Ans,ans[N],flag;
void Link(int r,int c)
{
    int t;
    t=(r-1)*m+c+base;
    //行列头的
    C[t]=c;H[t]=m+r;/*是m+r*/S[c]++;
    //t的
    U[t]=c_T[c];D[t]=c;
    L[t]=h_T[m+r];/*是m+r*/R[t]=m+r;
    //上一个t的
    U[c]=t;
    D[c_T[c]]=t;
    L[m+r]=t;
    R[h_T[m+r]]=t;
    //next
    c_T[c]=t;
    h_T[m+r]=t;
}
void remove(int c)
{
    int i,j;
    R[L[c]]=R[c];
    L[R[c]]=L[c];
    for(i=D[c];i!=c;i=D[i])
    {
        for(j=R[i];j!=i;j=R[j])
        {
            D[U[j]]=D[j];
            U[D[j]]=U[j];
            S[C[j]]--;
        }
    }
}
void resume(int c)
{
    int i,j;
    for(i=U[c];i!=c;i=U[i])
    {
        for(j=L[i];j!=i;j=L[j])
        {
            U[D[j]]=j;
            D[U[j]]=j;
            S[C[j]]++;
        }
    }
    R[L[c]]=c;
    L[R[c]]=c;
}
void Dance(int k)
{
    int i,l,j;
    int c,t=111111111;
    if(flag)    return ;
    if(R[head]==head)
    {
        flag=1;
        return ;
    }
    for(i=R[head];i!=head;i=R[i])
    {
        if(S[i]<t)
        {
            c=i;
            t=S[i];
        }
    }
    remove(c);
    for(i=D[c];i!=c;i=D[i])
    {
        ans[k]=H[i];Ans=k+1;
        for(j=R[i];j!=i;j=R[j])
        {
            if(m<j && j<=m+n)   continue;
            remove(C[j]);
        }
        Dance(k+1);
        if(flag)    return ;
        for(j=L[i];j!=i;j=L[j])
        {
            if(m<j && j<=m+n)   continue;
            resume(C[j]);
        }
    }
    resume(c);
}
int main()
{
    int i,l;
    int c,t,temp;
    while(scanf("%d%d",&n,&m)!=-1)
    {
        memset(S,0,sizeof(S));
        
        head=0;
        L[head]=m;R[head]=1;
        U[head]=m+n;D[head]=m+1;
        for(i=1;i<=m;i++)
        {
            L[i]=i-1;
            R[i-1]=i;
            C[i]=U[i]=D[i]=i;
            c_T[i]=i;
        }
        R[m]=0;
        for(i=m+2;i<=m+n;i++)
        {
            U[i]=i-1;
            D[i-1]=i;
            H[i]=L[i]=i;R[i]=i;
            h_T[i]=i;
        }
        U[m+1]=0;h_T[m+1]=m+1;
        D[m+n]=0;
        base=n+m;

        for(i=1;i<=n;i++)
        {
            scanf("%d",&c);
            while(c--)
            {
                scanf("%d",&t);
                Link(i,t);
            }
        }

        flag=0;
        Ans=0;
        Dance(0);

        if(!flag)   printf("NO\n");
        else
        {
            printf("%d",Ans);
            for(i=0;i<Ans;i++)  printf(" %d",ans[i]-m);
            printf("\n");
        }
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值