[HDU3639]Hawk-and-Chicken（tarjan+dfs）

题目描述

传送门
题意：求有最多的能到达它的点的点集，第一行输出最多的能到达它的点，第二行输出所有满足最大值的点集。

题解

首先Tarjan缩点，之后反向建图，可以知道拥有最大值的点一定在入度为0的点里。
那么可以从每个入度为0的点出发dfs。
刚开始想的是正向建图然后递推到它的儿子，发现貌似有些混乱，于是改姿势了。
而且这题还卡格式。

代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
using namespace std;

const int max_n=5e3+5;
const int max_m=3e4+5;

int T,n,m,x,y,temp,N,cnt,Max,ansnum;
int tot,point[max_n],nxt[max_m],v[max_m];
int Dfn[max_n],Low[max_n],belong[max_n],num[max_n],in[max_n],f[max_n],strack[max_n],ans[max_n];
bool vis[max_n],pd[max_n];
struct hp{int x,y;}edge[max_m];
queue <int> q;

inline void clear(){
    x=y=temp=N=cnt=Max=ansnum=0;
    for (int i=1;i<=n;++i) Dfn[i]=Low[i]=belong[i]=num[i]=in[i]=f[i]=ans[i]=vis[i]=pd[i]=strack[i]=0;
}
inline void addedge(int x,int y){++tot;nxt[tot]=point[x];point[x]=tot;v[tot]=y;}
inline void tarjan(int now){
    Dfn[now]=Low[now]=++N; vis[now]=true; strack[++temp]=now;
    for (int i=point[now];i;i=nxt[i])
      if (!Dfn[v[i]]){
        tarjan(v[i]);
        Low[now]=min(Low[now],Low[v[i]]);
      }
      else if (vis[v[i]]) Low[now]=min(Low[now],Dfn[v[i]]);
    if (Dfn[now]==Low[now]){
        ++cnt;
        while (strack[temp]!=now) belong[strack[temp]]=cnt,num[belong[strack[temp]]]++,vis[strack[temp]]=false,--temp;
        belong[strack[temp]]=cnt,num[belong[strack[temp]]]++,vis[strack[temp]]=false,--temp;
    }
}
inline int dfs(int now){
    if (vis[now]) return 0;
    else vis[now]=true;
    int ans=num[now];
    for (int i=point[now];i;i=nxt[i])
      ans+=dfs(v[i]);
    return ans;
}
int main(){
    scanf("%d",&T);
    for (int t=1;t<=T;++t){
        scanf("%d%d",&n,&m);
        clear();
        tot=0; memset(point,0,sizeof(point)); memset(nxt,0,sizeof(nxt)); memset(v,0,sizeof(v));
        for (int i=1;i<=m;++i){
            scanf("%d%d",&x,&y),addedge(++x,++y);
            edge[i].x=x,edge[i].y=y;
        }
        for (int i=1;i<=n;++i) if (!Dfn[i]) tarjan(i);
        tot=0; memset(point,0,sizeof(point)); memset(nxt,0,sizeof(nxt)); memset(v,0,sizeof(v));
        for (int i=1;i<=m;++i)
          if (belong[edge[i].x]!=belong[edge[i].y])
            addedge(belong[edge[i].y],belong[edge[i].x]),in[belong[edge[i].x]]++;
        for (int i=1;i<=cnt;++i)
          if (!in[i]){
            memset(vis,0,sizeof(vis));
            int k=dfs(i)-1;
            if (k>Max){
                Max=k; ans[ansnum=1]=i;
            }
            else if (k==Max) ans[++ansnum]=i;
          }
        for (int i=1;i<=ansnum;++i) pd[ans[i]]=true;
        printf("Case %d: %d\n",t,Max);
        bool flag=false;
        for (int i=1;i<=n;++i)
          if (pd[belong[i]]){
            if (!flag){
                flag=true;
                printf("%d",i-1);
            }
            else printf(" %d",i-1);
          }
        printf("\n");
    }
}