题目网址: http://soj.me/2500
问题描述:
N( N ≤ 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;
}