BZOJ 3511: 土地划分

题意：Y国有N座城市，并且有M条双向公路将这些城市连接起来，并且任意两个城市至少有一条路径可以互达。
Y国的国王去世之后，他的两个儿子A和B都想成为新的国王，但他们都想让这个国家更加安定，不会用武力解决问题。
于是他们想将这个国家分成两个小国家A国和B国。现在，A拥有1号城市，B拥有N号城市，其他的城市还尚未确定归属哪边（划分之后的国家内部城市可以不连通）。
由于大家都想让国家变得更好，而某些城市的人民愿意国王的A儿子作为他们的领袖，而某些城市更看好B，而为了交通的便捷，如果划分后的公路连接两个同一个国家的城市，那么更利于城市之间的交流。于是大臣们设计了一种对土地划分的评分机制，具体如下：
1. 对于城市i，如果它划分给A国，将得到VA[i]的得分；划分给B国，将得到VB[i]的得分。
2. 对于一条公路i，如果它连接两个A国的城市，将得到EA[i]的得分；连接两个B国的城市，将得到EB[i]的得分；否则，这条公路将失去意义，将扣除EC[i]的得分。
现请你找到最优的土地划分，使得这种它的评分最高。

很明显的二元关系的最小割，新建S,T，对于每个点，由S向I连VA[I],由I向T连VB[I].对于每条边，由S向x,y连EA/2,由T向x,y连EB/2.x,y互连EA/2+EB/2+EC，特别地，S向1连INF，N向T连INF.我们先将所有可能的收益加起来，网络流中的边权代表扣除的收益，这样跑最小割之后就是剩余的最大收益。
Tips:先将所有权值乘上2，最后除以2，这样避免实数运算。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<queue>
using namespace std;
const int maxn=10000+10;
const int INF=1000000000;
int n,m,s,t,st[maxn],ed[maxn],ans,cur[maxn],vis[maxn],dis[maxn];
struct edge
{
  int from,to,cap,flow;
};
vector<int> g[maxn];
vector<edge> edges;
void addedge(int from,int to,int cap)
{
  edges.push_back((edge){from,to,cap,0});
  edges.push_back((edge){to,from,0,0});
  int m1=edges.size();
  g[from].push_back(m1-2);
  g[to].push_back(m1-1);
}
bool bfs()
{
  memset(vis,0,sizeof(vis));
  queue<int> Q;Q.push(s);vis[s]=1;dis[s]=0;
  while(!Q.empty())
  {
    int x=Q.front();Q.pop();
    for(int i=0;i<g[x].size();i++)
    {
      edge e=edges[g[x][i]];
      if(e.cap>e.flow&&!vis[e.to])
      {
        vis[e.to]=1;
        dis[e.to]=dis[x]+1;
        Q.push(e.to);
      }
    }
  }
  return vis[t];
}
int dfs(int x,int a)
{
  if(x==t||a==0) return a;
  int flow=0,f;
  for(int &i=cur[x];i<g[x].size();i++)
  {
    edge &e=edges[g[x][i]];
    if(dis[e.to]==dis[x]+1&&(f=dfs(e.to,min(a,e.cap-e.flow)))>0)
    {
      flow+=f;
      e.flow+=f;
      edges[g[x][i]^1].flow-=f;
      a-=f;
      if(a==0) break;
    }
  }
  return flow;
}
int maxflow()
{
  int flow=0;
  while(bfs())
  {
    memset(cur,0,sizeof(cur));
    flow+=dfs(s,INF);
  }
  return flow;
}
int main()
{
  //freopen("3511.in","r",stdin);
  //freopen("3511.out","w",stdout);
  scanf("%d%d",&n,&m);s=0,t=n+1;st[1]=INF;ed[n]=INF;
  for(int i=2;i<=n-1;i++)
  {
    int x;scanf("%d",&x);
    st[i]+=2*x;
    ans+=2*x;
  }
  for(int i=2;i<=n-1;i++)
  {
    int x;scanf("%d",&x);
    ed[i]+=2*x;ans+=2*x;
  }
  for(int i=1;i<=m;i++)
  {
    int x,y,ea,eb,ec;
    scanf("%d%d%d%d%d",&x,&y,&ea,&eb,&ec);
    ans+=(ea+eb)*2;st[x]+=ea;st[y]+=ea;ed[x]+=eb;ed[y]+=eb;
    addedge(x,y,ea+eb+2*ec);addedge(y,x,ea+eb+2*ec);
  }
  for(int i=1;i<=n;i++) 
  {
    addedge(s,i,st[i]);
    addedge(i,t,ed[i]);
  }
  printf("%d\n",(ans-maxflow())/2);
  return 0;
}