bzoj2127(还没过+最小割)

这题对我来说还是太难了

利用最小割考虑。

对于原图中所有相邻的两个人A，B，我们建边：

s->A:cost[A文]+c[文][A][B]/2,s->B:cost[B文]+c[文][A][B]/2;

A->t:cost[A理]+c[理][A][B]/2,B->t:costB[理]+c[理][A][B]/2;

A<–>B:c[文][A][B]/2+c[理][A][B]/2

这样会出现两种割，分别对应两种相同，一种选文一种选理。

2016.8.26更新

首先这样的模型显然是最小割（二者选其一问题），实际上只要是做过在看还是可以产生思路的。

#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<queue>
#define rep(n,m) for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) 
#define FOR for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) 
using namespace std;
const int inf=0x3f3f3f3f;
int n,m,sum;
int a[105][105],b[105][105];
int id[105][105],S,T,cnt;
int head[105*105],tot;
struct aa
{
 int to,pre,flow,cap;
}edge[300005];
void addedge(int u,int v,int c)
{
 edge[++tot].to=v;edge[tot].pre=head[u];edge[tot].cap=c;head[u]=tot;
 edge[++tot].to=u;edge[tot].pre=head[v];edge[tot].cap=0;head[v]=tot;
}
int lev[105*105],cur[105*105];
bool bfs()
{
 memset(lev,0,sizeof(lev));lev[S]=1;
 queue<int> q;q.push(S);
 while (!q.empty())
 {
  int u=q.front();q.pop();
  for (int i=head[u];i;i=edge[i].pre)
  if (edge[i].cap>edge[i].flow&&lev[edge[i].to]==0)
  {
   lev[edge[i].to]=lev[u]+1;
   if (edge[i].to==T) return true;
   q.push(edge[i].to);
  }
 }
 return false;
}
int dfs(int u,int maxflow)
{
 if (u==T||maxflow==0) return maxflow;
 int ans=0; 
 for (int &i=cur[u];i;i=edge[i].pre)
 if (lev[edge[i].to]==lev[u]+1)
 {
  int flow=dfs(edge[i].to,min(edge[i].cap-edge[i].flow,maxflow));
  ans+=flow;
  maxflow-=flow;
  edge[i].flow+=flow;
  edge[((i-1)^1)+1].flow-=flow;
  if (maxflow==0) return ans;
 } 
 return ans;
}
 
int dinic()
{
 int ans=0;
 while (bfs())
 {
  for (int i=1;i<=cnt;i++) cur[i]=head[i];
  ans+=dfs(S,inf);
 }
 return ans;
}

void build()
{
 int x;
 rep(n-1,m)
 {
  scanf("%d",&x);
  a[i][j]+=x;sum+=x;a[i+1][j]+=x;
  addedge(id[i][j],id[i+1][j],x);
  addedge(id[i+1][j],id[i][j],x);
 }
 rep(n-1,m)
 {
  scanf("%d",&x);
  b[i][j]+=x;sum+=x;b[i+1][j]+=x;
  addedge(id[i][j],id[i+1][j],x);
  addedge(id[i+1][j],id[i][j],x);
 }
 rep(n,m-1)
 {
  scanf("%d",&x);
  a[i][j]+=x;sum+=x;a[i][j+1]+=x;
  addedge(id[i][j],id[i][j+1],x);
  addedge(id[i][j+1],id[i][j],x);
 }
 rep(n,m-1)
 {
  scanf("%d",&x);
  b[i][j]+=x;sum+=x;b[i][j+1]+=x;
  addedge(id[i][j],id[i][j+1],x);
  addedge(id[i][j+1],id[i][j],x);
 }
 rep(n,m)
 {
  addedge(S,id[i][j],a[i][j]);
  addedge(id[i][j],T,b[i][j]);
 }
}
int main()
{
 scanf("%d%d",&n,&m);
 FOR scanf("%d",&a[i][j]),sum+=a[i][j],a[i][j]<<=1;
 FOR scanf("%d",&b[i][j]),sum+=b[i][j],b[i][j]<<=1;
 FOR id[i][j]=++cnt;S=++cnt,T=++cnt;
 build();
 printf("%d",sum-(dinic()>>1));
 return 0;
}

总结：

1：最小割问题所代表的意义，是不像最大流问题的流那么形象化，代表一种方案，这道题中出现了/2的情况，并不平常（做多了就见怪不怪了，还是我见的少），就是通过割的思想来处理的并不是用流的方法。学长的经验就是最小割做题是绝对不能带着流的思维的！，要用割的思维！

2：这道题的实现技巧也很巧妙，因为出现了/2的情况，但是/2就会有单双数的问题要处理，如果转化为浮点数，真的是没法做，所以，解决办法就是一上来把所有的数据都乘以2，最后答案再除以2，以此让答案更加精确。