最长上升序列（Lis）

Description

A numeric sequence of ai is ordered if a1 < a2 < ... < aN. Let the subsequence of the given numeric sequence (a1, a2, ..., aN) be any sequence (ai1, ai2, ..., aiK), where 1 <= i1 < i2 < ... < iK <= N. For example, the sequence (1, 7, 3, 5, 9, 4, 8) has ordered subsequences, e.g., (1, 7), (3, 4, 8) and many others. All longest ordered subsequences of this sequence are of length 4, e.g., (1, 3, 5, 8).
Your program, when given the numeric sequence, must find the length of its longest ordered subsequence.
给出一个长度为N的数字串，找出一个严格上升的数字序列来.

Input

The first line of input contains the length of sequence N (1 <= N <= 1000). The second line contains the elements of sequence - N integers in the range from 0 to 10000 each, separated by spaces.

Output

Output must contain a single integer - the length of the longest ordered subsequence of the given sequence.

Sample Input

7
1 7 3 5 9 4 8

Sample Output

这道题对于现在的我来说已经算是水题了，但是当年我还是看了半个小时才看懂。

f[i]表示到i最长的上升序列

主要的思路就是枚举每个点，然后再与后面的点比较，加入后面的点比他大，长度就+1，即f[j]>f[i],f[j]=f[i]+1;(当然前提是f[i]+1要比原本的f[j]要大)

代码：

 1 #include<cstdio>
 2 #include<algorithm>
 3 using namespace std;
 4 int s[10002],n,a[10002],ans;
 5 int main()
 6 {
 7     scanf("%d",&n);
 8     for(int i=1;i<=n;i++)
 9     {
10         scanf("%d",&a[i]);
11         s[i]=1;
12     }
13     for(int i=1;i<=n;i++)
14         for(int j=i+1;j<=n;j++)
15             if(a[i]<a[j]&&s[i]+1>=s[j])s[j]=s[i]+1,ans=max(ans,s[j]);
16     printf("%d",ans);
17 }

这是O(N^2)的算法，当数据大一些的时候就不行了。

众所皆知，大多数O(N^2)的算法可以用二分优化到O(N log(N))。

没错，就是我们可以用一个数组来存，但是这个数组存的并不是答案，只是当前形成的上升序列，每进来一个数，都用二分查找第一个比他小的数，然后取而代之，如果没有比他小的数，就放到最后面，数组最后的元素个数就是答案。

代码：

 1 #include<bits/stdc++.h>
 2 using namespace std;
 3 int f[40001];
 4 int t,m,n;
 5 int main()
 6 {
 7     int m;
 8     cin>>m;
 9     for(int i=1;i<=m;i++)
10     {
11         int ans=0;
12         cin>>n;
13         for(int i=1;i<=n;i++)
14         {
15  
16             cin>>t;
17             if(i==1) f[++ans]=t;
18             else
19             {
20                 if(t>f[ans]) f[++ans]=t;
21                 else
22                 {
23                     int x=lower_bound(f+1,f+ans,t)-f;
24                     f[x]=t;
25                 }
26             }
27         }
28         cout<<ans<<endl;
29     }
30     return 0;
31 }

其实还有树状数组的做法，这里不给出了