//有向圖強連通分支,dfs鄰接陣形式,O(n^2)
//返回分支數,id返回1..分支數的值
//傳入圖的大小n和鄰接陣mat,不相鄰點邊權0
#define MAXN 100

void search(int n,int mat[][MAXN],int* dfn,int* low,int now,int& cnt,int& tag,int* id,int* st,int& sp){
	int i,j;
	dfn[st[sp++]=now]=low[now]=++cnt;
	for (i=0;i<n;i++)
		if (mat[now][i]){
			if (!dfn[i]){
				ssearch(n,mat,dfn,low,i,cnt,tag,id,st,sp);
				if (low[i]<low[now])
					low[now]=low[i];
			}
			else if (dfn[i]<dfn[now]){
				for (j=0;j<sp&&st[j]!=i;j++);
				if (j<cnt&&dfn[i]<low[now])
					low[now]=dfn[i];
			}
		}
	if (low[now]==dfn[now])
		for (tag++;st[sp]!=now;id[st[--sp]]=tag);
}

int find_components(int n,int mat[][MAXN],int* id){
	int ret=0,i,cnt,sp,st[MAXN],dfn[MAXN],low[MAXN];
	for (i=0;i<n;dfn[i++]=0);
	for (sp=cnt=i=0;i<n;i++)
		if (!dfn[i])
			search(n,mat,dfn,low,i,cnt,ret,id,st,sp);
	return ret;
}