BZOJ 4248 AAQQZ

最新推荐文章于 2019-07-16 15:22:00 发布

PoPoQQQ

最新推荐文章于 2019-07-16 15:22:00 发布

阅读量1.8k

点赞数

CC 4.0 BY-SA版权

分类专栏： BZOJ 文章标签： BZOJ BZOJ4248

本文链接：https://blog.youkuaiyun.com/PoPoQQQ/article/details/48225889

BZOJ 专栏收录该内容

667 篇文章

订阅专栏

博主历经5天成功解答BZOJ 4248题，题目要求对字符串区间排序以找到最长回文子串。限制条件包括n≤3000和字符集大小≤3000。博主分享了题解链接，并提到题中诸多细节，建议学习日语以提高解题能力，愿意通过QQ讨论解题技巧。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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

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

题解戳这里

多学日语有益于身体健康，实在想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;
}