bzoj3280: 小R的煩惱（費用流）

題目傳送門
。

解法：
一開始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;
}