【BZOJ 2324】[ZJOI2011]营救皮卡丘 费用流

本人实行诱骗拐卖（利用自然分层与实际意义），正解拼接补充（充分利用最大流限制（不浪费任何一个走出去的机会而不是不浪费任何一个已有的流）与问题转换）

#include <cstdio>
#include <cstring>
#include <algorithm>
const int N=;
const int M=;
const int K=;
const int Inf=0x3f3f3f3f;
const int treat=;
const int P=N*;
const int E=N*N*;
int n,m,k;
int via[N][N];
struct V{
  int to,next,f,c;
}c[E];
int head[P],t=;
inline void add(int x,int y,int z,int _){
  c[++t].to=y,c[t].next=head[x],c[t].f=z,c[t].c=_,head[x]=t;
}
int dis[P],anc[P];
int q[P],front,back;
bool in[P];
int S,T;
inline bool spfa(){
  memset(dis,0x3f,sizeof(dis));
  in[S]=true,q[back++]=S,dis[S]=;
  if(back==P)back=;
  while(front!=back){
    int x=q[front++];in[x]=false;
    if(front==P)front=;
    for(int i=head[x];i;i=c[i].next)
      if(c[i].f&&dis[x]+c[i].c<dis[c[i].to]){
        dis[c[i].to]=dis[x]+c[i].c,anc[c[i].to]=i;
        if(!in[c[i].to]){
          q[back++]=c[i].to,in[c[i].to]=true;
          if(back==P)back=;
        }
      }
  }
  return dis[T]!=Inf;
}
inline int shoot(){
  int f=Inf;
  for(int i=anc[T];i;i=anc[c[i^].to])f=std::min(f,c[i].f);
  for(int i=anc[T];i;i=anc[c[i^].to])c[i].f-=f,c[i^].f+=f;
  return f*dis[T];
}
int main(){
  memset(via,0x3f,sizeof(via));
  scanf("%d%d%d",&n,&m,&k);
  ++n;
  for(int i=,x,y,z;i<=m;++i){
    scanf("%d%d%d",&x,&y,&z);
    ++x,++y;
    via[x][y]=via[y][x]=std::min(via[x][y],z);
  }
  S=n+n+,T=n+n+;
  for(int i=;i<=n;++i){
    via[i][i]=;
    add(i,i+n,,-treat),add(i+n,i,,treat);
    add(i,i+n,k,),add(i+n,i,,);
    add(i+n,T,k,),add(T,i+n,,);
  }
  for(int a=;a<=n;++a)
    for(int b=;b<=n;++b)
      via[a][b]=std::min(via[a][b],via[a][]+via[][b]);
  for(int i=;i<=n;++i){
    for(int a=;a<=n;++a)
      for(int b=;b<=n;++b)
        via[a][b]=std::min(via[a][b],via[a][i]+via[i][b]);
    for(int j=;j<i;++j){
      add(j+n,i,k,via[j][i]);
      add(i,j+n,,-via[i][j]);
    }
  }
  add(S,,k,);
  add(,S,,);
  int ans=;
  while(spfa())ans+=shoot();
  ans+=n*treat;
  printf("%d",ans);
  return ;
}