bzoj 2879: [Noi2012]美食节

最新推荐文章于 2019-01-05 09:57:00 发布

weixin_33704591

最新推荐文章于 2019-01-05 09:57:00 发布

阅读量48

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/xydddd/p/5309015.html

本文介绍了一种基于最小费用流算法的实现方法，该方法能够在动态加点的情况下求解最短路径问题。通过定义一系列的数据结构和操作流程，如节点连接、宽度优先搜索等，最终实现了从源点到汇点的最小费用流计算。

  1 #include<iostream>
  2 #include<cstring>
  3 #include<cstdio>
  4 #include<queue>
  5 #define M 100005
  6 #define S 100001
  7 using namespace std;
  8 int n,m,head[M],next[30*M],u[30*M],v[30*M],c[30*M],qi[30*M],a[M],tot,t[45][106],cnt=1,ans,d[M],f[M],from[M];
  9 void jia(int a1,int a2,int a3,int a4)
 10 {
 11     cnt++;
 12     qi[cnt]=a1;
 13     u[cnt]=a2;
 14     v[cnt]=a3;
 15     c[cnt]=a4;
 16     next[cnt]=head[a1];
 17     head[a1]=cnt;
 18 }
 19 bool bfs()
 20 {
 21     for(int i=1;i<=S;i++)
 22       d[i]=0x7fffffff;
 23     queue<int>q;
 24     q.push(0);
 25     f[0]=1;
 26     for(;!q.empty();)
 27       {
 28         int p=q.front();
 29         q.pop();
 30         f[p]=0;
 31         for(int i=head[p];i;i=next[i])
 32           if(v[i]&&d[u[i]]>d[p]+c[i])
 33             {
 34                 d[u[i]]=d[p]+c[i];
 35                 from[u[i]]=i;
 36                 if(!f[u[i]])
 37                   {
 38                     q.push(u[i]);
 39                     f[u[i]]=1;
 40                   }
 41             }
 42       }
 43     if(d[S]==0x7fffffff)
 44       return 0;
 45     else
 46       return 1;
 47 }
 48 void wang()
 49 {
 50     int ma=0x7fffffff,x,y;
 51     for(int i=from[S];i;i=from[qi[i]])
 52       {
 53         ma=min(ma,v[i]);
 54         if(qi[i]==0)
 55           {
 56             int a1=u[i];
 57             x=(a1-1)/tot+1;
 58             y=a1%tot+1;
 59           }
 60       }
 61     for(int i=from[S];i;i=from[qi[i]])
 62       {
 63         v[i]-=ma;
 64         v[i^1]+=ma;
 65         ans+=c[i]*ma;
 66       }
 67     for(int i=1;i<=n;i++)
 68       {
 69         jia((x-1)*tot+y,m*tot+i,1,y*t[i][x]);
 70         jia(m*tot+i,(x-1)*tot+y,0,-y*t[i][x]);
 71       }
 72     return;
 73 }
 74 int main()
 75 {
 76     scanf("%d%d",&n,&m);
 77     for(int i=1;i<=n;i++)
 78       {
 79         scanf("%d",&a[i]);
 80         tot+=a[i];
 81       }
 82     for(int i=1;i<=n;i++)
 83       for(int j=1;j<=m;j++)
 84         scanf("%d",&t[i][j]);
 85     for(int i=1;i<=m*tot;i++)
 86       {
 87         jia(0,i,1,0);
 88         jia(i,0,0,0);
 89       }
 90     for(int i=1;i<=n;i++)
 91       {
 92         jia(m*tot+i,S,a[i],0);
 93         jia(S,m*tot+i,0,0);
 94       }
 95     for(int i=1;i<=m;i++)
 96       for(int j=1;j<=n;j++)
 97         {
 98           jia((i-1)*tot+1,m*tot+j,1,t[j][i]);
 99           jia(m*tot+j,(i-1)*tot+1,0,-t[j][i]);
100         }
101     for(;bfs();)
102       wang();
103     printf("%d",ans);
104     return 0;
105 }