poj 1236 Network of Schools(Tarjan)

依旧找强连通分量的题来做，YM的是，这题开始没看懂 = =。。。看了别人的思路才看懂了，YM。

 

Task A ： 在已经存在的图中，最少需要给多少个点，才能把这个图给遍历完。就是求入度为0的点。

 

Task B ： 最少需要连接多少条边才可以使图为一个强连通图。因为两个强连通合并必然是出度为0的连接入度为0的点，所以要解决掉入度为0，和出度为0的点，所以答案是这两个的最大值（我说的点指缩点）。

 

#include <stdio.h> #include <stdlib.h> #include <iostream> #include <string.h> #include <stack> #define MAX 110 using namespace std; typedef struct MAP{ int to; MAP *next; }MAP; MAP *map[MAX],node[MAX*MAX]; int n,cou,ind,scc,sccf[MAX],low[MAX],numed[MAX]; int ins[MAX]; stack <int> s; void init() { int i; for(i=0; i<=n; i++) numed[i] = i; cou = scc = 0; ind = 1; memset(ins,0,sizeof(ins)); memset(low,0,sizeof(low)); memset(node,'/0',sizeof(node)); memset(map,'/0',sizeof(map)); memset(numed,0,sizeof(numed)); memset(sccf,0,sizeof(sccf)); } void input() { int i,x; scanf("%d",&n); for(i=1; i<=n; i++) { while( scanf("%d",&x) && x ) { node[cou].to = x; node[cou].next = map[i]; map[i] = &node[cou++]; } } } void Tarjan(int u) { int v; low[u] = numed[u] = ind++; s.push(u); ins[u] = 1; MAP *head = map[u]; while( head != NULL ) { v = head->to; if( !numed[v] ) { Tarjan(v); if( low[v] < low[u] ) low[u] = low[v]; } else if( ins[v] && numed[v] < low[u] ) low[u] = numed[v]; head = head->next; } if( numed[u] == low[u] ) { scc++; do { v = s.top(); s.pop(); ins[v] = 0; sccf[v] = scc; }while( u != v ); } } int Atask() // 找入度为0的点。 { int indegree[MAX],ans = 0,i; MAP *head; memset(indegree,0,sizeof(indegree)); for(i=1; i<=n; i++) { head = map[i]; while( head != NULL ) { if( sccf[i] != sccf[head->to] ) indegree[sccf[head->to]]++; head = head->next; } } for(i=1; i<=scc; i++) if( indegree[i] == 0 ) ans++; return ans; } int Btask()//找出度为0的点。 { int outd[MAX],ans = 0,i; MAP *head; memset(outd,0,sizeof(outd)); for(i=1; i<=n; i++) { head = map[i]; while( head != NULL ) { if( sccf[i] != sccf[head->to] ) outd[sccf[i]]++; head = head->next; } } for(i=1; i<=scc; i++) if( outd[i] == 0 ) ans++; return ans; } int main() { int i,ansa,ansb; init(); input(); for(i=1; i<=n; i++) if( !numed[i] ) Tarjan(i); if( scc == 1 ) // 特殊处理。 printf("1/n0/n"); else { ansa = Atask(); ansb = Btask(); printf("%d/n%d/n",ansa,max(ansa,ansb)); } return 0; }