路径

在橙子国里,有 n+1 座城市,其中城市 n+1 为首都,其中有 m 条有向公路各个城市连接起
来。从每个城市 不一定能 走到首都,其中不同的定义是:每条边可以走多次,如 果走边的顺
序有一条不同即称两方式不同。
橙子国的居民想知道到达首都方式最多的城市到达首都的方法数,以及有多少个城市有这么
多方式,按照顺序输出城市编号。
如果最多不同方式超过了 36500 那么视做全部相等,方法数输出 zawsze。
输入格式
第 1 行,两个正整数,n 和 m。(n,m<=1000000)
第 2 至 m+1 行,两个正整数 u 和 v,代表存在一条城市 u 通向城市 v 的道路。
输出格式
第 1 行,一个正整数,代表到达首都方式最多的城市到达首都的方法数。
第 2 行,一个正整数 k,代表有多少个城市有这么多方式。
第 3 行,k 个正整数,代表拥有最多达到首都方式的城市编号,从小到大。
输入样例
35
12
13
23
34
34
输出样例
4
1
1

题解:

将边反过来,去掉不能到达n+1的点,然后统计入度,跑拓扑排序,注意存在环的情况(存在入度没有清0的点)

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
const int N=1000005;
int gi(){
    int str=0;char ch=getchar();
    while(ch>'9' || ch<'0')ch=getchar();
    while(ch>='0' && ch<='9')str=(str<<1)+(str<<3)+ch-48,ch=getchar();
    return str;
}
int n,m,q[N],num=0,head[N],lim=36501;
struct Lin{
    int next,to;
}a[N];
void init(int x,int y){
    a[++num].next=head[x];
    a[num].to=y;
    head[x]=num;
}
bool vis[N],mark[N];int f[N],st[N],top=0,mx=0,du[N];
void hfs(){
    int t=0,x,sum=1,u;q[1]=n+1;f[n+1]=1;
    while(t!=sum){
        t++;if(t==N)t-=N;x=q[t];
        for(int i=head[x];i;i=a[i].next){
            u=a[i].to;
            f[u]+=f[x];du[u]--;
            if(f[u]>lim)f[u]=lim;
            if(f[u]>mx){
                mx=f[u];
            }
            if(!du[u]){
                sum++;if(sum==N)sum-=N;q[sum]=u;
            }
        }
    }
}
void prework(){
    int t=0,x,sum=1,u;q[1]=n+1;mark[n+1]=true;
    while(t!=sum){
        t++;x=q[t];
        for(int i=head[x];i;i=a[i].next){
            u=a[i].to;
            if(mark[u])continue;
            mark[u]=true;q[++sum]=u;
        }
    }
    for(int i=1;i<=n+1;i++){
        if(!mark[i])continue;
        for(int j=head[i];j;j=a[j].next){
            du[a[j].to]++;
        }
    }
}
void work(){
    n=gi();m=gi();
    int x,y;
    for(int i=1;i<=m;i++){
        x=gi();y=gi();
        init(y,x);
    }
    prework();
    hfs();
    for(int i=1;i<=n+1;i++){
        if(du[i])f[i]=lim,mx=lim;
    }
    if(mx>=lim)printf("zawsze
");
    else printf("%d
",mx);
    for(int i=1;i<=n;i++){
        if(f[i]==mx)st[++top]=i;
    }
    printf("%d
",top);
    for(int i=1;i<=top;i++)printf("%d ",st[i]);
}
int main()
{
    freopen("path.in","r",stdin);
    freopen("path.out","w",stdout);
    work();
    return 0;
}