CodeForces - 1253F Cheap Robot (多源Dijkstra+最小生成树)

题意是有一张无向连通图,边带权,每走过一条边需要消耗一定的电量。有k个关键点可以充满电,问从一个关键点走到另一个关键点,需要的电池容量至少是多少。

首先可以分析一下极端情况:

1.k=2:问题转化成求两点间最短路

2.k=n:问题转化成求最小生成树

所以要把两者综合起来考虑。

首先对所有关键点求一个多源最短路,得到每个点到最近的关键点的距离d,之后会发现:

1.假设走到某个点时的电量小于d,那么这个点就废掉了(永远不可能走到其他关键点)
2.假设电池容量为c,在不废掉的前提下,到达每个点时的电量必然是c-d(如果小于这个值,就可以走到最近的关键点充满电再走回来,前提是电量不小于d)

所以,对于每条边(u,v,w),点u可通过这条边走到点v的条件是c-d[u]-w>=d[v],即d[u]+d[v]+w<=c,这个性质关于u和v是对称的。(这个结论成立的前提是起点和终点都必须是关键点,否则就不可做了)

因此可以把每条边的边权改为(u,v,d[u]+d[v]+w),用kruskal算法求个全局最小生成树(按秩合并,无路径压缩),之后可以$O(logn)$询问。总复杂度$O(mlogm+qlogn)$

(ps:如果对最小生成树再来个树剖,这样复杂度就成了$O(mlogm+qloglogn)$)

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 const int N=3e5+10,mod=1e9+7;
 5 ll inf=0x3f3f3f3f3f3f3f3f;
 6 int n,m,k,Q,fa[N],hd[N],ne,mxd[N];
 7 ll d[N],val[N];
 8 struct E {int v,c,nxt;} e[N<<1];
 9 struct E2 {
10     int u,v;
11     ll c;
12     bool operator<(const E2& b)const {return c<b.c;}
13 } e2[N];
14 void link(int u,int v,int c) {e[ne]= {v,c,hd[u]},hd[u]=ne++;}
15 struct D {
16     int u;
17     ll g;
18     bool operator<(const D& b)const {return g>b.g;}
19 };
20 void Dij() {
21     priority_queue<D> q;
22     memset(d,inf,sizeof d);
23     for(int i=1; i<=k; ++i)q.push({i,d[i]=0});
24     while(q.size()) {
25         int u=q.top().u;
26         ll g=q.top().g;
27         q.pop();
28         if(d[u]!=g)continue;
29         for(int i=hd[u]; ~i; i=e[i].nxt) {
30             int v=e[i].v,c=e[i].c;
31             if(d[v]>d[u]+c)d[v]=d[u]+c,q.push({v,d[v]});
32         }
33     }
34 }
35 int fd(int x) {return fa[x]?fd(fa[x]):x;}
36 void mg(int x,int y,ll c) {
37     int fx=fd(x),fy=fd(y);
38     if(fx==fy)return;
39     if(mxd[fx]>mxd[fy])swap(fx,fy);
40     fa[fx]=fy,val[fx]=c,mxd[fy]=max(mxd[fy],mxd[fx]+1);
41 }
42 void Kruskal() {
43     sort(e2,e2+m);
44     for(int i=1; i<=n; ++i)mxd[i]=1,fa[i]=0;
45     for(int i=0; i<m; ++i) {
46         int u=e2[i].u,v=e2[i].v;
47         ll c=e2[i].c;
48         mg(u,v,c);
49     }
50 }
51 ll qry(int u,int v) {
52     ll ret=0;
53     for(; u!=v; u=fa[u]) {
54         if(mxd[u]>mxd[v])swap(u,v);
55         ret=max(ret,val[u]);
56     }
57     return ret;
58 }
59 int main() {
60     memset(hd,-1,sizeof hd),ne=0;
61     scanf("%d%d%d%d",&n,&m,&k,&Q);
62     for(int i=0; i<m; ++i) {
63         int u,v,c;
64         scanf("%d%d%d",&u,&v,&c);
65         link(u,v,c);
66         link(v,u,c);
67         e2[i]= {u,v,c};
68     }
69     Dij();
70     for(int i=0; i<m; ++i)e2[i].c+=d[e2[i].u]+d[e2[i].v];
71     Kruskal();
72     while(Q--) {
73         int u,v;
74         scanf("%d%d",&u,&v);
75         printf("%lld\n",qry(u,v));
76     }
77     return 0;
78 }

 

posted @ 2020-03-05 10:46  jrltx  阅读(360)  评论(0编辑  收藏  举报