zoukankan      html  css  js  c++  java
  • 【BZOJ】【3280】小R的烦恼

    网络流/费用流


      和软件开发那题基本相同,只是多加了一个“雇佣研究生”的限制:不同价格的研究生有不同的数量……

      那么只需加一个附加源点,对每一种研究生连边 S->ss 容量为l[i],费用为p[i];再改由附加源点向每天的流入点(i+n)连边即可。

      1 /**************************************************************
      2     Problem: 3280
      3     User: Tunix
      4     Language: C++
      5     Result: Accepted
      6     Time:276 ms
      7     Memory:5972 kb
      8 ****************************************************************/
      9  
     10 //BZOJ 3280
     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,k,ans,tot,flow;
     32 struct edge{int from,to,v,c;};
     33 struct Net{
     34     edge E[M];
     35     int head[N],next[M],cnt;
     36     void ins(int x,int y,int z,int c){
     37         E[++cnt]=(edge){x,y,z,c};
     38         next[cnt]=head[x]; head[x]=cnt;
     39     }
     40     void add(int x,int y,int z,int c){
     41         ins(x,y,z,c); ins(y,x,0,-c);
     42     }
     43     int from[N],Q[M],d[N],S,T,ss;
     44     bool inq[N];
     45     bool spfa(){
     46         int l=0,r=-1;
     47         F(i,1,T) d[i]=INF;
     48         d[S]=0; Q[++r]=S; inq[S]=1;
     49         while(l<=r){
     50             int x=Q[l++];
     51             inq[x]=0;
     52             for(int i=head[x];i;i=next[i])
     53                 if(E[i].v>0 && d[x]+E[i].c<d[E[i].to]){
     54                     d[E[i].to]=d[x]+E[i].c;
     55                     from[E[i].to]=i;
     56                     if (!inq[E[i].to]){
     57                         Q[++r]=E[i].to;
     58                         inq[E[i].to]=1;
     59                     }
     60                 }
     61         }
     62         return d[T]!=INF;
     63     }
     64     void mcf(){
     65         int x=INF;
     66         for(int i=from[T];i;i=from[E[i].from])
     67             x=min(x,E[i].v);
     68         for(int i=from[T];i;i=from[E[i].from]){
     69             E[i].v-=x;
     70             E[i^1].v+=x;
     71         }
     72         flow+=x;
     73         ans+=x*d[T];
     74     }
     75     void init(){
     76         n=getint(); m=getint(); k=getint();
     77         flow=tot=ans=0; cnt=1; memset(head,0,sizeof head);
     78         S=0; ss=n+n+1; T=n+n+2;
     79         int x,y,z;
     80         F(i,1,n){
     81             x=getint(); tot+=x;
     82             add(S,i,x,0);
     83             if (i<n) add(i,i+1,INF,0);
     84             add(ss,i+n,INF,0);
     85             add(i+n,T,x,0);
     86         }
     87         F(i,1,m){
     88             x=getint(); y=getint();
     89             add(S,ss,x,y);
     90         }
     91         F(i,1,k){
     92             x=getint(); y=getint();
     93             F(j,1,n) if(j+x+1<=n) add(j,j+n+x+1,INF,y);
     94             else break;
     95         }
     96         while(spfa()) mcf();
     97         if (flow==tot) printf("%d
    ",ans);
     98         else puts("impossible");
     99     }
    100 }G1;
    101  
    102 int main(){
    103 #ifndef ONLINE_JUDGE
    104     freopen("3280.in","r",stdin);
    105     freopen("3280.out","w",stdout);
    106 #endif
    107     int T=getint();
    108     F(i,1,T){
    109         printf("Case %d: ",i);
    110         G1.init();
    111     }
    112     return 0;
    113 }
    View Code

    3280: 小R的烦恼

    Time Limit: 10 Sec  Memory Limit: 512 MB
    Submit: 210  Solved: 100
    [Submit][Status][Discuss]

    Description

    小R最近遇上了大麻烦,他的程序设计挂科了。于是他只好找程设老师求情。善良的程设老师答应不挂他,但是要求小R帮助他一起解决一个难题。
    问题是这样的,程设 老师最近要进行一项邪恶的实验来证明P=NP,这个实验一共持续n天,第i天需要a[i]个研究生来给他搬砖。研究生毕竟也是人,所以雇佣研究生是需要钱 的,机智的程设老师已经联系好了m所大学,第j所大学共有l[j]个研究生,同时雇佣这所大学的一个研究生需要p[j]元钱。
    本来程设老师满心欢 喜的以为,这样捡最便宜的max{a[i]}个研究生雇来,就可以完成实验;结果没想到,由于他要求硕士生们每天工作25个小时不许吃饭睡觉上厕所喝水说 话咳嗽打喷嚏呼吸空气,因此一天下来给他搬砖的所有研究生都会进入濒死状态。濒死状态的研究生,毫无疑问,就不能再进行工作了。但是机智的老师早早联系好 了k家医院,第i家医院医治一个濒死的研究生需要d[i]天,并且需要q[i]元钱。

    现在,程设老师想要知道,最少花多少钱,能够在这n天中满足每天的需要呢?若无法满足,则请输出”impossible”。注意,由于程设老师良心大大的坏,所以他是可以不把濒死的研究生送去医院的!

    Input

           本题包含多组数据;第一行是一个数T(T<=11),表示数据组数,以下T组数据。
    对于每一组数据,第一行三个数,n,m,k;
    以下一行n个数,表示a[1]…a[n]
    接着一行2m个数,表示l[1],p[1]…l[n],p[n]
    接着一行2k个数,表示d[1],q[1]…d[n],q[n]

    Output

           对于每组数据以样例的格式输出一行,两个数分别表示第几组数据和最少钱数。

    Sample Input

    2
    3 2 1
    10 20 30
    40 90 15 100
    1 5
    3 2 1
    10 20 30
    40 90 15 100
    2 5

    Sample Output

    Case 1: 4650
    Case 2: impossible

    HINT

    样例解释:买下90块钱的那40个研究生,另外再买10个100块钱的。这样,第一天用完的10个人全部送到医院,那么他们在第三天可以继续使用;同时,第二天和第三天都用新的研究生来弥补,这样一共需要花费40*90 + 10*100 + 5*10 = 4650元。

    数据规模:

    对于30%的数据中的每组数据,

    满足n<=5,m,k<=2,其余数均小于等于100或者

    n<=10,m,k<=10,其余数均小于等于20.

    对于100%的数据

    n,m,k<=50,其余数均小于等于100.

    Source

    [Submit][Status][Discuss]
  • 相关阅读:
    《算法导论》读书笔记
    【原创】POI操作Excel导入导出工具类ExcelUtil
    10-JMM
    09-字节码执行引擎
    08-类加载机制
    07-前端编译与优化(待补充)
    06-字节码指令
    05-类文件结构
    04-垃圾回收(2)
    03-垃圾回收(1)
  • 原文地址:https://www.cnblogs.com/Tunix/p/4352205.html
Copyright © 2011-2022 走看看