当前位置:网站首页>bzoj3280: 小R的烦恼(费用流)

bzoj3280: 小R的烦恼(费用流)

2020-12-07 15:21:27 osc_s6zj0021

题目传送门

解法:
一开始yy了个建图。
完全没毛病。
结果不知道怎么判impossible了。。
很尴尬。



每天分成两个点。
一个点x表示今天收到多少个活人。。
一个点y表示多少个濒死的人。

今天雇的人一定会用。所以变成死人。大学连y,容量无限费用q[i]
每天最多a[i]个死人。st连y容量a[i]费用0
每天最多a[i]个活人。x连ed容量a[i]费用0
死人救活了。懒得写怎么连了。
今天的活人可以留到明天。x连x+1容量无限费用0



最小费用最大流

代码实现:

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<queue>
using namespace std;
struct node {
  
  int x,y,c,d,next,other;}a[1100000];int len,last[110000];
void ins(int x,int y,int c,int d) {
    int k1,k2;
    len++;k1=len;a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
    a[len].next=last[x];last[x]=len;
    len++;k2=len;a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
    a[len].next=last[y];last[y]=len;
    a[k1].other=k2;a[k2].other=k1;
}
int list[11000],head,tail,st,ed,d[11000],fbian[11000];bool v[11000];
bool spfa() {
    head=1;tail=2;list[1]=st;
    for(int i=1;i<=ed;i++)d[i]=999999999;d[st]=0;
    memset(v,false,sizeof(v));v[st]=true;memset(fbian,0,sizeof(fbian));
    while(head!=tail) {
        int x=list[head];
        for(int k=last[x];k;k=a[k].next) {
            int y=a[k].y;
            if(d[y]>d[x]+a[k].d&&a[k].c>0) {
                d[y]=d[x]+a[k].d;fbian[y]=k;
                if(v[y]==false) {
                    v[y]=true;list[tail++]=y;
                    if(tail==ed+1)tail=1;
                }
            }
        }v[x]=false;head++;if(head==ed+1)head=1;
    }return d[ed]==999999999?false:true;
}
int F;
int wh() {
    int x=ed,flow=999999999;
    while(x!=st) {
  
  int k=fbian[x];flow=min(flow,a[k].c);x=a[k].x;}
    x=ed;F+=flow;
    while(x!=st) {
  
  int k=fbian[x];a[k].c-=flow;a[a[k].other].c+=flow;x=a[k].x;}
    return flow*d[ed];
}
const int inf=999999999;
int main() {
    int T;scanf("%d",&T);
    for(int tt=1;tt<=T;tt++) {
        int n,m,k;scanf("%d%d%d",&n,&m,&k);
        int sum=0;st=n*2+m+1;ed=st+1;
        len=0;memset(last,0,sizeof(last));
        for(int i=1;i<=n;i++) {
            int x;scanf("%d",&x);sum+=x;
            ins(st,i+m,x,0);ins(i+n+m,ed,x,0);
            if(i!=n)ins(i+n+m,i+n+m+1,inf,0);
        }
        for(int i=1;i<=m;i++) {
            int x,y;scanf("%d%d",&x,&y);
            ins(st,i,x,0);
            for(int j=1;j<=n;j++) ins(i,j+n+m,inf,y);
        }
        while(k--) {
            int x,y;scanf("%d%d",&x,&y);
            for(int i=1;i<=n-x-1;i++) ins(i+m,i+m+n+x+1,inf,y);
        }F=0;
        int ans=0;
        while(spfa()==true) ans+=wh();
        printf("Case %d: ",tt);
        if(F<sum)printf("impossible\n");
        else printf("%d\n",ans);
    }
    return 0;
}

版权声明
本文为[osc_s6zj0021]所创,转载请带上原文链接,感谢
https://my.oschina.net/u/4356138/blog/4778677