Cod sursa(job #479386)

Utilizator irene_mFMI Irina Iancu irene_m Data 23 august 2010 20:53:12
Problema PScPld Scor 80
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.69 kb
#include <cstdio>
#include <string>
#define MaxN 1000024
#define infile "pscpld.in"
#define outfile "pscpld.out"

long long N,sol[MaxN],solp[MaxN],rez;
char a[MaxN];

void read()
{
      int i;

      scanf("%s",&a);
      N=strlen(a);
}

long long minim(long long x,long long y)
{
      if(x<y)
            return x;
      return y;
}

void solve()
{
      long long lasti,lastp,l,i;

      lasti=lastp=1;
      if(a[1]==a[0])
            solp[0]=1, rez++;

      for(i=1;i<N-1;i++)
      {
            if(a[i+1]==a[i])
                  rez++;

            if(i<lasti+sol[lasti])
                  sol[i]=minim(sol[lasti-(i-lasti)],lasti+sol[lasti]-i);

            l=sol[i];
            if(lasti+sol[lasti]<=i+sol[i])
            {
                  lasti=i;
                  while(i-l>=0 && i+l<N && a[i-l]==a[i+l])
                        l++;
                  if(l>=1)
                        sol[i]=l-1;
            }
            rez+=sol[i];


            //
            if(i<lastp+solp[lastp])
                  solp[i]=minim(solp[lastp-(i-lastp)], lastp+solp[lastp]-i);

            l=solp[i];
            if(lastp+solp[lastp]<=i+solp[i])
            {
                  lastp=i;
                  while(i-l>=0 && i+l+1<N && a[i-l]==a[i+l+1])
                        l++;
                  if(l>=1)
                        solp[i]=l-1;
            }
            rez+=solp[i];
      }
}

void write()
{
      rez+=N;
      printf("%lld\n",rez);
}

int main()
{
      freopen(infile,"r",stdin);
      freopen(outfile,"w",stdout);

      read();
      solve();
      write();

      fclose(stdin);
      fclose(stdout);
      return 0;
}