【Luogu】 P4649 [IOI2007] training 训练路径

题目链接

点击打开链接

题目解法

好神仙的题啊!!!

首先一个合法的选路径方案等价于没有偶环出现

我们先判掉和树边能组成偶环的非树边

然后考虑一个结论是:如果有一条边被两个偶环都经过了一次,那么这个方案不合法

为什么?考虑把这两条路径的交去掉这两条路径的并,剩下的是一个偶环

考虑把删边变为加边,需要加上权值和最大的边使得方案合法

可以发现每个点的度数很小,于是考虑状压

令 f u , S f_{u,S} fu,S 为在 u u u 的子树中, u u u 的儿子集合 S S S 不考虑在内的方案数

这样只需要枚举每一条 l c a ( x , y ) = u lca(x,y)=u lca(x,y)=u 的非树边,然后转移即可

时间复杂度 O ( 2 10 m ) O(2^{10}m) O(210m)

cpp 复制代码
#include <bits/stdc++.h>
using namespace std;
const int N=1100,M=5100;
int n,m,dp[N][1<<10],depth[N],fa[N];
int e[N<<1],ne[N<<1],h[N],idx;
int id[N][N],rv[N][20];
struct Node{ int x,y,z;}E[M];
vector<Node> qry[N];
inline int read(){
	int FF=0,RR=1;
	char ch=getchar();
	for(;!isdigit(ch);ch=getchar()) if(ch=='-') RR=-1;
	for(;isdigit(ch);ch=getchar()) FF=(FF<<1)+(FF<<3)+ch-48;
	return FF*RR;
}
void add(int x,int y){ e[idx]=y,ne[idx]=h[x],h[x]=idx++;}
void dfs(int u){
	depth[u]=depth[fa[u]]+1;
	for(int i=h[u];~i;i=ne[i]) if(e[i]!=fa[u]) fa[e[i]]=u,dfs(e[i]);
}
int get_lca(int x,int y){
	if(depth[x]>depth[y]) swap(x,y);
	while(depth[y]>depth[x]) y=fa[y];
	while(x!=y) x=fa[x],y=fa[y];
	return x;
}
void dfs2(int u){
	int cnt=0;
	for(int i=h[u];~i;i=ne[i]) if(e[i]!=fa[u]){
		id[u][e[i]]=1<<cnt,rv[u][cnt]=e[i],cnt++;
		dfs2(e[i]);
	}
	for(int S=0;S<1<<cnt;S++)
		for(int i=0;i<cnt;i++) if(!(S>>i&1)) dp[u][S]+=dp[rv[u][i]][0];
	for(Node t:qry[u]){
		int x=t.x,y=t.y,res=t.z,exc=0;
		if(x!=u){
            res+=dp[x][0];
			while(fa[x]!=u) res+=dp[fa[x]][id[fa[x]][x]],x=fa[x];
			exc|=id[fa[x]][x];
		}
		if(y!=u){
            res+=dp[y][0];
			while(fa[y]!=u) res+=dp[fa[y]][id[fa[y]][y]],y=fa[y];
			exc|=id[fa[y]][y];
		}
		for(int S=0;S<1<<cnt;S++) if((S|exc)==S+exc) dp[u][S]=max(dp[u][S],dp[u][S+exc]+res);
	}
}
int main(){
	n=read(),m=read();
	memset(h,-1,sizeof(h));
    int cnt=0;
	for(int i=1;i<=m;i++){
		int x=read(),y=read(),z=read();
		if(!z) add(x,y),add(y,x);
		else E[++cnt]={x,y,z};
	}
	dfs(1);
	int tot=0;
	for(int i=1;i<=cnt;i++){
		int lca=get_lca(E[i].x,E[i].y);tot+=E[i].z;
        if(~(depth[E[i].x]+depth[E[i].y]-2*depth[lca])&1) qry[lca].push_back(E[i]);
	}
	dfs2(1);
	printf("%d\n",tot-dp[1][0]);
	fprintf(stderr,"%d ms\n",int(1e3*clock()/CLOCKS_PER_SEC));
	return 0;
}
相关推荐
三毛的二哥9 小时前
BEV:典型BEV算法总结
人工智能·算法·计算机视觉·3d
南宫萧幕10 小时前
自控PID+MATLAB仿真+混动P0/P1/P2/P3/P4构型
算法·机器学习·matlab·simulink·控制·pid
故事和你9111 小时前
洛谷-数据结构1-4-图的基本应用1
开发语言·数据结构·算法·深度优先·动态规划·图论
我叫黑大帅11 小时前
为什么map查找时间复杂度是O(1)?
后端·算法·面试
炽烈小老头12 小时前
【每天学习一点算法 2026/04/20】除自身以外数组的乘积
学习·算法
skilllite作者12 小时前
AI agent 的 Assistant Auto LLM Routing 规划的思考
网络·人工智能·算法·rust·openclaw·agentskills
py有趣14 小时前
力扣热门100题之不同路径
算法·leetcode
_日拱一卒14 小时前
LeetCode:25K个一组翻转链表
算法·leetcode·链表
啊哦呃咦唔鱼14 小时前
LeetCodehot100-394 字符串解码
算法
小欣加油14 小时前
leetcode2078 两栋颜色不同且距离最远的房子
数据结构·c++·算法·leetcode·职场和发展