求最长上升子序列的长度

今天做到一个题需要用到最大上升子序列，然后就学了下这个算法，用到了二分查找复杂度为O(nlogn)，这里有这个算法的思想，写得挺好，http://www.slyar.com/blog/longest-ordered-subsequence.html，然后我自己整理了下对代码原理的理解，看了思想再看下我的整理，应该效果不错

#include<stdio.h>
#include<string.h>
#define MAXN 10010
int a[MAXN];
int dp[MAXN];
int main()
{
	int n,i;
	while(scanf("%d",&n)!=EOF)
	{
		memset(dp,0,sizeof(dp));
		for(i=1;i<=n;i++)
			scanf("%d",&a[i]);
		int ans=1;
		//最长递减
		//dp[1]=a[n];
		//for(i=n;i>=1;i--)
		dp[1]=a[1];
		for(i=1;i<=n;i++)
		{
			int l=1,r=ans;
			while(l<=r)
			{
				int m=(l+r)>>1;
				if(a[i]<=dp[m])//这个二分查找的目的是找到第一个大于或等于a[i]的数
							   //dp[]本身是递增的，如果a[i]比dp的顶部元素还大，那么循环终止后会有l==ans+1
							   //如果a[i]不比顶部元素大，那么肯定能找到第一个>=a[i]的元素，
							   //循环终止后这个元素就是dp[l],  r为l-1
					r=m-1;
				else
					l=m+1;
			}
			dp[l]=a[i];//如果是第一种情况，那么会在末尾添加a[i]
						//第二种情况会找到l满足dp[l]>=a[i],在dp[l]处用a[i]覆盖之前的数字
			if(l>ans)//只会在第一种情况出现l==ans+1，相当于当前的上升子序列的长度增加了1
				ans++;
		}
		printf("%d\n",ans);
	}
	return 0;
}