最大流-最大利益

来源:互联网 发布:淘宝商品怎么推广 编辑:程序博客网 时间:2024/04/27 20:41

。。起初以为基本的最大流最少也能过个4 ~5组数据。。可是不然。。。
就拿NOI2006-Day2的profit来说吧。。。
【问题描述】
新的技术正冲击着手机通讯市场,对于各大运营商来说,这既是机遇,更是挑战。THU集团旗下的CS&T通讯公司在新一代通讯技术血战的前夜,需要做太多的准备工作,仅就站址选择一项,就需要完成前期市场研究、站址勘测、最优化等项目。
在前期市场调查和站址勘测之后,公司得到了一共N个可以作为通讯信号中转站的地址,而由于这些地址的地理位置差异,在不同的地方建造通讯中转站需要投入的成本也是不一样的,所幸在前期调查之后这些都是已知数据:建立第i个通讯中转站需要的成本为Pi(1≤i≤N)。
另外公司调查得出了所有期望中的用户群,一共M个。关于第i个用户群的信息概括为Ai, Bi和Ci:这些用户会使用中转站Ai和中转站Bi进行通讯,公司可以获益Ci。(1≤i≤M, 1≤Ai, Bi≤N)
THU集团的CS&T公司可以有选择的建立一些中转站(投入成本),为一些用户提供服务并获得收益(获益之和)。那么如何选择最终建立的中转站才能让公司的净获利最大呢?(净获利 = 获益之和 - 投入成本之和)
【输入格式】
输入文件中第一行有两个正整数N和M 。
第二行中有N个整数描述每一个通讯中转站的建立成本,依次为P1, P2, …, PN 。
以下M行,第(i + 2)行的三个数Ai, Bi和Ci描述第i个用户群的信息。
所有变量的含义可以参见题目描述。
【输出格式】
你的程序只要向输出文件输出一个整数,表示公司可以得到的最大净获利。
【输入样例】
5 5
1 2 3 4 5
1 2 3
2 3 4
1 3 3
1 4 2
4 5 3
【输出样例】
4
【样例说明】
选择建立1、2、3号中转站,则需要投入成本6,获利为10,因此得到最大收益4。
【评分方法】
本题没有部分分,你的程序的输出只有和我们的答案完全一致才能获得满分,否则不得分。
【数据规模和约定】
80%的数据中:N≤200,M≤1 000。
100%的数据中:N≤5 000,M≤50 000,0≤Ci≤100,0≤Pi≤100。

最基本的最大流用邻接表只能过2组。。。不用邻接表。。1组
我起初做没用邻接表。。只过了一组

这是我宽搜加邻接表后的程序

program as;var   i,j,k,m,n,sum,min,big : longint;   q1          : boolean;   a,less      : array[0..2500,0..2500] of integer;   b,c,s  : array[0..2000] of integer;   bool          : array[0..2000] of boolean;procedure findway(h :integer);begin   min:=s[h];   repeat      dec(a[b[c[h]],b[h]],min);      inc(a[b[h],b[c[h]]],min);      h:=c[h];   until h=1;   sum:=sum+min;end; { findway }procedure dps;var   i,j,k,q,h,min :integer;begin   q:=1;   h:=2;   b[1]:=1;   c[1]:=1;   s[1]:=10000;   while q<>h do   begin      for j:=1 to less[b[q],0] do        begin         i:=less[b[q],j];  if (a[b[q],i]<>0) and (bool[i]) then     begin        b[h]:=i;        c[h]:=q;        bool[i]:=false;        if a[b[q],b[h]]<=s[q] then    s[h]:=a[b[q],b[h]]                  else                    s[h]:=s[q];        if i=n then    begin       findway(h);       exit;    end;        inc(h);     end;          end;      inc(q);   end;   q1:=false;end; { shensou }procedure duru;var   i,j,k,x,y,z,s : integer;   b:array[1..201] of integer;begin   readln(n,m);   for i:=1 to n do     read(b[i]);   q1:=true;   fillchar(a,sizeof(a),0);   for i:=1 to m do     begin       read(x,y,z);       a[i+1,x+m+1]:=100;       a[i+1,y+m+1]:=100;       inc(less[i+1,0]);       a[1,i+1]:=z;       big:=big+z;     end;    s:=0;   for i:=m+2 to m+n+1 do     begin      inc(s);      a[i,m+n+2]:=b[s];     end;   n:=n+m+2;   for i:=1 to n do     for j:=1 to n do       if (a[i,j]<>0)or(a[j,i]<>0) then         begin           inc(less[i,0]);           less[i,less[i,0]]:=j;         end;end; { duru }begin   assign(input,'profit.in');   reset(input);   assign(output,'profit.ans');   rewrite(output);   duru;   repeat      fillchar(bool,sizeof(bool),true);      fillchar(b,sizeof(b),0);      fillchar(s,sizeof(s),0);      fillchar(c,sizeof(c),0);      bool[1]:=false;      dps;   until q1=false;   writeln(big-sum);   close(input);   close(output);end.

过了8组,后两组超内存了。。。。 不过我已经满足了,币以往的2,3组来说我已经很满足了

0 0