BZOJ-[Noi2010]Plane 航空管制2（贪心+并查集

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=2535

事实上这个就是航空管制，只是前面那个没有SPJ没法AC而已额。。。

首先对于第一个问题，对于一对限制a,b，要求a比b先起飞，在一个有向图中连边b->a，这样之后进行一次拓扑排序，尽量把起飞 顺序放后，可以发现在有解情况下，这样一定可以出解；

对于第二个问题，考虑到其他航班的影响，我们希望当前航班除外的都尽量靠后，但是可以由当前航班代表的节点到达的节点除外，删去这些点之后再做一次第一问，那么剩下的所有空位中最后一个就是答案了。

代码：

#include <cstdio>

#include <algorithm>

#include <cstring>

#include <bitset>

 

using namespace std ;

 

#define travel( x ) for ( edge *p = head[ x ] ; p ; p = p -> next )

#define rep( i , x ) for ( int i = 0 ; i ++ < x ; )

#define DOWN( i , r , l ) for ( int i = r ; i >= l ; -- i )

 

const int maxn = 2010 , maxm = 10100 ;

 

struct edge {

    edge *next ;

    int t ;

} E[ maxm ] ;

 

edge *pt = E , *head[ maxn ] ;

 

inline void addedge( int s , int t ) {

    pt -> t = t , pt -> next = head[ s ] ; head[ s ] = pt ++ ;

}

 

int father[ maxn ] , n , m , ti[ maxn ] ;

 

inline int getfa( int now ) {

    int i = now ; for ( ; father[ i ] ; i = father[ i ] ) ;

    for ( int j = now , k ; father[ j ] ; k = father[ j ] , father[ j ] = i , j = k ) ;

    return i ;

}

 

int d[ maxn ] , sta[ maxn ] , tp , dp[ maxn ] , ans[ maxn ] ;

 

inline void Dp( int del ) {

    rep( i , n ) d[ i ] = ans[ i ] = father[ i ] = 0 , dp[ i ] = min( n , ti[ i ] ) ;

    rep( i , n ) travel( i ) ++ d[ p -> t ] ;

    tp = 0 ;

    rep( i , n ) if ( ! d[ i ] ) sta[ ++ tp ] = i ;

    int now , pos ;

    while ( tp ) {

        if ( ( now = sta[ tp -- ] ) == del ) continue ;

        pos = getfa( dp[ now ] ) ; ans[ pos ] = now , father[ pos ] = pos - 1 ;

        travel( now ) {

            dp[ p -> t ] = min( dp[ p -> t ] , dp[ now ] ) ;

            if ( ! ( -- d[ p -> t ] ) ) sta[ ++ tp ] = p -> t ;

        }

    }

}

 

int main(  ) {

    memset( head , 0 , sizeof( head ) ) ;

    scanf( "%d%d" , &n , &m ) ;

    rep( i , n ) scanf( "%d" , ti + i ) ;

    while ( m -- ) {

        int s , t ; scanf( "%d%d" , &t , &s ) ;

        addedge( s , t ) ;

    }

    Dp( 0 ) ;

    rep( i , n ) printf( "%d " , ans[ i ] ) ; printf( "\n" ) ;

    rep( i , n ) {

        Dp( i ) ;

        DOWN( i , n , 1 ) if ( ! ans[ i ] ) {

            printf( "%d " , i ) ; break ;

        }

    }

    printf( "\n" ) ;

    return 0 ;

}