[LA 6393]Self-Assembly(拓扑排序)

题目链接

http://acm.hust.edu.cn/vjudge/problem/viewProblem.action?id=46836

题目大意

以下中文题面转自清澄OJ
自动化学制造（Automatic Chemical Manufacturing，简称ACM）正在对一个叫自组装（self-assembly）的过程进行实验。在这个过程中，有着天然相互吸引力的分子被混合在溶液中，任由它们聚集组合成更大的结构。但是有一个问题随之出现：有时候，分子们会把自身组合成一个无限大的结构体，以至于把容器撑爆。
　　你需要写一个程序来判断一个给定的分子集合是否可能组合成一个无限大的结构体。为了使问题简化，你可以作以下两个假设：
　　1. 问题被限制在二维平面上。
　　2. 分子集合中的每个分子都被表示成一个正方形。其中正方形的四条边分别代表分子间相连接的四个表面。
　　你将从数据中得到每种分子的描述。每种分子有四个连接标识来分别表示每条边能与另外分子的哪种边相连。连接标识有两种：
　　● 一个大写字母(A,…,Z)加上一个 ‘＋’ 号或一个 ‘－’ 号。两条边能并在一起当且仅当两者的字母相等且符号相反。比方说，‘A＋’ 与 ‘A－’ 兼容，但与 ‘A＋’ 或 ‘B－’ 不兼容
　　● 两个零 ‘00’。这条边将不和任意一条边兼容（包括‘00’）。
　　假设每种分子都有无限个，并且每个分子都可以旋转和翻转。当分子将自身组成一个结构体时，相互贴合的边必须能够相互兼容，当然，无论边的连接标识是什么，它都可以不与另外边贴合。

图1 是一个由三种分子组成的一个有限的结构体（它们也有可能组成另外的有限结构体）。

思路

首先很容易将此题与图论问题联想起来。我们可以尝试将每个方格看成点，将方格与方格可以挨着的邻接边用无向边相连。若这个图存在一个环，则表明该结构体是无限的。

然而这个做法是错误的，有可能会出现一个情况：找出的环的链接并不是合法的，可能出现一个方块的一条边链接了多个方块的情况。问题就在于这个模型没有细化方块与方块之间是通过哪两条边链接的，只标明了哪些方块之间链接了。因此我们可以进一步改进模型：将每种非”0 0”的方块边均看成点，方块与方块的链接方式看成有向边，这样由于找出的环上不可能有一条边出现多次，因此不会出现一个方块的一条边链接了多个方块的情况。具体连边方式就是对于每个方块上的每个可链接的边，如“A+”，从图中的点“A+”向其相反点“A-”连有向边，然后再做拓扑排序

代码

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <algorithm>
#include <queue>

#define MAXV 51000
#define MAXE 5100000

using namespace std;

int n;

struct edge
{
    int u,v,next;
}edges[MAXE];

int head[MAXV],nCount=0;
int inDegree[MAXV];

void AddEdge(int U,int V)
{
    inDegree[V]++;
    edges[++nCount].u=U;
    edges[nCount].v=V;
    edges[nCount].next=head[U];
    head[U]=nCount;
}

int anti(int x) //相反点
{
    if(x<=26) x+=26;
    else x-=26;
    return x;
}

queue<int>q;

bool toposort() //返回true表示有环
{
    for(int i=1;i<=52;i++)
        if(inDegree[i]==0)
            q.push(i);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int p=head[u];p!=-1;p=edges[p].next)
        {
            int v=edges[p].v;
            inDegree[v]--;
            if(!inDegree[v]) q.push(v);
        }
    }
    for(int i=1;i<=52;i++) if(inDegree[i]) return true;
    return false;
}

char s[20];

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        nCount=0;
        memset(head,-1,sizeof(head));
        memset(inDegree,0,sizeof(inDegree));
        for(int i=1;i<=n;i++)
        {
            int stack[10],top=0;
            scanf("%s",s+1);
            if(s[1]!='0')
            {
                stack[++top]=s[1]-'A'+1;
                if(s[2]=='-') stack[top]+=26;
            }
            if(s[3]!='0')
            {
                stack[++top]=s[3]-'A'+1;
                if(s[4]=='-') stack[top]+=26;
            }
            if(s[5]!='0')
            {
                stack[++top]=s[5]-'A'+1;
                if(s[6]=='-') stack[top]+=26;
            }
            if(s[7]!='0')
            {
                stack[++top]=s[7]-'A'+1;
                if(s[8]=='-') stack[top]+=26;
            }
            for(int k=1;k<=top;k++)
                for(int j=1;j<=top;j++)
                    if(k!=j)
                        AddEdge(stack[k],anti(stack[j]));
        }
        if(toposort()) printf("unbounded\n");
        else printf("bounded\n");
    }
    return 0;
}