标签：des style blog http color io os ar for
传送门：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_JUDGE 299 freopen("in.txt", "r", stdin); 300 // freopen("out.txt", "w", stdout); 301 #endif 302 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
标签：des style blog http color io os ar for
原文地址：http://www.cnblogs.com/vinceliang/p/3991313.html
