最大流-最大利益

。。起初以为基本的最大流最少也能过个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
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值