码迷,mamicode.com
首页 > 其他好文 > 详细

【BZOJ】【2245】【SDOI2011】工作安排

时间:2015-03-19 21:39:57      阅读:127      评论:0      收藏:0      [点我收藏+]

标签:

网络流/费用流


  裸题吧……直接建模就好了……所谓的“分段函数”就是吓唬你的,其实就是对于每个人分开建几条流量不同、费用不同的弧而已。

  对每种产品,连S->i ,(c[i],0);对每个工作人员 j ,连多条 j+n->T,流量分别为s[k]-s[k-1],对应的费用为w[j][k],至于那个矩阵……如果工作人员 j 可以做产品 i 就连边 i->j+n (INF,0);反正这题只有人有花费= =其他的都设置费用为0就行了啊……

  注意所有数据都在$10^5$以内,所以答案是会超过int的,需要使用long long

技术分享
  1 /**************************************************************
  2     Problem: 2245
  3     User: Tunix
  4     Language: C++
  5     Result: Accepted
  6     Time:1964 ms
  7     Memory:6216 kb
  8 ****************************************************************/
  9  
 10 //BZOJ 2245
 11 #include<vector>
 12 #include<cstdio>
 13 #include<cstring>
 14 #include<cstdlib>
 15 #include<iostream>
 16 #include<algorithm>
 17 #define rep(i,n) for(int i=0;i<n;++i)
 18 #define F(i,j,n) for(int i=j;i<=n;++i)
 19 #define D(i,j,n) for(int i=j;i>=n;--i)
 20 #define pb push_back
 21 using namespace std;
 22 inline int getint(){
 23     int v=0,sign=1; char ch=getchar();
 24     while(ch<0||ch>9){ if (ch==-) sign=-1; ch=getchar();}
 25     while(ch>=0&&ch<=9){ v=v*10+ch-0; ch=getchar();}
 26     return v*sign;
 27 }
 28 const int N=510,M=200000,INF=~0u>>2;
 29 typedef long long LL;
 30 /******************tamplate*********************/
 31 int n,m,c[N],a[251][251],s[6];
 32 LL ans;
 33 struct edge{int from,to,v,c;};
 34 struct Net{
 35     edge E[M];
 36     int head[N],next[M],cnt;
 37     void ins(int x,int y,int z,int c){
 38         E[++cnt]=(edge){x,y,z,c};
 39         next[cnt]=head[x]; head[x]=cnt;
 40     }
 41     void add(int x,int y,int z,int c){
 42         ins(x,y,z,c); ins(y,x,0,-c);
 43     }
 44     int from[N],Q[M],d[N],S,T,ed;
 45     bool inq[N],sign;
 46     bool spfa(){
 47         int l=0,r=-1;
 48         F(i,1,T) d[i]=INF;
 49         d[S]=0; Q[++r]=S; inq[S]=1;
 50         while(l<=r){
 51             int x=Q[l++];
 52             inq[x]=0;
 53             for(int i=head[x];i;i=next[i])
 54                 if(E[i].v>0 && d[x]+E[i].c<d[E[i].to]){
 55                     d[E[i].to]=d[x]+E[i].c;
 56                     from[E[i].to]=i;
 57                     if (!inq[E[i].to]){
 58                         Q[++r]=E[i].to;
 59                         inq[E[i].to]=1;
 60                     }
 61                 }
 62         }
 63         return d[T]!=INF;
 64     }
 65     void mcf(){
 66         int x=INF;
 67         for(int i=from[T];i;i=from[E[i].from])
 68             x=min(x,E[i].v);
 69         for(int i=from[T];i;i=from[E[i].from]){
 70             E[i].v-=x;
 71             E[i^1].v+=x;
 72         }
 73         ans+=(LL)x*d[T];
 74     }
 75     void init(){
 76         m=getint(); n=getint(); cnt=1;
 77         S=0; T=n+m+1;
 78         int x,y,z;
 79         F(i,1,n){
 80             x=getint();
 81             add(S,i,x,0);
 82         }
 83         F(j,1,m) F(i,1,n){
 84             x=getint();
 85             if(x) add(i,j+n,INF,0);
 86         }
 87         F(i,1,m){
 88             x=getint();
 89             if (x) F(j,1,x) s[j]=getint();
 90             s[0]=0; s[x+1]=INF;
 91             F(j,1,x+1){
 92                 y=getint();
 93                 add(i+n,T,s[j]-s[j-1],y);
 94             }
 95         }
 96         while(spfa()) mcf();
 97         printf("%lld\n",ans);
 98     }
 99 }G1;
100  
101 int main(){
102 #ifndef ONLINE_JUDGE
103     freopen("2245.in","r",stdin);
104     freopen("2245.out","w",stdout);
105 #endif
106     G1.init();
107     return 0;
108 }
View Code

 

【BZOJ】【2245】【SDOI2011】工作安排

标签:

原文地址:http://www.cnblogs.com/Tunix/p/4351792.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!