题目: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 ;
}