Sicily 2500---看风景

  
        题目网址: http://soj.me/2500
         问题描述:
        NN ≤ 1000)个士兵站成一排,每个士兵都希望能看到路边美好的风景。一个士兵能够看到左边(右边)风景的条件是在他左边(右边)的士兵都比他矮。由于某些原因,一些士兵不能看到风景(“看到”指看到左边或者看到右边,当然两边都能看到也是允许的),输入N和N个士兵的身高,求最少出列多少个士兵,使得剩下的士兵都至少能够看到一边的风景。
       思路:
        转化为最长递增子序列问题,将士兵分成两半,左半部分递增,右半部分递减,然后枚举这个分割点,求出最长的序列。不在序列中的士兵就是要出列的。
        怎么求解最长递增子序列呢?
        设A[i]表示序列中的第i个数,F[i]表示从1到i这一段中以A[i]结尾的最长递增子序列的长度,初始时设F[i] = 0(i = 1, 2, ..., len(A))。则有动态规划方程:F[i] = max{1, F[j] + 1} (j = 1, 2, ..., i - 1, 且A[j] < A[i])。
        但是这样每次都要从1遍历到i-1,不理想,可以优化下。 用一个变量len记录到目前为止所找出来的最长递增序列的长度,另外准备一个数组b[],用这个数组表示长度为j的递增序列中最后一个元素的值。在这里长度为j的递增序列不止一个,我们所要保存是那个最小的。为什么呢?因为最后一个元素越小,那么这个递增序列在往后被延长的机会越大。初始化b[0] = -1;len = 0;遍历a[i],如果 a[i]( 1 <= i <= n)比长为len的序列的最大值大(len在变动,b数组是递增的),则把这个元素直接添加到b数组的后面;如果这个元素比b数组的第一个元素还要小则把这个元素赋给b数组的第一个值;否则进行二分查找,在b数组里面找到第一个刚好大于或等于a[i]的数,赋值为a[i],最后得出的len就是最长递增子序列的长度。
       代码:
       
#include <iostream>
#include <cstring>
using namespace std;

int n;
int a[1003];

int find(int *a,int len,int val)//修改后的二分查找,若返回值为x,则a[x]>=val
{
       int left=0,right=len,mid=(left+right)/2;
       while(left<=right)
       {
              if(val>a[mid]) left=mid+1;
              else if(val<a[mid]) right=mid-1;
              else return mid;
              mid=(left+right)/2;
       }
       return left;
}

int leftmaxlength(int pos)  //0到pos之间的最长子序列
{

	  int b[1003],c[1003];
	  memset(b,0,sizeof(b));
	  memset(c,0,sizeof(c));
	  int len,j,i;
      b[0]=1;//b[i]表示以a[i]结尾的最长递增子序列的长度
      c[0]=-1;//c[i]表示长度为i的最长递增子序列的最小末尾元素
      c[1]=a[0];
      len=1;//此时只有c[1]求出来,最长递增子序列的长度为1.
      for(i=1;i<=pos;i++)
      {
             j=find(c,len,a[i]);//在c(递增)中找比a[i]大的最小下标
             c[j]=a[i];
             if(j>len)//要更新len,另外补充一点:由二分查找可知j只可能比len大1
                    len=j;//更新len
      }
	  return len;
}

int rightmaxlength(int pos)//pos+1  到 n-1之间的最长子序列
{
	  if(pos+1>=n) return 0;
	  int b[1003],c[1003],j;
	  memset(b,0,sizeof(b));
	  memset(c,0,sizeof(c));
	  int len,i;
      b[0]=1;//b[i]表示以a[i]结尾的最长递增子序列的长度
      c[0]=-1;//c[i]表示长度为i的最长递增子序列的最小末尾元素
      c[1]=a[n-1];
      len=1;//此时只有c[1]求出来,最长递增子序列的长度为1.
      for(i=n-2;i>pos;i--)
      {
             j=find(c,len,a[i]);//在c(递增)中找比a[i]大的最小下标
             c[j]=a[i];
             if(j>len)//要更新len,另外补充一点:由二分查找可知j只可能比len大1
                    len=j;//更新len
      }
	  return len;
}

int main()
{
       
		   cin>>n;

		   for(int i=0;i<n;i++)
                     cin>>a[i];
		   if(n<=2) {cout<<0<<endl;return 0;}
		   int maxm=-1;
		   for(int pos=0;pos<n;pos++)
		   {
			   int temp=leftmaxlength(pos)+rightmaxlength(pos);
			   if(temp>maxm) maxm=temp;
		   }
		   cout<<(n-maxm)<<endl;
              


       return 0;
}

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值