Cod sursa(job #736234)
#include <vector>
#include <fstream>
#include <cstdlib>
#include <iterator>
#include <algorithm>
#define N_MAX 100011
using namespace std;
ofstream out("scmax.out");
vector<int> v, w;
vector<int>::iterator it, iend;
int index[N_MAX], l[N_MAX], f[N_MAX], aib[N_MAX];
inline void Update(int& N, int xIndex, int& xValue)
{
for(; xIndex <= N; xIndex+=(xIndex & -xIndex) )
if(l[aib[xIndex]] < l[xValue])
aib[xIndex]=xValue;
}
inline int Query(int xIndex)
{
int m=0;
for(; xIndex; xIndex-=(xIndex & -xIndex) )
if(l[aib[xIndex]] > l[m])
m=aib[xIndex];
return m;
}
int main()
{
int N, i, j, maxIndex=0;
ifstream in("scmax.in");
in>>N;
v.push_back(0);
copy(istream_iterator<int>(in), istream_iterator<int>(), back_inserter(v));
w=v;
it=w.begin(), iend=w.end();
sort(it, iend);
iend=unique(it, iend);
for(i=1; i <= N; ++i)
index[i]=lower_bound(it, iend, v[i])-it;
for(i=1; i <= N; ++i)
{
f[i]=Query(index[i]-1);
l[i]=l[f[i]]+1;
Update(N, index[i], i);
if(l[i] > l[maxIndex])
maxIndex=i;
}
out<<l[maxIndex]<<'\n';
for(i=maxIndex, j=l[maxIndex]; i; i=f[i], --j)
w[j]=v[i];
copy(w.begin()+1, w.begin()+l[maxIndex]+1, ostream_iterator<int>(out, " "));
return EXIT_SUCCESS;
}