HDU
来源:互联网 发布:程序员示爱网页源代码 编辑:程序博客网 时间:2024/05/20 07:57
题解:
这题不分在最短路里面,我以为是个数学题呢,最后是个图论题。转自大佬的http://www.cnblogs.com/kuangbin/archive/2012/08/17/2644557.html
显然,题目给的是一个0/1规划模型。
解题的关键在于如何看出这个模型的本质。
3个条件明显在刻画未知数之间的关系,从图论的角度思考问题,容易得到下面3个结论:
1.X12+X13+...X1n=1 于是1号节点的出度为1
2..X1n+X2n+...Xn-1n=1 于是n号节点的入度为1
3.∑Xki =∑Xij 于是2~n-1号节点的入度必须等于出度
于是3个条件等价于一条从1号节点到n号节点的路径,故Xij=1表示需要经过边(i,j),代价为Cij。Xij=0表示不经过边(i,j)。注意到Cij非负且题目要求总代价最小,因此最优答案的路径一定可以对应一条简单路径。
最终,我们直接读入边权的邻接矩阵,跑一次1到n的最短路即可,记最短路为path。
以上情况设为A
非常非常非常非常非常非常非常非常抱歉,简单路径只是充分条件,但不必要。(对造成困扰的队伍深表歉意)
漏了如下的情况B:
从1出发,走一个环(至少经过1个点,即不能是自环),回到1;从n出发,走一个环(同理),回到n。
容易验证,这是符合题目条件的。且A || B为该题要求的充要条件。
由于边权非负,于是两个环对应着两个简单环。
因此我们可以从1出发,找一个最小花费环,记代价为c1,再从n出发,找一个最小花费环,记代价为c2。(只需在最短路算法更新权值时多加一条记录即可:if(i==S) cir=min(cir,dis[u]+g[u][i]))
故最终答案为min(path,c1+c2)
/*HDU 4370 0 or 1转换思维的题啊,由一道让人不知如何下手的题,转换为了最短路基本思路就是把矩阵看做一个图,图中有n个点,1号点出度为1,n号点入度为1,其它点出度和入度相等,路径长度都是非负数,等价于一条从1号节点到n号节点的路径,故Xij=1表示需要经过边(i,j),代价为Cij。Xij=0表示不经过边(i,j)。注意到Cij非负且题目要求总代价最小,因此最优答案的路径一定可以对应一条简单路径。最终,我们直接读入边权的邻接矩阵,跑一次1到n的最短路即可,记最短路为path。漏了如下的情况B:从1出发,走一个环(至少经过1个点,即不能是自环),回到1;从n出发,走一个环(同理),回到n。也就是1和n点的出度和入度都为1,其它点的出度和入度为0.由于边权非负,于是两个环对应着两个简单环。因此我们可以从1出发,找一个最小花费环,记代价为c1,再从n出发,找一个最小花费环,记代价为c2。(只需在最短路算法更新权值时多加一条记录即可:if(i==S) cir=min(cir,dis[u]+g[u][i]))故最终答案为min(path,c1+c2)*//*本程序用SPFA来完成最短路。但是由于要计算从出发点出发的闭环的路径长度。所以要在普通SPFA的基础上做点变化。就是把dist[start]设为INF。同时一开始并不是让出发点入队,而是让出发点能够到达的点入队。*/#include<stdio.h>#include<iostream>#include<string.h>#include<algorithm>using namespace std;const int INF=0x3f3f3f3f;const int MAXN=330;int cost[MAXN][MAXN];//保存路径长度的邻接矩阵int dist[MAXN];int que[MAXN];//注意队列的循环利用,建成循环队列bool vis[MAXN];//是否在队列中标记void SPFA(int start,int n){ int front=0,rear=0; for(int v=1;v<=n;v++)//初始化 { if(v==start)//由于要找start的闭环,所以dist[start]设为INF,且不入队 { dist[v]=INF; vis[v]=false; } else if(cost[start][v]!=INF) { dist[v]=cost[start][v]; que[rear++]=v; vis[v]=true; } else//即dist[start][v]==INF情况,对本题没有这种情况 { dist[v]=INF; vis[v]=false; } } while(front!=rear)//注意这个条件是不等,因为是循环队列 { int u=que[front++]; for(int v=1;v<=n;v++) { if(dist[v]>dist[u]+cost[u][v]) { dist[v]=dist[u]+cost[u][v]; if(!vis[v])//不在队列 { vis[v]=true; que[rear++]=v; if(rear>=MAXN) rear=0;//循环队列 } } } vis[u]=false; if(front>=MAXN)front=0; }}int main(){ //freopen("in.txt","r",stdin); //freopen("out.txt","w",stdout); int n; while(scanf("%d",&n)!=EOF) { for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&cost[i][j]); SPFA(1,n); int ans=dist[n];//1到n的最短路 int loop1=dist[1];//1的闭环长度 SPFA(n,n); int loopn=dist[n];//n的闭环长度 ans=min(ans,loop1+loopn); printf("%d\n",ans); } return 0;}
/*用堆栈实现SPFA,有时候比队列快*/#include<stdio.h>#include<iostream>#include<string.h>#include<algorithm>using namespace std;const int MAXN=330;const int INF=0x3f3f3f3f;int cost[MAXN][MAXN];int dist[MAXN];int Q[MAXN];bool vis[MAXN];void SPFA(int start,int n){//堆栈实现,有时候比队列快 int top=0; for(int v=1;v<=n;v++) { if(v==start) { dist[v]=INF; vis[v]=false; } else { dist[v]=cost[start][v]; vis[v]=true; Q[top++]=v; } } while(top!=0) { int u=Q[--top]; for(int v=1;v<=n;v++) { if(dist[v]>dist[u]+cost[u][v]) { dist[v]=dist[u]+cost[u][v]; if(!vis[v]) { vis[v]=true; Q[top++]=v; } } } vis[u]=false; }}int main(){ //freopen("in.txt","r",stdin); //freopen("out.txt","w",stdout); int n; while(scanf("%d",&n)!=EOF) { for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&cost[i][j]); SPFA(1,n); int ans=dist[n];//1到n的最短路 int loop1=dist[1];//1的闭环长度 SPFA(n,n); int loopn=dist[n];//n的闭环长度 ans=min(ans,loop1+loopn); printf("%d\n",ans); } return 0;}
- hdu
- hdu
- HDU
- hdu ()
- hdu
- hdu
- HDU
- HDU
- hdu
- hdu
- HDU
- Hdu
- hdu
- hdu-
- hdu
- hdu
- hdu
- HDU
- SVN
- Linux Centos7安装Oracle12c第二版本
- Html5前端图片压缩
- iOS 解决NSTimer在UITableView滑动中暂停问题
- 获取json数据和ImageLoader图片显示到listview
- HDU
- 《一个人的朝圣》读后感
- 阿里云Centos下安装mysql找不到mysql-sever安装包的解决方案
- 网狐IOS端微信登录遇到的问题
- 【51nod 1085】背包问题(0-1背包模板)
- 如何同步本地代码到SVN服务器
- 1、java容易忘记的基础语法知识
- stm32:error: #268: declaration may not appear after executable statement in block
- GNU GPL介绍