hdu5195 二分+线段树+拓扑序

这题说的给了n个点m条边要求保证是一个有向无环图，可以删除至多k条边使得这个图的拓扑序的字典序最大，我们知道如果我们要排一个点的时候一定要考虑比他大的点是否可以、通过拆边马上拆出来，如果可以拆当然是拆，肯定保证字典序最大，如果不能拆，就不拆留着以后拆，当初这个比他大的点度数小于k的，最大是多少，这个方法我一直想不出，后来看了题解，二分加线段树，可以做到，线段树维护每个点的d[i],然后通过二分找出小于k的最大点是多少。

#include <iostream>
#include <algorithm>
#include <string.h>
#include <vector>
#include <cstdio>
#include<queue>
using namespace std;
const int maxn =100005;
struct Edge{
   int u,v;
   Edge(int uu=0, int vv=0){
       u= uu; v =vv;
   }
};
vector<Edge> E;
vector<int> G[maxn],REG[maxn];
void add_edg(int u, int v){
    E.push_back(Edge(u,v));
    G[u].push_back(E.size()-1);
    REG[v].push_back(E.size()-1);
}
int ind[maxn];
int cL, cR,cans;
struct Itree{
    int id[maxn*4];
    void build(int O, int L, int R){
        if(L==R){
              id[O] = ind[L] ; return ;
        }
        int mid = ( L + R ) >> 1 ;
        build(O*2,L,mid);
        build(O*2+1,mid+1,R);
        id[O]=min(id[O*2],id[O*2+1]);
    }
    void update(int o, int L, int R){
        if(  L == R ){
              id[o] = cans;  return ;
        }
        int mid=(L+R)>>1;
        if(cL<=mid){
             update(o*2,L,mid);
        }else update(o*2+1,mid+1,R);
        id[o]=min(id[o*2],id[o*2+1]);
    }
    void query(int o , int L, int R){
        if(cL<=L && cR>=R){
             cans =min(cans,id[o]);return ;
        }
        int mid = (L+R)>>1;
        if(cL <= mid) query(o*2,L,mid);
        if(cR>mid) query(o*2+1,mid+1,R);
    }
}T;
bool use[maxn];
void init(int n){
    for(int i=0; i<n; i++){
         G[i].clear();
         REG[i].clear();
    }
    E.clear();
    memset(ind,0,sizeof(ind));
    memset(use,true,sizeof(use));
}
int jud(int k, int n){
   int L=1,R=n,ans=-1;
   while(L<=R){
        int mid = (L+R)>>1;
        cans = k+1;
        cL=mid; cR=R;
        T.query(1,1,n);
        if(cans<=k){
             ans=mid; L=mid+1;
        }
        else R=mid-1;
   }
   return ans;
}
priority_queue<int> Q;
void solve1(int nk,int k,int n){

     for(int i=0; i<REG[nk].size(); i++){
          use[ REG[nk][i] ] = false;
     }
     cans = k+1;
     cL=cR=nk;
     T.update(1,1,n);
}
void solve2(int nk,int k,int n){
    for(int i=0; i<G[nk].size(); i++){
        int numedg = G[nk][i];
        if(use[numedg]){
             use[numedg]=false;
             Edge q = E[numedg];
             ind[q.v]--;
             if(ind[q.v]==0){
                Q.push(q.v); ind[q.v]=k+1;
             }
             cL = cR = q.v;
             cans =ind[q.v];
             T.update(1,1,n);
        }
    }
}
int ans[maxn];
int main()
{
     int n,m,k;

     while(scanf("%d%d%d",&n,&m,&k)==3){
         init(n);
         for(int i=0; i<m; i++){
             int u,v;
             scanf("%d%d",&u,&v);
             add_edg(u,v);
             ind[v]++;
         }
         while(!Q.empty())Q.pop();
         for(int i =1; i<=n; i++){
             if(ind[i]==0){
                 Q.push(i);
                 ind[i]=k+1;
             }
         }
         T.build(1,1,n);
         int st=0;
         while(!Q.empty()){
             int top = Q.top();
             int nk = jud(k,n);
             if(nk>top){
                 k-=ind[nk];
                 solve1(nk,k,n);
                 Q.push(nk);
             }else{
                 Q.pop();
                 ans[st++]=top;
                 solve2(top,k,n);
             }
         }
         for(int i=0; i<n-1; i++) printf("%d ",ans[i]);
         printf("%d
",ans[n-1]);
     }

    return 0;
}