Cod sursa(job #642301)

Utilizator a_h1926Heidelbacher Andrei a_h1926 Data 30 noiembrie 2011 21:56:58
Problema Cutii Scor 40
Compilator cpp Status done
Runda Arhiva de probleme Marime 3.17 kb
#include <fstream>
#include <algorithm>

#define NMax 3505

using namespace std;

ifstream fin ("cutii.in");
ofstream fout ("cutii.out");

typedef struct
{
    int X, Y, Z;
}
Box;
int N, AIB[NMax][NMax], S;
Box B[NMax];

inline int Max (int a, int b)
{
    if (a>b)
    {
        return a;
    }
    return b;
}

inline bool Compare (Box a, Box b)
{
    if (a.X<b.X)
    {
        return true;
    }
    return false;
}

void Read ()
{
    for (int i=1; i<=N; ++i)
    {
        fin >> B[i].X >> B[i].Y >> B[i].Z;
    }
    sort (B+1, B+N+1, Compare);
}

inline int LSB (int X)
{
    return X&(-X);
}

void Update (int X, int Y, int V)
{
    for (int i=X; i<=N; i+=LSB (i))
    {
        for (int j=Y; j<=N; j+=LSB (j))
        {
            AIB[i][j]=Max (AIB[i][j], V);
            if (V==0)
            {
                AIB[i][j]=0;
            }
        }
    }
}

int Query (int X, int Y)
{
    int Q=0;
    for (int i=X; i>0; i-=LSB (i))
    {
        for (int j=Y; j>0; j-=LSB (j))
        {
            Q=Max (Q, AIB[i][j]);
        }
    }
    return Q;
}

void LIS ()
{
    S=0;
    for (int i=1; i<=N; ++i)
    {
        int DP=1+Query (B[i].Y-1, B[i].Z-1);
        Update (B[i].Y, B[i].Z, DP);
        S=Max (S, DP);
    }
}

void ClearAIB ()
{
    for (int i=1; i<=N; ++i)
    {
        Update (B[i].Y, B[i].Z, 0);
    }
}

void Print ()
{
    fout << S << "\n";
}

int main()
{
    int T;
    fin >> N >> T;
    for (; T>0; --T)
    {
        Read ();
        LIS ();
        Print ();
        ClearAIB ();
    }
    fin.close ();
    fout.close ();
    return 0;
}

/*
#include<stdio.h>
#include<algorithm>

#define maxn 3505

using namespace std;

FILE*f=fopen("cutii.in","r");
FILE*g=fopen("cutii.out","w");

int n,t,i,j,best,crt;
int aib[maxn][maxn];

struct cutie{
	int x,y,z;
}A[maxn];

struct cmp{
	inline bool operator () ( const cutie &a , const cutie &b ){
		return a.x < b.x;
	}
};

inline int lsb ( int x ){
	return x & -x;
}

inline void update_aib2D ( int x , int y , int val ){

	for ( int i = x ; i <= n ; i += lsb(i) ){
		for ( int j = y ; j <= n ; j += lsb(j) ){
			if ( val > aib[i][j] )
				aib[i][j] = val;
		}
	}
}

inline int query_aib2D ( int x , int y ){
	int bst = 0;

	for ( int i = x ; i > 0 ; i -= lsb(i) ){
		for ( int j = y ; j > 0 ; j -= lsb(j) ){
			if ( aib[i][j] > bst )
				bst = aib[i][j];
		}
	}

	return bst;
}

inline void set0_aib2D ( int x , int y ){

	for ( int i = x ; i <= n ; i += lsb(i) ){
		for ( int j = y ; j <= n ; j += lsb(j) ){
			aib[i][j] = 0;
		}
	}
}

int solve () {

	for ( i = 1 ; i <= n ; ++i ){
		fscanf(f,"%d %d %d",&A[i].x,&A[i].y,&A[i].z);
	}
	sort(A+1,A+n+1,cmp());

	best = 0;
	for ( i = 1 ; i <= n ; ++i ){
		crt = query_aib2D(A[i].y-1,A[i].z-1) + 1;
		if ( crt > best )	best = crt;
		update_aib2D(A[i].y,A[i].z,crt);
	}

	for ( i = 1 ; i <= n ; ++i ){
		set0_aib2D(A[i].y,A[i].z);
	}

	return best;
}

int main () {

	fscanf(f,"%d %d",&n,&t);
	for ( int i = 1 ; i <= t ; ++i ){
		fprintf(g,"%d\n",solve());
	}

	fclose(f);
	fclose(g);

	return 0;
}*/