POJ 1236 Network of Schools

题目大意:N（2<N<100）各学校之间有单向的网络，每个学校得到一套软件后，可以通过单向网络向周边的学校传输，问题1：初始至少需要向多少个学校发放软件，使得网络内所有的学校最终都能得到软件。2，至少需要添加几条传输线路（边），使任意向一个学校发放软件后，经过若干次传送，网络内所有的学校最终都能得到软件。

 

分析要点：

1.如果网络中存在强连通分量G，从G强连通分量中任意一个点出发，都能遍历G这个连通分量，而任意通过G以外的点进入G，都可以通过G中的任意一个出度大于1的点离开G。这样，可以将网络中的强连通分量看成一个点，也就是传说中的“缩点”（详见《 缩点》），这样图就得到了简化。寻找强连通分量的方法用kosaraju算法。

2.新生成的图是一个有向无环图，如果将图中的无向环，也就是点从u出发经过多条分支重新汇到一点u，从u到v的环都可以再次收缩成一条从u到v的边。完成了这两部工作，最后的到的图将是一个森林。

3.从森林里任意一个点出发都可以遍历它的子树，因此，为了遍历所有的点，所有点的公共祖先就是问题1所求点，这些点的特征应该是入度为0

4.问题2的要求是添加最少边，使图完全连通，现在转化成将一个森林变化成完全联通，很显然，只需要将一棵树的叶子结点轮流着连接到“相邻的”树形结构的祖先点，具体需要连接的数量视总的祖先点和叶子点的较大值，叶子节点的特征是出度为0。必须注意的是有一种特殊情况，如果采取记录度的方式判断，需要注意当最终生成图中只有一个点（也就是说原图中只存在一个强联通分量），不需要将出度为零的点与入度为零的点相连（那样得到答案2为1），正确答案是0

 

 #include <iostream>
#include <fstream>
#include <memory.h>
using namespace std;

#define init(l) memset(l,0,sizeof(l))
#define N 101
struct edge
{
 int dest;
 edge *next;
};
edge *G1[N],*G2[N],*G[N];
int used[N],path[N],part[N],mark[N][N];
int m,n;
void readdata();
void dfs1(int s);
void dfs2(int s);
void kosaraju();
void addedge(edge *G[],int i, int j);
void shrink();
void work();

int main()
{
 readdata();
 kosaraju();
 shrink();
 work();
 return 0;
}

void work()
{
 int i,j,A,B,ind[N],outd[N];
 init(ind);
 init(outd);
 edge *L;
 for (m=part[0],i=1;i<=m;i++)
  for (L=G[i];L!=NULL;L=L->next)
  {
  outd[i]++;
  ind[L->dest]++;
  }
 for (A=B=0,i=1;i<=m;i++)
 {
  if (!ind[i])
  A++;
  if (!outd[i])
  B++;
 }
 B=A>B?A:B;
 if (m==1)
  B=0;
 cout<<A<<endl<<B<<endl;
}

void addedge(edge *T[],int i, int j)
{
 edge *L;
 L=new edge;
 L->dest=j;
 L->next=T[i];
 T[i]=L;
}

void shrink()
{
 int i,j,k;
 edge *L;
 init(mark);
 for (k=1;k<=n;k++)
  for (L=G1[k],i=part[k];L!=NULL;L=L->next)
  {
  j=part[L->dest];
  if (i!=j&&!mark[i][j])
  {
  mark[i][j]=1;
  addedge(G,i,j);
  }
  }
}

void dfs1(int s)
{
 if(!used[s])
 {
  used[s]=1;
  for (edge *l=G1[s];l!=NULL;l=l->next)
  dfs1(l->dest);
  path[0]++;
  path[path[0]]=s;
 }
}

void dfs2(int s)
{
 if (!used[s])
 {
  used[s]=1;
  for (edge *l=G2[s];l!=NULL;l=l->next)
  dfs2(l->dest);
  part[s]=part[0];
 }
}

void kosaraju()
{
 int i;
 init(used);
 path[0]=part[0]=0;
 for (i=1;i<=n;i++)
  dfs1(i);
 init(used);
 for (i=n;i>=1;i--)
  if (!used[path[i]])
  {
  part[0]++;
  dfs2(path[i]);
  }
}

void readdata()
{
// ifstream cin("data.txt");
 int i,j;
 cin>>n;
 for (i=1;i<=n;i++)
  G1[i]=G2[i]=G[i]=NULL;
 for (i=1;i<=n;i++)
  while (cin>>j&&j)
  {
  addedge(G1,i,j);
  addedge(G2,j,i);
  }
}