ETO的公开赛T5《猎杀蓝色空间号》题解

这道题别看题面这么长，其实题意很简单

就是让你求从起点开始的最长合法区间

合法的要求有两个：兜圈子和直飞

且这两个条件相互独立

（也就是说兜圈子的末尾不会对下面可能出现的直飞造成影响）

举个例子：

1 2 3 2 1 5 4 3 8 9

这个序列他的合法长度是8

因为直飞是 5 4 3 8 9

1是兜圈子的末尾，对直飞无影响

这样看来，兜圈子比直飞优秀的多

因为如果直飞的某段属于兜圈子

那么把这一段归于兜圈子后对序列的合法性无影响

但如果兜圈子的某段属于直飞，那归于直飞后

剩下的这部分可能就不是兜圈子了

所以我们先用manacher跑出回文区间

用差分数组维护一下回文区间

如果在回文区间内，则不考虑

否则判断是否直飞即可

只要确保完全符合，就记下来这个区间和这个区间的末尾值

后面的就是贪心了

最优情况自然是你选的那条正好是飞船所在

飞船所在正好是最短的

最坏情况是你最后一次选到的是最短的

且这个位置是飞船所在

几率我们会惊讶的发现其实是一样的

排个序第二个也很好求

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define rii register int i
#define rij register int j
#define mod 19260817
using namespace std;
char s[],cp[];
int p[],cd,t,n,ans,cnt,jl[],cf[],lnt;
char zcq1[];
int zcq2[],m;
void ycl()
{
    int k=;
    cp[k]=;
    k++;
    for(rii=;i<=cd-;i++)
    {
        cp[k]='#';
        k++;
        cp[k]=s[i];
        k++;
    }
    cp[k]='#';
    k++;
    cd=k;
}
void manacher()
{
    ycl();
    int mx=;
    int last=;
    for(rii=;i<=cd-;i++)
    {
        if(mx>i)
        {
            p[i]=min(p[*last-i],mx-i);
        }
        else
        {
            p[i]=;
        }
        while(cp[i+p[i]]==cp[i-p[i]])
        {
            p[i]++;
        }
        if(p[i]+i>mx)
        {
            mx=p[i]+i;
            last=i;
        }
    }
}
void cfs()
{
    for(rii=;i<=cnt;i++)
    {
        if(zcq2[i]>)
        {
            cf[i-(zcq2[i])/]++;
            cf[i+(zcq2[i])/]--;
        }
    }
}
int yz()
{
    int sl=,pre=,pd=;
    for(rii=;i<=cnt;i++)
    {
        sl+=cf[i];
        if(sl==&&cf[i]>=)
        {
            if(pre==)
            {
                pre=zcq1[i];
                continue;
            }
            if(pd==)
            {
                if(zcq1[i]>pre)
                {
                    pd=;
                }
                else
                {
                    if(zcq1[i]==pre)
                    {
                        pd=;
                    }
                    else
                    {
                        pd=;
                    }
                }
                pre=zcq1[i];
                continue;
            }
            if(pd==)
            {
                if(pre>zcq1[i])
                {
                    return i;
                }
                else
                {
                    pre=zcq1[i];
                }
            }
            if(pd==)
            {
                if(pre<zcq1[i])
                {
                    return i;
                }
                else
                {
                    pre=zcq1[i];
                }
            }
        }
        else
        {
            pre=;
        }
    }
    return cnt+;
}
long long pw(long long ds,int ms)
{
    long long res=;
    while(ms!=)
    {
        if(ms%==)
        {
            ms--;
            res*=ds;
            res%=mod;
        }
        else
        {
            ds*=ds;
            ds%=mod;
            ms/=;
        }
    }
    return res;
}
bool cmp(int kl,int lk)
{
    return kl<lk;
}
void sr()
{
    memset(s,,sizeof(s));
    memset(p,,sizeof(p));
    memset(cf,,sizeof(cf));
    for(rij=;j<t;j++)
    {
        s[j]=getchar();
    }
    getchar();
    getchar();
}
void fz()
{
    for(rij=;j<=cd-;j+=)
    {
        cnt++;
        zcq1[cnt]=cp[j];
        zcq2[cnt]=p[j]-;
    }
}
int main()
{
    scanf("%d %d",&n,&t);
    scanf("%d ",&m);
    for(rii=;i<=n;i++)
    {
           sr();
        cd=t;
        manacher();
        cnt=;
        fz();
        cfs();
        int bnt=cnt;
        int ltt=yz();
        if(ltt==m+)
        {
            lnt++;
            jl[lnt]=pw(zcq1[bnt],zcq1[bnt]);
        }
//        printf("%d\n",ans-1);
    }
    if(lnt==)
    {
        cout<<"-1";
        return ;
    }
    if(lnt==)
    {
        cout<<jl[]<<" "<<jl[]<<" 1.00000000 1.00000000";
        return ;
    }
    sort(jl+,jl+lnt+,cmp);
    int minx=jl[];
//    jl[1];
    int ant=;
    long long js=;
    for(rii=;i<=lnt;i++)
    {
        if(jl[i]==minx)
        {
            ant++;
        }
        else
        {
            break;
        }
    }
    for(rii=;i<=lnt;i++)
    {
        js+=jl[i];
    }
    js*=;
    js+=minx;
    cout<<minx<<" "<<js<<" ";
    double kn=((1.0)*ant)/lnt;
    kn/=lnt;
    printf("%.8lf %.8lf",kn,kn);
}