2013黑龙江省赛G Fly（缩点+最长路）

有n个城市，从起点出发最多可以到几个城市。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;
#define clr(x)memset(x,0,sizeof(x))
#define INF 0x1f1f1f1f
#define maxn 100010
struct node{
    int to,next,w;
}e[150010],e2[150010];
int tot,tot2;
int head[100010];
int head2[100010];
void add(int s,int u){
    e[tot].to=u;
    e[tot].next=head[s];
    head[s]=tot++;
}
void add2(int s,int u,int w){
    e2[tot2].to=u;
    e2[tot2].w=w;
    e2[tot2].next=head2[s];
    head2[s]=tot2++;
}
int sn,ti,top;
bool ins[maxn];
int dfn[maxn];
int low[maxn];
int stack[maxn];
int col[maxn];
void tarjan(int u){
    dfn[u]=low[u]=++ti;
    stack[++top]=u;
    ins[u]=true;
    int i,k;
    for(i=head[u];i;i=e[i].next){
        k=e[i].to;
        if(dfn[k]==0){
            tarjan(k);
            if(low[k]<low[u])low[u]=low[k];
        }
        else if(ins[k]&&dfn[k]<low[u])
            low[u]=dfn[k];
    }
    if(dfn[u]==low[u]){
        sn++;
        do{
            k=stack[top--];
            ins[k]=false;
            col[k]=sn;
        }while(k!=u);
    }
}
int val[maxn];
void init(){
    clr(col);
    clr(low);
    clr(dfn);
    clr(ins);
    clr(head);
    clr(head2);
    clr(val);
    tot=1;
    tot2=1;
    sn=top=ti=0;
}
int bfs(int src,int n)
{
    int d[maxn];
    int vis[maxn];
    for(int i=1;i<=n;i++){vis[i]=0,d[i]=-INF;}
    queue<int>q;
    while(!q.empty())q.pop();
    q.push(src);
    vis[src]=1;
    d[src]=0;
    while(!q.empty()){
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head2[u];i;i=e2[i].next){
            int v =e2[i].to;
            int w =e2[i].w;
            if(d[u]+w>d[v]){
                d[v]=d[u] + w;
                if(!vis[v]){
                    q.push(v);
                    vis[v]=1;
                }
            }
        }
    }
    int res=-INF;
    for(int i=1;i<=n;i++)res=max(res,d[i]);
    return res;
}

int main(){
    int a,b,n,m,i,j,k,res,st;
    while( scanf("%d%d%d",&n,&m,&st)!=EOF){
        init();
        while(m--){
            scanf("%d%d",&a,&b);
            add(a,b);
        }
        for(i=1;i<=n;i++)
            if(dfn[i]==0)tarjan(i);
        for(i=1;i<=n;i++){
            val[col[i]]++;
        }
        for(i=1;i<=n;i++)
            for(j=head[i];j;j=e[j].next){
                k=e[j].to;
                if(col[k]!=col[i]){
                    add2(col[i],col[k],val[col[k]]);
                }
            }
        int ans=bfs(col[st],sn);
        printf("%d\n",ans+val[col[st]]);
    }
    return 0;
}