Codeforces 1253 F. Cheap Robot —— 想法,dijkstra+kruskal重构树+并查集+树链剖分+lca,有丶东西

本文主要是介绍Codeforces 1253 F. Cheap Robot —— 想法,dijkstra+kruskal重构树+并查集+树链剖分+lca,有丶东西,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

This way

题意:

给你一张大小为n的图,并且点1~k是充电站。有q个询问,每次询问你从a走到b最少需要的电池大小为多少。
假设你当前电量为c,当走过一个权值为w的边时,你的点会变成c-w,如果这个点是充电站,那么你的电量会便会变回电池容量。
问你你每次需要的电池容量最小是多少。

题解:

有一说一,这道题很厉害。
首先先用dijkstra找到每个点距离它最近的充电站的距离,然后的话我们可以知道,如果一条边要用到,当且仅当 d i s [ x ] + d i s [ y ] + e [ i ] . w ≤ c dis[x]+dis[y]+e[i].w\le c dis[x]+dis[y]+e[i].wc,因为它可以去最近的充电站再往终点走。那么我们要求的就是将每条边的代价算出来,然后从a到b就是这个路径上的边权的最大值。那么我们要使最大值最小,就用kruskal就行了。
之后我们建出最小生成树之后将边权存到儿子的点权上,然后查询的话就查a~b的路径上的最大值即可。注意lca的点权不能算进去。
代码非常长,好像还有别的做法更容易点。时间复杂度一般。

在这里插入图片描述

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=3e5+5;
const ll inf=1e18;
struct edge{int to,next;ll w;
}e[N*2];
int cnt,head[N];
void add(int x,int y,ll w){e[cnt].to=y;e[cnt].w=w;e[cnt].next=head[x];head[x]=cnt++;
}
struct node{int id;ll w;bool operator< (const node& a)const {return w>a.w;}
};
struct Kruskal{int x,y;ll w;bool operator< (const Kruskal& a)const {return w<a.w;}
}edg[N];
priority_queue<node>Q;
ll dis[N],mx[N*4];
void dij(){while(!Q.empty()){node u=Q.top();Q.pop();if(u.w!=dis[u.id])continue;for(int i=head[u.id];~i;i=e[i].next){int ne=e[i].to;if(dis[ne]>dis[u.id]+e[i].w){dis[ne]=dis[u.id]+e[i].w;Q.push({ne,dis[ne]});}}}
}
int fa[N];
int finds(int x){return x==fa[x]?fa[x]:fa[x]=finds(fa[x]);}
int f[N][25],dep[N],siz[N],son[N],tim,pos[N],top[N];
int n,m,k,q;
void dfs1(int u)
{siz[u]=1;for (int i=head[u];~i;i=e[i].next){int v=e[i].to;if (v!=f[u][0]){f[v][0]=u;dep[v]=dep[u]+1;dfs1(v);siz[u]+=siz[v];if (son[u]==0||siz[v]>siz[son[u]])son[u]=v;}}
}
void dfs2(int u,int ff)
{tim++;top[u]=ff;pos[u]=tim;if (son[u]) dfs2(son[u],ff);for (int i=head[u];~i;i=e[i].next){int v=e[i].to;if(v!=f[u][0]&&v!=son[u])dfs2(v,v);}
}void deal(){for(int j=1;(1<<j)<=n;j++)for(int i=1;i<=n;i++)f[i][j]=f[f[i][j-1]][j-1];
}
int lca(int x,int y){if(dep[x]<dep[y])swap(x,y);for(int i=20;i>=0;i--)if(dep[f[x][i]]>=dep[y])x=f[x][i];if(x==y)return x;for(int i=20;i>=0;i--)if(f[x][i]!=f[y][i])x=f[x][i],y=f[y][i];return f[x][0];
}void add(int l,int r,int root,int p,ll val)
{if(l==r){mx[root]=val;return ;}int mid=l+r>>1;if(mid>=p)add(l,mid,root<<1,p,val);elseadd(mid+1,r,root<<1|1,p,val);mx[root]=max(mx[root<<1],mx[root<<1|1]);
}
ll ask(int l,int r,int root,int ql,int qr)
{if(l>=ql&&r<=qr)return mx[root];int mid=l+r>>1;ll ans=0;if(mid>=ql)ans=ask(l,mid,root<<1,ql,qr);if(mid<qr)ans=max(ask(mid+1,r,root<<1|1,ql,qr),ans);return ans;
}
ll query(int x,int y)
{ll ans=0;while (top[x]!=top[y]){if (dep[top[x]]<dep[top[y]]) swap(x,y);ans=max(ans,ask(1,n,1,pos[top[x]],pos[x]));x=f[top[x]][0];}if (pos[x]>pos[y]) swap(x,y);ans=max(ans,ask(1,n,1,pos[x],pos[y]));return ans;
}
int in[N];
int main()
{memset(head,-1,sizeof(head));scanf("%d%d%d%d",&n,&m,&k,&q);int x,y;ll w;for(int i=1;i<=m;i++)scanf("%d%d%lld",&edg[i].x,&edg[i].y,&edg[i].w),add(edg[i].x,edg[i].y,edg[i].w),add(edg[i].y,edg[i].x,edg[i].w);for(int i=1;i<=n;i++)dis[i]=inf,fa[i]=i;for(int i=1;i<=k;i++)Q.push({i,0}),dis[i]=0;dij();for(int i=1;i<=m;i++)edg[i].w+=dis[edg[i].x]+dis[edg[i].y];sort(edg+1,edg+1+m);cnt=0;memset(head,-1,sizeof(head));for(int i=1;i<=m;i++){int fax=finds(edg[i].x),fay=finds(edg[i].y);if(fax==fay)continue;add(edg[i].x,edg[i].y,edg[i].w),add(edg[i].y,edg[i].x,edg[i].w);fa[fay]=fax;in[i]=1;}dep[1]=1;dfs1(1);dfs2(1,1);deal();for(int i=1;i<=m;i++){if(!in[i])continue;x=edg[i].x,y=edg[i].y;if(dep[x]>dep[y])add(1,n,1,pos[x],edg[i].w);elseadd(1,n,1,pos[y],edg[i].w);}while(q--){int x,y;scanf("%d%d",&x,&y);int l=lca(x,y),firstson;ll ans=0;if(y!=l){firstson=y;for(int i=20;i>=0;i--)if(dep[f[firstson][i]]>dep[l])firstson=f[firstson][i];ans=query(firstson,y);}if(x!=l){firstson=x;for(int i=20;i>=0;i--)if(dep[f[firstson][i]]>dep[l])firstson=f[firstson][i];ans=max(ans,query(firstson,x));}printf("%lld\n",ans);}return 0;
}

这篇关于Codeforces 1253 F. Cheap Robot —— 想法,dijkstra+kruskal重构树+并查集+树链剖分+lca,有丶东西的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/881138

相关文章

poj1330(LCA最近公共祖先)

题意:求最近公共祖先 思路:之前学习了树链剖分,然后我就用树链剖分的一小部分知识就可以解这个题目了,记录每个结点的fa和depth。然后查找时,每次将depth大的结点往上走直到x = y。 代码如下: #include<iostream>#include<algorithm>#include<stdio.h>#include<math.h>#include<cstring>

poj 1502 MPI Maelstrom(单源最短路dijkstra)

题目真是长得头疼,好多生词,给跪。 没啥好说的,英语大水逼。 借助字典尝试翻译了一下,水逼直译求不喷 Description: BIT他们的超级计算机最近交货了。(定语秀了一堆词汇那就省略吧再见) Valentine McKee的研究顾问Jack Swigert,要她来测试一下这个系统。 Valentine告诉Swigert:“因为阿波罗是一个分布式共享内存的机器,所以它的内存访问

uva 10801(乘电梯dijkstra)

题意: 给几个电梯,电梯0 ~ n-1分别可以到达很多层楼。 换乘电梯需要60s时间。 问从0层到target层最小的时间。 解析: 将进入第0层的电梯60s也算上,最后减。 坑点是如果target为0输出0。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algori

hdu 3790 (单源最短路dijkstra)

题意: 每条边都有长度d 和花费p,给你起点s 终点t,要求输出起点到终点的最短距离及其花费,如果最短距离有多条路线,则输出花费最少的。 解析: 考察对dijkstra的理解。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstrin

poj 3255 次短路(第k短路) A* + spfa 或 dijkstra

题意: 给一张无向图,求从1到n的次短路。 解析: A* + spfa 或者 dijkstra。 详解见上一题:http://blog.csdn.net/u013508213/article/details/46400189 本题,spfa中,stack超时,queue的效率最高,priority_queue次之。 代码: #include <iostream>#i

poj 1182 并查集 食物链类

题意: 有n只动物,分别编号1....n。所有动物都属于A,B,C中的一种,已知A吃B,B吃C,C吃A。 按顺序给出下面两种共K条信息: 1. x 和 y 属于同一类。 2. x 吃 y 。 然而这些信息可能会出错,有可能有的信息和之前给出的信息矛盾,也有的信息可能给出的 x 和 y 不在n的范围内。 求k条信息中有多少条是不正确的。 解析: 对于每只动物,创建3个元素 i

Codeforces Round #240 (Div. 2) E分治算法探究1

Codeforces Round #240 (Div. 2) E  http://codeforces.com/contest/415/problem/E 2^n个数,每次操作将其分成2^q份,对于每一份内部的数进行翻转(逆序),每次操作完后输出操作后新序列的逆序对数。 图一:  划分子问题。 图二: 分而治之,=>  合并 。 图三: 回溯:

Codeforces Round #261 (Div. 2)小记

A  XX注意最后输出满足条件,我也不知道为什么写的这么长。 #define X first#define Y secondvector<pair<int , int> > a ;int can(pair<int , int> c){return -1000 <= c.X && c.X <= 1000&& -1000 <= c.Y && c.Y <= 1000 ;}int m

Codeforces Beta Round #47 C凸包 (最终写法)

题意慢慢看。 typedef long long LL ;int cmp(double x){if(fabs(x) < 1e-8) return 0 ;return x > 0 ? 1 : -1 ;}struct point{double x , y ;point(){}point(double _x , double _y):x(_x) , y(_y){}point op

Codeforces Round #113 (Div. 2) B 判断多边形是否在凸包内

题目点击打开链接 凸多边形A, 多边形B, 判断B是否严格在A内。  注意AB有重点 。  将A,B上的点合在一起求凸包,如果凸包上的点是B的某个点,则B肯定不在A内。 或者说B上的某点在凸包的边上则也说明B不严格在A里面。 这个处理有个巧妙的方法,只需在求凸包的时候, <=  改成< 也就是说凸包一条边上的所有点都重复点都记录在凸包里面了。 另外不能去重点。 int