BZOJ 4248 AAQQZ

历时5days我终于TM把这道傻逼题A掉了

题目大意：给定一个字符串，要求选择一个区间sort一下，使得到的新字符串的最长回文子串最长
n≤3000 ，字符集大小 ≤3000

题解戳这里

多学日语有益于身体健康，实在想AC这道题的加我QQ私聊吧，这题乱七八糟的细节实在太多。。。
我正在准备埋了出题人……

(翻译：这样就真的绝对能AC了！要是AC不了，把我拖到树底下埋了也没关系哦~)

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 3030
using namespace std;
int n,c,ans;
int a[M],cnt[M],f[M][M];
void This_Is_Not_Manacher()
{
    int i;
    for(i=2;i<=n<<1;i++)
    {
        int pos1=i>>1,pos2=i+1>>1;
        while( pos1>=1 && pos2<=n && a[pos1]==a[pos2] )
            --pos1,++pos2;
        ans=max(ans,pos2-pos1-1);
    }
}
void Dynamic_Programming()
{
    int i,j;
    for(i=1;i<=n;i++)
        for(j=n;j>=i;j--)
            f[i][j]=a[i]==a[j]?f[i-1][j+1]+1:0;
}
void Calculate1()
{
    static int cnt[M],sum[M];
    int i,j;
    for(i=1;i<=n<<1;i++)
    {
        int pos1=i>>1,pos2=i+1>>1;
        memset(cnt,0,sizeof cnt);
        memset(sum,0,sizeof sum);
        while( pos1>=1 && pos2<=n && a[pos1]==a[pos2] )
            --pos1,++pos2;
        for(j=pos1;j;j--)
        {
            if( j==pos1 || a[j]>=a[j+1] )
                cnt[a[j]]++;
            else
                break;
        }
        for(j=1;j<=c+1;j++)
            sum[j]=sum[j-1]+cnt[j];
        int p1=1,p2=c+1;
        while(cnt[p1]<=0&&p1<=c)
            ++p1;
        for(j=pos2;j<=n;j++)
        {
            if(--cnt[a[j]]<0)
                p2=min(p2,a[j]);
            while(cnt[p1]<=0&&p1<=c)
                ++p1;
            int fix=min(sum[p1]-cnt[p1],sum[p2]);
            ans=max(ans,pos2-pos1-1+fix*2+(j-pos2+1==fix?f[pos1-fix][pos2+fix]*2:0));
        }
    }
}
void Calculate2()
{
    static int cnt[M],sum[M];
    int i,j;
    for(i=1;i<n;i++)
    {
        int pos1=i,pos2=i+1;
        memset(cnt,0,sizeof cnt);
        memset(sum,0,sizeof sum);
        for(j=pos1;j;j--)
        {
            if( j==pos1 || a[j]>=a[j+1] )
                cnt[a[j]]++;
            else
                break;
        }
        for(j=1;j<=c+1;j++)
            sum[j]=sum[j-1]+cnt[j];
        int p1=1,p2=c+1;
        int mid=-1,mid_cnt=0;
        while(cnt[p1]<=0&&p1<=c)
            ++p1;
        for(j=pos2;j<=n;j++)
        {
            if(--cnt[a[j]]<0)
            {
                if( a[j]<=a[pos1] && mid==-1 || mid==a[j] )
                    mid=a[j],++mid_cnt;
                else
                    p2=min(p2,a[j]);
            }
            while(cnt[p1]<=0&&p1<=c)
                ++p1;
            int fix=min(sum[p1]-cnt[p1],sum[p2]);
            ans=max(ans,mid_cnt+fix*2+(j-pos2+1==mid_cnt+fix?f[pos1-fix][pos2+mid_cnt+fix]*2:0));
        }
    }
}
int main()
{
    //freopen("aaqqz.in","r",stdin);
    //freopen("aaqqz.out","w",stdout);
    int i;
    cin>>n>>c;
    for(i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
        cnt[a[i]]++;
    }
    for(i=1;i<=c;i++)
        ans=max(ans,cnt[i]);
    This_Is_Not_Manacher();
    Dynamic_Programming();
    Calculate1();
    Calculate2();
    for(i=1;i<n+1-i;i++)
        swap(a[i],a[n+1-i]);
    for(i=1;i<=n;i++)
        a[i]=c+1-a[i];
    Dynamic_Programming();
    Calculate1();
    Calculate2();
    cout<<ans<<endl;
    return 0;
}