树是指由n个点,n-1条边构成的联通无向图。如果有一棵树,它的每一条边(u,v)都有一个权值l(u,v),我们把这样的树称作带权树。
我们知道对于树上的任意两个点,他们之间的路径是唯一的。对于两个点u,v来说,我们可以算出u与v之间的路径上的所有边权之和,将其称作u与v之间路径的长度,记作d(u,v)。
你的任务是计算:∑d(u,v) (u!=v) 【输入文件】
输入文件的第一行为一个正整数n。
接下来的n-1行,每行三个非负整数x,y,w表示点x与点y之间有一条权值为w的一条边。
保证输入的图示一棵树。
【输出文件】
输出文件仅包含一行一个数,即答案。因为结果可能很大,请将答案模100 000 007输出
【输入样例】
4 1 2 4 1 3 4 1 4 4
【输出样例】
72
【样例解释】
d(1,2)+d(1,3)+d(1,4)=4+4+4=12 d(2,1)+d(2,3)+d(2,4)=4+8+8=20 d(3,1)+d(3,2)+d(3,4)=4+8+8=20 d(4,1)+d(4,2)+d(4,3)=4+8+8=20 ans=12+20+20+20=72 【数据范围】
n≤3×10^5 0≤w<1000000007

【分析】:

看n的范围肯定不能枚举……那么我们就换个角度,考虑每条边对答案的贡献!

首先随便选一个点作根,dfs出每个子树的size大小与deep深度,那么可以发现对于deep[v]=deep[u]+1的两点,他们之间的边l(u,v)在答案中会被选size[v]*(n-size[v])次,枚举每条边看看即可。最优因为点对是无序的,但是我们只是算了从u那一边选v,所以我们还要×2算从v那一边选u。

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <string>
#include <cmath>
#include <ctime>
#include <algorithm>
#define ll long long
#define N 1000001
#define mod 100000007
using namespace std;
ll n,ans,root,cnt,first[N],size[N],depth[N];
bool vis[N];
struct edge{int next,from,to;ll len;}e[N];
int read()
{
int x=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
while(ch>='0'&&ch<='9') {x=x*10+ch-'0';ch=getchar();}
return x*f;
}
void add(int x,int y,int w)
{
e[++cnt].next=first[x];
e[cnt].from=x;
e[cnt].to=y;
e[cnt].len=w;
first[x]=cnt;
}
void dfs(int now)
{
size[now]=1,vis[now]=true;
for(int i=first[now];i;i=e[i].next)
{
int next=e[i].to;
if(!vis[next])
{
depth[next]=depth[now]+1;
dfs(next);
size[now]+=size[next];
}
}
}
int main()
{
int sizee = 50 << 20; //开50MB系统栈
char *p = (char*)malloc(sizee) + sizee;
__asm__("movl %0, %%esp\n" :: "r"(p));
//freopen("tree1.in","r",stdin);
//freopen("tree.out","w",stdout);
n=read();
for(int i=1;i<n;i++)
{
int x=read(),y=read(),w=read();
add(x,y,w),add(y,x,w);
}
dfs(1);
for(int i=1;i<=cnt;i++)
{
int x=e[i].from,y=e[i].to;ll w=e[i].len;
if(depth[y]==depth[x]+1)
{
ans+=(ll)(size[y]%mod*(n-size[y])%mod*(ll)w%mod);
ans%=mod;
}
}
cout<<ans*2ll%mod<<endl;
return 0;
}

概念

树的dfs遍历就是对于一个树上每个点 root,它向下的多个分支,选择一个分支一直走下去,直至走完并**回溯**到root再走其他分支形成的遍历。

Code
```
void dfs(int u,int fa)
{
vis[u]=1;
//邻接表枚举i的每个相邻节点
for(int i=link[u]; i; i=e[i].next)
{
int v = e[i].u;
if(v!=fa)
dfs(v,u);
}
}
```

求树的深度

每个节点x的深度用deep[x]表示。
代码只要在dfs向下搜索前deep[e[i].y]=deep[x]+1;

求子树大小

以x节点为根节点的子树大小size[x];
在dfs向下搜索之后size[x]+=size[e[i].y];

树的重心

对于一个节点x,如果把它从树中删除,原来的一棵树可能会分成若干个不相连的部分,每部分都是一颗子树。
设 Maxp(x)表示在删除节点x后分成的子树中,包含节点最多的那颗子树的节点数。
对于所有点来说,如果使得Maxp(x)最小的节点x就成为整颗树的重心。

int pos;//记录重心的编号
void dfs(int x,int fa)
{
v[x]=1;
sz[x]=1;
int Maxp=0;
//邻接表枚举i的每个相邻节点
for(int i=link[x]; i; i=e[i].next)
{
int y=e[i].y;
if(y!=fa)
{
dfs(y,x);
sz[x]+=sz[y];
Maxp=max(Maxp,sz[y]);
}
}
Maxp=max(Maxp,n-sz[x]);
if(Maxp<ans)
{
ans=Maxp;
pos=x;
}
}

树的DFS序就是在对树进行DFS的时候,对树的节点进行重新编号;
DFS序有一个很强的性质: 一颗子树的所有节点在DFS序内是连续的一段, 利用这个性质我们可以解决很多问题。
```
void DFS(int u, int fa)
{
L[u] = ++dfs_clock;
for(int k = head[u]; ~k; k = E[k].next)
{
int v = E[k].to;
if(v == fa) continue;
DFS(v, u);
}
R[u] = dfs_clock;
}
```

树上各种DFS姿势算法笔记的更多相关文章

  1. 【bzoj4940】[Ynoi2016]这是我自己的发明 DFS序+树上倍增+莫队算法

    题目描述 给一个树,n 个点,有点权,初始根是 1. m 个操作,每次操作: 1. 将树根换为 x. 2. 给出两个点 x,y,从 x 的子树中选每一个点,y 的子树中选每一个点,如果两个点点权相等, ...

  2. 算法笔记--sg函数详解及其模板

    算法笔记 参考资料:https://wenku.baidu.com/view/25540742a8956bec0975e3a8.html sg函数大神详解:http://blog.csdn.net/l ...

  3. 算法笔记--数位dp

    算法笔记 这个博客写的不错:http://blog.csdn.net/wust_zzwh/article/details/52100392 数位dp的精髓是不同情况下sta变量的设置. 模板: ]; ...

  4. 算法笔记--lca倍增算法

    算法笔记 模板: vector<int>g[N]; vector<int>edge[N]; ][N]; int deep[N]; int h[N]; void dfs(int ...

  5. 算法笔记_023:拓扑排序(Java)

    目录 1 问题描述 2 解决方案 2.1 基于减治法实现 2.2 基于深度优先查找实现 1 问题描述 给定一个有向图,求取此图的拓扑排序序列. 那么,何为拓扑排序? 定义:将有向图中的顶点以线性方式进 ...

  6. 学习Java 以及对几大基本排序算法(对算法笔记书的研究)的一些学习总结(Java对算法的实现持续更新中)

    Java排序一,冒泡排序! 刚刚开始学习Java,但是比较有兴趣研究算法.最近看了一本算法笔记,刚开始只是打算随便看看,但是发现这本书非常不错,尤其是对排序算法,以及哈希函数的一些解释,让我非常的感兴 ...

  7. BZOJ_4238_电压_树上差分+dfs树

    BZOJ_4238_电压_树上差分+dfs树 Description 你知道Just Odd Inventions社吗?这个公司的业务是“只不过是奇妙的发明(Just Odd Inventions)” ...

  8. 算法笔记--STL中的各种遍历及查找(待增)

    算法笔记 map: map<string,int> m; map<string,int>::iterator it;//auto it it = m.begin(); whil ...

  9. 算法笔记--priority_queue

    算法笔记 priority_queue<int>que;//默认大顶堆 或者写作:priority_queue<int,vector<int>,less<int&g ...

随机推荐

  1. Appium与python自动测试环境及demo详解

    App--UI自动化这种高端的名词已经被越来越多的人所高呼,可是从实际角度来讲,个人觉得还是有点鸡肋,不如接口自动化敏捷度高,工作量 也是接口自动化的好几倍.但是,[划重点了]  在技术时代中,作为测 ...

  2. ASP.NET Core 2.1 源码学习之 Options[2]:IOptions 【转】

    原文链接:https://www.cnblogs.com/RainingNight/p/strongly-typed-options-ioptions-in-asp-net-core.html 在 上 ...

  3. DataSource的设置

    1.Centos和redhat,Fedora等版本无须设置,直接在cloud.cfg指定,默认是EC2 datasource_list: ['ConfigDrive','OpenStack'] 2.u ...

  4. [ecmagent][redis学习][1初识redis] redis安装+redis快速教程+python操作redis

    # redis安装 # redis安装教程 -- 服务器(ubuntu)安装redis服务 sudo apt-get install redis-server -- 源码安装 -- $ wget ht ...

  5. A - 最长上升子序列

    A - 最长上升子序列 Time Limit: 1000/1000MS (C++/Others) Memory Limit: 65536/65536KB (C++/Others) Problem De ...

  6. hexo站点地图

    将网站链接提交到百度 百度搜索引擎提交入口 有三种验证方式,我选择Html标签验证,在themes\next\layout\_partials\head.swing中添加验证代码: <meta ...

  7. NO JSP Support for /, did not find org.eclipse.jetty.jsp.JettyJspServlet

    springMVC 内嵌jetty时,出现了NO JSP Support for /, did not find org.eclipse.jetty.jsp.JettyJspServlet: 添加了 ...

  8. axis2实践(二)Restful入门示例

    1. 实例说明 本示例直接参照了RESTful Web Services with Apache Axis2,本示例基本就是沿用的原示例,就是一个对学生信息(包括姓名,年龄,课程)的管理的例子,提供如 ...

  9. 关于Yarn源码那些事-前传之ResourceManager篇(一)初始化

    在关于Yarn那些事的博客里,介绍的主要是针对任务提交的一个动态流程说明,而其中牵涉到的一些细节问题,必须通过Resourcemanager的启动和NodeManager的启动,来更好的说明. 而本系 ...

  10. [bzoj] 3343 教主的魔法 || 带修改分块

    原题 长度为n的序列,有两种操作: 1.[l,r]区间每个数+w 2.询问[l,r]区间有多少个数>c 记录lazy数组即可. #include<cstdio> #include&l ...