HDU3231Box Relations【拓扑排序】

题意：有很多立方体，给出他们的位置关系，I代表i，j有相交部分，X代表 Xi 完全小于 Xj（Y、Z同理），输出满足情况的一种。

思路：相当于在三条坐标轴上分别拓扑排序，把一个立方体在一条坐标轴方向的两个面看成两个点，相交的时候没说i在左，还有j在左，找到这两种情况的共性，x1 < x2' && x2 < x1',对于每个立方体的两个面x1 < x1'

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<stdlib.h>
#include<math.h>
#include<vector>
#include<list>
#include<map>
#include<stack>
#include<queue>
#include<algorithm>
#include<numeric>
#include<functional>
using namespace std;
typedef long long ll;
const int maxn = 2005;

int counts[maxn][5],n;
vector<int> g[maxn][5];
int d[maxn][5];
queue<int> q;

void init(int n)
{
    memset(counts,0,sizeof(counts));
    memset(d,0,sizeof d);
    for(int i = 0; i <= 2*n;i++)
	    for(int j = 0; j <= 3; j++)
	        g[i][j].clear();
}

void add(int u,int v,int w)
{
    g[u][w].push_back(v);
    counts[v][w]++;
}

bool topsort(int w)
{
	int sum = 2*n;  //n 节点总数
	bool flag=false;
	while(!q.empty())
        q.pop();
	for(int i = 1; i <= 2*n; i++)  //入度为0入栈
	{
		if(counts[i][w]==0)
			q.push(i);
	}
	while(!q.empty())
	{	
		sum--;
		int j=q.front();
		q.pop();
		for(int k=0;k<g[j][w].size();k++)
		{
			d[ g[j][w][k] ][w] = max(  d[ g[j][w][k] ][w] , d[j][w] + 1  );
			if(--counts[g[j][w][k]][w]==0)
				q.push(g[j][w][k]);
		}
	}
	if(sum>0)
		return false;  //表明存在回路
	return true;        //OK
}


int main(void)
{
	int m,kase = 1,a,b;
	char s[50];
	while(scanf("%d%d",&n,&m)!=EOF && n+m)
	{
		init(n);
		for(int i = 1; i <= n; i++)
		{
			for(int j = 0; j < 3; j++)
				add(2*i-1,2*i,j);
		}
		while(m--)
		{
			scanf("%s%d%d",s,&a,&b);
			if(s[0] == 'I')
			{
				for(int i = 0; i < 3; i++) // 三维
				{
					add(2*b-1,2*a,i);
					add(2*a-1,2*b,i);
				}
			}
			else if(s[0] == 'X')
				add(2*a,2*b-1,0);
			else if(s[0] == 'Y')
				add(2*a,2*b-1,1);
			else if(s[0] == 'Z')
				add(2*a,2*b-1,2);	
		}
		int flag = 1;
		for(int i = 0; i < 3; i++)
		{
			if(!topsort(i))
			{
				flag = 0;
				break;
			}
		}
		if(flag)
		{
			printf("Case %d: POSSIBLE\n",kase++);
			for(int i = 1; i <= n; i++)
			{
				printf("%d %d %d ",d[2*i-1][0],d[2*i-1][1],d[2*i-1][2]);
				printf("%d %d %d\n",d[2*i][0],d[2*i][1],d[2*i][2]);
			}
		}
		else
			printf("Case %d: IMPOSSIBLE\n",kase++);
		putchar('\n');
	}
	return 0;
}