标签:
Time Limit: 12000/4000 MS (Java/Others) Memory Limit: 104857/104857 K (Java/Others)
Total Submission(s): 3568 Accepted Submission(s): 1235
#include<stdio.h> #include<string.h> #include<algorithm> #include<stack> #include<vector> #define MAX 50010 #define INF 0x3f3f3f using namespace std; int n,m; int ans,head[MAX]; int low[MAX],dfn[MAX]; int instack[MAX],sccno[MAX]; vector<int>newmap[MAX]; vector<int>scc[MAX]; int scccnt,dfsclock; int in[MAX],out[MAX]; stack<int>s; struct node { int beg,end,next; }edge[MAX]; void init() { ans=0; memset(head,-1,sizeof(head)); } void add(int beg,int end) { edge[ans].beg=beg; edge[ans].end=end; edge[ans].next=head[beg]; head[beg]=ans++; } void getmap() { int i,a,b; while(m--) { scanf("%d%d",&a,&b); add(a,b); } } void tarjan(int u) { int v,i,j; s.push(u); instack[u]=1; low[u]=dfn[u]=++dfsclock; for(i=head[u];i!=-1;i=edge[i].next) { v=edge[i].end; if(!dfn[v]) { 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]) { scccnt++; while(1) { v=s.top(); s.pop(); instack[v]=0; sccno[v]=scccnt; if(v==u) break; } } } void find(int l,int r) { memset(low,0,sizeof(low)); memset(dfn,0,sizeof(dfn)); memset(instack,0,sizeof(instack)); memset(sccno,0,sizeof(sccno)); dfsclock=scccnt=0; for(int i=l;i<=r;i++) { if(!dfn[i]) tarjan(i); } } void suodian() { int i; for(i=1;i<=scccnt;i++) { newmap[i].clear(); in[i]=0;out[i]=0; } for(i=0;i<ans;i++) { int u=sccno[edge[i].beg]; int v=sccno[edge[i].end]; if(u!=v) { newmap[u].push_back(v); in[v]++;out[u]++; } } } void solve() { int i,j; if(scccnt==1) { printf("0\n"); return ; } else { int minn=0; int maxx=0; for(i=1;i<=scccnt;i++) { if(!in[i]) minn++; if(!out[i]) maxx++; } printf("%d\n",max(minn,maxx)); } } int main() { while(scanf("%d%d",&n,&m)!=EOF) { init(); getmap(); find(1,n); suodian(); solve(); } return 0; }
hdoj 3836 Equivalent Sets【scc&&缩点】【求最少加多少条边使图强连通】
标签:
原文地址:http://www.cnblogs.com/tonghao/p/4820724.html