POJ 2449 Remmarguts' Date

题目大意：

        好男人决不能错过和女生的约会。

        从前有一个UDF（联合三角洲自由王国）王国，一天王国的邻居国家的公主Uyuw被派往UDF完成一项外交使命，而在这之前Uyuw已经通知过UDF的王子Remmarguts她将会到UDF首都的大厅中与王子商讨商业事宜，但是商讨能开展的条件就是王子必须通过k短路径到达大厅（公主并不想见王子，故意为其设置障碍）。出于对公主的倾心以及对国家商贸发展的重视王子决定前往。

        UDF首都共有N个驿站（1 ≤ N ≤ 1,000）和M条泥泞小路（都是单向的，1 ≤ M ≤ 100,000），S为王子的出发地点，T是目标地点（即大厅），为了走k短路径同一个驿站可能会经过两次或者两次以上，包括S和T。

        现有多个测例（测例数无上限），每个测例中都会给出N和M，驿站按1-N编号，每条路都会给出起点、终点和路上所花的时间t（1 ≤ t ≤ 100），最后给出起点S、终点T和第k短路径的k（1 ≤ k ≤ 1,000），对于每个测例都需要打印k短路径的长度（即总时间），如果k短路径不存在则打印-1即可。

题目链接

注释代码：

/*                    * Problem ID : POJ 2449 Remmarguts' Date * Author     : Lirx.t.Una                    * Language   : G++        * Run Time   : 657 ms                    * Run Memory : 8808 KB                   */ #pragma G++ optimize("O2")#include <iostream>#include <cstring>#include <cstdio>#include <queue>#defineNIL0//infinity//这里指k短路径的最大时间//INF ≥ MAXSIDEWAYN × MAXT（单条路径的最大时间是100）#defineINF1000000000//maximum number of stations//站的最大数量#defineMAXSTATIONN1001//maximum number of sideways//有向路径的最大数量#defineMAXSIDEWAYN100001using namespace std;structArc {intu;intv;intt;//time，时间权值Arc(void) {}Arc( int xu, int xv, int xt ) :u(xu), v(xv), t(xt) {}};Arcarc[MAXSIDEWAYN];inte;inthead[MAXSTATIONN];intrhead[MAXSTATIONN];//reverse，逆图intnext[MAXSIDEWAYN];intrnext[MAXSIDEWAYN];intstk[MAXSTATIONN];boolin[MAXSTATIONN];inttop;intdp[MAXSTATIONN];//dp[i]，i点到目标点的最短时间structNode {intith;inttot;//total time，从起始到该点所经过的总时间intf;//f为估价函数，等于tot加上该店到目标点的单源最短路径Node(void) {}Node( int xith, int pst, int dt ) :ith(xith), tot( pst + dt ),f( tot + dp[xith] ) {}booloperator<(const Node &oth)const {return f > oth.f;}};voidaddarc( int u, int v, int t ) {arc[e]   = Arc( u, v, t );next[e] = head[u];head[u] = e;//初始化逆图rnext[e] = rhead[v];rhead[v] = e++;}voidinit( int n, int m ) {intu, v;intt;inti;memset( head + 1, NIL, n * sizeof(int));memset( rhead + 1, NIL, n * sizeof(int));memset( in + 1, false, n * sizeof(bool));for ( i = 1; i <= n; i++ )dp[i] = INF;e = 1;while ( m-- ) {scanf("%d%d%d", &u, &v, &t);addarc( u, v, t );}}boolrelax( int u, int v, int dt ) {//SPFA弹性判断intttmp;ttmp = dp[u] + dt;if ( ttmp < dp[v] ) {dp[v] = ttmp;return true;}return false;}voidspfa(int end) {//求各点到目标点的最短时间intu, v;inti;dp[end] = 0;top= 1;stk[1]= end;in[end]= true;while ( top ) {u   = stk[top--];in[u] = false;for ( i = rhead[u]; i; i = rnext[i] ) {v = arc[i].u;if ( relax( u, v, arc[i].t ) && !in[v] ) {stk[++top] = v;in[v]   = true;}}}}intastar( int bgn, int end, int k ) {//A*启发式搜索//第k次搜索到目标结点表示搜索成功//begin，end，k，起始结点，终止结点，和第k短路径priority_queue<Node>heap;Nodenode;//临时结点intpst;//past time，从bgn到父结点所经历的时间intu;//临时结点的标号inti;//计数变量if ( bgn == end )//自己到自己不计，因此需要减去//等价就是将k加1k++;heap.push( Node( bgn, 0, 0 ) );while ( !heap.empty() ) {node = heap.top();u = node.ith;pst  = node.tot;if ( u == end && !(--k) )//第k次搜索到终点时便得到答案//直接将k短路径的长度返回即可return pst;heap.pop();for ( i = head[u]; i; i = next[i] )heap.push( Node( arc[i].v, pst, arc[i].t ) );}return -1;//有可能除了最短路径之外，从其它点经过时会走到死路//比如      1//       2  3  4//     5    6    7//          8//5和7都是死路（起始是1，终点是8）}intmain() {intn, m;intbgn, end;intk;while ( ~scanf("%d%d", &n, &m) ) {init( n, m );scanf("%d%d%d", &bgn, &end, &k);spfa(end);if ( INF == dp[bgn] ) {//连最短路径都是无穷就意味着//bgn到end不可达puts("-1");continue;}printf("%d\n", astar( bgn, end, k ));}return 0;}

无注释代码：

#pragma G++ optimize("O2")#include <iostream>#include <cstring>#include <cstdio>#include <queue>#defineNIL0#defineINF1000000000#defineMAXSTATIONN1001#defineMAXSIDEWAYN100001using namespace std;structArc {intu;intv;intt;Arc(void) {}Arc( int xu, int xv, int xt ) :u(xu), v(xv), t(xt) {}};Arcarc[MAXSIDEWAYN];inte;inthead[MAXSTATIONN];intrhead[MAXSTATIONN];intnext[MAXSIDEWAYN];intrnext[MAXSIDEWAYN];intstk[MAXSTATIONN];boolin[MAXSTATIONN];inttop;intdp[MAXSTATIONN];structNode {intith;inttot;intf;Node(void) {}Node( int xith, int pst, int dt ) :ith(xith), tot( pst + dt ),f( tot + dp[xith] ) {}booloperator<(const Node &oth)const {return f > oth.f;}};voidaddarc( int u, int v, int t ) {arc[e]   = Arc( u, v, t );next[e] = head[u];head[u] = e;rnext[e] = rhead[v];rhead[v] = e++;}voidinit( int n, int m ) {intu, v;intt;inti;memset( head + 1, NIL, n * sizeof(int));memset( rhead + 1, NIL, n * sizeof(int));memset( in + 1, false, n * sizeof(bool));for ( i = 1; i <= n; i++ )dp[i] = INF;e = 1;while ( m-- ) {scanf("%d%d%d", &u, &v, &t);addarc( u, v, t );}}boolrelax( int u, int v, int dt ) {intttmp;ttmp = dp[u] + dt;if ( ttmp < dp[v] ) {dp[v] = ttmp;return true;}return false;}voidspfa(int end) {intu, v;inti;dp[end] = 0;top= 1;stk[1]= end;in[end]= true;while ( top ) {u   = stk[top--];in[u] = false;for ( i = rhead[u]; i; i = rnext[i] ) {v = arc[i].u;if ( relax( u, v, arc[i].t ) && !in[v] ) {stk[++top] = v;in[v]   = true;}}}}intastar( int bgn, int end, int k ) {priority_queue<Node>heap;Nodenode;intpst;intu;inti;if ( bgn == end )k++;heap.push( Node( bgn, 0, 0 ) );while ( !heap.empty() ) {node = heap.top();u = node.ith;pst  = node.tot;if ( u == end && !(--k) )return pst;heap.pop();for ( i = head[u]; i; i = next[i] )heap.push( Node( arc[i].v, pst, arc[i].t ) );}return -1;}intmain() {intn, m;intbgn, end;intk;while ( ~scanf("%d%d", &n, &m) ) {init( n, m );scanf("%d%d%d", &bgn, &end, &k);spfa(end);if ( INF == dp[bgn] ) {puts("-1");continue;}printf("%d\n", astar( bgn, end, k ));}return 0;}

单词解释：

date：n, 约会

appointment：n, 约定，预约

mandarin：n, 普通话，国语，官话

softly：adv, 温柔的

prince：n, 王子

delta：n, 三角洲

princess：n, 公主

diplomatic：adj, 外交的

mission：n, 使命，任务

erenow：adv, 在此之前

inform：vt, 通知

hall：n, 大厅

commercial：adj, 商业的

via：prep, 通过，经由

enamored：adj, 迷恋的，倾心的

prime minister：n, 首相

muddy：adj, 泥泞的

sideway：n, 人行道，小巷

directed：adj, 定向的

disparate：adj, 不同的，相异的

station：n, 车站，驻地

quote：n, 引用，引号