Cod sursa(job #795208)

Utilizator radustn92Radu Stancu radustn92 Data 7 octombrie 2012 19:47:25
Problema Subsir crescator maximal Scor 70
Compilator cpp Status done
Runda Arhiva educationala Marime 1.36 kb
#include <stdio.h>
#include <algorithm>
#define NMAX 100005
using namespace std;
int n,A[NMAX],C[NMAX],ord[NMAX],aib[NMAX],Vmax,best[NMAX],rez,rez_poz;
void read()
{
	scanf("%d",&n);
	int i;
	for (i=1; i<=n; i++)
	{
		scanf("%d",&A[i]);
		ord[i]=i;
	}
}
bool comp(int x,int y)
{
	return A[x]<A[y];
}
void normal()
{
	int i;
	sort(ord+1,ord+n+1,comp);
	C[ord[1]]=1;
	for (i=2; i<=n; i++)
		if (A[ord[i]]==A[ord[i-1]])
			C[ord[i]]=C[ord[i-1]];
		else
			C[ord[i]]=C[ord[i-1]]+1;
	Vmax=C[ord[n]];
}
inline int max(int x,int y)
{
	return x>y ? x : y;
}
int lsb(int x)
{
	return x & -x;
}
int find_max(int poz)
{
	int i,s=0;
	for (i=poz; i>0; i-=lsb(i))
		s=max(s,aib[i]);
	return s;
}
void update(int poz,int val)
{
	int i;
	for (i=poz; i<=Vmax; i+=lsb(i))
		aib[i]=max(aib[i],val);
}
void solve()
{
	int i;
	for (i=1; i<=n; i++)
	{
		best[i]=find_max(C[i]-1)+1;
		update(C[i],best[i]);
		if (best[i]>rez)
			rez=best[i],rez_poz=i;
	}
	printf("%d\n",rez);
}
void recons(int poz,int left)
{
	if (left==0)
		return ;
	int i;
	for (i=poz-1; i ; i--)
		if (best[i]==left-1 && C[i]<C[poz])
		{
			recons(i,left-1);
			break ;
		}
	printf("%d ",A[poz]);
}
int main()
{
	freopen("scmax.in","r",stdin);
	freopen("scmax.out","w",stdout);
	read();
	normal();
	solve();
	recons(rez_poz,rez);
	printf("\n");
	return 0;
}