[BZOJ]1922 [SDOI2010] 大陆争霸 最短路 + 堆优化

1922: [Sdoi2010]大陆争霸

Time Limit: 10 Sec   Memory Limit: 64 MB
Submit: 1978   Solved: 878
[ Submit][ Status][ Discuss]

Description

在一个遥远的世界里有两个国家：位于大陆西端的杰森国和位于大陆东端的 克里斯国。两个国家的人民分别信仰两个对立的神：杰森国信仰象征黑暗和毁灭 的神曾·布拉泽，而克里斯国信仰象征光明和永恒的神斯普林·布拉泽。 幻想历 8012年 1月，杰森国正式宣布曾·布拉泽是他们唯一信仰的神，同 时开始迫害在杰森国的信仰斯普林·布拉泽的克里斯国教徒。 幻想历 8012年 3月2日，位于杰森国东部小镇神谕镇的克里斯国教徒发动 起义。 幻想历 8012年 3月7日，神谕镇的起义被杰森国大军以残酷手段镇压。 幻想历 8012年 3月8日，克里斯国对杰森国宣战。由数十万大军组成的克 里斯军团开至两国边境，与杰森军团对峙。 幻想历 8012年 4月，克里斯军团攻破杰森军团防线进入神谕镇，该镇幸存 的克里斯国教徒得到解放。 战争随后进入胶着状态，旷日持久。战况惨烈，一时间枪林弹雨，硝烟弥漫， 民不聊生。 幻想历 8012年 5月12日深夜，斯普林·布拉泽降下神谕：“Trust me, earn eternal life.”克里斯军团士气大增。作为克里斯军团的主帅，你决定利用这一机 会发动奇袭，一举击败杰森国。具体地说，杰森国有 N 个城市，由 M条单向道 路连接。神谕镇是城市 1而杰森国的首都是城市 N。你只需摧毁位于杰森国首都 的曾·布拉泽大神殿，杰森国的信仰，军队还有一切就都会土崩瓦解，灰飞烟灭。 为了尽量减小己方的消耗，你决定使用自爆机器人完成这一任务。唯一的困 难是，杰森国的一部分城市有结界保护，不破坏掉结界就无法进入城市。而每个 城市的结界都是由分布在其他城市中的一些结界发生器维持的，如果想进入某个 城市，你就必须破坏掉维持这个城市结界的所有结界发生器。 现在你有无限多的自爆机器人，一旦进入了某个城市，自爆机器人可以瞬间 引爆，破坏一个目标（结界发生器，或是杰森国大神殿），当然机器人本身也会 一起被破坏。你需要知道：摧毁杰森国所需的最短时间。

Input

第一行两个正整数 N, M。 接下来 M行，每行三个正整数 ui, vi, wi，表示有一条从城市ui到城市 vi的单 向道路，自爆机器人通过这条道路需要 wi的时间。 之后 N 行，每行描述一个城市。首先是一个正整数 li，维持这个城市结界所 使用的结界发生器数目。之后li个1~N 之间的城市编号，表示每个结界发生器的 位置。如果 Li = 0，则说明该城市没有结界保护，保证L1 = 0 。

Output

仅包含一个正整数 ，击败杰森国所需的最短时间。

Sample Input

6 6
1 2 1
1 4 3
2 3 1
2 5 2
4 6 2
5 3 2
0
0
0
1 3
0
2 3 5

Sample Output

5

HINT

对于 20%的数据，满足 N≤15，M≤50； 
对于 50%的数据，满足 N≤500，M≤6,000； 
对于 100%的数据，满足 N≤3,000，M≤70,000，1≤wi≤108
。 
输入数据保证一定有解，且不会存在维持某个城市结界的结界发生器在这个
城市内部。 
连接两个城市的道路可能不止一条， 也可能存在一个城市自己到自己的道路。

Source

第一轮Day1

[ Submit][ Status][ Discuss]

HOME Back

  对基础算法有长见识了...这道题进入到某个城市是有条件的 -- 必须要先把相应的结界器给破坏掉. 这就有点拓扑的感觉， 一个城市如果能够到达了的话， 把这个城市破坏掉， 相应的这个城市保护的其他城市的保护就少了一个. 若某个城市对应的结界器全部破坏， 那么就可以到达了. 因为要求最短时间， 所以想到用最短路去求解.

设两个数组dis 和 act， 一个是可以到达的时间dis， 一个是允许进入的时间act(这个点结界器破坏完之时就是允许进入之时). 两者取max就是城市真实能到达的时间. 因为你能够到达某个城市， 但是可能他的结界还没有开放. 如果某个城市很早开放的话， 你也有可能需要花很长的时间才能到达那里.  那么max(dis[n], act[n])就是答案.

考虑怎么维护这两个数组.  我们用dijkstra跑最短路.

对于dis数组.在跑最短路的时候， 用取出的当前点的真实值(取max的那个值)去更新相连的其他点的dis值， 为什么只是dis值?因为其他店可能还没有开放， 所以你当前维护的只能是dis值， 如果你连的那个点结界被破坏掉了的话， 那么你就把他push到堆里， 不过push进去的距离值应该是那个点的真实值. 为什么不是dis值?因为比如说当前这个点更新其他点的话， 我们要得到其他点的dis值的话， 必须要保证你当前这个点的值是真实值才能去更新. 因为dis虽然来说只是可以到达的值， 但是你起码要保证可以到达啊， 所以你要用这个点到达的时间(即真实值)去更新相连的其他点， 才能保证其他点是真正可以到达的(虽然那些点可能还没有开放).

对于act数组. 你从堆里取出来的点相当于就是被爆破掉了， 只需要把当前点保护的其他城市所拥有的结界器--， 并更新他的允许进入时间. 如果那些点拥有结界器全没了， 就一样的把他push进堆里.

只能用Dijkstra， 因为SPFA的话会不断松弛不断入队， 当前取出的点的dis值并非尘埃落定， 以后可能还会被松弛. 那既然你当前点不是真实的， 这个点也就不能去更新其他点的dis值， act值. 要求真实. 对于类似Dijkstra这种性质的运用也可以去刘汝佳大神写的算法艺术与信息学竞赛里了解.

#include<stdio.h>
#include<cstring>
#include<queue>
#include<algorithm>
#define P pair<int, int>
#define fufil(a) memset(a, 60, sizeof(a))
using namespace std;
const int maxn = 3005;
priority_queue<P, vector<P>, greater<P> > q;
int n, m, num;
bool vis[maxn];
int h[maxn], a[maxn][maxn], dis[maxn], act[maxn], d[maxn], pro[maxn];
inline const int read(){
	register int x = 0;
	register char ch = getchar();
	while(ch < '0' || ch > '9') ch = getchar();
	while(ch >= '0' && ch <= '9') x = (x << 3) + (x << 1) + ch - '0', ch = getchar();
	return x;
}
struct edge{ int nxt, v, w;}e[70005];
inline void add(int u, int v, int w){ e[++num].v = v, e[num].nxt = h[u], e[num].w = w, h[u] = num;}

inline void Dijkstra(){
	fufil(dis);
	dis[1] = 0, q.push(P(0, 1));
	while(!q.empty()){
		int u = q.top().second, mx = q.top().first;
		q.pop();
		if(vis[u]) continue; vis[u] = true;
		for(int i = h[u]; i; i = e[i].nxt)
			if(dis[e[i].v] > mx + e[i].w){
				dis[e[i].v] = mx + e[i].w;
				if(!d[e[i].v]) q.push(P(max(dis[e[i].v], act[e[i].v]), e[i].v));
			}
		for(int i = 1; i <= pro[u]; ++i){
			int v = a[u][i];
			d[v]--, act[v] = max(act[v], mx);
			if(!d[v]) q.push(P(max(dis[v], act[v]), v));
		}
	}
}
int main(){
	n = read(), m = read();
	int u, v, w;
	for(int i = 1; i <= m; ++i){
		u = read(), v  = read(), w = read();
		if(u ^ v) add(u, v, w);
	}
	for(int i = 1; i <= n; ++i){
		d[i] = read();
		for(int j = 1; j <= d[i]; ++j) u = read(), a[u][++pro[u]] = i;
	}
	Dijkstra();
	printf("%d\n", max(dis[n], act[n]));
}