【动态规划：防卫导弹】：离散化+DP；15ms;

防卫导弹
时限：1000ms 内存限制：10000K  总时限：3000ms
描述：

一种新型的防卫导弹可截击多个攻击导弹。它可以向前飞行，也可以用很快的速度向下飞行，可以毫无损伤地截击进攻导弹，但不可以向后或向上飞行。但有一个缺点，尽管它发射时可以达到任意高度，但它只能截击比它上次截击导弹时所处高度低或者高度相同的导弹。现对这种新型防卫导弹进行测试，在每一次测试中，发射一系列的测试导弹（这些导弹发射的间隔时间固定，飞行速度相同），该防卫导弹所能获得的信息包括各进攻导弹的高度，以及它们发射次序。现要求编一程序，求在每次测试中，该防卫导弹最多能截击的进攻导弹数量，一个导弹能被截击应满足下列两个条件之一：
a)它是该次测试中第一个被防卫导弹截击的导弹；
b)它是在上一次被截击导弹的发射后发射，且高度不大于上一次被截击导弹的高度的导弹。

输入：

多个测例。
每个测例第一行是一个整数n（n不超过100），第二行n个整数表示导弹的高度（数字的顺序即发射的顺序）。
n＝0表示输入结束。

输出：

每个测例在单独的一行内输出截击导弹的最大数目。

输入样例：

5
5 6 100 6 61

输出样例：

2

#include <iostream> using namespace std; int *pos; int *col; int **mar; int c_p; int n; int compar (const void *a ,const void *b) { int *aa=(int *)a; int *bb=(int *)b; if(*aa>*bb) return 1; if(*aa==*bb) return 0; if(*aa<*bb) return -1; return 0; } int maxa(int a,int b) { return a>b? a:b; } int dp() { int max; int **get=new int*[n]; for(int f=0;f<n;++f) { get[f]=new int[n]; } get[0][n-1]=0; for(int i=0;i<c_p;++i) { for(int j=n;j>=0;--j) { int temp1=0,temp2=0; if(i-1>=0) { temp1=get[i-1][j]; } if(j+1<n) { temp2=get[i][j+1]; } get[i][j]=mar[i][j]+maxa(temp1,temp2); } } max=get[0][0]; for(int h=1;h<c_p;++h) { if(get[h][0]>max) { max=get[h][0]; } } return max; } int main() { cin>>n; while(n!=0) { col=new int[n]; pos=new int[n]; mar=new int*[n]; c_p=0; for(int i=0;i<n;++i) { cin>>col[i]; bool flag=true; for(int j=0;j<c_p;++j) { if(pos[j]==col[i]) { flag=false; break; } } if(flag==true) { pos[c_p++]=col[i]; } mar[i]=new int[n]; } qsort(pos,c_p,sizeof(int),compar); for(int k=0;k<n;++k) { for(int l=0;l<c_p;++l) { mar[l][k]=0; if(pos[l]==col[k]) { mar[l][k]=1; } } } cout<<dp()<<endl; cin>>n; } return 0; }

 

解题思路:  首先将问题抽象成模型,在这里就是一个矩阵表示的炮弹分布图.  有的炮弹一样高,有的炮弹不一样高.

 

所以,矩阵的高应该是不同的高度的数量,矩阵的长应该是炮弹的数量n.

 

这里用到了离散化数据的技巧,不能理所当然的将10000米的炮弹设置为矩阵10000高,这样太费空间了.

 

所以,必须将不同的高度作为一个点处理,将不同高度的点压缩到一个很小的矩形上,只描述它们之间的次序,而不描述它们的具体高度.

 

例如下图分布: 有5个炮弹, 高度是100，9999，20，100，700，离散化后的矩阵应该是形如这样的：

 

       9999

                                700

100                  100

                 20

 

所以矩阵只需要高4,长5,就足以表示炮弹之间的位置关系了.

 

我采用O(n2)的时间复杂度进行离散.   方法如下:

 

利用col[n]存储输入的每个炮弹的高度, 利用pos[n]（如果高度有重复，则用不了n个）存储炮弹出现过的不重复高度,

例如: 5  100 9999 20 100 700， 那么输入数据结束时,col[5]={100,9999,20,100,700} ,pos[4]={100,9999,20,700}

，这个记录过程是在输入数据过程中完成的. 无非一个col记录所有炮弹的高度，一个pos记录有哪几种高度。

 

现在,需要生成上面炮弹分布图M了,依据的就是col,pos数组.,  首先把pos按高度自小到大排序 ， 从col[0]开始到col[4]，对每个col，去pos里找它所在的位置，例如先找100所在高度,发现它是pos里第1矮的，那么就标记M[1][0]=1,表示这里有炮弹。 然后9999，发现它是pos里第3矮的，那么标记M[3][1]=1.  对于20则是 M[0][2]=1, 对于100，M[1][3]=1,对于700，M[2][4]=1

 

现在炮弹分布图生成了，除了标记1的，其他都标记0 ，现在我们就有了一个非常直观的0，1矩阵. 1表示有炮弹,0表示没炮弹.

现在用动态规划O(n2)就可以解决问题了，sum[i][j]表示从(i,j)打出的炮弹能打到的最多数量,最终要求的是max{sum[0 to 3][0]}，就是我们想要的结果.

 

由于出发点(i,j)不同,那么它能打到的炮弹的路线也不同，所以要求max{}。另外，DP转移式为sum[i][j]=max{sum[i-1][j],sum[i][j+1]+M[i][j]}, M[i][j] 就是那个炮弹分布图, 表示当前格子有没有炮弹,有炮弹当然也要算上1个。