[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 ...
随机推荐
- 封装document.ready方法
function $(fn){ if(document.addEventListener){ //W3C document.addEventListener('DOMContentLoaded',fu ...
- DevExpress后置代码中初始化SQL数据源的方法
//初始化SQL数据源的提供者和连接字符串 函数 OK public virtual void InitSqlDataSource_ConStr(SqlDataSource sql_ds) { Con ...
- 并行编程之CountdownEvent的用法
教程:http://blog.gkarch.com/threading/part5.html#the-parallel-class http://www.cnblogs.com/huangxinche ...
- FileUpload控件「批次上传 / 多档案同时上传」的范例--以「流水号」产生「变量名称」
原文出處 http://www.dotblogs.com.tw/mis2000lab/archive/2013/08/19/multiple_fileupload_asp_net_20130819. ...
- C# 平时碰见的问题【3】
今天发现一个问题纳闷了半个小时, 需求是处理project文件里边的数据内容,其中需要判断任务名称不存在重复; 在测试的时候弄了两行一样的任务,如预想: 任务[xxx]重复 然后删掉重复的任务行,继续 ...
- 第五节:什么导致Finalize方法被调用
Finalize方法在垃圾回收结束时被调用,下面有5种事件会导致开始垃圾回收 1.第0代已满 第0代已满,垃圾回收会自动开始.该事件是目前导致Finalize方法被调用的最常见的一种方式,因为虽 ...
- C# 页面抓取获取快递信息
通过页面抓取信息可以获得很多我们想要的信息,比如现在常会用到的快递查询,主要抓取的网站为http://www.kuaidi100.com/ 通过IE的网络分析我们可以得到下面信息 通过对这个网站的分析 ...
- Activity(一)
一个应用程序中至少包含一个Activity,Activity启动流程:当启动一个应用程序时,android操作系统会访问该应用程序的AndroidManifest.xml文件(该文件中说明了应用程序使 ...
- C# WinForm自定义控件响应键盘事件
自己定义的winform控件,用其他键盘事件都无法响应,只有用ProcessCmdKey事件可以达到目的(别忘了主窗体的KeyPreview属性要设置为true),写法如下: prot ...
- hdu 2066 一个人的旅行
题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=2066 一个人的旅行 Description 虽然草儿是个路痴(就是在杭电待了一年多,居然还会在校园里迷 ...