2324: [ZJOI2011]营救皮卡丘
Description
皮卡丘被火箭队用邪恶的计谋抢走了!这三个坏家伙还给小智留下了赤果果的挑衅!为了皮卡丘,也为了正义,小智和他的朋友们义不容辞的踏上了营救皮卡丘的道路。
火箭队一共有 N N 个据点,据点之间存在条双向道路。据点分别从1到 N N 标号。小智一行人从真新镇出发,营救被困在 N N 号据点的皮卡丘。为了方便起见,我们将真新镇视为0号据点,一开始个人都在0号点。
由于火箭队的重重布防,要想摧毁K号据点,必须按照顺序先摧毁1到 K K -1号据点,并且,如果-1号据点没有被摧毁,由于防御的连锁性,小智一行任何一个人进入据点 K K ,都会被发现,并产生严重后果。因此,在-1号据点被摧毁之前,任何人是不能够经过 K K 号据点的。
为了简化问题,我们忽略战斗环节,小智一行任何一个人经过号据点即认为 K K 号据点被摧毁。被摧毁的据点依然是可以被经过的。
个人是可以分头行动的,只要有任何一个人在 K K -1号据点被摧毁之后,经过号据点, K K 号据点就被摧毁了。显然的,只要号据点被摧毁,皮卡丘就得救了。
野外的道路是不安全的,因此小智一行希望在摧毁 N N 号据点救出皮卡丘的同时,使得个人所经过的道路的长度总和最少。
请你帮助小智设计一个最佳的营救方案吧!
Input
第一行包含三个正整数 N N ,, K K 。表示一共有+1个据点,分别从0到 N N 编号,以及条无向边。一开始小智一行共 K K 个人均位于0号点。
接下来行,每行三个非负整数,第 i i 行的整数为, Bi B i , Li L i 。表示存在一条从 Ai A i 号据点到 Bi B i 号据点的长度为 Li L i 的道路。
Output
仅包含一个整数S,为营救皮卡丘所需要经过的最小的道路总和。
Sample Input
3 4 2
0 1 1
1 2 1
2 3 100
0 3 1
Sample Output
3
样例说明
小智和小霞一起前去营救皮卡丘。在最优方案中,小智先从真新镇前往1号点,接着前往2号据点。当小智成功摧毁2号据点之后,小霞从真新镇出发直接前往3号据点,救出皮卡丘。
HINT
对于100%的数据满足 N N ≤ 150,≤ 20 000, 1 ≤ K K ≤ 10, ≤ 10 000, 保证小智一行一定能够救出皮卡丘。至于为什么 K K ≤ 10,你可以认为最终在小智的号召下,小智,小霞,小刚,小建,小遥,小胜,小光,艾莉丝,天桐,还有去日本旅游的黑猫警长,一同前去大战火箭队。
解:
一个很重要的条件就是要按顺序摧毁所有节点。
那我们单看一个人的话他不可能从编号大的节点走向编号小的,那么我们就把这个图转化成拓扑图了。先floyd预处理一下,从到
j
j
走的编号不大于或
j
j
的最短路。
我们现在需要找条路径使其覆盖所有点并且经过的边权和最小。上费用流模板,一个点拆两个,有路径连边,源点向
i
i
连流量为1的边,向汇点连流量为1的边。因为最后一定会满流,不用担心会不合法。特殊的,源点向0号点连流量为
k
k
的边,限制有个人。因为有传递闭包,不会有人走重复路径。
floyd那里的技巧:从小到大加点,如果小于 i i 或才更新。
code:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
struct lxy{
int to,flow,w,next;
}eg[30005];
int dis[155][155];
int n,m,k,s,t,a1,a2,a3,ans,cnt=-1;
int head[365],layer[365];
bool vis[365];
queue <int> d;
void add(int op,int ed,int flow,int w){
eg[++cnt].next=head[op];
eg[cnt].to=ed;
eg[cnt].flow=flow;
eg[cnt].w=w;
head[op]=cnt;
}
bool bfs(){
memset(vis,0,sizeof(vis));
memset(layer,0x3f,sizeof(layer));
d.push(s);layer[s]=0;vis[s]=1;
while(!d.empty()){
int now=d.front();d.pop();vis[now]=0;
for(int i=head[now];i!=-1;i=eg[i].next)
if(eg[i].flow!=0&&layer[eg[i].to]>layer[now]+eg[i].w){
layer[eg[i].to]=layer[now]+eg[i].w;
if(vis[eg[i].to]==0){
vis[eg[i].to]=1;
d.push(eg[i].to);
}
}
}
if(layer[t]==0x3f3f3f3f) return false;
return true;
}
int dfs(int u,int a){
if(u==t||a==0) return a;
vis[u]=1;
int f,flow=0;
for(int i=head[u];i!=-1;i=eg[i].next)
if(eg[i].flow!=0&&vis[eg[i].to]==0&&layer[eg[i].to]==layer[u]+eg[i].w){
f=dfs(eg[i].to,min(a,eg[i].flow));
eg[i].flow-=f;eg[i^1].flow+=f;
a-=f;flow+=f;
ans+=f*eg[i].w;
if(a==0) break;
}
return flow;
}
int dinic(){
int ret=0;
while(bfs()) ret+=dfs(s,0x7f7f7f7f);
return ret;
}
int main()
{
memset(head,-1,sizeof(head));
memset(dis,0x3f,sizeof(dis));
scanf("%d%d%d",&n,&m,&k);n++;
s=0,t=2*n+1;
for(int i=1;i<=m;i++){
scanf("%d%d%d",&a1,&a2,&a3);a1++,a2++;
dis[a1][a2]=min(dis[a1][a2],a3);
dis[a2][a1]=min(dis[a2][a1],a3);
}
for(int p=1;p<=n;p++)
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(p<=j||p<=i)
dis[i][j]=min(dis[i][j],dis[i][p]+dis[p][j]);
for(int i=1;i<n;i++)
for(int j=i+1;j<=n;j++)
if(dis[i][j]!=0x3f3f3f3f){
add(i+n,j,0x3f3f3f3f,dis[i][j]);
add(j,i+n,0,-dis[i][j]);
}
add(s,1+n,k,0);add(1+n,s,0,0);
for(int i=2;i<=n;i++){
add(s,i+n,1,0);add(i+n,s,0,0);
add(i,t,1,0);add(t,i,0,0);
}
dinic();
printf("%d",ans);
}