题目描述
传送门
题意:求有最多的能到达它的点的点集,第一行输出最多的能到达它的点,第二行输出所有满足最大值的点集。
题解
首先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");
}
}