标签:path char AC pat amp 平均值 turn tar complex
要求最后所有仓库的库存相同,那么就可以想到求出这个值,也就是仓库库存量的平均值delta。
如果某一个仓库的库存量x大于这个平均值,那么就从源点向这个店连一条流量为x-delta的边,表示这个仓库可以向外搬运的货物数量,费用为0。
反之,如果某一个仓库的存货量小于x这个平均值,就从这个店向汇点连一条流量为delta-x的边,表示这个仓库需要运来的货物数量,费用为0。
再把相邻的仓库互相连边,流量为INF(如果这个仓库有多余的货物,就可以一直往外搬),单位费用为1。
最后跑费用流就可以了。
#include<complex> #include<cstdio> #include<queue> #include<cstring> using namespace std; const int INF=0x3f3f3f3f; const int N=201,M=1001; struct node{ int v,f,w,nxt; }e[M<<1]; int n,s,t,Enum=1; long long num; int front[N],dis[N],path[N],a[N],incf[N]; bool vis[N]; queue<int>q; int qread() { int x=0; char ch=getchar(); while(ch<‘0‘ || ch>‘9‘)ch=getchar(); while(ch>=‘0‘ && ch<=‘9‘){x=x*10+ch-‘0‘;ch=getchar();} return x; } void Insert(int u,int v,int f,int w) { e[++Enum].v=v;e[Enum].f=f;e[Enum].w=w;e[Enum].nxt=front[u];front[u]=Enum; e[++Enum].v=u;e[Enum].f=0;e[Enum].w=-w;e[Enum].nxt=front[v];front[v]=Enum; } bool SPFA() { memset(dis,0x3f,sizeof(dis)); memset(vis,0,sizeof(vis)); dis[s]=0;q.push(s); incf[s]=INF; int u,v; while(!q.empty()) { u=q.front();q.pop(); vis[u]=0; for(int i=front[u];i;i=e[i].nxt) { v=e[i].v; if(e[i].f && dis[v]>dis[u]+e[i].w) { dis[v]=dis[u]+e[i].w; incf[v]=min(incf[u],e[i].f); path[v]=i; if(!vis[v]) { vis[v]=1; q.push(v); } } } } return dis[t]<INF; } int CostFlow() { int res=0; while(SPFA()) { int x=t; while(x!=s) { int i=path[x]; e[i].f-=incf[t]; e[i^1].f+=incf[t]; x=e[i^1].v; } res+=dis[t]*incf[t]; } return res; } int main() { scanf("%d",&n); s=0;t=n+1; int x; for(int i=1;i<=n;i++) { a[i]=qread(); num+=a[i]; } num/=n; for(int i=1;i<=n;i++) if(a[i]-num>0)Insert(s,i,a[i]-num,0); else if(a[i]-num<0)Insert(i,t,num-a[i],0); Insert(1,2,INF,1);Insert(1,n,INF,1); Insert(n,n-1,INF,1);Insert(n,1,INF,1); for(int i=2;i<n;i++) { Insert(i,i+1,INF,1); Insert(i,i-1,INF,1); } printf("%d\n",CostFlow()); return 0; }
标签:path char AC pat amp 平均值 turn tar complex
原文地址:https://www.cnblogs.com/LeTri/p/9033515.html