http://poj.org/problem?id=2186
题意:有向图求能被所有其他点到达的点的个数。
题解:强连通缩点,求出度为0的分量,如果只有一个输出分量大小,如果大于一个,则这两个分量不连通,不存在这样的点。
代码:
#include <iostream>
#include<cstdio>
#include<stack>
#include<cstring>
using namespace std;
const int maxn = 10005;
const int maxm = 50005;
int n,m;
struct node{
int to,next;
}e1[maxm],e2[maxm],e3[maxm];
//bool mp[maxn][maxn];
int head1[maxn],head2[maxn],head3[maxn];
int belong[maxn];
bool vis1[maxn],vis2[maxn];
int tot1,tot2,tot3;
int cnt1,cnt2;
int st[maxn];
int num;
int setnum[maxn];
void add(int u,int v){
e1[tot1].to=v;
e1[tot1].next=head1[u];
head1[u]=tot1++;
e2[tot2].to=u;
e2[tot2].next=head2[v];
head2[v]=tot2++;
}
void add2(int u,int v){
e3[tot3].to=v;
e3[tot3].next=head3[u];
head3[u]=tot3++;
}
void dfs1(int u){
vis1[u]=true;
int i;
for(i=head1[u];i!=-1;i=e1[i].next){
if(!vis1[e1[i].to]){
dfs1(e1[i].to);
}
}
st[cnt1++]=u;
}
void dfs2(int u){
vis2[u]=true;
num++;
belong[u]=cnt2;
int i;
for(i=head2[u];i!=-1;i=e2[i].next){
if(!vis2[e2[i].to]){
dfs2(e2[i].to);
}
}
}
void solve(){
int i;
cnt1=cnt2=0;
memset(vis1,false,sizeof(vis1));
memset(vis2,false,sizeof(vis2));
for(i=1;i<=n;i++)
if(!vis1[i])dfs1(i);
for(i=cnt1-1;i>=0;i--){
if(!vis2[st[i]]){
num=0;
dfs2(st[i]);
setnum[cnt2++]=num;
}
}
}
int out[maxn];
void init(){
//memset(mp,false,sizeof(mp));
memset(head1,-1,sizeof(head1));
memset(head2,-1,sizeof(head2));
memset(head3,-1,sizeof(head1));
memset(out,0,sizeof(out));
tot1=tot2=tot3=0;
}
int main()
{
int a,b;
init();
scanf("%d%d",&n,&m);
int i,j,k;
for(i=1;i<=m;i++){
scanf("%d%d",&a,&b);
add(a,b);
}
solve();
for(i=1;i<=n;i++){
for(j=head1[i];j!=-1;j=e1[j].next){
if(belong[i]==belong[e1[j].to])continue;
int u=belong[i];
int v=belong[e1[j].to];
for(k=head3[u];k!=-1;k=e3[k].next){
if(e3[k].to==v)break;
}
add2(u,v);
out[u]++;
}
}
int tmp=-1;
for(i=0;i<cnt2;i++){
if(out[i]==0){
if(tmp==-1)tmp=setnum[i];
else tmp=-2;
}
}
if(tmp>0)printf("%d\n",tmp);
else printf("0\n");
return 0;
}