luoguP4383 [八省联考2018]林克卡特树(树上dp,wqs二分)

uoguP4383 [八省联考2018]林克卡特树(树上dp,wqs二分)

Luogu

题解时间

$ k $ 条边权为 $ 0 $ 的边。

是的,边权为零。

转化成选正好 $ k+1 $ 条链。

$ k \le 100 $ 的部分。

毫无疑问是树上打背包dp。

但具体设计还要注意一下。

一个问题是单点成链,这个要特判。

之后由于选择的都是链,所以每个点的度数不会超过2.

这样方程就出来了。

$ k \le n $ 的部分。

很明显不能背包了。

但“选正好k个求最大权值和”这个要求如果熟悉的话可能想到wqs二分。

打表试一下发现确实是上凸函数。

之后就按wqs二分的套路来,二分加权mid,求 $ dp_{x,k,deg} - mid * k $ 就完事了。

#include<bits/stdc++.h>
using namespace std;
typedef long long lint;
template<typename TP>inline void read(TP &tar)
{
	TP ret=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){ret=ret*10+(ch-'0');ch=getchar();}
	tar=ret*f;
}
namespace RKK
{
const int N=300011;
const int inf=0x3f3f3f3f;
const lint linf=0x3f3f3f3f3f3f3f;
struct sumireko{int to,ne;lint w;}e[N<<1];int he[N],ecnt;
void addline(int f,int t,lint w){e[++ecnt].to=t,e[ecnt].w=w;e[ecnt].ne=he[f],he[f]=ecnt;}
struct pat
{
	lint v,k;pat(lint v=-linf,lint k=inf):v(v),k(k){}
	bool operator<(const pat &p)const{return v==p.v?k>p.k:v<p.v;}
	pat friend operator+(const pat &pa,const pat &pb){return pat(pa.v+pb.v,pa.k+pb.k);}
}dp[N][3],dg[3],dtmp;

int n,kap;lint km;
void dfs(int x,int f)
{
	for(int i=he[x],t=e[i].to,w=e[i].w;i;i=e[i].ne,t=e[i].to,w=e[i].w)if(t!=f)
	{
		dfs(t,x);for(int p=0;p<3;p++) dg[p]=pat();
		for(int j=0;j<3;j++)for(int p=0;p<3;p++) dg[j]=max(dg[j],dp[x][j]+dp[t][p]);
		dg[1]=max(dg[1],dp[x][0]+max(dp[t][0]+pat(w-km,1),dp[t][1]+pat(w,0)));
		dg[2]=max(dg[2],dp[x][1]+max(dp[t][0]+pat(w,0),dp[t][1]+pat(w+km,-1)));
		memcpy(dp[x],dg,sizeof(dg));
	}
}
void dpclr(){for(int i=1;i<=n;i++) dp[i][0]=pat(0,0),dp[i][1]=pat(),dp[i][2]=pat(-km,1);}

int main()
{
	#ifdef RDEBUG
	freopen("sample.in","r",stdin);
	#endif
	read(n),read(kap),kap++;for(int i=2,x,y,w;i<=n;i++) read(x),read(y),read(w),addline(x,y,w),addline(y,x,w);
	lint ql=-1e8,qr=1e8,qm,qa;
	while(ql<=qr)
	{
		qm=ql+qr>>1,km=qm;
		dpclr(),dfs(1,0),dtmp=max(dp[1][0],max(dp[1][1],dp[1][2]));
		if(dtmp.k==kap){printf("%lld\n",dtmp.v+km*kap);return 0;}
		else if(dtmp.k>kap) ql=qm+1;
		else qr=qm-1,qa=qm;
	}
	km=qa;
	dpclr(),dfs(1,0),dtmp=max(dp[1][0],max(dp[1][1],dp[1][2]));
	printf("%lld\n",dtmp.v+km*kap);
	return 0;
}
}
int main(){return RKK::main();}
上一篇:【机器学习系列】MCMC第三讲:理解MCMC前必先弄懂这两点


下一篇:ADO连接数据库