HDU 3231 Box Relations

http://acm.hdu.edu.cn/showproblem.php?pid=3231

很巧妙地将问题转化为拓扑排序

#include <cstdio>
#include <cstring>
using namespace std;

const int N=2010;
bool g[3][N][N];
int ind[3][N],n;
int topo[3][N];
int x[N],y[N],z[N];
bool toposort(int k)
{
    memset(ind[k],0,sizeof(ind[k]));
    for(int v=1;v<=n;v++)
        for(int u=1;u<=n;u++)
            if(g[k][u][v]) ind[k][v]++;
    for(int i=0;i<n;i++)
    {
        int u;
        for(u=1;u<=n;u++)
            if(ind[k][u]==0) break;
        if(u>n) return false;
        topo[k][i]=u; ind[k][u]--;
        for(int v=1;v<=n;v++)
            if(g[k][u][v]) ind[k][v]--;
    }
    return true;
}
int main()
{
    int m,C=0;
    while(scanf("%d%d",&n,&m),n||m)
    {
        memset(g,0,sizeof(g));
        for(int i=0;i<3;i++)
            for(int u=1;u<=n;u++) g[i][u][u+n]=1;
        while(m--)
        {
            char op[2];
            int u,v;
            scanf("%s%d%d",op,&u,&v);
            if(op[0]=='I')
            {
                for(int i=0;i<3;i++)
                    g[i][v][u+n]=g[i][u][v+n]=1;
            }
            else
            {
                int k=op[0]-'X';
                g[k][u+n][v]=1;
            }
        }
        n*=2;
        bool flag=true;
        for(int i=0;flag&&i<3;i++)
            if(!toposort(i)) flag=false;
        n/=2;
        if(flag)
        {
            printf("Case %d: POSSIBLE\n",++C);
            for(int i=0;i<2*n;i++)
                x[topo[0][i]]=y[topo[1][i]]=z[topo[2][i]]=i;
            for(int i=1;i<=n;i++)
                printf("%d %d %d %d %d %d\n",x[i],y[i],z[i],x[i+n],y[i+n],z[i+n]);
        }
        else printf("Case %d: IMPOSSIBLE\n",++C);
        printf("\n");
    }
    return 0;
}