NYOJ--第214题单调递增数列（二）

最新推荐文章于 2021-03-27 14:28:21 发布

iteye_1261

最新推荐文章于 2021-03-27 14:28:21 发布

阅读量149

点赞数

这个题如果单纯的按照以前的方法来解决的话肯定会超时，所以需要我们对代码进行优化。优化用到了二分法，把以前那种（n^2）变成了（nlogn）其实我也是第一次接触二分法，之前对二分法并不是特别熟悉，所以就去百度了一下基本的概念。我觉得二分法可以单独整理一下。呵呵，废话不多说了，下面是代码和证明（证明是看的别人的，我也不是很明白，不过再思考一下应该可以明白）

原题地址：点击打开链接。

代码如下：#include<stdio.h> #include<string.h> #include<iostream> using namespace std; int a[100002],d[100002]; int erfen(int num,int len)//二分函数 { int left,right,mid; left=1; right=len; mid=(left+right)/2; while(left<=right) { if(d[mid]<num) left=mid+1; else if(d[mid]>num) right=mid-1; else return mid; mid=(left+right)/2; }return left; } int main() { int n,len,i,j; while(scanf("%d",&n)!=EOF) { for(i=0;i<n;i++) scanf("%d",&a[i]); memset(d,0,sizeof(0)); len=1; d[1]=a[0]; a[0]=-1000000; for(i=1;i<n;i++) { j=erfen(a[i],len); d[j]=a[i]; if(j>len) len=j; }printf("%d\n",len); }return 0; }

现在来证明这个算法为什么是正确的。要使算法正确只须证如下命题：

命题1：每一次循环结束数组B中元素总是按递增顺序排列的。

证明：用数学归纳法，对循环次数ｉ进行归纳。

　　当i=0时，即程序还没进入循环时，命题显然成立。

设i<k时命题成立，当i=k时，假设存在j1<j2,B[j₁]>B[j₂]，因为第i次循环之前数组B是递增的，因此第i次循环时B[j₁]或B[j₂]必有一个更新，假设B[j₁]被更新为元素a_i+1，由于a_i+1=B[j₁]> B[j₂]，按算法a_i+1应更新B[j₂]才对，因此产生矛盾；假设B[j₂]被更新，设更新前的元素为s,更新后的元素为a_i+1，则由算法可知第i次循环前有B[j₂]＝s< a_i+1< B[j₁]，这与归纳假设矛盾。命题得证。

命题2：B[c]中存储的元素是当前所有最长递增子序列长度为c的序列中，最小的最末元素，即设当前循环次数为i，有B[c]={a_j|f(k)=f(j)=c∧k,j≤i+1→a_j≤a_k}(f(i)为与第二种算法中的f(i)含义相同)。

证明：程序中每次用元素a_i更新B[c]时(c=f(i))，设B[c]原来的值为s，则必有a_i<s，不然a_i就能接在s的后面形成长度为c+1的最长递增子序列，而更新B[c+1]而不是B[c]了。所有B[c]中存放的总是当前长度为ｃ的最长递增子序列中，最小的最末元素。

命题3：设第i次循环后得到的p为p(i+1)，那么p(i)为以元素a_i为最末元素的最长递增子序列的长度。

证明：只须证p(i)等于第二种算法中的f(i)。显然一定有p(i)<＝f(i)。假设p(i)<f(i)，那么有两种情况，第一种情况是由二分查找法找到的p(i)不是数组B中能让a_i接在后面成为新的最长递增子序列的最大的元素，由命题1和二分查找的方法可知，这是不可能的；第二种情况是能让a_i接在后面形成长于p(i)的最长递增子序列的元素不在数组B中，由命题2可知，这是不可能的，因为B[c]中存放的是最末元素最小的长度为c的最长递增子序列的最末元素，若a_i能接在长度为L(L> p(i))的最长递增子序列后面，就应该能接在B[L]后面，那么就应该有p(i)=L,与L> p(i)矛盾。因此一定有p(i)＝f(i)，命题得证。