HDU6165(tarjan算法 强连通缩点,拓扑排序)

这里说明一下，tarjan是一个人，他发明的算法都以他命名，所以，tarjan后面要加后缀，不然不知道是哪个算法。

一道比较标准的模板题，这里有篇很好的讲强 tarjan求强连通分量 的博客点击打开链接

我这里直接把新图构建出来了，然后拓扑排序即可，只要那一层出现两个入度为0的点即不能到达。

#include<iostream>
#include<vector>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
int T,n,m;
int sccn,dfni;
int scc[1005],dfn[1005],low[1005];
vector<int> mp[1005],newmp[1005];
stack<int> stac;
int instac[1005];

void init() {
	sccn=0;
	dfni=1;
	memset(low,0,sizeof(low));
	memset(scc,0,sizeof(scc));
	memset(instac,0,sizeof(instac));
	memset(dfn,0,sizeof(dfn));
	for(int i=1;i<=n;i++) {
		mp[i].clear();
	} 
}

void tarjan(int p) {
	dfn[p]=low[p]=dfni;
	dfni++;
	instac[p]=1;
	stac.push(p);
	for(int i=0;i<mp[p].size();i++) {
		int v=mp[p][i];
		if(dfn[v]==0) {
			tarjan(v);
			if(low[v]<low[p])
			    low[p]=low[v];
		}else if(instac[v]&&dfn[v]<low[p]) {
			low[p]=dfn[v];
		}
	}
	if(dfn[p]==low[p]) {
		sccn++;
		int te=0;
		while(te!=p) {
			te=stac.top();
			scc[te]=sccn;
			instac[te]=0;
			stac.pop();
		}
	}
}

bool visit[1005][1005];//新图以scc为单位 
void newmap() {
	memset(visit,0,sizeof(visit));
	for(int i=1;i<=sccn;i++) {
		newmp[i].clear();
	}
	for(int i=1;i<=n;i++) {
		for(int j=0;j<mp[i].size();j++) {
			if(scc[i]!=scc[mp[i][j]]&&visit[scc[i]][scc[mp[i][j]]]==0) {
				newmp[scc[i]].push_back(scc[mp[i][j]]);
				visit[scc[i]][scc[mp[i][j]]]=1;
			} 
		}
	}
}

int degree[1005],pi,pnum;
int toposort() {
	memset(degree,0,sizeof(degree));
	for(int i=1;i<=sccn;i++) {
		for(int j=0;j<newmp[i].size();j++) {
		    degree[newmp[i][j]]++;
		}
	}
	while(1) {
		pi=0,pnum=0;
		for(int i=1;i<=sccn;i++) {
			if(degree[i]==0) {
				pnum++;
				pi=i;
			}
	    }
	    if(pnum==0) {
	    	return 1;
		}else if(pnum>=2) {
			return 0;
		}else if(pnum==1) {
			degree[pi]=-1;
			for(int i=0;i<newmp[pi].size();i++) {
				degree[newmp[pi][i]]--;
			}	
		}
	}
	return 1;	
}

int main() {
	//freopen("fuck.txt","r",stdin);
	scanf("%d",&T);
	while(T--) {
		scanf("%d%d",&n,&m);
		init();
		for(int i=1;i<=m;i++) {
			int u,v;
			scanf("%d%d",&u,&v);
			mp[u].push_back(v);
		}
		for(int i=1;i<=n;i++) {
			if(scc[i]==0) {
				tarjan(i);
			}
		}
		newmap();
		if(toposort()) {
			printf("I love you my love and our love save us!\n");
		}else {
			printf("Light my fire!\n");
		}
	}
	return 0;
}