Cod sursa(job #642014)

Utilizator a_h1926Heidelbacher Andrei a_h1926 Data 30 noiembrie 2011 13:31:25
Problema Cutii Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.47 kb
#include <cstdio>
#include <algorithm>

#define NMax 3505

using namespace std;

typedef struct
{
    int X, Y, Z;
}
Box;
int N, AI[3*NMax][3*NMax];
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)
    {
        scanf ("%d %d %d", &B[i].X, &B[i].Y, &B[i].Z);
    }
    sort (B+1, B+N+1, Compare);
}

void UpdateY (int KX, int KY, int L, int R, int Y, int V)
{
    int Mid=(L+R)/2;
    if (R<=Y)
    {
        AI[KX][KY]=Max (AI[KX][KY], V);
        if (V==0)
        {
            AI[KX][KY]=0;
        }
        return;
    }
    if (Y<=Mid)
    {
        UpdateY (KX, 2*KY, L, Mid, Y, V);
    }
    else
    {
        UpdateY (KX, 2*KY+1, Mid+1, R, Y, V);
    }
    AI[KX][KY]=Max (AI[KX][2*KY], AI[KX][2*KY+1]);
}

void UpdateX (int K, int L, int R, int X, int Y, int V)
{
    int Mid=(L+R)/2;
    if (R<=X)
    {
        UpdateY (K, 1, 0, N, Y, V);
        return;
    }
    if (X<=Mid)
    {
        UpdateX (2*K, L, Mid, X, Y, V);
    }
    else
    {
        UpdateX (2*K+1, Mid+1, R, X, Y, V);
    }
    UpdateY (K, 1, 0, N, Y, V);
}

int QueryY (int KX, int KY, int L, int R, int Y)
{
    int Mid=(L+R)/2;
    if (R<=Y)
    {
        return AI[KX][KY];
    }
    if (Y<=Mid)
    {
        return QueryY (KX, 2*KY, L, Mid, Y);
    }
    return Max (QueryY (KX, 2*KY, L, Mid, Y), QueryY (KX, 2*KY+1, Mid+1, R, Y));
}

int QueryX (int K, int L, int R, int X, int Y)
{
    int Mid=(L+R)/2;
    if (R<=X)
    {
        return QueryY (K, 1, 0, N, Y);
    }
    if (X<=Mid)
    {
        return QueryX (2*K, L, Mid, X, Y);
    }
    return Max (QueryX (2*K, L, Mid, X, Y), QueryX (2*K+1, Mid+1, R, X, Y));
}

void LIS ()
{
    for (int i=1; i<=N; ++i)
    {
        int DP=1+QueryX (1, 0, N, B[i].X-1, B[i].Y-1);
        UpdateX (1, 0, N, B[i].X, B[i].Y, DP);
    }
}

void Clear ()
{
    for (int i=1; i<=N; ++i)
    {
        UpdateX (1, 1, N, B[i].X, B[i].Y, 0);
    }
}

void Print ()
{
    printf ("%d\n", AI[1][1]);
}

int main()
{
    freopen ("cutii.in", "r", stdin);
    freopen ("cutii.out", "w", stdout);
    int T;
    scanf ("%d %d", &N, &T);
    for (; T>0; --T)
    {
        Read ();
        LIS ();
        Print ();
        Clear ();
    }
    return 0;
}