首页 > 代码库 > CodeForces 449B - Jzzhu and Cities

CodeForces 449B - Jzzhu and Cities

传送门:Jzzhu and Cities

 

题意:

题意一定要看仔细,这里说一个国家,有N个城市,城市与城市之间有M条路,距离为W,编号为1的城市是首都,1和K个城市有火车道,距离为X

现在国家要节省开支,想把其中部分铁路关闭,但是又想关闭之后从首都到每一个城市本来的最短距离不变,求最多能关闭多少条铁路。

 

分析:

注意最重要的几点:

1、求关闭最多的铁路

2、关闭之后,从首都到每一个城市的最短距离不变。

3、从首都到某个城市的铁路可能有多条

 

算法:

很明显,求最短距离不变,那么就是最短路径问题,DIJK算法用上。

从首都到某个城市,有两种方式

1、走路

2、坐火车

但是为了节省,如果走路的距离==坐火车,那么这条火车道就干掉。

那么好了,知道这点,我们就必须先求走路的距离,再求一下坐火车的距离,判断一下是否相等就可以了。

如果火车比走路更远,干掉。

但是要注意一点:

可能从首都有多条铁路到某一个城市,那么一定一定要想到,我们只需要最短那条铁道来判断,其他比他远的,留来干嘛!!!!(我这里没判重,坑了好久)。

接下来就是第二点:

从首都到某个城市的距离,我们中途是可以坐火车的!!!!那么这个距离可能跟你从首都直接坐火车到这个城市的距离相等,这条铁路也得干掉!!!!(第二个坑)

 

好了,现在就可以求最多多少条了。

 

首先,最短路径,先去掉所有铁路求一次,那么就得到从首都走路到其他城市的所有最短距离

 

再把铁路全部加上,那么就求到(可能)加上坐火车到其他城市的最短距离。

 

之后就是判断,这条铁路是否需要就可以了。

 

判断的依据是,在DIJK算法里面加一个pre[v] = u ,记录v是从u这个城市过来的。

 

转移的时候,如果u是首都1,并且存在一个u‘,跟从首都1直接坐火车到v的距离相等,那么我们就抛弃从首都坐火车这条,改为从首都到 u‘再走路过来v这条路,那么就又少一条铁路了。

 

之后把这些情况加上,就是最多能删多少条铁路了。

 

 

 

  1 #ifndef ONLINE_JUDGE  2 #define __DEBUG__ 0  3 #endif  4 /*  5 * =======================================  6 *    FileName: code.cpp  7 *        Desc:   8 *      Author: vinceliang  9 *       Email: liangguoqiu@gmail.com 10 *    HomePage:  11 *     Version: 0.0.1 12 *  LastChange: 2014-07-15 12:16:49 13 *     History: 14 *======================================== 15 */ 16 #ifdef __DEBUG__ 17 #define Log(...) printf(__VA_ARGS__) 18 #else 19 #define Log(...) // 20 #endif 21 #include <algorithm> 22 #include <string> 23 #include <complex> 24 #include <cassert> 25 #include <memory> 26 #include <set> 27 #include <stack> 28 #include <map> 29 #include <list> 30 #include <deque> 31 #include <numeric> 32 #include <cctype> 33 #include <cstddef> 34 #include <vector> 35 #include <queue> 36 #include <iostream> 37 #include <iomanip> 38 #include <iterator> 39 #include <cmath> 40 #include <cstdio> 41 #include <cstdlib> 42 #include <sstream> 43 #include <fstream> 44 #include <ctime> 45 #include <cstring> 46 #include <functional> 47 #include <bitset> 48 using namespace std; 49 #if defined(_MSC_VER) || defined(__BORLANDC__) 50 typedef unsigned __int64 uint64; 51 typedef unsigned __int64 ULL; 52 typedef signed __int64 int64; 53 typedef signed __int64 LL; 54 #else 55 typedef unsigned long long uint64; 56 typedef unsigned long long ULL; 57 typedef signed long long int64; 58 typedef signed long long LL; 59 #endif 60 typedef vector<int> VI; 61 typedef vector<string> VS; 62 typedef pair<int, int> PII; 63 typedef pair<int64, int64> PLL; 64 typedef vector<int64> VL; 65 #define pb push_back 66 #define pob pop_back 67 #define mp make_pair 68 #define fi first 69 #define se second 70 #define PII pair< int,int > 71 #define PSS pair< string,string > 72 #define PDD pair< double,double > 73 #define PLL pair< LL, LL > 74 #define FOR( i, a, b ) for ( int i=(a), _n=(b); i < _n; ++i ) 75 #define FORE( i, a, b ) for ( int i=(a), _n=(b); i <= _n; ++i ) 76 #define FORD( i, a, b ) for ( int i=(a), _n=(b); i >= _n; --i ) 77 #define REP( i, n ) FOR( i, 0, (n) ) 78 #define REPE( i, n ) FORE( i, 1, (n) ) 79 #define ALL( c ) (c).begin(), (c).end() 80 #define SORT( c ) sort( ALL( c ) ) 81 #define LEN( s ) (int)( (s).size() ) 82 #define CLS( a ) memset( (a),0,sizeof(a) ) 83 #define IOS ios::sync_with_stdio( false )  84 const double PI = 3.1415926535897932384626433832795028841971; 85 const double EPS = 1E-9; 86 const int64 INF64 = (int64)1E18; 87 const int INF = 1 << 30; 88 const int maxn = 1000005; 89  90  91 int eh[maxn], tot; 92 LL dist[maxn]; 93 LL dist2[maxn]; 94 bool vist[maxn]; 95 int n, m, k; 96  97  98 struct Edge { 99     int u, v;100     LL cost;101     int next;102     Edge() {};103     Edge( int a, LL b ) { v = a, cost = b; }104     Edge( int a, int b, LL c, int d ) : u(a), v(b), cost(c), next(d) {}105     bool operator < ( const Edge &x ) const {106         return cost > x.cost;107     }108 };109 110 111 Edge edge[maxn];112 int pre[maxn];113 114 115 int Dijk( int s )116 {117     for ( int i = 1; i <= n; ++i ) 118     {119         vist[i] = false, dist[i] = INF64;120         pre[i] = 1;121     }122     dist[s] = 0;123     priority_queue< Edge > que;124     que.push( Edge( s, 0 ) );125     int ans = 0;126     while ( !que.empty() )127     {128         Edge tmp = que.top();129         que.pop();130         int u = tmp.v, cost = tmp.cost;131         if ( vist[u] ) continue;132         ans += cost, vist[u] = true;133         for (  int i = eh[u]; i != -1; i = edge[i].next )134         {135             int v = edge[i].v;136             if ( !vist[v] && ( dist[v] > dist[u] + edge[i].cost || (pre[v] == s && dist[v] == dist[u] + edge[i].cost) ) )137             {138                 dist[v] = dist[u] + edge[i].cost;139                 pre[v] = u;140                 que.push( Edge( v, dist[v] ) );141             }142         }143 144 145     }146     return ans;147 }148 149 150 void addedge( int a, int b, LL c )151 {152     Edge e = Edge( a, b, c, eh[a] );153     edge[tot] = e;154     eh[a] = tot++;155 }156 157 158 void init()159 {160     tot = 0;161     memset( eh, -1, sizeof( eh ) );162 }163 164 165 int point[maxn];166 LL value[maxn];167 int pcnt;168 169 170 int Point1[maxn];171 int Point2[maxn];172 LL Value[maxn];173 LL MinVal[maxn];174 bool hascheck[maxn];175 176 177 void run()178 {179     //IOS;180     while( cin >> n >> m >> k ) 181     {182         pcnt = 1;183         init();184         for ( int i = 1; i <= m; ++i )185         {186             int u, v;187             LL w;188             cin >> u >> v >> w;189             addedge( u, v, w );190             addedge( v, u, w );191             Point1[i] = u;192             Point2[i] = v;193             Value[i] = w;194         }195 196 197         Dijk( 1 );198 199 200         for ( int i = 1; i <= n; ++i )201         {202             MinVal[i] = -1;203             hascheck[i] = 0;204             dist2[i] = dist[i];205         }206 207 208         init();209 210 211         for ( int i = 1; i <= m; ++i )212         {213             if ( Point1[0] == 89457 )214             {215                 if ( Point1[i] == 1 || Point2[i] == 1 || Value[i] == 0 )            216                 {217                     cout << Point1[i] << " " << Point2[i] << " " << Value[i] << endl;218                 }219             }220             addedge( Point1[i], Point2[i], Value[i] );221             addedge( Point2[i], Point1[i], Value[i] );222         }223 224 225         for ( int i = 1; i <= k; ++i )226         {227             int s;228             LL x;229             cin >> s >> x;230             if ( MinVal[s] == -1 )231                 MinVal[s] = x;232             else if ( MinVal[s] > x )233                 MinVal[s] = x;234 235 236             //if ( n == 100000 )237             //    cout << s << " " << x << endl;238             point[pcnt] = s;239             value[pcnt] = x;240             pcnt++;241 242 243             addedge( 1, s, x );244             addedge( s, 1, x );245 246 247             if ( x == 0 )248                 cout << s << " " << x << endl;249 250 251         }252 253 254         Dijk( 1 );255 256 257         int cnt = 0;258 259 260         for ( int i = 1; i <= k; ++i )261         {262             int s;263             LL x;264             s = point[i];265             x = value[i];266 267 268             if ( MinVal[s] != x || hascheck[s] == 1 )269             {270                 cnt++;271                 continue;272             }273 274 275             if ( dist2[s] <= x )276                 cnt++;277             else if ( dist[s] != x )278                 cnt++;279             else if ( pre[s] != 1 )280                 cnt++;281 282 283 284 285             hascheck[s] = 1;286         }287 288 289         cout << cnt << endl;290     }291 }292 293 294 int main( int argc, char **argv )295 {296 297 298 #ifndef ONLINE_JUDGE299        freopen("in.txt", "r", stdin);300 //      freopen("out.txt", "w", stdout);301 #endif302        time_t st = clock();303        run();304        Log("\n=============\n");305        Log("Time: %.2lf sec\n", (clock() - st) / double(CLOCKS_PER_SEC));306        return 0;307 }

 

CodeForces 449B - Jzzhu and Cities