拦截导弹(文件名:missile.c/cpp/pas)[O(n*log(n)),二分优化]

   拦截导弹(文件名:missile.c/cpp/pas)
 
某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷:
虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度。某天,雷
达捕捉到敌国的导弹来袭。由于该系统还在试用阶段,所以只有一套系统,因此有可能不能拦截所有
的导弹。
  
输入数据:
第一行为一个整数 N,表示飞来的导弹个数,N<=100000
第二行为 N 个整数,依次表示导弹飞来的高度,高度数据为不大于 30000 的正整数。
 
输出数据:
第一行,输出计算这套系统最多能拦截多少导弹
第二行,输出要拦截所有导弹最少要配备多少套这种导弹拦截系统。  
 
样例
输入文件:missile.in
8
389 207 155 300 299 170 158 65  
 
输出文件:missile.out
6  
2

===================

1问:

开一个g[i]表记录可拦截i个导弹的最小高度...

可证明g[i]表中的数是单调递增的

  因为:拦截1个导弹的最小高度永远<=拦截2个导弹的最小高度...

二分所能打的个数..

--------------------------------

2问:最长上升序列...

证明:利用第一问分析模拟打导弹的情况...

就可以得出结论..

===========================

var
  n:longint;
  high,g:array[1..100000]of longint;
//  g_len:longint;
  ans:longint;
procedure init;
begin
  assign(input,'missile.in');
  assign(output,'missile.out');
  reset(input); rewrite(output);
end;

procedure terminate;
begin
  close(input); close(output);
  halt;
end;

procedure divide_two1(h:longint);
var
  l,r,m:longint;
begin
  l:=1; r:=n;
  ans:=0;
  while l<=r do
    begin
      m:=(l+r) shr 1;
      if h<=g[m] then begin ans:=m; l:=m+1 end
                 else r:=m-1;
    end;
end;

procedure divide_two2(h:longint);
var
  l,r,m:longint;
begin
  ans:=0;
  l:=1; r:=n;
  while l<=r do
    begin
      m:=(l+r) shr 1;
      if h>=g[m] then begin ans:=m; l:=m+1 end
                 else r:=m-1;
    end;
end;

procedure main;
var
  i:longint;
  ans1,ans2:longint;
begin
  readln(n);
  for i:=1 to n do
    begin
      read(high[i]);
      //f[i]:=1;
    end;
    
  fillchar(g,sizeof(g),0);
  g[1]:=high[1];        //表示拦截此次数为i的最大高度...
  //g_h[1]:=high[1];
  //g_len:=1;
  ans1:=0;
  
  for i:=2 to n do
    begin
      divide_two1(high[i]);
      if high[i]>g[ans+1] then g[ans+1]:=high[i];
      if ans1<ans+1 then ans1:=ans+1;
    end;

  for i:=2 to n do
    g[i]:=maxlongint;
  ans2:=0;
  g[1]:=high[1];
  for i:=2 to n do
    begin
      divide_two2(high[i]);
      if high[i]<g[ans+1] then g[ans+1]:=high[i];
      if ans2<ans+1 then ans2:=ans+1;
    end;
  writeln(ans1);
  writeln(ans2);
end;

begin
  init;
  main;
  terminate;
end.


 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值