[poj 1741]Tree 点分治
题意
求树上距离不超过k的点对数,边权<=1000
题解
点分治。
点分治的思想就是取一个树的重心,这种路径只有两种情况,就是经过和不经过这个重心,如果不经过重心就把树剖开递归处理,经过就把两边的点瞎那啥统计一下,因为会有完全在子树内的路径,还要容斥算算。
点分治是O(logn)的,但是每次操作是O(nlogn)那么总时间就是
#include<map>
#include<stack>
#include<queue>
#include<cstdio>
#include<string>
#include<vector>
#include<cstring>
#include<complex>
#include<iostream>
#include<assert.h>
#include<algorithm>
using namespace std;
#define inf 1001001001
#define infll 1001001001001001001LL
#define ll long long
#define dbg(vari) cerr<<#vari<<" = "<<(vari)<<endl
#define gmax(a,b) (a)=max((a),(b))
#define gmin(a,b) (a)=min((a),(b))
#define Ri register int
#define gc getchar()
#define il inline
il int read(){
bool f=true;Ri x=0;char ch;while(!isdigit(ch=gc))if(ch=='-')f=false;while(isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=gc;}return f?x:-x;
}
#define gi read()
#define FO(x) freopen(#x".in","r",stdin),freopen(#x".out","w",stdout);
struct edge{
int to,next,v;
} e[23333];
int last[20000],cnt;
void link(int a,int b,int c){
e[++cnt]=(edge){b,last[a],c};last[a]=cnt;
e[++cnt]=(edge){a,last[b],c};last[b]=cnt;
}
int siz[23333],f[23333],vis[23333],deep[23333],d[23333],root,sum,ans,n,k;
void zy(int x,int fa=-1){
siz[x]=1;f[x]=0;
for(int i=last[x];i;i=e[i].next){
if(e[i].to!=fa&&!vis[e[i].to]){
zy(e[i].to,x);
siz[x]=siz[x]+siz[e[i].to];
f[x]=max(f[x],siz[e[i].to]);
}
}
f[x]=max(f[x],sum-siz[x]);
if(f[x]<f[root])root=x;
}
void dfs(int x,int fa=-1){
deep[++deep[0]]=d[x];
for(int i=last[x];i;i=e[i].next){
if(e[i].to!=fa&&!vis[e[i].to]){
d[e[i].to]=d[x]+e[i].v;
dfs(e[i].to,x);
}
}
}
int js(int x,int now){
d[x]=now;deep[0]=0;
dfs(x);
sort(deep+1,deep+deep[0]+1);
int zzyy=0,l,r;
for(int l=1,r=deep[0];l<r;){
if(deep[l]+deep[r]<=k)zzyy=zzyy+r-l,l++;
else --r; }//f**k;
return zzyy;
}
void dfz(int x){
ans+=js(x,0);
vis[x]=1;
for(int i=last[x];i;i=e[i].next){
if(!vis[e[i].to]){
ans-=js(e[i].to,e[i].v);
sum=siz[e[i].to];
root=0;
zy(e[i].to,root);
dfz(root);
}
}
}
int main(){
while(n=gi,k=gi,n&&k){
ans=0,cnt=0,root=0;
memset(vis,0,sizeof(vis));
memset(last,0,sizeof(last));
for(int i=1;i<n;i++){
int a,b,c;
a=gi;b=gi;c=gi;
link(a,b,c);
}
sum=n;f[0]=inf;
zy(1);
dfz(root);
printf("%d\n",ans);
}
return 0;
}
[poj 1741]Tree 点分治的更多相关文章
- POJ 1741.Tree 树分治 树形dp 树上点对
Tree Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 24258 Accepted: 8062 Description ...
- POJ 1741 Tree 树分治
Tree Description Give a tree with n vertices,each edge has a length(positive integer less than 1 ...
- [bzoj 1468][poj 1741]Tree [点分治]
Description Give a tree with n vertices,each edge has a length(positive integer less than 1001). Def ...
- POJ 1741 Tree(点分治点对<=k)
Description Give a tree with n vertices,each edge has a length(positive integer less than 1001). Def ...
- POJ 1741 Tree ——点分治
[题目分析] 这貌似是做过第三道以Tree命名的题目了. 听说树分治的代码都很长,一直吓得不敢写,有生之年终于切掉这题. 点分治模板题目.自己YY了好久才写出来. 然后1A了,开心o(* ̄▽ ̄*)ブ ...
- POJ - 1741 - Tree - 点分治 模板
POJ-1741 题意: 对于带权的一棵树,求树中距离不超过k的点的对数. 思路: 点分治的裸题. 将这棵树分成很多小的树,分治求解. #include <algorithm> #incl ...
- poj 1741 Tree(树的点分治)
poj 1741 Tree(树的点分治) 给出一个n个结点的树和一个整数k,问有多少个距离不超过k的点对. 首先对于一个树中的点对,要么经过根结点,要么不经过.所以我们可以把经过根节点的符合点对统计出 ...
- POJ 1741.Tree and 洛谷 P4178 Tree-树分治(点分治,容斥版) +二分 模板题-区间点对最短距离<=K的点对数量
POJ 1741. Tree Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 34141 Accepted: 11420 ...
- POJ 1741 Tree 求树上路径小于k的点对个数)
POJ 174 ...
随机推荐
- ThinkPHP实现联动菜单;
联动菜单,首先给你看看前端是怎么写的:
- web app 开发
去除手机浏览器标签默认高亮边框 -webkit-tap-highlight-color 属性 属性描述:这个属性可以指设置透明度.如果未设置透明度,iOS上的Safari会给予颜色一个默认的透明度.把 ...
- 用Python作GIS之五:从示例入手—example函数
进入STARS后,最简单的学习方法就是演示示例数据.对于源码的分析也可以从这里入手. 以下为出发菜单项“Example Project”的函数example:def example(se ...
- 应聘复习基础笔记1:网络编程之TCP与UDP的优缺点,TCP三次握手、四次挥手、传输窗口控制、存在问题
重要性:必考 一.TCP与UDP的优缺点 ①TCP---传输控制协议,提供的是面向连接.可靠的字节流服务.当客户和服务器彼此交换数据前,必须先在双方之间建立一个TCP连接,之后才能传输数据.TCP提供 ...
- hdu 5280 Senior's Array
题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=5280 Senior's Array Description One day, Xuejiejie ge ...
- iOS学习之Object-C语言内存管理
一.内存管理的方式 1.iOS应用程序出现Crash(闪退),90%的原因是因为内存问题. 2.内存问题: 1)野指针异常:访问没有所有权的内存,如果想要安全的访问,必须 ...
- 如何在Window上使用Git
开始的时候同事只给了一个地址,类似这样:git@111.111.1.1:ABCDEF (1)如何在Windows上使用Git 有一篇博客不错:http://www.tuicool.com/articl ...
- 有关ZxMiddleTier构想
终于决定动手写一个自己的东西,没有其它想法,人已经30多了,感觉原来学的东西都是零散的,想通过这样一个中间层的项目,串联起原来所学的东西,形成一个体系,也算是对自己这么多年的程序员的生涯做出一个交代, ...
- UIAlertView(已经过时) UIActionView swift
// // ViewController.swift // UILabelTest // // Created by mac on 15/6/23. // Copyright (c) 2015年 fa ...
- ActiveX控件的基本操作方法以及如何在VS2010下使用控件
在此,小编就介绍下ActiveX控件的基本操作方法以及如何在VS2010下使用控件,我们以一个程序为例, (1) 打开VS2010编译器(右键以管理员身份运行,因为ActiveX需要注册), ...