bzoj3832 [Poi2014]Rally（拓扑排序+堆）

非常厉害的思路.
设f[i],g[i]分别表示从起点到i的最长路径和从点i开始的最长路径.
这个显然可以用拓扑排序处理.

对于一条边(u,v),它能产生的最长路径是f[u]+g[v]+1
我们用一个可删除堆来维护当前所有合法边产生的最长路径的最大值.

首先把g[i]数组加入线段树.

按拓扑序依次枚举每一个点.
把指向它的边产生的最长路径和从这个点开始的最长路径从堆中删除.
当前全局最大值就是删掉这个点的最长路径.
然后把这个点连出去的边产生的最长路径加入堆
也要把到这个点的最长路径加入堆

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,m,h[N],num=0,h1[N],du[N],q[N],qh=1,qt=0,f[N],g[N],ans,id;
struct edge{
    int to,next;
}data[N<<1],data1[N<<1];
struct Heap{
    priority_queue<int>q1,q2;
    inline void push(int x){q1.push(x);}
    inline void del(int x){q2.push(x);}
    inline int top(){while(!q2.empty()&&q1.top()==q2.top()) q1.pop(),q2.pop();return q1.top();}
}Q;
int main(){
//  freopen("a.in","r",stdin);
    n=read();m=read();ans=n;
    while(m--){
        int x=read(),y=read();du[y]++;
        data[++num].to=y;data[num].next=h[x];h[x]=num;
        data1[num].to=x;data1[num].next=h1[y];h1[y]=num;
    }for(int i=1;i<=n;++i) if(!du[i]) q[++qt]=i;
    while(qh<=qt){
        int x=q[qh++];
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(--du[y]==0) q[++qt]=y;
        }
    }for(int ii=n;ii>=1;--ii){
        int x=q[ii];
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;g[x]=max(g[x],g[y]+1);
        }
    }for(int ii=1;ii<=n;++ii){
        int x=q[ii];
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;f[y]=max(f[y],f[x]+1);
        }
    }for(int i=1;i<=n;++i) Q.push(g[i]);
    for(int i=1;i<=n;++i){
        int x=q[i];
        for(int i=h1[x];i;i=data1[i].next){
            int y=data1[i].to;Q.del(f[y]+g[x]+1);
        }Q.del(g[x]);int res=Q.top();
        if(res<ans) ans=res,id=x;
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;Q.push(f[x]+g[y]+1);
        }Q.push(f[x]);
    }printf("%d %d\n",id,ans);
    return 0;
}