HDU 3231 Box Relations

题目大意： 给定一些正方体的关系，要求一组符合这些关系的正方体坐标，如果不存在符合条件的正方体坐标，IMPOSSIBLE。(Special Judge)

实力还是太弱了，完全不会……

#include <cstdio>
#include <stdlib>
#include <iostream>
#define MAXN 2010
#define MAXR 500000
#define MAX 999999

typedef struct edges{
    int v,w,next;
}edge;

int N, R;
edge edge_X[MAXR], edge_Y[MAXR], edge_Z[MAXR];
int s_X[MAXN], s_Y[MAXN], s_Z[MAXN];
int index_X[MAXN], index_Y[MAXN], index_Z[MAXN];
int degree_X[MAXN], degree_Y[MAXN], degree_Z[MAXN];

void Add(struct edges e[], int start, int end, int index[], int degree[], int &tot){
    e[tot].v = end;
    e[tot].next = index[start];
    index[start] = tot++;
    degree[end]++;
}

bool init(){
    int i, totx, toty, totz;
    scanf("%d%d", &N, &R);
    if (N == 0 && R == 0) return false;
    memset(index_X, 255, sizeof(*index_X)*MAXN); //-1 indicates end of link
    memset(index_Y, 255, sizeof(*index_Y)*MAXN);
    memset(index_Z, 255, sizeof(*index_Z)*MAXN);
    memset(degree_X, 0, sizeof(*degree_X)*MAXN);
    memset(degree_Y, 0, sizeof(*degree_Y)*MAXN);
    memset(degree_Z, 0, sizeof(*degree_Z)*MAXN);
    totx = toty = totz = 0;

    for (i=1; i<=N; i++){
        // node 0是超级汇点，最大的点。值为0
        Add(edge_X, 0, i, index_X, degree_X, totx);         //x1 of n-th Cube
        Add(edge_Y, 0, i, index_Y, degree_Y, toty);
        Add(edge_Z, 0, i, index_Z, degree_Z, totz);
        Add(edge_X, i, i+N, index_X, degree_X, totx);
        Add(edge_Y, i, i+N, index_Y, degree_Y, toty);
        Add(edge_Z, i, i+N, index_Z, degree_Z, totz);
    }
    for (i=1; i<=R; i++){
        char t[10];
        int  A, B;
        scanf("%s%d%d", t, &A, &B);
        if (t[0] == 'I'){
            Add(edge_X, A, B+N, index_X, degree_X, totx);        //x1(A) < x1(B)
            Add(edge_X, B, A+N, index_X, degree_X, totx);        //x2(A) > x1(B) or x1(B) < x2(A)
            Add(edge_Y, A, B+N, index_Y, degree_Y, toty);
            Add(edge_Y, B, A+N, index_Y, degree_Y, toty);
            Add(edge_Z, A, B+N, index_Z, degree_Z, totz);
            Add(edge_Z, B, A+N, index_Z, degree_Z, totz);
        }
        if (t[0] == 'X')
            Add(edge_X, A+N, B, index_X, degree_X, totx);    //x2(A) < x1(B)
        if (t[0] == 'Y')
            Add(edge_Y, A+N, B, index_Y, degree_Y, toty);    //y2(A) < y1(B)
        if (t[0] == 'Z')
            Add(edge_Z, A+N, B, index_Z, degree_Z, totz);    //z2(A) < z1(B)
    }
    return true;
}

bool NonPreFirstTopSort(edge e[], int index[], int degree[], int n, int s[]){
    int i, count(0), head(0), tail(1);
    int Q[MAXN];

    Q[0] = 0;

    while (head != tail){
        s[Q[head]] = count++;
        i = index[Q[head]];
        while (i != -1){
            if (--degree[e[i].v] == 0) Q[tail++] = e[i].v;
            i = e[i].next;
        }
        ++head;
    }

    if (count < n) return false;
    else return true;
}

int main(){
    int Cases = 0, i;
    while (init() && ++Cases){
        bool flag;
        flag = NonPreFirstTopSort(edge_X, index_X, degree_X, N*2 + 1, s_X);
        if (flag) flag = NonPreFirstTopSort(edge_Y, index_Y, degree_Y, N*2 + 1, s_Y);
        if (flag) flag = NonPreFirstTopSort(edge_Z, index_Z, degree_Z, N*2 + 1, s_Z);
        if (flag){
            printf("Case %d: POSSIBLE\n", Cases);
            for (i=1; i<=N; i++)
                printf("%d %d %d %d %d %d\n", s_X[i], s_Y[i], s_Z[i], s_X[i+N],  s_Y[i+N], s_Z[i+N]);
            printf("\n");

        }
        else
            printf("Case %d: IMPOSSIBLE\n\n", Cases);
    }
    //system("pause");
    return 0;
}