LightOJ1086 Jogging Trails(欧拉回路+中国邮递员问题+SPFA)
2016-01-18 21:19
459 查看
题目求从某点出发回到该点经过所有边至少一次的最短行程。
这个问题我在《图论算法理论、实现及应用》中看过,是一个经典的问题——中国邮递员问题(CPP, chinese postman problem)也称为中国邮路问题,是我国数学家管梅谷教授于1962年首次提出的,引起了世界不少数学家的关注。例如1973年匈牙利数学家Edmonds和Johnsom对中国邮路问题提出了一种有效算法。
解决的方法就是通过添加边,使其成为欧拉回路。
而这一题的问题就可以转化给这个无向图加最短的边使其所有点的度位偶数(为什么不说欧拉回路,因为题目没说图一定连通。。
这道题在LightOJ分类是DP。。我想大概就是15个点二进制压缩的状压DP。也就是把各点度为偶数奇数分别看作二进制各位0和1,然后每次可以通过在任意两个奇度或偶度点加一条边来转移。。
不过显然这是不行的,这样的状态不满足DP的无后效性,状态是点转移是边那就可以是个有向有环图——有向有环图——于是我就想到最短路,显然可以用最短路做,不过时间复杂度不好说。。具体来说:
偶度点记0奇度记1压缩成的二进制就是要求最短路的图中的顶点,这样顶点最多215个;
而如果两个顶点有且仅有i和j这两位的二进制不同且原图中i到j有路径,那么这两个顶点间有一条边且权为原图中i到j中最短的那条路径;
最短路的起点就是原图度序列压缩表示的顶点,终点就是0。
如此构图,跑个SPFA,然后就AC了。
这个问题我在《图论算法理论、实现及应用》中看过,是一个经典的问题——中国邮递员问题(CPP, chinese postman problem)也称为中国邮路问题,是我国数学家管梅谷教授于1962年首次提出的,引起了世界不少数学家的关注。例如1973年匈牙利数学家Edmonds和Johnsom对中国邮路问题提出了一种有效算法。
解决的方法就是通过添加边,使其成为欧拉回路。
而这一题的问题就可以转化给这个无向图加最短的边使其所有点的度位偶数(为什么不说欧拉回路,因为题目没说图一定连通。。
这道题在LightOJ分类是DP。。我想大概就是15个点二进制压缩的状压DP。也就是把各点度为偶数奇数分别看作二进制各位0和1,然后每次可以通过在任意两个奇度或偶度点加一条边来转移。。
不过显然这是不行的,这样的状态不满足DP的无后效性,状态是点转移是边那就可以是个有向有环图——有向有环图——于是我就想到最短路,显然可以用最短路做,不过时间复杂度不好说。。具体来说:
偶度点记0奇度记1压缩成的二进制就是要求最短路的图中的顶点,这样顶点最多215个;
而如果两个顶点有且仅有i和j这两位的二进制不同且原图中i到j有路径,那么这两个顶点间有一条边且权为原图中i到j中最短的那条路径;
最短路的起点就是原图度序列压缩表示的顶点,终点就是0。
如此构图,跑个SPFA,然后就AC了。
#include<cstdio> #include<cstring> #include<queue> using namespace std; #define INF (1<<30) int n,map[15][15]; int d[1<<15]; bool vis[1<<15]; int SPFA(int vs){ for(int i=0; i<(1<<15); ++i){ d[i]=INF; vis[i]=0; } d[vs]=0; vis[vs]=1; queue<int> que; que.push(vs); while(!que.empty()){ int u=que.front(); que.pop(); for(int i=0; i<n; ++i){ for(int j=i+1; j<n; ++j){ if(map[i][j]==-1) continue; int v=u; v^=(1<<i); v^=(1<<j); if(d[v]>d[u]+map[i][j]){ d[v]=d[u]+map[i][j]; if(!vis[v]){ vis[v]=1; que.push(v); } } } } vis[u]=0; } return d[0]; } int main(){ int t,m,a,b,c; scanf("%d",&t); for(int cse=1; cse<=t; ++cse){ scanf("%d%d",&n,&m); memset(map,-1,sizeof(map)); int deg[15]={0},res=0; while(m--){ scanf("%d%d%d",&a,&b,&c); --a; --b; ++deg[a]; ++deg[b]; res+=c; if(map[a][b]==-1 || map[a][b]>c) map[a][b]=map[b][a]=c; } int vs=0; for(int i=0; i<n; ++i){ if(deg[i]&1) vs|=(1<<i); } res+=SPFA(vs); printf("Case %d: %d\n",cse,res); } return 0; }
相关文章推荐
- 关于Cookie中的path、domain、expire、secure属性详述
- Factorial Trailing Zeroes
- Gitlab RPC failed; result=22, HTTP code = 413 解决方法
- ftk学习记(waitbox篇)
- DataOutputStream、DataInputStream
- Closed kinematic chain - how to control it
- Aix6.1操作系统安装教程
- 7.4 决策树
- 3.4wait的条件发生变化导致逻辑错误
- raid详细介绍
- (assign,retain,copy):
- 美国末日AI之状态机(The Last Of Us AI)
- GPG key generation: Not enough random bytes available.
- Project Euler 102:Triangle containment 包含原点的三角形
- Syntax error, parameterized types are only available if source level is 1.5 or greater
- RAID5单块和多块硬盘故障如何恢复
- 3.3wait(long)
- GTP协议的分析 http://blog.csdn.net/stephen_yin/article/details/6951258
- <蜡烛图精解Candlestick Charting Explained>-- 笔记1
- GTP(GPRS Tunnelling Protocol)协议http://blog.csdn.net/stephen_yin/article/details/6951237