SPOJ 1825 Free Tour | 终极之树分治
求树上最长路径使得经过的拥挤节点个数不超过K
//欢迎访问这个博客!http://www.cnblogs.com/luyouqi233/p/8036828.html #include<cstdio>
#include<algorithm>
#include<cstring>
#include<vector>
#define N 200010
using namespace std;
typedef long long ll;
typedef pair<int,int>pii;
int head[N],ecnt,dis[N],fa[N],son[N],sze[N],ans,n,m,k,que[N],t[N],nump[N],maxn[N];
//dis[i] i到根节点的距离,fa[i] 父亲,son[i] i的最大子树大小,sze[i] BFS时以i为根的子树大小(son[]和sze[])用来求重心
//t[i] 当前子树节点到子树根且经过至多i个黑点的路径的最大权值和
//nump[i] 统计i到根节点的路径上的黑点数
//maxn[i] 以根节点为起点,终点在前面已经搜过的子树中,且该路径至多经过i个黑点,这样的路径的最大权值和。
bool vis[N],iscrowd[N];//iscrowd[i] 是否是拥挤节点
inline int read()
{
int ret=0,neg=1;
char j=getchar();
for (;j>'9' || j<'0';j=getchar())
if (j=='-') neg=-1;
for (;j>='0' && j<='9';j=getchar())
ret=ret*10+j-'0';
return ret*neg;
}
struct edge
{
int nxt,v,w;
}e[2*N];
vector <pair<int,int> > s;
void add(int u,int v,int w)
{
e[++ecnt].v=v,e[ecnt].w=w,e[ecnt].nxt=head[u],head[u]=ecnt;
e[++ecnt].v=u,e[ecnt].w=w,e[ecnt].nxt=head[v],head[v]=ecnt;
}
int calcG(int sv)
{
int u,v,mx=n,G,qn;
que[qn=1]=sv,fa[sv]=0;
for (int ql=1;ql<=qn;ql++)
{
sze[u=que[ql]]=1,son[u]=0;
for (int i=head[u];i;i=e[i].nxt)
{
if (vis[v=e[i].v] || v==fa[u]) continue;
fa[v]=u,que[++qn]=v;
}
}
for (int ql=qn;ql>=1;ql--)
{
u=que[ql],v=fa[u];
if (qn-sze[u]>son[u]) son[u]=qn-sze[u];
if (son[u]<mx) G=u,mx=son[u];
if (!v) break;
sze[v]+=sze[u];
if (sze[u]>son[v]) son[v]=sze[u];
}
return G;
}
inline int getmaxp(int st,ll L)
{
int qn=0,maxp=0;
que[++qn]=st;
nump[st]=iscrowd[st];
dis[st]=L;
fa[st]=0;
maxp=max(maxp,nump[st]);
for (int ql=1,u;ql<=qn;ql++)
for (int i=head[u=que[ql]],v;i;i=e[i].nxt)
{
if (vis[v=e[i].v] || v==fa[u]) continue;
fa[v]=u;
dis[v]=dis[u]+e[i].w;
nump[v]=nump[u]+iscrowd[v];
maxp=max(maxp,nump[v]);
que[++qn]=v;
}
return maxp;
}
inline void getmaxt(int st)
{
int qn=0;
que[++qn]=st;
t[nump[st]]=max(t[nump[st]],dis[st]);
for (int ql=1,u;ql<=qn;ql++)
for (int i=head[u=que[ql]],v;i;i=e[i].nxt)
{
if (vis[v=e[i].v] || v==fa[u]) continue;
t[nump[v]]=max(t[nump[v]],dis[v]);
que[++qn]=v;
}
return ;
}
void solve(int x)
{
int G=calcG(x);
vis[G]=1;s.clear();
for (int i=head[G],v;i;i=e[i].nxt)
if (!vis[v=e[i].v])
s.push_back(pii(getmaxp(v,e[i].w),v));
sort(s.begin(),s.end());
if (iscrowd[G]) k--;
for (int i=0;i<s.size();i++)
{
getmaxt(s[i].second);
int now=0;
if (i)
{
for (int j=s[i].first;j>=0;j--)
{
while (now+j<k && now<s[i-1].first)
now++,maxn[now]=max(maxn[now],maxn[now-1]);
if(now+j<=k)ans=max(ans,maxn[now]+t[j]);
}
}
if (i!=s.size()-1)
for (int j=0;j<=s[i].first;j++)
maxn[j]=max(maxn[j],t[j]),t[j]=0;
else
for (int j=0;j<=s[i].first;j++)
{
if (j<=k)
ans=max(ans,max(t[j],maxn[j]));
maxn[j]=t[j]=0;
}
}
if (iscrowd[G]) k++;
for (int i=head[G],v;i;i=e[i].nxt)
if (!vis[v=e[i].v])
solve(v);
return ;
}
int main()
{
scanf("%d%d%d",&n,&k,&m);
for (int i=1,x;i<=m;i++)
iscrowd[read()]=1;
for (int i=1,x,y,z;i<n;i++)
x=read(),y=read(),z=read(),add(x,y,z);
solve(1);
printf("%d",ans);
return 0;
}
SPOJ 1825 Free Tour | 终极之树分治的更多相关文章
- SPOJ 1825 Free tour II (树的点分治)
题目链接 Free tour II 题意:有$N$个顶点的树,节点间有权值, 节点分为黑点和白点. 找一条最长路径使得 路径上黑点数量不超过K个 这是树的点分治比较基本的题,涉及树上启发式合并……仰望 ...
- SPOJ:Free tour II (树分治+启发式合并)
After the success of 2nd anniversary (take a look at problem FTOUR for more details), this 3rd year, ...
- SPOJ 1825 Free tour II 树分治
题意: 给出一颗边带权的数,树上的点有黑色和白色.求一条长度最大且黑色节点不超过k个的最长路径,输出最长的长度. 分析: 说一下题目的坑点: 定义递归函数的前面要加inline,否则会RE.不知道这是 ...
- 【SPOJ】1825. Free tour II(点分治)
http://www.spoj.com/problems/FTOUR2/ 先前看了一会题解就自己yy出来了...对拍过后交tle.................. 自己造了下大数据........t ...
- spoj 1825 Free tour II
http://www.spoj.com/problems/FTOUR2/ After the success of 2nd anniversary (take a look at problem FT ...
- [spoj] FTOUR2 FREE TOUR II || 树分治
原题 给出一颗有n个点的树,其中有M个点是拥挤的,请选出一条最多包含k个拥挤的点的路径使得经过的权值和最大. 正常树分治,每次处理路径,更新答案. 计算每棵子树的deep(本题以经过拥挤节点个数作为d ...
- 算法笔记--树的直径 && 树形dp && 虚树 && 树分治 && 树上差分 && 树链剖分
树的直径: 利用了树的直径的一个性质:距某个点最远的叶子节点一定是树的某一条直径的端点. 先从任意一顶点a出发,bfs找到离它最远的一个叶子顶点b,然后再从b出发bfs找到离b最远的顶点c,那么b和c ...
- hdu-5977 Garden of Eden(树分治)
题目链接: Garden of Eden Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/ ...
- 【BZOJ-1468】Tree 树分治
1468: Tree Time Limit: 10 Sec Memory Limit: 64 MBSubmit: 1025 Solved: 534[Submit][Status][Discuss] ...
随机推荐
- LVM(扩展)
LVM(扩展)======================== [root@aminglinux newdir]# fdisk -l /dev/sdb 磁盘 /dev/sdb:10.7 GB, 107 ...
- 神经网络系列学习笔记(四)——神经网络之RNN学习笔记
不同于传统的FNNs(Feed-forward Neural Networks,前向反馈神经网络),RNNs引入了定向循环,能够处理那些输入之间前后关联的问题. RNNs的目的是用来处理序列数据. 具 ...
- macOs 使用Homebrew升级到MySQL 8系列之后,php无法连接解决方法
当前时间2018-9-28 在使用brew install mysql 默认安装为 MySQL 8,但是使用php连接到数据库之后,出现了这种错误 (Unexpected server respose ...
- c#字符显示转换{0:d} string.Format()
这一篇实际和前几个月写的没什么本质上的区别.但是这篇更明确一点,学起来easy c#字符显示转换{0:d} C#:String.Format数字格式化输出 : int a = 12345678; // ...
- 在庫購買管理(MM)
■購買管理■ [購買依頼]ME51N: 登録ME52N: 変更ME53N: 照会 [購買発注]ME21N: 登録ME22N: 変更ME23N: 照会 [見積依頼]ME41: 登録ME42: 変更ME4 ...
- urllib使用二
urlopen方法返回一个html 对html使用info()方法返回HTTPMessage对象实例 import urllib def print_list(lists): for i in lis ...
- 手把手写代码学习C语言基础
- 9,Linux下的python3,virtualenv,Mysql、nginx、redis安装配置
常用服务安装部署 学了前面的Linux基础,想必童鞋们是不是更感兴趣了?接下来就学习常用服务部署吧! 安装环境: centos7 + vmware + xshell MYSQL(mariadb) ...
- 对于STM32别名区的理解 (转载)
1. 什么是位段.位带别名区? 2. 它有什么好处? 答1: 是这样的,记得MCS51吗? MCS51就是有位操作,以一位(BIT)为数据对象的操作, MCS51可以简单的将P1口的第2位 ...
- 非阻塞IO模板
服务端 from socke import * server = socket(AF_INET, SOCK_STREAM) server.bind(('127.0.0.1',8083)) server ...