【bzoj2559】【IOI2011】【Race】【点分治】

2024-02-20 15:58
文章标签 分治 race bzoj2559 ioi2011

本文主要是介绍【bzoj2559】【IOI2011】【Race】【点分治】,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Description

给一棵树,每条边有权.求一条路径,权值和等于K,且边的数量最小.

Input

第一行 两个整数 n, k
第二..n行 每行三个整数 表示一条无向边的两端和权值 (注意点的编号从0开始)

Output

一个整数 表示最小边数量 如果不存在这样的路径 输出-1

Sample Input

4 3
0 1 1
1 2 2
1 3 4

Sample Output

2
题解:直接上点分治。
对每次计算的子树搞一个f数组,f[i]表示在当前子树中边权为i的路径的最小长度。
然后我们再维护一个dis数组和d数组,分别表示在当前子树中点i到根的权值和和深度。
显然对于当前计算的子树的每个点x,f[dis[x]]=min(f[dis[x]],d[x]);
ans=min(ans,f[k-dis[x]]+d[x]);
在dfs中维护着两个量即可。
注意每次做完一颗子树要再dfs一边把f数组清空。
代码:
#include<iostream>
#include<cstdio>
#define N 200010
using namespace std;
int point[N],n,k,cnt,next[N*2],x,y,s[N],h[N],ans,f[N*10],dis[N],d[N],r,temp,w;
bool visit[N];
struct use{int st,en,v;}e[N*2];
void add(int x,int y,int v){next[++cnt]=point[x];point[x]=cnt;e[cnt].st=x;e[cnt].en=y;e[cnt].v=v;
}
void dfs(int x,int fa){s[x]=1;h[x]=0;for (int i=point[x];i;i=next[i]){if (e[i].en!=fa&&!visit[e[i].en]){dfs(e[i].en,x);s[x]+=s[e[i].en];h[x]=max(h[x],s[e[i].en]);}}h[x]=max(h[x],temp-s[x]);if (h[x]<h[r]) r=x;
}
void calc(int x,int fa){if (dis[x]<=k) ans=min(ans,f[k-dis[x]]+d[x]);for (int i=point[x];i;i=next[i])if (e[i].en!=fa&&!visit[e[i].en]){d[e[i].en]=d[x]+1;dis[e[i].en]=dis[x]+e[i].v;calc(e[i].en,x);}             
}
void updata(int x,int fa,bool p){if(dis[x]<=k){if(p)f[dis[x]]=min(f[dis[x]],d[x]);else f[dis[x]]=99999999;}for(int i=point[x];i;i=next[i])if(e[i].en!=fa&&!visit[e[i].en])updata(e[i].en,x,p);
}
void solve(int x){visit[x]=1;f[0]=0;for (int i=point[x];i;i=next[i]){if (!visit[e[i].en]){d[e[i].en]=1;dis[e[i].en]=e[i].v;calc(e[i].en,0);updata(e[i].en,0,1);}}for (int i=point[x];i;i=next[i])if (!visit[e[i].en]) updata(e[i].en,0,0);for (int i=point[x];i;i=next[i]){if (!visit[e[i].en]){r=0;temp=s[e[i].en];dfs(e[i].en,0);solve(r);}}
}
int main(){scanf("%d%d",&n,&k);for (int i=1;i<n;i++){scanf("%d%d%d",&x,&y,&w);x++;y++;add(x,y,w);add(y,x,w);}for (int i=1;i<=k;i++) f[i]=n;temp=ans=h[0]=n;dfs(1,0);solve(r);if (ans!=n) cout<<ans<<endl;else cout<<"-1"<<endl;
}




这篇关于【bzoj2559】【IOI2011】【Race】【点分治】的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

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

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

分治算法与凸包问题

1. 什么是凸包问题? 凸包问题是计算几何中的经典问题。给定二维平面上的点集,凸包是一个最小的凸多边形,它包含了点集中所有的点。你可以把凸包想象成一根松紧带将所有点紧紧包裹住的样子,凸包的边缘仅沿着最外面的点延伸。 2. 分治法简介 分治算法是解决复杂问题的强大策略,它的思想是将问题分解为多个子问题,分别解决这些子问题后再合并得到最终解。凸包问题可以通过分治算法高效地解决,时间复杂度可以达到

【COGS】577 蝗灾 cdq分治

传送门:【COGS】577 蝗灾 题目分析:cdq分治入门题= =。。。。用差分思想将矩阵分成四块来计算。。排序一维,另一维用树状数组解决。 代码如下: #include <cstdio>#include <vector>#include <cstring>#include <algorithm>using namespace std ;#define REP(

【ACdream】1157 Segments cdq分治

传送门:【ACdream】1157 Segments 题目分析:第一题cdq(陈丹琦)分治!cdq_____Orz! 听说cdq分治可以写,就去学了cdq分治了。。 在我们平常使用的分治中,每一个子问题只解决它本身(可以说是封闭的)。 而在cdq分治中,对于划分出来的两个子问题,前一个子问题用来解决后一个子问题而不是它本身。 具体算法流程如下: 1.将整个操作序列分为两个长

【HDU】4871 Shortest-path tree 最短路+点分治

传送门:【HDU】4871 Shortest-path tree 题目分析: 学了点分治后来看这道题,简直就是水题。。。 但是我竟然花了将近一个晚上才写出来。。。就因为一个地方写漏了T U T。。 首先根据题意求一棵树,最短路一下,然后最小字典序就按照编号顺序遍历邻接表给节点标记。 剩下的就是树分治的事了。 在以重心X为根的子树中,按照X的子节点v的子树中最长路径包含节点数升序遍

【HDU】4812 D Tree 点分治

传送门:【HDU】4812 D Tree 题目分析:点分治搞之。乘积等于K的路径。 首先我们定义一个path[ i ]用以记录从根结点x在子树x内的第 i 条路径的值(乘积)。然后每次我们搞完当前重心rt的一棵子树以后,我们用判断K*逆元[ path[ i ] * val[ rt ] %MOD ] % MOD 是否存在来确定乘积为K的路径是否存在,然后再用这个path[ i ]去更

【SPOJ】1825 Free tour II 点分治

传送门:【SPOJ】1825 Free tour II 题目分析:敲了两遍。。。 本题是论文题,具体见漆子超论文《分治算法在树的路径问题中的应用》。 在以root为根的第 i 棵子树上,我们用G[ i ,j ]表示root的第 i 棵子树的路径上严格有 j 个黑点的路径的最长长度。用F[ i ,j ]表示在root为根的第 i 棵子树的路径上不超过 j 个黑点的路径的最长长度。因

【HDU】4670 Cube number on a tree 点分治

传送门:【HDU】4670 Cube number on a tree 题目分析:首先因为至多30个素数,3^30在long long以内,如果一条路径上的数的乘积是个立方数,则这条路径上每个素数因子的个数都应该是3的倍数,于是我们用三进制表示含有素数的状态,当且仅当状态为0(即所有素数的个数都是3的倍数)时这条路径上数的乘积为完全立方数。考虑树分治,每层分治,求出当前重心的一个儿子的一个

【BZOJ】2152: 聪聪可可 点分治

传送门:【BZOJ】2152: 聪聪可可 题目分析:记录权值和%3的路径的个数。。。然后去重。。没了。。 代码如下: #include <vector>#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>using namespace std ;typedef lo