【BZOJ 1070】[SCOI2007]修车 费用流

就是拆个点限制一下(两点一排一大片），这道题让我注意到了限制这个重要的词。我们跑网络流跑出来的图都是有一定意义的，一般这个意义就对应了问题的一种方案，一般情况下跑一个不知道对不对的方案是相对容易的我们要做的就是尽可能去限制他，使他正确，有的时候（特别是费用流）我们也要使他计算答案更加简便。

#include <cstdio>
#include <cstring>
#include <algorithm>
#define car(a) (a)
#define tech(a,b) ((a)*n+(b))
const int N=;
const int M=;
const int P=N*M;
const int E=N*N*M*;
const int Inf=0x3f3f3f3f;
struct V{
  int to,next,c,f;
}c[E];
int head[P],t=;
inline void add(int x,int y,int z,int _){
  c[++t].to=y,c[t].f=z,c[t].c=_,c[t].next=head[x],head[x]=t;
}
int dis[P],anc[P],S,T,n,m;
int q[P],front,back;
bool in[P];
int ans;
inline bool spfa(){
  memset(dis,0x3f,sizeof(dis));
  dis[S]=,in[S]=true,q[back++]=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&&c[i].c+dis[x]<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(){
  scanf("%d%d",&m,&n);
  S=n*m+n+,T=S+;
  for(int i=;i<=n;++i){
    add(S,car(i),,);
    add(car(i),S,,);
  }
  for(int i=,x;i<=n;++i)
    for(int j=;j<=m;++j){
      scanf("%d",&x);
      for(int k=;k<=n;++k){
        add(car(i),tech(j,k),,k*x);
        add(tech(j,k),car(i),,-k*x);
      }
      add(tech(j,i),T,,);
      add(T,tech(j,i),,);
    }
  while(spfa())ans+=shoot();
  printf("%.2f",(double)ans/(double)n);
  return ;
}