CF100015C
主要找到环上任意一条边,有比较dis(u,v),dis(u,a)+w+dis(b,v),dis(u,b)+w+dis(a,u)
然后,然后没了
lca求dis(u,v);dis(u,v)=dis[u]+dis[v]-2*dis[lca(u,v)];
我加油慢慢写
咳咳,我可真是写了半天才写过了样例,lca经测试无误。
可是啊

哭了...
还是附上代码
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
struct node{
int u,v,w,nxt;
}e[1000101];
int fir[1000101],cnt=0,top[1000101];
void add(int x,int y,int z){
e[++cnt].nxt=fir[x];e[cnt].u=x;e[cnt].v=y;e[cnt].w=z;fir[x]=cnt;
}
int dis[1000101]={0};
int vis[1000101]={0},id[100101],tot=0,p=0,a=0,b=0,ww,siz[1000101],dep[1000101],son[1000101]={0},fa[1000101];
void dfs(int x,int f){
for(int i=fir[x];i;i=e[i].nxt){
int vv=e[i].v;//cout<<vv<<' '<<x<<' '<<i<<endl;
if(vv==f)continue;
if(vis[vv]==1&&a==0&&b==0){
p=i;a=e[i].u;b=e[i].v;ww=e[i].w;e[i].w=-1;return;
}
if(!vis[vv]){
vis[vv]=1; dfs(vv,x);
}
}
}
void dfs1(int x,int f){
siz[x]=1;dep[x]=dep[f]+1;fa[x]=f;
for(int i=fir[x];i;i=e[i].nxt){
int vv=e[i].v;//cout<<vv<<' '<<x<<' '<<i<<endl;
if((a==x&&vv==b)||(vv==a&&x==b))e[i].w=-1;
if(vv==f||(a==x&&vv==b)||(vv==a&&x==b))continue;
if(!vis[vv]&&e[i].w!=-1){
vis[vv]=1;dis[vv]=dis[x]+e[i].w;dfs1(vv,x);
}
siz[x]+=siz[vv];
if(siz[vv]>siz[son[x]])son[x]=vv;
}
}
void dfs2(int x,int topp){
top[x]=topp;
if(son[x])dfs2(son[x],topp);
for(int i=fir[x];i;i=e[i].nxt){
int vv=e[i].v;
if(vv==fa[x]||vv==son[x]||(a==x&&vv==b)||(vv==a&&x==b))continue;
dfs2(vv,vv);
}
}
int lca(int x,int y){
while(top[x]!=top[y]){
if(dep[top[x]]>=dep[top[y]])x=fa[top[x]];
else y=fa[top[y]];
}
return dep[x]<dep[y]?x:y;
}
int diss(int x,int y,int zz){
return dis[x]+dis[y]-2*dis[zz];
}
int n,m;
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++){
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);add(y,x,z);
}
dfs(0,0);
//cout<<a<<' '<<b<<endl;
memset(vis,0,sizeof(vis));siz[0]=0;dep[0]=0;
dfs1(0,0);
dfs2(0,0);
int q;
cin>>q;
//for(int i=1;i<=n;i++)cout<<dis[i]<<' ';
while(q--){
int x,y;
scanf("%d%d",&x,&y);
int zz=lca(x,y),tt=lca(x,a),kk=lca(b,y);
int ans=min(diss(x,y,zz),diss(x,a,tt)+ww+diss(b,y,kk));
tt=lca(x,b);kk=lca(a,y);
ans=min(ans,diss(x,b,tt)+ww+diss(a,y,kk));
printf("%d\n",ans);
}
int lla;
cin>>lla;
return 0;
}
我觉得,照理,树剖的lca是跑得草鸡快的,至少比神马倍增跑得快
可是啊...不想说了
CF100015C的更多相关文章
随机推荐
- Python time.md
time模块 Comparing Clocks time.clock():在Unix 上,返回当前的处理器时间,以浮点数秒数表示. time.monotonic():返回一个单调时钟的值(在分秒内), ...
- 谷歌希望让 Swift 成为安卓的优先选择,以取代由 Oracle 开发的 Java 程序语言。
http://news.coolban.com/Web/Index/land/app/2/id/405239
- ResultJsonInfo<T>
using System; using System.Collections.Generic; using System.Linq; using System.Web; namespace QY.We ...
- 最长公共子序列&最长公共子串
首先区别最长公共子串和最长公共子序列 LCS(计算机科学算法:最长公共子序列)_百度百科 最长公共子串,这个子串要求在原字符串中是连续的.而最长公共子序列则并不要求连续. 最长公共子序列: http ...
- ethereumjs/ethereumjs-common-3-test
查看test能够让你更好滴了解其API文档的使用 ethereumjs-common/tests/chains.js const tape = require('tape') const Common ...
- js之radio应用实例
radio和checkbox还有select,可谓是前后端常用三剑客啊!特别是checkbox和select,关于这两个今天不讲,因为在下面这几篇文章,我已经比较详细的讲解了. SpringMVC之a ...
- Hadoop源码分析(mapreduce.lib.partition/reduce/output)
Map的结果,会通过partition分发到Reducer上.Reducer做完Reduce操作后,通过OutputFormat,进行输出.以下我们就来分析參与这个过程的类. Mapper的结果, ...
- springboot mybatis 后台框架平台 shiro 权限 集成代码生成器
1.代码生成器: [正反双向](单表.主表.明细表.树形表,快速开发利器)+快速表单构建器freemaker模版技术 ,0个代码不用写,生成完整的一个模块,带页面.建表sql脚本.处理类.servic ...
- scala_类的继承
Scala继承一个基类跟Java很相似, 但我们需要注意以下几点: 重写一个非抽象方法必须使用override修饰符,以及重写父类属性也必须使用override修饰符. 只有主构造函数才可以往基类的构 ...
- Scala中class、object、case class、case object区别
/** class.object.case class.case object区别 * * class 类似Java中的class: * object Scala不能定义静态成员,用定义单例对象代之: ...