Cod sursa(job #636745)

Utilizator a_h1926Heidelbacher Andrei a_h1926 Data 19 noiembrie 2011 23:02:54
Problema PalM Scor 30
Compilator cpp Status done
Runda .com 2011 Marime 1.68 kb
#include <cstdio>
#include <cstring>

#define NMax 505

using namespace std;

char String[NMax];
int N, DP[NMax][NMax], MaxS;
int AI[3*NMax][255];

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

void Read ()
{
    freopen ("palm.in", "r", stdin);
    scanf ("%s", String);
    N=strlen (String);
}

void Update (int K, int L, int R, int P, int C, int V)
{
    int Mid=(L+R)/2;
    if (L==R)
    {
        for (int c=C; c>=(int)'a'; --c)
        {
            AI[K][c]=Max (AI[K][c], V);
        }
        return;
    }
    if (P<=Mid)
    {
        Update (2*K, L, Mid, P, C, V);
    }
    else
    {
        Update (2*K+1, Mid+1, R, P, C, V);
    }
}

int Query (int K, int L, int R, int Q, int V)
{
    int Mid=(L+R)/2;
    if (L==R)
    {
        return AI[K][V];
    }
    if (Q<=Mid)
    {
        return Query (2*K, L, Mid, Q, V);
    }
    return Query (2*K+1, Mid+1, R, Q, V);
}

void Solve ()
{
    for (int i=N-1; i>=0; --i)
    {
        DP[i][i]=1;
        for (int j=i+1; j<N; ++j)
        {
            if (String[i]!=String[j])
            {
                continue;
            }
            DP[i][j]=2+Query (1, 0, N-1, j-1, (int)String[i]);
            MaxS=Max (MaxS, DP[i][j]);
        }
        for (int j=i; j<N; ++j)
        {
            if (String[i]!=String[j])
            {
                continue;
            }
            Update (1, 0, N-1, j, (int)String[j], DP[i][j]);
        }
    }
}

void Print ()
{
    freopen ("palm.out", "w", stdout);
    printf ("%d\n", MaxS);
}

int main()
{
    Read ();
    Solve ();
    Print ();
    return 0;
}