spfa + slf优化

最近在练习费用流 ， 不是要用spfa吗 ，我们教练说：ns学生写朴素的spfa说出去都让人笑 。

QwQ，所以就去学了一下优化 。

slf优化就是双向队列优化一下，本来想用lll优化，可是优化后我tm居然t了（那道题特地卡spfa），所以lll优化太迷了 ，还是只用slf优化好 。

 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <deque>
 const int inf =  <<  , maxn =  +  , M =  +   ;
 using namespace std ;//
 int n , m , head[maxn]  , dis[maxn] , cnt , sum , tot ;
 bool mark[maxn] ;
 struct id
 {
     int nxt ,to , val ;
 } edge[M] ;
 deque < int > Q ;

 inline void Init ( )
 {
     freopen( "NSOOJ#10719.in" , "r" , stdin  ) ;
     freopen( "NSOOJ#10719.out" , "w" , stdout ) ;
 }

 int read( )
 {
     char ch = getchar( ) ; int k =  , ret =  ;
     while( ch < '' || ch > '' ) { if( ch == '-' ) k = - ; ch = getchar( ) ; }
     while( ch >= '' && ch <= '' ) ret = ret *  + ch - '' , ch = getchar( ) ;
     return k * ret ;
 }

 void add( int u , int v , int va )
 {
     edge[++cnt].nxt = head[u] , edge[cnt].to = v ;
     edge[cnt].val = va , head[u] = cnt ;
 }

 void input(  )
 {
     n = read()  , m = read( ) ;
     int u ,v , c ;
     memset( head , - , sizeof(head)) ;
     for( int x =  ; x <= m ; ++x )
     {
         u = read( ) , v = read( ) , c = read( ) ;
         add( u ,v , c )  ;
     }
 }

 void spfa( )
 {
     memset( dis , / , sizeof(dis) ) ;
     dis[] =  , mark[] = true ;
     Q.push_back(  ) ;
     while( !Q.empty( ) )
     {
         int u = Q.front( ) ; Q.pop_front( ) ; mark[u] = false ;

         for( int i = head[u] ; ~i ; i = edge[i].nxt )
         {
             int v = edge[i].to ;
             if( dis[v] > dis[u] + edge[i].val )
             {
                 dis[v] = dis[u] + edge[i].val ;
                 if( !mark[v] )
                 {
                     mark[v] = true ;
                     if( Q.empty( ) || dis[v] > dis[Q.front( )]  ) Q.push_back( v ) ;
                     else Q.push_front( v ) ;

                 }

             }
         }
     }
     if( dis[n] ==  ) printf( "%d\n" , - ) ;
     else printf( "%d\n" , dis[n] ) ;
 }

 int main( )
 {
 //    Init( ) ;
     input( ) ;
     spfa( ) ;
 //    fclose( stdin ) ;
 //       fclose( stdout ) ;
     return  ;
 }

纠正，其实lll优化一般效率还是挺高的，只是学长说如果要卡的话可以卡到指数级（=_=）,所以下面是lll优化的代码 。

 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <deque>
 const int inf =  <<  , maxn =  +  , M =  +   ;
 using namespace std ;//
 int n , m , head[maxn]  , dis[maxn] , cnt , sum , tot ;
 bool mark[maxn] ;
 struct id
 {
     int nxt ,to , val ;
 } edge[M] ;
 deque < int > Q ;

 inline void Init ( )

 {
     freopen( "NSOOJ#10719.in" , "r" , stdin  ) ;
     freopen( "NSOOJ#10719.out" , "w" , stdout ) ;
 }

 int read( )
 {
     char ch = getchar( ) ; int k =  , ret =  ;
     while( ch < '' || ch > '' ) { if( ch == '-' ) k = - ; ch = getchar( ) ; }
     while( ch >= '' && ch <= '' ) ret = ret *  + ch - '' , ch = getchar( ) ;
     return k * ret ;
 }

 void add( int u , int v , int va )
 {
     edge[++cnt].nxt = head[u] , edge[cnt].to = v ;
     edge[cnt].val = va , head[u] = cnt ;
 }

 void input(  )
 {
     n = read()  , m = read( ) ;
     int u ,v , c ;
     memset( head , - , sizeof(head)) ;
     for( int x =  ; x <= m ; ++x )
     {
         u = read( ) , v = read( ) , c = read( ) ;
         add( u ,v , c )  ;
     }
 }

 void spfa( )
 {
     memset( dis , / , sizeof(dis) ) ;
     dis[] =  , mark[] = true ;
     Q.push_back(  ) ; tot =  ;
     while( !Q.empty( ) )
     {
         int u = Q.front( ) ; Q.pop_front( ) ; mark[u] = false ;
         tot-- ; sum -= dis[u] ;
         for( int i = head[u] ; ~i ; i = edge[i].nxt )
         {
             int v = edge[i].to ;
             if( dis[v] > dis[u] + edge[i].val )
             {
                 dis[v] = dis[u] + edge[i].val ;
                 if( !mark[v] )
                 {
                     mark[v] = true ;
                     if( Q.empty( ) || dis[v] * tot <= sum ) Q.push_back( v ) ;
                     else Q.push_front( v ) ;
                     tot++ ; sum += dis[v] ;
                 }

             }
         }
     }
     if( dis[n] ==  ) printf( "%d\n" , - ) ;
     else printf( "%d\n" , dis[n] ) ;
 }

 int main( )
 {
 //    Init( ) ;
     input( ) ;
     spfa( ) ;
 //    fclose( stdin ) ;
 //       fclose( stdout ) ;
     return  ;
 }

好的，下面是两个都加的代码 。

 #include <iostream>
 #include <cstdlib>
 #include <cstring>
 #include <cstdio>
 #include <deque>
 const int inf =  <<  , maxn =  +  , M =  +   ;
 using namespace std ;//
 int n , m , head[maxn]  , dis[maxn] , cnt , sum , tot ;
 bool mark[maxn] ;
 struct id
 {
     int nxt ,to , val ;
 } edge[M] ;
 deque < int > Q ;

 inline void Init ( )
 {
     freopen( "NSOOJ#10719.in" , "r" , stdin  ) ;
     freopen( "NSOOJ#10719.out" , "w" , stdout ) ;
 }

 int read( )
 {
     char ch = getchar( ) ; int k =  , ret =  ;
     while( ch < '' || ch > '' ) { if( ch == '-' ) k = - ; ch = getchar( ) ; }
     while( ch >= '' && ch <= '' ) ret = ret *  + ch - '' , ch = getchar( ) ;
     return k * ret ;
 }

 void add( int u , int v , int va )
 {
     edge[++cnt].nxt = head[u] , edge[cnt].to = v ;
     edge[cnt].val = va , head[u] = cnt ;
 }

 void input(  )
 {
     n = read()  , m = read( ) ;
     int u ,v , c ;
     memset( head , - , sizeof(head)) ;
     for( int x =  ; x <= m ; ++x )
     {
         u = read( ) , v = read( ) , c = read( ) ;
         add( u ,v , c )  ;
     }
 }

 void spfa( )
 {
     memset( dis , / , sizeof(dis) ) ;
     dis[] =  , mark[] = true ;
     Q.push_back(  ) ; tot =  ;
     while( !Q.empty( ) )
     {
         int u = Q.front( ) ; Q.pop_front( ) ; mark[u] = false ;
         tot-- ; sum -= dis[u] ;
         for( int i = head[u] ; ~i ; i = edge[i].nxt )
         {
             int v = edge[i].to ;
             if( dis[v] > dis[u] + edge[i].val )
             {
                 dis[v] = dis[u] + edge[i].val ;
                 if( !mark[v] )
                 {
                     mark[v] = true ;
                     if( Q.empty( ) || dis[v] > dis[Q.front( )] || dis[v] * tot <= sum ) Q.push_back( v ) ;
                     else Q.push_front( v ) ;
                     tot++ ; sum += dis[v] ;
                 }

             }
         }
     }
     if( dis[n] ==  ) printf( "%d\n" , - ) ;
     else printf( "%d\n" , dis[n] ) ;
 }

 int main( )
 {
 //    Init( ) ;
     input( ) ;
     spfa( ) ;
 //    fclose( stdin ) ;
 //       fclose( stdout ) ;
     return  ;
 }