博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
【bzoj2324】[ZJOI2011]营救皮卡丘 最短路-Floyd+有上下界费用流
阅读量:5278 次
发布时间:2019-06-14

本文共 2823 字,大约阅读时间需要 9 分钟。

原文地址:http://www.cnblogs.com/GXZlegend/p/6832504.html


题目描述

皮卡丘被火箭队用邪恶的计谋抢走了!这三个坏家伙还给小智留下了赤果果的挑衅!为了皮卡丘,也为了正义,小智和他的朋友们义不容辞的踏上了营救皮卡丘的道路。

火箭队一共有N个据点,据点之间存在M条双向道路。据点分别从1到N标号。小智一行K人从真新镇出发,营救被困在N号据点的皮卡丘。为了方便起见,我们将真新镇视为0号据点,一开始K个人都在0号点。

由于火箭队的重重布防,要想摧毁K号据点,必须按照顺序先摧毁1到K-1号据点,并且,如果K-1号据点没有被摧毁,由于防御的连锁性,小智一行任何一个人进入据点K,都会被发现,并产生严重后果。因此,在K-1号据点被摧毁之前,任何人是不能够经过K号据点的。

为了简化问题,我们忽略战斗环节,小智一行任何一个人经过K号据点即认为K号据点被摧毁。被摧毁的据点依然是可以被经过的。

K个人是可以分头行动的,只要有任何一个人在K-1号据点被摧毁之后,经过K号据点,K号据点就被摧毁了。显然的,只要N号据点被摧毁,皮卡丘就得救了。

野外的道路是不安全的,因此小智一行希望在摧毁N号据点救出皮卡丘的同时,使得K个人所经过的道路的长度总和最少。

请你帮助小智设计一个最佳的营救方案吧!

输入

第一行包含三个正整数N,M,K。表示一共有N+1个据点,分别从0到N编号,以及M条无向边。一开始小智一行共K个人均位于0号点。 

接下来M行,每行三个非负整数,第i行的整数为Ai,Bi,Li。表示存在一条从Ai号据点到Bi号据点的长度为Li的道路。

输出

仅包含一个整数S,为营救皮卡丘所需要经过的最小的道路总和。

样例输入

3 4 2

0 1 1
1 2 1
2 3 100
0 3 1

样例输出

3


题解

最短路-Floyd+有上下界费用流

先用Floyd求出任意两点间距离,注意这里的路径是带有条件的,若为i与j之间的距离,则中间枚举点k必须满足k<=i或k<=j。

因为必须在编号小的点都被摧毁的条件下才能算编号大的点的路径。

这样就求出了题目条件下的两点最短路。

然后就转化为类似于 的问题,拆点费用流即可。

#include 
#include
#include
#define inf 0x3f3f3f3fusing namespace std;queue
q;int map[160][160] , head[400] , to[100000] , val[100000] , cost[100000] , next[100000] , cnt = 1 , s , t , dis[400] , from[400] , pre[400];void add(int x , int y , int v , int c){ to[++cnt] = y , val[cnt] = v , cost[cnt] = c , next[cnt] = head[x] , head[x] = cnt; to[++cnt] = x , val[cnt] = 0 , cost[cnt] = -c , next[cnt] = head[y] , head[y] = cnt;}bool spfa(){ int x , i; memset(dis , 0x3f , sizeof(dis)); memset(from , -1 , sizeof(from)); dis[s] = 0 , q.push(s); while(!q.empty()) { x = q.front() , q.pop(); for(i = head[x] ; i ; i = next[i]) if(val[i] && dis[to[i]] > dis[x] + cost[i]) dis[to[i]] = dis[x] + cost[i] , from[to[i]] = x , pre[to[i]] = i , q.push(to[i]); } return ~from[t];}int mincost(){ int ans = 0 , i , k; while(spfa()) { k = inf; for(i = t ; i != s ; i = from[i]) k = min(k , val[pre[i]]); ans += dis[t] * k; for(i = t ; i != s ; i = from[i]) val[pre[i]] -= k , val[pre[i] ^ 1] += k; } return ans;}int main(){ int n , m , p , i , j , k , x , y , z; scanf("%d%d%d" , &n , &m , &p); memset(map , 0x3f , sizeof(map)); while(m -- ) scanf("%d%d%d" , &x , &y , &z) , map[x][y] = map[y][x] = min(map[x][y] , z); for(k = 0 ; k <= n ; k ++ ) for(i = 0 ; i <= n ; i ++ ) for(j = 0 ; j <= n ; j ++ ) if((k <= i || k <= j) && map[i][j] > map[i][k] + map[k][j]) map[i][j] = map[i][k] + map[k][j]; s = 2 * n + 2 , t = 2 * n + 3 , add(2 * n + 1 , 0 , p , 0); for(i = 1 ; i <= n ; i ++ ) add(0 , i , inf , map[0][i]); for(i = 1 ; i <= n ; i ++ ) { add(s , i + n , 1 , 0) , add(i , t , 1 , 0) , add(i + n , 2 * n + 1 , inf , 0); for(j = i + 1 ; j <= n ; j ++ ) if(map[i][j] != inf) add(i + n , j , inf , map[i][j]); } printf("%d\n" , mincost()); return 0;}

 

转载于:https://www.cnblogs.com/GXZlegend/p/6832504.html

你可能感兴趣的文章
TYVJ.1864.[Poetize I]守卫者的挑战(概率DP)
查看>>
0925 韩顺平java视频
查看>>
iOS-程序启动原理和UIApplication
查看>>
mysql 8.0 zip包安装
查看>>
awk 统计
查看>>
模板设计模式的应用
查看>>
实训第五天
查看>>
平台维护流程
查看>>
2012暑期川西旅游之总结
查看>>
12010 解密QQ号(队列)
查看>>
2014年辛星完全解读Javascript第一节
查看>>
装配SpringBean(一)--依赖注入
查看>>
java选择文件时提供图像缩略图[转]
查看>>
方维分享系统二次开发, 给评论、主题、回复、活动 加审核的功能
查看>>
Matlab parfor-loop并行运算
查看>>
string与stringbuilder的区别
查看>>
2012-01-12 16:01 hibernate注解以及简单实例
查看>>
iOS8统一的系统提示控件——UIAlertController
查看>>
PAT甲级——1101 Quick Sort (快速排序)
查看>>
python创建进程的两种方式
查看>>