poj 2553 （tarjan的应用）

趁着刚学这个算法，就连续刷刷。。。掌握才是关键。。。

这个题刚开始怎么也看不懂，后来在discuss里看到网友们有说“ if for every node w in G that is reachable from v, v is also reachable from w. ”这句很重要的时候 ，豁然开朗呀。。。只怪自己没怎么认真看题呀

 

/* * 关于sink：如果v可以到达w，那么w一定要能到达v；如果到不了，也是可以的。 * 这就是那句话if for every node w in G that is reachable from v, v is also reachable from w. */ #include <iostream> #include<cstdlib> using namespace std; const int maxl=5010; struct node{ int v,next; }edge[maxl*3]; int edgehead[maxl],k=0,n; int low[maxl],dfn[maxl],out[maxl],belong[maxl],inde=0,cnt=0; int stack[maxl],head=0; bool instack[maxl]; void addedge(int u,int v){ edge[k].v=v;edge[k].next=edgehead[u]; edgehead[u]=k++; } void tarjan(int u){ int v; dfn[u]=low[u]=inde++; instack[u]=1; stack[head++]=u; for(int i=edgehead[u];i!=-1;i=edge[i].next){ v=edge[i].v; if(dfn[v]==-1) { tarjan(v); low[u]=min(low[u],low[v]); } else if(instack[v]) low[u]=min(low[u],dfn[v]); } if(low[u]==dfn[u]) { cnt++; do{ v=stack[--head]; instack[v]=0; belong[v]=cnt; }while(v!=u); } } int main(int argc, char** argv) { int m,i,j,u,v; while(scanf("%d",&n) && n){ memset(edgehead,-1,sizeof(edgehead)); memset(instack,0,sizeof(instack)); memset(dfn,-1,sizeof(dfn)); memset(out,0,sizeof(out)); //刚开始没初始化，一直错 scanf("%d",&m); head=cnt=k=0; inde=1; while(m--){ scanf("%d%d",&u,&v); addedge(u,v); } for(i=1;i<=n;i++) if(dfn[i]==-1) tarjan(i); for(i=1;i<=n;i++) for(j=edgehead[i];j!=-1;j=edge[j].next) if(belong[i]!=belong[edge[j].v]){ out[belong[i]]++; break; } for(i=1;i<=n;i++) if(!out[belong[i]]) printf("%d ",i); printf("/n"); } return 0; }