POJ2553解题报告.pptVIP

POJ2553解题报告 ——唐波（CS3） 类型：有向图的强连通分支 题目大意：求一个图的“底”，就是一个图中所有的那些如果从一点可以到达另外一个点，那么另外的这个点也可以到达这个点。 bottom(G)={v∈V|∀w∈V:(v→w)⇒(w→v)}. （就是可以相互到达，但不能是其中一定可以到达另外一个，而另外一个不能到达这个点）。 解题思路：很简单，就是先求出强连通分支图Gscc，然后在Gscc图中找到所有的出度为0的强连通分支。这里面的点就是要求的点。 代码 #include stdio.h#include string.h#include stdlib.h#define NV 5001#define NEstruct{int v,next;}e[NE*2];int v1[NV],v2[NV],ne,n,m;  int v1[NV],v2[NV],ne,n,m; void init(){memset(v1,-1,4*n);memset(v2,-1,4*n);ne = 0;} void insert(int x,int y,int v[]) //建图{int temp;e[ne].v = y;temp = v[x];v[x] = ne;e[ne].next = temp;//ne++;} int id[NV],cnt,scnt,out[NV],post[NV];void dfs1(int s,int v[]) {int i;id[s] = 1; //第一次用来标记此点是否被访问过for(i = v[s];i != -1; i = e[i].next)   if(id[e[i].v] == -1)    dfs1(e[i].v,v);post[cnt++] = s; //按DFS的访问结束时间排序，第一次DFS没用 } void dfs2(int s,int v[]){int i;id[s] = scnt; //第二次标记改点是否被访问过，还有标记每一个点的所在的强连通分支for(i = v[s];i != -1; i = e[i].next)   if(id[e[i].v] == -1)    dfs2(e[i].v,v);} void Gscc() //求强连通分支 {int i,j;memset(id,-1,4*n);cnt = 0;for(i = 0;i n; i++)   if(id[i] == -1)dfs1(i,v2);memset(id,-1,4*n);scnt = 0;for(i = cnt - 1; i = 0; i--)   if(id[post[i]] == -1){dfs2(post[i],v1);scnt++;}memset(out,0,4*n);for(i = 0;i n; i++)    for(j = v1[i];j != -1; j = e[j].next)    if(id[i] != id[e[j].v]){out[id[i]]++;}} int cmp(void const *a,void const *b){return *((int *)a) - *((int *)b);} int main(){int i,j,k,x,y,temp,p[NV],m,anw[NV];bool hash[NV];while(scanf(%d,n) != EOFn){   init();   scanf(%d,m);   for( i = 0 ;i m ;i++)   {    scanf(%d%d,x,y);    x--;y--;    insert(x,y,v1);    insert(y,x,v2);   }   Gscc();   memset(hash,0,scnt);   for(i = 0; i scnt; i++)   {    if(out[i] == 0)hash[i] = 1;   }   for(i = 0 ,x = 0;i n ;i++)    if(hash[id[i]] == 1)anw[x++] = i;   qsort(anw,x,sizeof(int),cmp);   for(i = 0;i x;i ++)    printf(%d ,anw[i]+1);   printf(\n);}return 1;}