bzoj3280

bzoj1221的变形版本，弄懂1221，这题还是很简单的

 const inf=;
 type node=record
        point,next,flow,cost:longint;
      end;

 var edge:array[..] of node;
     q:array[..] of longint;
     cur,pre,d,p,a:array[..] of longint;
     v:array[..] of boolean;
     i,j,x,y,e,n,m,k,t,te,len,sf:longint;

 function min(a,b:longint):longint;
   begin
     if a>b then exit(b) else exit(a);
   end;

 procedure add(x,y,f,c:longint);
   begin
     inc(len);
     edge[len].point:=y;
     edge[len].flow:=f;
     edge[len].cost:=c;
     edge[len].next:=p[x];
     p[x]:=len;
   end;

 function spfa:boolean;
   var f,r,i,x,y:longint;
   begin
     for i:= to t do
       d[i]:=inf;
     d[]:=;
     f:=;
     r:=;
     fillchar(v,sizeof(v),false);
     v[]:=true;
     while f<=r do
     begin
       x:=q[f];
       i:=p[x];
       v[x]:=false;
       while i<>- do
       begin
         y:=edge[i].point;
         if edge[i].flow> then
           if d[y]>d[x]+edge[i].cost then
           begin
             d[y]:=d[x]+edge[i].cost;
             pre[y]:=x;
             cur[y]:=i;
             if not v[y] then
             begin
               inc(r);
               q[r]:=y;
               v[y]:=true;
             end;
           end;
         i:=edge[i].next;
       end;
       inc(f);
     end;
     if d[t]=inf then exit(false) else exit(true);
   end;

 procedure mincost;
   var i,j,neck,f,c:longint;
   begin
     f:=;
     c:=;
     while spfa do
     begin
       neck:=inf;
       i:=t;
       while i<> do
       begin
         j:=cur[i];
         neck:=min(neck,edge[j].flow);
         i:=pre[i];
       end;
       i:=t;
       while i<> do
       begin
         j:=cur[i];
         dec(edge[j].flow,neck);
         inc(edge[j xor ].flow,neck);
         i:=pre[i];
       end;
       f:=f+neck;
       c:=c+d[t]*neck;
   //    writeln(f,' ',c);
     end;
     if f=sf then writeln(c) else writeln('impossible');
   end;

 begin
   readln(te);
   for e:= to te do
   begin
     len:=-;
     fillchar(p,sizeof(p),);
     readln(n,m,k);
     t:=n*+;
     sf:=;
     for i:= to n do
     begin
       read(a[i]);
       add(i,t,a[i],);
       add(t,i,,);
       sf:=sf+a[i];
     end;
     readln;
     for i:= to m do
     begin
       read(x,y);
       add(,,x,y);
       add(,,,-y);
     end;
     readln;
     for i:= to k do
     begin
       read(x,y);
       for j:= to n do
         if j++x<=n then
         begin
           add(j+n,j++x,inf,y);
           add(j++x,j+n,,-y);
         end
         else break;
     end;
     for i:= to n do
     begin
       add(,i+n,a[i],);
       add(i+n,,,);
       if i<>n then
       begin
         add(i,i+,inf,);  //今天没有用的人可以明天用
         add(i+,i,,);
         add(i+n,i+n+,inf,);
         add(i+n+,i+n,,);
       end;
     end;
     write('Case ',e,': ');
     mincost;
   end;
 end.