【dp】导弹拦截

拦截导弹

(missile.pas/c/cpp)

来源:NOIP1999(提高组) 第一题

【问题描述】

    某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度。某天,雷达捕捉到敌国的导弹来袭。由于该系统还在试用阶段,所以只有一套系统,因此有可能不能拦截所有的导弹。

   输入导弹依次飞来的高度(雷达给出的高度数据是不大于30000的正整数),计算这套系统最多能拦截多少导弹,如果要拦截所有导弹最少要配备多少套这种导弹拦截系统。

【输入文件】missile.in

  单独一行列出导弹依次飞来的高度。

【输出文件】missile.out

  两行,分别是最多能拦截的导弹数,要拦截所有导弹最少要配备的系统数

【输入样例】

389 207 155 300 299 170 158 65

【输出样例】

6

2

【问题分析】

有经验的选手不难看出这是一个求最长非升子序列问题,显然标准算法是动态规划。

以导弹依次飞来的顺序为阶段,设计状态opt[i]表示前i个导弹中拦截了导弹i可以拦截最多能拦截到的导弹的个数。

状态转移方程:

opt[i]=max(opt[j])+1  (h[i]>=h[j],0=<j<i)   {h[i]存,第i个导弹的高度}

最大的opt[i]就是最终的解。

这只解决了第一问,对于第二问最直观的方法就是求完一次opt[i]后把刚才要打的导弹去掉,在求一次opt[i]直到打完所有的导弹,但这样做就错了。

不难举出反例: 6 1 7 3 2       

错解: 6 3 2/1/7   正解:6 1/7 3 2

其实认真分析一下题就回发现:每一个导弹最终的结果都是要被打的,如果它后面有一个比它高的导弹,那打它的这个装置无论如何也不能打那个导弹了,经过这么一分析,这个问题便抽象成在已知序列里找最长上升序列的问题。

求最长上升序列和上面说的求最长非升序列是一样的,这里就不多说了。

复杂度:时间复杂度为O(N2),空间复杂度为O(N)。

 


 

 

题意:一种导弹拦截系统的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能高于前一发的高度。某天,雷达捕捉

到敌国的导弹来袭。由于该系统还在试用阶段,所以只有一套系统,因此有可能不能拦截所有的导弹。输入导弹依次飞来的高

度,计算这套系统最多能拦截多少导弹,如果要拦截所有导弹最少要配备多少套这种导弹拦截系统?

 

 

分析:

第一问很简单,就是求最长不上升子序列。对于第二问不能用贪心的方法来做,因为有反例:7 5 4 1 6 3 2.

我们把第二问的问题抽象出来,那就是:把一个数列划分成最少的最长不升子序列,这里我们要介绍一个很优美的定理。

 

Dilworth定理:对于一个偏序集,最少链划分等于最长反链长度。

Dilworth定理的对偶定理:对于一个偏序集,其最少反链划分数等于其最长链的长度。

 

也就是说把一个数列划分成最少的最长不升子序列的数目就等于这个数列的最长上升子序列的长度。

 


下面来说说这个定理是怎么来的:

 

偏序集的定义:偏序是在集合X上的二元关系≤(这只是个抽象符号,不是“小于或等于”,它满足自反性、反对称性和传递

性)。即,对于X中的任意元素a,b和c,有:

 

(1)自反性:a≤a;

(2)反对称性:如果a≤b且b≤a,则有a=b;

(3)传递性:如果a≤b且b≤c,则a≤c 。

 

带有偏序关系的集合称为偏序集。

 

 


令(X,≤)是一个偏序集,对于集合中的两个元素a、b,如果有a≤b或者b≤a,则称a和b是可比的,否则a和b不可比。

在这个例子(反链)中元素Ri<=Rj是指(i<=j) and (ai>=aj)

 

一个反链A是X的一个子集,它的任意两个元素都不能进行比较。

一个链C是X的一个子集,它的任意两个元素都可比。

 

【定理】

在X中,对于元素a,如果任意元素b,都有a≤b,则称a为极小元。

定理1:令(X,≤)是一个有限偏序集,并令r是其最大链的大小。则X可以被划分成r个但不能再少的反链。

 

其对偶定理称为Dilworth定理:

令(X,≤)是一个有限偏序集,并令m是反链的最大的大小。则X可以被划分成m个但不能再少的链。

虽然这两个定理内容相似,但第一个定理证明要简单一些。此处就只证明定理1。

 

证明:设p为最少反链个数

(1)先证明X不能划分成小于r个反链。由于r是最大链C的大小,C中任两个元素都可比,因此C中任两个元素都不能属于同一反

链。所以p>=r。

(2)设X1=X,A1是X1中的极小元的集合。从X1中删除A1得到X2。注意到对于X2中任意元素a2,必存在X1中的元素a1,使得

a1<=a2。令A2是X2中极小元的集合,从X2中删除A2得到X3……,最终会有一个Xk非空而Xk+1为空。于是A1,A2,…,Ak就是X的

反链的划分,同时存在链a1<=a2<=…<=ak,其中ai在Ai内。由于r是最长链大小,因此r>=k。由于X被划分成了k个反链,因此

r>=k>=p。

(3)因此r=p,定理1得证。

 

【解决】

要求最少的覆盖,按照Dilworth定理

最少链划分 = 最长反链长度

所以最少系统 = 最长导弹高度上升序列长度。

 


 

 

 1 #include<bits/stdc++.h>
 2 
 3 int arr[1000010];
 4 int dp[1000010];
 5 
 6 int main(){
 7     int n;
 8     while(~scanf("%d", &n)){
 9         for(int i = 0; i < n; i++) scanf("%d", &arr[i]);
10 
11         for(int i = 0; i < n; i++) dp[i] = 1;
12         for(int i = 1; i < n; i++){
13             for(int j = 0; j < i; j++){
14                 if(arr[i] <= arr[j] && dp[j]+1 > dp[i]) dp[i] = dp[j]+1;
15             }
16         }
17         int ans = 0;
18         for(int i = 0; i < n; i++) if(ans < dp[i]) ans = dp[i];
19         printf("%d\n", ans);
20 
21         for(int i = 0; i < n; i++) dp[i] = 1;
22         for(int i = 1; i < n; i++){
23             for(int j = 0; j < i; j++){
24                 if(arr[i] > arr[j] && dp[j]+1 > dp[i]) dp[i] = dp[j]+1;
25             }
26         }
27         ans = 0;
28         for(int i = 0; i < n; i++) if(ans < dp[i]) ans = dp[i];
29         printf("%d\n", ans);
30     }
31 
32     return 0;
33 }

 

转载于:https://www.cnblogs.com/miaowTracy/p/5927744.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/356703.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

安装tron_具有Tron效果的JavaFX 2 Form

安装tron这是一个具有TRON效果的简单JavaFX登录表单。 在此示例中&#xff0c;我使用CSS设置TextField和Button的样式。 这是CSS和Effect代码的片段&#xff1a; .text-field{-fx-background-color: transparent;-fx-border-color: #00CCFF;-fx-text-fill: white; }.password-f…

linux版本photoshop,在linux上使用ps(转载)

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼用Photoshop的&#xff0c;虽说GIMP也是不错&#xff0c;但是用起来总有点不顺的感觉&#xff0c;决定搞一个运行在Linux下的Photoshop&#xff0c;找了好多资料都没有虚拟CS的&#xff0c;不过倒是了解了不少虚拟Photoshop 这方面…

优美三角剖分_递归

问题 D: 优美三角剖分 时间限制: 1 Sec 内存限制: 64 MB提交: 13 解决: 4[提交][状态][讨论版]题目描述 小X同学为了搞好和小C同学的关系&#xff0c;特意寻找了一些优美的图像作为礼物。这是一些由无穷无尽三角形组成的极为优美的图形&#xff0c;小X同学很想实现这些极富美…

在Spring 4.2中更简单地处理异步事务绑定事件

介绍 如您可能已经知道的&#xff08;例如&#xff0c;从我以前的博客文章中 &#xff09;&#xff0c;不再需要创建一个单独的类&#xff0c;该类使用onApplicationEvent方法实现ApplicationListener以便能够对应用程序事件做出响应&#xff08;包括来自Spring Framework本身和…

linux chattr 无权限,从零开始学习Linux(二十八):文件权限之chattr权限

1、chattr命令命令格式&#xff1a; chattr [-] [选项] 文件名或者目录名&#xff1b;参数说明&#xff1a;&#xff1a;增加权限&#xff1b;-&#xff1a;删除权限&#xff1b; 等于某权限&#xff1b;选项说明&#xff1a;i&#xff1a;如果对文件设置i属性&#xff0c;则不…

初等数论及其应用——中国剩余定理

在线性代数中&#xff0c;我们用高斯消元解决多元的线性方程组&#xff0c;而在数论中&#xff0c;面对一元变量的线性模方程组&#xff0c;我们利用中国剩余定理去求解x。 转载于:https://www.cnblogs.com/rhythmic/p/5928483.html

linux下编译openjdk7,Ubuntu 14.04.3 64位环境下OpenJDK7编译

系统环境&#xff1a;Ubuntu 14.04.3 -amd64 in VMWare1.安装openjdk7依赖sudo apt-get build-dep openjdk-7sudo apt-get install openjdk-7-jdksudo apt-get install build-essential gawk m4 libasound2-dev libcups2-dev libxrender-dev xorg-dev xutils-dev x11proto-prin…

jQuery选择器的效率问题

jQuery提供了功能强大&#xff0c;并兼容多种css版本的选择器&#xff0c;不过发现很多同学在使用选择器时并未注重效率的问题。 a) 尽量使用Id选择器&#xff0c;jQuery选择器使用的API都是基于getElementById或getElementsByTagName&#xff0c;因此可以知道效率最高的是Id选…

HTTP请求和标头参数的CDI拦截器-简单示例

在使用Java EE开发和发展REST API的过程中&#xff0c;有些情况下您希望以更细粒度的方式&#xff08;视情况而定&#xff09;对传入的HTTP请求&#xff08;特别是标头参数&#xff09;进行“观察” Servlet过滤器或更具体的ContainerRequestFilters 。 我发现在某些情况下非常…

csgo服务器搭建 linux,请问搭建csgo服务器切换创意工坊地图报错怎么办

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼Unknown command "sv_maxcmdrate"Unknown command "cash_team_elimination_hostage_map"Unknown command "cash_team_win_by_time_running_out"Unknown command "sv_vote_creation_time"W…

1405 树的距离之和

1405 树的距离之和 基准时间限制&#xff1a;1 秒 空间限制&#xff1a;131072 KB给定一棵无根树&#xff0c;假设它有n个节点&#xff0c;节点编号从1到n, 求任意两点之间的距离&#xff08;最短路径&#xff09;之和。Input第一行包含一个正整数n (n < 100000)&#xff0c…

linux c 11 运行库,11.1.3 运行库与I/O

11.1.3 运行库与I/O在了解了glibc和MSVC的入口函数的基本思路之后&#xff0c;让我们来深入了解各个初始化部分的具体实现。但在具体了解初始化之前&#xff0c;我们要先了解一个重要的概念&#xff1a;I/O。IO(或I/O)的全称是Input/Output&#xff0c;即输入和输出。对于计算…

设置Elasticsearch N-Gram Word分析器的过程

说n-gram是一个大话题&#xff0c;这是轻描淡写的。 进行快速搜索&#xff0c;您会发现自己盯着语言学和语言模型&#xff0c;数据挖掘或特定蛋白质的分解对首次亮相文化的衰落所隐含的信息量。 好吧&#xff0c;我在跟那最后一个开玩笑。 但是&#xff0c;如果您是开发人员&a…

linux windows 丢失,Win10预览版9879硬盘丢失的Linux解决方案

IT之家讯 12月3日消息&#xff0c;最近IT之家论坛网友九仙仙总结了Win10预览版9879硬盘问题的解决方法&#xff0c;并在论坛中发布出来。经测试&#xff0c;此为快速有效的解决方法&#xff0c;故公之于众供朋友们参考。以下为作者原文。开头说明两点&#xff1a;1、这是我个人…

linux单网卡多拨Adsl,ROS单线多拨pppoe

#设置内网IP地址/ ip add add interfaceether1address192.168.0.254/255.255.255.0#设置共享上网/ ip firewall nat add actionmasquerade chainsrcnatsrc-address192.168.0.0/255.255.255.0#设置ADSL拔号/ interface pppoe-client add name"pppoe-out1" interfaceet…

Android开发——Android系统启动以及APK安装、启动过程

0. 前言 从Android手机打开开关&#xff0c;到我们可以使用其中的app时&#xff0c;这个启动过程到底是怎么样的&#xff1f;1. 系统上电当给Android系统上电&#xff0c;在电源接通的瞬间&#xff0c;CPU内的寄存器和各引脚均会被置为初始状态&#xff0c;CPU复位之后&#…

linux生成文件清单,Linux 获取文件名称生成列表 txt - create_filelist

Linux 获取文件名称生成列表 txt - create_filelist1. find/home/strong/MOTChallenge/MOT16/MOT16/train/MOT16-04/img1/ 文件夹下所有 *.jpg (000001.jpg - 001050.jpg) 的路径 图片名信息写入 txt 文件。1.1 llstrongforeverstrong:~/MOTChallenge/MOT16/MOT16/train/MOT16…

jvm线程分析命令_JVM:如何分析线程转储

jvm线程分析命令本文将教您如何分析JVM线程转储&#xff0c;并查明问题的根本原因。 以我的观点&#xff0c;线程转储分析是掌握Java EE生产支持的任何个人最重要的技能。 您可以从线程转储快照中获取的信息量通常远远超出您的想象。 我的目标是与您分享我在过去10年中积累的有…

[na]华为acl(traffic-filter)和dhcp管理

这个是财务网络的一个问题, 要求财务的某台机器能访问其他部门区的打印机. 其他部门是不能访问到财务网络的. 华为alc配置实例:-traffic-filter # 在VLAN100上配置基于ACL的报文过滤&#xff0c;允许源IP地址为192.168.0.2/32的报文通过&#xff0c;丢弃其他报文。 <HUAWEI&…

svn 1.6 linux 下载,LINUX下Subversion1.6.17 部署

1、服务器版本和安装软件服务器版本&#xff1a;Redhat Enterprise 5.5 32bit安装软件&#xff1a;apr-1.3.12.tar.gzapr-util-1.3.12.tar.gzhttpd-2.2.9.tar.gzsqlite-amalgamation-3.7.3.tarsubversion-1.6.17.tar.bz2subversion-deps-1.6.17.tar.bz22、安装检查操作系统http…