LIS最长不下降（上升或持平）子序列 n * logn版

输入多个数字，直到EOF，计算其LIS

我们暂且把n * logn算法叫做“歪门邪道”，其机理如下：
“初始时dp[0] = s[0]， 从i = 1时遍历原数列， 将每个数与dp数列的末尾的数进行比较， 如果大于末尾的数， 就把s[i]放在dp数列的最后， 如果小于末尾的数，那么就用s[i]替换掉dp中第一个比s[i]大的数。最后dp数列的长度就是LIS的长度，但是要注意dp中的数不是LIS的元素。”
嗯，之所以可以这样做，是因为可以保证“歪门邪道”子序列末元素是LIS的元素（你写个例子体会一下），必定存在那么一个序列，以“歪门邪道”子序列末元素为末尾，且符合“不下降”规则…你说琢磨这个算法的人他是咋想出来的…

#include <cstdio>
#include <iostream>
#include <cmath>
#include <algorithm>
#include <cstring>
using namespace std;

int main()
{
    int s[25];
    int dp[25];
    memset(s, 0, sizeof(s));
    memset(dp, 0, sizeof(dp));
    int k = 0;
    int a;
    while(cin >> a)
    {
        s[k] = a;
        k++;
//        char b = getchar();
//        if(b == '\n' )
//            break;
    }
    dp[0] = s[0];
    int len = 0;
    for(int i = 1; i < k; ++i)
    {
        if(s[i] >= dp[len])
            dp[++len] = s[i];
        else
            *lower_bound(dp, dp + len + 1, s[i]) = s[i];
        if(s[i] <= dp[len])
            dp[++len] = s[i];
//        else
//        {
//            for(int j = 0; j < len + 1; ++j)
//            {
//                if(dp[j] < s[i])
//                {
//                    dp[j] = s[i];
//                    break;
//                }
//            }
//        }
    }
//    for(int i = 0; i < len + 1; ++i)
//    {
//        cout << dp[i] << " ";
//    }
//    cout << '\n';
    cout << len + 1 << '\n';
    return 0;
}