[八省联考2018]林克卡特树lct
题面在这里
description
一个\(N\)个点的\(Tree\),每条边有一个整数边权\(v_i\),表示走这条边会获得\(v_i\)的收益;
小\(L\)需要控制主角\(Link\),\(Cut\)掉树上的恰好\(K\)条边,然后再连接\(K\)条边权为$ 0\(的边,得到一棵新的树。
接着,他会选择树上的两个点\)p,q,\(并沿着树上连接这两点的简单路径从\)p\(走到\)q\(,并为经过的每条边获取相应收益。
请求出\)Link$能得到的最大总收益。
data range
对于\(10%\)的数据,\(k = 0\);
对于另外\(10%\)的数据,\(k = 1\);
对于另外\(15%\)的数据,\(k = 2\);
对于另外\(25%\)的数据,\(k \le 100\);
对于全部的测试数据,
保证有\(1\le N\le 3\times 10^5,1 \le x_i, y_i \le N; |v_i| \le 10^6\)。
Hint
题目并不难。
solution
上面的提示并不能让我这种蒟蒻得到什么有用的信息
不过分析一波之后还是能发现题意求的是树上\(k\)条点不相交直径的最大边权和
60 point
使用树形DP就好啦
设\(f[x][k][d=0,1,2]\)表示在\(x\)的子树中选取\(k\)调不相交的链,
且\(x\)的度数为\(d\)时的最大链长和
合并的时候需要注意讨论一下,由于不是重点所以没有写得那么详细,
具体实现请看代码
100 point
当树的形态一定时,选取\(k\)条链时的最优解\(best_k\)是上凸的
真是神一般的性质啊
于是我们想到用直线去切这个凸包
那么用一条斜率为\(k\)的直线切这个凸包\(best_x\)的意义是什么呢?
本蒟蒻是这么理解的(图画得不好,看看就行了):
当我们要用一条直线\(y=kx\)切凸包的时候,这个切点满足\(best_x-kx\)最大。
那么我们可以转而在树上选任意条数的链,但每条链需要花费\(k\)的代价,
之后树形DP求出满足这个条件的链条数和代价(这是可以\(O(n)\)做的),
就是我们直线要切的点\((x,best_x-kx)\)
二分斜率,于是我们就做完了这道题?
复杂度为\(O(nlog(\sum w))\)
code
这里是60分的DP(未经在线测试)
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#define FILE "lct"
#define RG register
#define il inline
using namespace std;
typedef long long ll;
typedef double dd;
typedef vector<int> VI;
const int N=300010;
const ll INF=1e18;
il ll read(){
RG ll d=0,w=1;char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')w=-1,ch=getchar();
while(ch<='9'&&ch>='0')d=d*10+ch-48,ch=getchar();
return d*w;
}
il void file(){
freopen(FILE".in","r",stdin);
freopen(FILE".out","w",stdout);
}
int n,k;
int head[N],nxt[N<<1],to[N<<1],val[N<<1],cnt;
il void add(int u,int v,int w){
to[++cnt]=v;
nxt[cnt]=head[u];
val[cnt]=w;
head[u]=cnt;
}
ll f[N][105][3],ans;
il void update(ll &a,ll b){a=a>b?a:b;}
#define fs (max(f[son][j][0],max(f[son][j][3],f[son][j][2])))
il void DP(int son,int fa,int e){
for(RG int i=k;i;i--)
for(RG int j=0;i-j+1>=0;j++){
update(f[fa][i][2],f[fa][i-j+1][4]+f[son][j][5]+val[e]);
if(i-j>=0)update(f[fa][i][2],f[fa][i-j][2]+fs);
}
for(RG int i=k;i;i--)
for(RG int j=0;i-j>=0;j++){
if(i-j-1>=0)
update(f[fa][i][6],f[fa][i-j-1][0]+f[son][j][0]);
update(f[fa][i][7],f[fa][i-j][0]+f[son][j][8]+val[e]);
update(f[fa][i][9],f[fa][i-j][10]+fs);
}
for(RG int i=k;i;i--)
for(RG int j=0;i-j>=0;j++)
update(f[fa][i][0],f[fa][i-j][0]+fs);
}
void dfs(int u,int fa){
f[u][0][0]=f[u][11][1]=0;
for(RG int i=head[u];i;i=nxt[i]){
RG int v=to[i];if(v==fa)continue;
dfs(v,u);DP(v,u,i);
}
}
int main()
{
n=read();k=read()+1;
for(RG int i=1,u,v,w;i<n;i++){
u=read();v=read();w=read();
add(u,v,w);add(v,u,w);
}
memset(f,-127/3,sizeof(f));
dfs(1,0);
for(RG int i=1;i<=n;i++)
ans=max(max(ans,f[i][k][0]),max(f[i][k][12],f[i][k][2]));
printf("%lld\n",ans);
return 0;
}
这里是AC代码,搞了个pair方便统计答案
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<queue>
#include<vector>
#define FILE "lct"
#define RG register
#define il inline
#define line first
#define val second
using namespace std;
typedef long long ll;
typedef double dd;
typedef vector<int> VI;
typedef pair<ll,ll> PL;
const int N=300010;
const ll INF=1e18;
il ll read(){
RG ll d=0,w=1;char ch=getchar();
while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
if(ch=='-')w=-1,ch=getchar();
while(ch<='9'&&ch>='0')d=d*10+ch-48,ch=getchar();
return d*w;
}
il void file(){
freopen(FILE".in","r",stdin);
freopen(FILE".out","w",stdout);
}
int n,k;
int head[N],nxt[N<<1],to[N<<1],val[N<<1],cnt;
il void add(int u,int v,int w){
to[++cnt]=v;
nxt[cnt]=head[u];
val[cnt]=w;
head[u]=cnt;
}
PL f[N][3];
PL operator +(PL a,PL b){
return (PL){a.line+b.line,a.val+b.val};
}
PL operator +(PL a,ll v){a.val+=v;return a;}
il PL max(PL a,PL b){
if(a.val>b.val)return a;
else if(a.val==b.val&&a.line>b.line)return a;
return b;
}
il void update(PL &a,PL b,int l){b.line+=l;a=max(a,b);}
#define fs (max(f[son][0],max(f[son][1],f[son][2])))
il void DP(int son,int fa,int e,ll K){
update(f[fa][2],f[fa][2]+fs,0);
update(f[fa][2],f[fa][1]+f[son][1]+val[e]+K,-1);
update(f[fa][1],f[fa][1]+fs,0);
update(f[fa][1],f[fa][0]+f[son][0]+(-K),1);
update(f[fa][1],f[fa][0]+f[son][1]+val[e],0);
update(f[fa][0],f[fa][0]+fs,0);
}
void dfs(int u,int fa,ll K){
f[u][0]=(PL){0,0};f[u][1]=(PL){1,-K};
f[u][2]=(PL){0,-INF};
for(RG int i=head[u];i;i=nxt[i]){
RG int v=to[i];if(v==fa)continue;
dfs(v,u,K);DP(v,u,i,K);
}
}
PL work(ll K){
RG PL ans;ans=(PL){0,-INF};
dfs(1,0,K);
for(RG int i=1;i<=n;i++)
ans=max(max(ans,f[i][0]),max(f[i][1],f[i][2]));
return ans;
}
int main()
{
n=read();k=read()+1;RG ll L=0,R=0,M;RG PL ret;
for(RG int i=1,u,v,w;i<n;i++){
u=read();v=read();w=read();
add(u,v,w);add(v,u,w);R+=abs(w);L-=abs(w);
}
while(L<R){
M=(L+R)>>1;ret=work(M);
if(ret.line<=k)R=M;else L=M+1;
}
ret=work(R);
printf("%lld\n",ret.val+ret.line*R);
return 0;
}
诸位珍重,望来日以富贵相见。