bzoj 1927 星际竞速 —— 最小费用最大流

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1927

首先注意到这是个DAG；

考虑每个点从哪里来，可以是瞬移来的，也可以是从某个点走过来的，而从每个点走出去只能用一次；

所以拆点，i 表示从这个点走出去，n+i 表示来到这个点；

建图：

1.瞬移：S 向 n+i 连边权 a[i]，流量1的边

2.走过来：如果 i 能走到 j，那么 i 向 n+j 连边权 w，流量1的边

然后 S 向 i 连边权0，流量1的边，表示一个点只能走出去一次；

n+i 向 T 连边权0，流量1的边，表示一个点只能到达一次；

因为是DAG，所以不会有环，而且一定会有一个点瞬移开始；

建图真妙啊！

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
int const xn=,xm=4e4+,inf=0x3f3f3f3f;
int n,m,S,T,hd[xn],ct=,to[xm],nxt[xm],w[xm],c[xm],dis[xn],inc[xn],pre[xn];
bool vis[xn];
queue<int>q;
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return f?ret:-ret;
}
void ade(int x,int y,int z,int f){to[++ct]=y; nxt[ct]=hd[x]; hd[x]=ct; w[ct]=z; c[ct]=f;}
void add(int x,int y,int z,int f){ade(x,y,z,f); ade(y,x,-z,);}
bool bfs()
{
  memset(dis,0x3f,sizeof dis);
  dis[S]=; q.push(S); vis[S]=; inc[S]=inf;//inc
  while(q.size())
    {
      int x=q.front(); q.pop(); vis[x]=;
      for(int i=hd[x],u;i;i=nxt[i])
    if(dis[u=to[i]]>dis[x]+w[i]&&c[i])
      {
        dis[u]=dis[x]+w[i]; pre[u]=i;
        inc[u]=min(inc[x],c[i]);
        if(!vis[u])vis[u]=,q.push(u);
      }
    }
  return dis[T]!=inf;
}
void upd()
{
  int x=T;
  while(x!=S)
    {
      int i=pre[x];
      c[i]-=inc[T]; c[i^]+=inc[T];
      x=to[i^];
    }
}
int main()
{
  n=rd(); m=rd(); S=; T=(n<<)+;
  for(int i=,x;i<=n;i++)
    x=rd(),add(S,i,,),add(S,n+i,x,),add(n+i,T,,);
  for(int i=,x,y,z;i<=m;i++)
    {
      x=rd(); y=rd(); z=rd();
      if(x>y)swap(x,y);
      add(x,n+y,z,);//n+y
    }
  int ans=;
  while(bfs())ans+=dis[T]*inc[T],upd();
  printf("%d\n",ans);
  return ;
}