球的移动(move)

有n个盒子(1<=N<=1000)围成一个圈,每个盒子有ai个球,所有盒子的球的总数小于等于n.每一次移动,可以把一个球移动到它的一个相邻的盒子内.
现在要使得每个盒子的球数<=1,求最少的移动次数

输入格式：

n
a1
a2
…
an

输出格式：

最少移动次数

样例输入：

12
0
0
2
4
3
1
0
0
0
0
0
1

样例输出：

19

网络流真乃xjb题的解题法宝

费用流。

先设超级源点0和超级汇点2*n+1，然后将每个点拆成源点和汇点，然后每个源点先与对应汇点连一条容量无限大，费用为0的边；然后与相邻点的汇点连容量无限大，费用为1的边，最后将各源点与超级源点连一条容量为a[i]费用为0的边，各汇点与超级汇点连一条容量为1费用为0的边，然后跑一遍费用流即可。（为什么？自己思考吧）

 

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
struct edge{
    int v,c,f;
}e[];
int fi[],ne[],la[],q[],d[],flow[],x,ansflow,anscost,tot,pe[],n,m,s,t,b[],pre[],u,v,w,f;
void add(int x,int y,int f,int c){
    e[++tot].v=y;e[tot].c=c;e[tot].f=f;
    !fi[x]?fi[x]=tot:ne[la[x]]=tot;la[x]=tot;
    e[++tot].v=x;e[tot].c=-c;e[tot].f=;
    !fi[y]?fi[y]=tot:ne[la[y]]=tot;la[y]=tot;
}
bool spfa(int s,int t){
    int H=,T=;
    memset(d,0x3f,sizeof(d));
    memset(flow,,sizeof(flow));
    q[++T]=s;b[s]=;d[s]=;pre[s]=;flow[s]=0x7fffffff;
    while(H<T){
        int x=q[++H];
        for(int i=fi[x];i;i=ne[i]){
            int y=e[i].v,z=e[i].f,c=e[i].c;
            if(z>&&d[y]>d[x]+c){
                d[y]=d[x]+c;
                pre[y]=x;
                pe[y]=i;
                flow[y]=min(flow[x],z);
                if(!b[y])b[y]=,q[++T]=y;
            }
        }

        b[x]=;
    }
    return d[t]!=0x3f3f3f3f;
}
void maxflow(){
    while(spfa(s,t)){
        int k=t;
        while(k!=s){
            e[pe[k]].f-=flow[t];
            e[pe[k]^].f+=flow[t];
            k=pre[k];
        }
        ansflow+=flow[t];
        anscost+=d[t]*flow[t];
    }
}
int main(){
    scanf("%d",&n);
    tot=;
    for(int i=;i<=n;i++){
        scanf("%d",&x);
        int pr=i==?n:i-,nx=i==n?:i+;
        add(,i,x,);
        add(i,nx+n,0x3f3f3f3f,);
        add(nx+n,i,0x3f3f3f3f,);
        add(i,pr+n,0x3f3f3f3f,);
        add(pr+n,i,0x3f3f3f3f,);
        add(i,i+n,0x3f3f3f3f,);
        add(i+n,i,0x3f3f3f3f,);
        add(i+n,*n+,,);
    }
    s=;t=*n+;
    maxflow();
    printf("%d",anscost);
}