[最小权路径覆盖 & 网络流] BZOJ2324 :[ZJOI2011] 营救皮卡丘

做了那道80人环游世界，感觉这题也是一个套路…

自己很naive地把每个点拆成三个点建图…

其实重构一下图就很简单了

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <queue>
#include <cstring>

using namespace std;

const int N=550,inf=1<<29;

int n,m,k,S,T,ss,tt,cnt=1,G[N];
struct edge{
  int t,nx,w,f;
}E[N*N*10];

inline void link(int x,int y,int f,int w){
  E[++cnt].t=y; E[cnt].nx=G[x]; G[x]=cnt; E[cnt].f=f; E[cnt].w=w;
  E[++cnt].t=x; E[cnt].nx=G[y]; G[y]=cnt; E[cnt].f=0; E[cnt].w=-w;
}

queue<int> Q;
int vis[N],dis[N],frm[N],ans;

inline bool spfa(){
  for(int i=0;i<=T;i++) dis[i]=inf,vis[i]=0;
  vis[S]=1; Q.push(S); dis[S]=0;
  while(!Q.empty()){
    int x=Q.front(); vis[x]=0; Q.pop();
    for(int i=G[x];i;i=E[i].nx)
      if(E[i].f && dis[E[i].t]>dis[x]+E[i].w){
    dis[E[i].t]=dis[x]+E[i].w; frm[E[i].t]=i;
    if(!vis[E[i].t])
      vis[E[i].t]=1,Q.push(E[i].t);
      }
  }
  if(dis[T]==inf) return false;
  int mf=inf;
  for(int i=T;i!=S;i=E[frm[i]^1].t) mf=min(mf,E[frm[i]].f);
  ans+=mf*dis[T];
  for(int i=T;i!=S;i=E[frm[i]^1].t) E[frm[i]].f-=mf,E[frm[i]^1].f+=mf;
  return true;
}

int dist[N][N];

int main(){
  scanf("%d%d%d",&n,&m,&k);
  ss=n+n+1; tt=n+n+2; S=tt+1; T=S+1;
  for(int i=0;i<=n;i++)
    for(int j=0;j<=n;j++)
      dist[i][j]=inf;
  for(int i=1,a,b,c;i<=m;i++){
    scanf("%d%d%d",&a,&b,&c);
    dist[a][b]=dist[b][a]=min(dist[a][b],c);
  }
  for(int t=0;t<=n;t++)
    for(int i=0;i<=n;i++)
      for(int j=0;j<=n;j++)
    if(t<i || t<j) dist[i][j]=min(dist[i][j],dist[i][t]+dist[t][j]);
  for(int i=1;i<=n;i++){
    link(ss,i,inf,dist[0][i]); link(i+n,tt,inf,0);
    link(S,i+n,1,0); link(i,T,1,0); link(i,i+n,inf,0);
  }
  for(int i=1;i<=n;i++)
    for(int j=i+1;j<=n;j++)
      if(dist[i][j]!=inf) link(i+n,j,inf,dist[i][j]);
  link(tt,ss,k,0);
  while(spfa());
  printf("%d\n",ans);
  return 0;
}