[ZJOI2011]营救皮卡丘（费用流）

2324: [ZJOI2011]营救皮卡丘

Description

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

火箭队一共有$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$个人所经过的道路的长度总和最少。

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

Input

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

接下来$M$行，每行三个非负整数，第$i$行的整数为$A_i$，$B_i$，$L_i$。表示存在一条从$A_i$号据点到$B_i$号据点的长度为$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$≤ 150,$M$≤ 20 000, 1 ≤$K$≤ 10, $L_i$ ≤ 10 000, 保证小智一行一定能够救出皮卡丘。至于为什么$K$≤ 10，你可以认为最终在小智的号召下，小智，小霞，小刚，小建，小遥，小胜，小光，艾莉丝，天桐，还有去日本旅游的黑猫警长，一同前去大战火箭队。

解：

一个很重要的条件就是要按顺序摧毁所有节点。
那我们单看一个人的话他不可能从编号大的节点走向编号小的，那么我们就把这个图转化成拓扑图了。先floyd预处理一下，从$i$到$j$走的编号不大于$i$或$j$的最短路。
我们现在需要找$k$条路径使其覆盖所有点并且经过的边权和最小。上费用流模板，一个点拆两个，有路径连边，源点向$i$连流量为1的边，$i'$向汇点连流量为1的边。因为最后一定会满流，不用担心会不合法。特殊的，源点向0号点连流量为$k$的边，限制有$k$个人。因为有传递闭包，不会有人走重复路径。

floyd那里的技巧：从小到大加点，如果小于$i$或$j$才更新。

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);
}