我一直都不會dij的堆優化,今天搞了一下。。。就是先弄一個優先隊列,存每個點的數據,然后這個題就加了一點不一樣的東西,每次的最短路算兩次,一次是自己的最短路,另一次是機關的最短路,兩者取最大值才是該點的真正的最短路。
dij堆優化鏈接
題干:
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
代碼:
#include<iostream> #include<cstdio> #include<cmath> #include<ctime> #include<queue> #include<algorithm> #include<cstring> #include<vector> using namespace std; #define duke(i,a,n) for(int i = a;i <= n;i++) #define lv(i,a,n) for(int i = a;i >= n;i--) #define clean(a) memset(a,0,sizeof(a)) const int INF = 1 << 30; typedef long long ll; typedef double db; template <class T> void read(T &x) {char c;bool op = 0;while(c = getchar(), c < '0' || c > '9')if(c == '-') op = 1;x = c - '0';while(c = getchar(), c >= '0' && c <= '9')x = x * 10 + c - '0';if(op) x = -x; } template <class T> void write(T x) {if(x < 0) putchar('-'), x = -x;if(x >= 10) write(x / 10);putchar('0' + x % 10); } struct node {int l,r,nxt,w;bool operator < (const node &other) const{return w < other.w;} } a[70005]; struct point {int x,dis;bool operator < (const point &other) const{return (other.dis < dis) || (x < other.x && dis == other.dis);} }; int len = 0,last[3005],m,n; int vis[3005]; int d1[3005],d2[3005],d[3005]; int st[3005][3005]; priority_queue <point> qu; void add(int x,int y,int w) {a[++len].l = x;a[len].r = y;a[len].w = w;a[len].nxt = last[x];last[x] = len; } void dij() {vis[1] = 0;d1[1] = 0;qu.push((point){1,max(d1[1],d2[1])});while(!qu.empty()){point g = qu.top();qu.pop();if(vis[g.x])continue;if(g.dis != max(d1[g.x],d2[g.x]))continue;for(int i = last[g.x]; i; i=a[i].nxt){int j = a[i].r;if(d1[j]>g.dis+a[i].w){d1[j] = g.dis + a[i].w;if(!d[j]) qu.push((point){j,max(d1[j],d2[j])});}}for(int i=1; i<=st[g.x][0]; ++i){d[st[g.x][i]]--;if(!d[st[g.x][i]]){d2[st[g.x][i]]=g.dis;qu.push((point){st[g.x][i],max(d1[st[g.x][i]],d2[st[g.x][i]])});}}vis[g.x] = 1; } } int main() {read(n);read(m);duke(i,1,m){int x,y,w;read(x);read(y);read(w);add(x,y,w);}memset(d1,127,sizeof(d1));duke(i,1,n){int p;read(d[i]);duke(j,1,d[i]){read(p);st[p][++st[p][0]] = i;}}dij();printf("%d\n",max(d1[n],d2[n]));return 0; } /* 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 */
?