拓扑排序的扩展

拓扑图：有向图中，在排完（拓扑）序以后，会有每个条边都指向后面的点

拓扑图等价于有向无环图

拓扑排序的步骤：

1.先将所有入度为0的点入队

2.宽搜

---

1.家谱树（拓扑排序）

信息学奥赛一本通（C++版）在线评测系统

本题要求父亲比儿子先输出 那么就是入度为0的点先输出 然后再遍历他的所有孩子输出 裸题求拓扑序

#include <bits/stdc++.h>
using namespace std;
const int N=110,M=N*N/2;
int n,m;
int h[N],e[M],ne[M],idx;
int q[N];
int d[N];
void add(int a,int b)
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void topsort()
{
    int hh=0,tt=-1;
    for(int i=1;i<=n;i++)
    {
        if(!d[i]) q[++tt]=i;
    }
    while(hh<=tt)
    {
        int t=q[hh++];
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(--d[j]==0)
            {
                q[++tt]=j;
            }
        }
    }
}
int main()
{
    cin>>n;
    memset(h,-1,sizeof h);
    for(int i=1;i<=n;i++)
    {
        int son;
        while(cin>>son,son)
        {
            add(i,son);
            d[son]++;
        }
    }
    topsort();
    for(int i=0;i<n;i++)
        cout<<q[i]<<" ";
    return 0;
}

2.奖金(拓扑排序+最长路)

信息学奥赛一本通（C++版）在线评测系统

简化版的差分约束求最小值 那么就是最长路

然后本题有一个特别特殊的性质：边权全都是1 在前面分析几节内容中 如果是不能满足差分约束条件的话就是最长路不存在正环

然后再强连通分量那里也有一个特殊的题 就是差分约束的图的规定的边权是0和1 那么在这种图中只有在同一个连通分量的点中的边的权值都是0

让这里就更加简化了 边权全为0（如果a比b多 那么就可以成a>=b+1 就是从b往a连一条长度为1的边） 那么这样一优化 就是一旦存在强连通分量（存在环）就一定无解了 那么这图是拓扑图就是有解的 那么就是在拓扑图上求最长路

总结：

 细节：由于每个员工的奖金至少为100 那么可以建一个虚拟源点 往所有边都加一条100的边

            然后实际上不用建虚拟源点 只需要把每个点的初值弄为100就行了

#include <bits/stdc++.h>
using namespace std;
const int N=10010,M=20010;
int n,m;
int h[N],e[M],ne[M],idx;
int q[N];
int d[N];
int dis[N];
void add(int a,int b)
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
bool topsort()
{
    int hh=0,tt=-1;
    for(int i=1;i<=n;i++)
    {
        if(!d[i])
            q[++tt]=i;
    }
    while(hh<=tt)
    {
        int t=q[hh++];
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(--d[j]==0)
            {
                q[++tt]=j;
            }
        }
    }
    return tt==n-1;
}
int main()
{
    cin>>n>>m;
    memset(h,-1,sizeof h);
    while(m--)
    {
        int a,b;
        cin>>a>>b;
        add(b,a);
        d[a]++;
    }
    if(!topsort()) puts("Poor Xed");
    else
    {
        for(int i=1;i<=n;i++)
            dis[i]=100;
        for(int i=0;i<n;i++)
        {
            int j=q[i];
            for(int k=h[j];~k;k=ne[k])
                dis[e[k]]=max(dis[e[k]],dis[j]+1);
        }
        int res=0;
        for(int i=1;i<=n;i++)
            res+=dis[i];
        cout<<res<<endl;
    }
    return 0;
}

3.可达性统计（拓扑排序+dp）

164. 可达性统计 - AcWing题库

拓扑图都能dp

按照拓扑排序的逆序去dp

分析一下时间复杂度(n*m/32) 假定bitset使用unsigned int 实现的

1.先说普通方法 如果是用一个数组st[i][j]存第i个点存了哪些数 那么我的 j 就要取到 N。定义为第i个点能不能取到第 j 个点 那么我每一次去找就是n的时间复杂度 然后遍历m条边就是n*m的复杂度

由于unsigned int是32位的 那么我要存N位的01串的话就需要N位 进可以对应为 N/32个unsigned int 那么我每次去找只需要N/32的时间复杂度(bitset内部就是这个时间复杂度) 然后遍历m条边 就ok了

压位详情

 

#include <bits/stdc++.h>
using namespace std;
const int N=30010,M=30010;
int n,m;
int h[N],e[M],ne[M],idx;
int d[N],q[N];
bitset<N> f[N];
void add(int a,int b)
{
    e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void topsort()
{
    int hh=0,tt=-1;
    for(int i=1;i<=n;i++)
        if(!d[i]) q[++tt]=i;
    while(hh<=tt)
    {
        int t=q[hh++];
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(--d[j]==0)
            {
                q[++tt]=j;
            }
        }
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    memset(h,-1,sizeof h);
    for(int i=0;i<m;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        add(a,b);
        d[b]++;
    }
    topsort();
    for(int i=n-1;i>=0;i--)
    {
        int j=q[i];
        f[j][j]=1;
        for(int k=h[j];~k;k=ne[k])
            f[j]|=f[e[k]];
    }
    for(int i=1;i<=n;i++)
        cout<<f[i].count()<<endl;
    return 0;
}

4.车站分级

456. 车站分级 - AcWing题库

还是简单的差分约束问题 与第二题一样 最长路+拓扑排序 但是由于每一次建边都会建n*n条边 所以最多的边有2.5*1e8条边 用领接表去存的话 也要循环所有的边和点 所以一定超时

优化:由于每次建边都是集合两边建边，假如直接建边就大概建很多条边，容易超内存，所有每个车次都建个虚拟中间点为n+i,则建 虚拟点>=a , b>=虚拟点+1就行 这样边数就会变成n+m(原来是n*m)

#include<bits/stdc++.h>
using namespace std;
const int N=2010,M=1000010;
int h[N],e[M],ne[M],w[M],idx;
int n,m;
int d[N],q[N];
int dist[N];
bool st[N];
void add(int a,int b,int c)
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
    d[b]++;
}
void topsort()
{
    int hh=0,tt=-1;
    for(int i=1;i<=n+m;i++)//因为增加了虚拟中间点
        if(!d[i])
             q[++tt]=i;
    while(hh<=tt)
    {
        int t=q[hh++];
        for(int i=h[t];~i;i=ne[i])
        {
            int j=e[i];
            if(--d[j]==0)
                q[++tt]=j;
        }
    }
}
int main()
{
  memset(h,-1,sizeof h);
  cin>>n>>m;
  for(int i=1;i<=m;i++)
  {
      memset(st,0,sizeof st);//清空上一层的车站
     int cnt;
     cin>>cnt;
     int start=n,end=1;//获取开始车站与结尾车站
     while(cnt--)
     {
         int stop;
         cin>>stop;
         start=min(start,stop);
         end=max(end,stop);
         st[stop]=true;//标记这个是车站
     }
     int ver=n+i;//虚拟中间点
     for(int k=start;k<=end;k++)//枚举经过的所有车站
        if(st[k]) add(ver,k,1);//假如是车站，则ver->k条1的边
        else add(k,ver,0);//反之连k->ver条0的边

  }
  topsort();//跑一遍拓扑排序
  for(int i=1;i<=n;i++) dist[i]=1;//初始化为最少为1
  for(int i=0;i<n+m;i++)//取出拓扑序
  {
      int j=q[i];
      for(int k=h[j];~k;k=ne[k])//枚举子节点
        dist[e[k]]=max(dist[e[k]],dist[j]+w[k]);//更新最长路
  }
  int res=0;
  for(int i=1;i<=n;i++) res=max(res,dist[i]);//更新最大值
  cout<<res<<endl;
    return 0;
}