二分匹配的灵活运用

3041还是比较好想的,考虑到横排/竖排射一枪就能搞定这一行/一列的所有点,

我们以行数为点集x,列数为点集y,在目标点(xi,yi)之间连一条边

这样最小射击次数=最小点覆盖(边两端点至少有一个点在覆盖集中)=最大匹配

poj2226是它的加强版

这里一块木板不能覆盖非泥泞点,也就是说一块木板不一定能搞定那一行所有的点

于是我们考虑到将连续的泥泞点标号,表示这些泥泞点是由哪个木板盖住(横竖分别标号)

这样,我们以横排木板为点集x,竖排木板为点集y,对于每个点,在盖住它的横竖排木板之间连边

(显然,每个点至多被一个横木板,一个竖木板覆盖)

于是根据poj3041,做匈牙利即可

 var a:array[..,..] of char;
    h1,h2:array[..,..] of longint;
    map:array[..,..] of boolean;
    cx,cy:array[..] of longint;
    v:array[..] of boolean;
    s1,s2,ans,i,j,n,m:longint; function find(x:longint):longint;   //增广路
  var i:longint;
  begin
    for i:= to s2 do
      if map[x,i] and not v[i] then
      begin
        v[i]:=true;
        if (cy[i]=-) or (find(cy[i])>) then
        begin
          cx[x]:=i;
          cy[i]:=x;
          exit();
        end;
      end;
    exit();
  end; begin
  readln(n,m);
  for i:= to n do
  begin
    for j:= to m do
    begin
      read(a[i,j]);
    end;
    readln;
  end;
  for i:= to n do         //横排标号
  begin
    j:=;
    while j<=m do
    begin
      if a[i,j]='*' then
      begin
        inc(s1);
        h1[i,j]:=s1;
        inc(j);
        while a[i,j]='*' do     
        begin
          h1[i,j]:=s1;
          inc(j);
        end;
      end
      else inc(j);
    end;
  end;
  for i:= to m do      //竖排标号
  begin
    j:=;
    while j<=n do
    begin
      if a[j,i]='*' then
      begin
        inc(s2);
        h2[j,i]:=s2;
        inc(j);
        while a[j,i]='*' do
        begin
          h2[j,i]:=s2;
          inc(j);
        end;
      end
      else inc(j);
    end;
  end;
  for i:= to n do
    for j:= to m do
      map[h1[i,j],h2[i,j]]:=true;    //连边
  fillchar(cx,sizeof(cx),);
  fillchar(cy,sizeof(cy),);
  for i:= to s1 do     //匈牙利
    if cx[i]=- then
    begin
      fillchar(v,sizeof(v),false);
      ans:=ans+find(i);
    end;
  writeln(ans);
end.

一定要注意问题建模的灵活运用;

poj3041,poj2226的更多相关文章

  1. POJ-1325 Machine Schedule,和3041有着异曲同工之妙,好题!

    Machine Schedule Time Limit: 1000MS   Memory Limit: 10000K       Description As we all know, machine ...

  2. ZOJ 1654 Place the Robots(最大匹配)

    Robert is a famous engineer. One day he was given a task by his boss. The background of the task was ...

  3. Pyhton开源框架(加强版)

    info:Djangourl:https://www.oschina.net/p/djangodetail: Django 是 Python 编程语言驱动的一个开源模型-视图-控制器(MVC)风格的 ...

  4. MPlayer

    名称   mplayer − 电影播放器 mencoder − 电影编解码器 概要   mplayer [选项] [文件|URL|播放列表|−] mplayer [选项] 文件1 [指定选项] [文件 ...

  5. python 爬取腾讯微博并生成词云

    本文以延参法师的腾讯微博为例进行爬取并分析 ,话不多说 直接附上源代码.其中有比较详细的注释. 需要用到的包有 BeautifulSoup WordCloud jieba # coding:utf-8 ...

  6. 面经 cisco

    1. 优先级反转问题及解决方法 (1)什么是优先级反转 简单从字面上来说,就是低优先级的任务先于高优先级的任务执行了,优先级搞反了.那在什么情况下会生这种情况呢? 假设三个任务准备执行,A,B,C,优 ...

  7. linux驱动(续)

    网络通信 --> IO多路复用之select.poll.epoll详解 IO多路复用之select.poll.epoll详解      目前支持I/O多路复用的系统调用有 select,psel ...

  8. HttpServletRequest对象(一)

    javaweb学习总结(十)——HttpServletRequest对象(一) 一.HttpServletRequest介绍 HttpServletRequest对象代表客户端的请求,当客户端通过HT ...

  9. 【poj2226】 Muddy Fields

    http://poj.org/problem?id=2226 (题目链接) 题意 给出一个只包含‘.’和‘*’的矩阵,用任意长度的宽为1的木板覆盖所有的‘*’而不覆盖‘.’,木板必须跟矩形的长或宽平行 ...

随机推荐

  1. aix 安装redis

    下载最新rpm安装包 http://www.perzl.org/aix/index.php?n=Main.Redis # uname -aAIX rhjf 1 6 00C5CC964C00# pwd/ ...

  2. SPI协议及其工作原理详解

    一.概述. SPI, Serial Perripheral Interface, 串行外围设备接口, 是 Motorola 公司推出的一种同步串行接口技术. SPI 总线在物理上是通过接在外围设备微控 ...

  3. QT定制有标题的扁平化下拉框控件

    关键字:QT,QComboBox,QLineEdit,QListView,QPushButton,QMenu,QWidgetAction,setStyleSheet OS:Windows 7 方法一: ...

  4. PL/SQL中如何执行DDL、SCL?

    PL/SQL程序中不能直接执行DDL语句.为什么? 假设我们在pl/sql程序中有这样的一条DDL语句—— drop table emp:在第一次解析pl/sql中的“drop table emp;” ...

  5. 【EF】疑难杂症

    用户和购物车数据  主从表 添加 ADO.NET实体数据类型 [Test] public void EntiyConnect() { var context = new projectDatabase ...

  6. Spring execution 表达式

    execution(modifiers-pattern? ret-type-pattern declaring-type-pattern? name-pattern(param-pattern) th ...

  7. 【BZOJ [1878】[SDOI2009]HH的项链

    Description HH有一串由各种漂亮的贝壳组成的项链.HH相信不同的贝壳会带来好运,所以每次散步 完后,他都会随意取出一段贝壳,思考它们所表达的含义.HH不断地收集新的贝壳,因此, 他的项链变 ...

  8. 第一个js库文件

    <!DOCTYPE html> <html xmlns=;         ;                     }                 };     })(); ...

  9. vbe6ext.olb不能被加载 宏内存溢出

    今天想玩一下PowerPoint的宏,却发现玩不起来!!! 另外,每次打开ppt时都会提示vbe6ext.olb不能加载. 网上说重新下载个vbe6ext.olb然后复制到相应的路径.我也试着下载,然 ...

  10. oracle 导出数据和导入数据

    导出数据 exp zl_gj/zlkj@gqxt  grants=y tables=(zl_gj.ckgj,zl_gj.gjlx,zl_gj.rkgj) file=c:\gj.dmp log=c:\g ...