最大递增子序列——POJ1836

     这道题目最简单的算法是枚举所有i(1<=i<=n),然后求h[1]到h[i]的最长递增子序列和h[i+1]到h[n]的最长递减子序列。算法的时间复杂度是O(n^2log(n))。

      还有一种算法，就是构造一个基准序列，基准序列是h数组中的所有元素的按递增顺序排序（去除重复元素）+h数组中的所有元素的按递减顺序排序（去除重复元素）所形成的序列。之后,我们只要求h数组和基准序列的最大公共子序列就行了。算法的时间复杂度是O(n^2)。

      还有一种时间复杂度为O(nlog(n))的算法。我们只需要构造出数组f1,f2，其中f1[i](0<=i<=n)保存h[1]到h[i]这个子序列的最长递增子序列长度（定义h[1]到h[0]的最大递增子序列长度为0）；f1[i](1<=i<=n+1)保存h[n]到h[i]这个子序列的最长递增子序列长度（定义h[n]到h[n+1]的最大递增子序列长度为0）。于是，最后只需要求f1[i]+f2[i+1](0<=i<=n)的最大值即可。

#include <iostream> using namespace std; const int MAXN=1010; double height[MAXN],mheight[MAXN]; int f1[MAXN],f2[MAXN]; int n; int bSearch(double v) { int lo=0,hi=n,mid,res=0; while(lo<=hi) { mid=(lo+hi)/2; if(mheight[mid]>=v)hi=mid-1; else { res=mid; lo=mid+1; } } return res; } int dp() { int i; /* 顺序最大递增序列 f1[i]保存从heght[1]到height[i]这个子序列中的最大递增子序列 */ for(i=0;i<=n;i++)mheight[i]=10; mheight[0]=0; for(i=1;i<=n;i++) { mheight[f1[i]=bSearch(height[i])+1]=height[i]; if(f1[i]<f1[i-1])f1[i]=f1[i-1]; } /* 逆序最大递增序列 f1[i]保存从heght[n]到height[i]这个子序列中的最大递增子序列 */ for(i=0;i<=n;i++)mheight[i]=10; mheight[0]=0; for(i=n;i>=1;i--) { mheight[f2[i]=bSearch(height[i])+1]=height[i]; if(f2[i]<f2[i+1])f2[i]=f2[i+1]; } int maxn=0; for(i=0;i<=n;i++) if(maxn<f1[i]+f2[i+1])maxn=f1[i]+f2[i+1]; return n-maxn; } int main() { scanf("%d",&n); for(int i=1;i<=n;i++)scanf("%lf",&height[i]); printf("%d/n",dp()); return 0; }