计蒜客拓扑排序题目 回收元件

#include<bits/stdc++.h>
using namespace std;
struct node{
	int x1,y1,x2,y2;
};
int n;
const int maxn=1e4+7;
const int MAX=1e6+7;
node it[maxn];
struct edge{
	int v,next;
}e[MAX<<3];
int cnt=0;
int head[maxn];
void add(int u,int v){
	e[cnt].v=v;
	e[cnt].next=head[u];
	head[u]=cnt++;
}
queue<int> q;
int degree[maxn];
void st(){
	for(int i=1;i<=n;i++) if(degree[i]==0) q.push(i);
	while(!q.empty()){
		int tp=q.front();
		q.pop();
		cout<<tp<<' ';
		for(int i=head[tp];i!=-1;i=e[i].next){
			int v=e[i].v;
			degree[v]--;
			if(degree[v]==0) q.push(v);
		}
	}
} 
int main(){
	cin>>n;
	memset(head,-1,sizeof head);
	for(int i=1;i<=n;i++){
		cin>>it[i].x1>>it[i].y1>>it[i].x2>>it[i].y2;
		if(it[i].x1>it[i].x2) swap(it[i].x1,it[i].x2),swap(it[i].y1,it[i].y2);
	} 
	for(int i=1;i<=n;i++){
		for(int j=1;j<=n;j++){
			if(i!=j&&it[i].x1>=it[j].x1&&it[i].x1<=it[j].x2){
				double y;
				if(it[j].x1==it[j].x2){
					if(it[i].x1==it[i].x2){
						if(max(it[j].y1,it[j].y2)<min(it[i].y1,it[i].y2)) add(j,i),degree[i]++;
						else add(i,j),degree[j]++;
					}
					else{
						if(it[i].y1<min(it[j].y1,it[j].y2)) add(i,j),degree[j]++;
						else add(j,i),degree[i]++;
					}
				}
				else{
					y=1.0*(it[i].x1-it[j].x1)*(it[j].y2-it[j].y1)*1.0/(it[j].x2-it[j].x1)+it[j].y1;
					if(y>it[i].y1) add(i,j),degree[j]++;
					else add(j,i),degree[i]++;
				}
			}
		}
	}
	st();
}