bzoj4486: [Jsoi2015]串分割

肉丝哥哥钦定好题 话说我的blog现在为什么到处都是肉丝哥哥

先来想一个弱化版，假如能够n整除K怎么做？

把每个数字看成一个字符串，按字典序排名，这个可以后缀数组解决，然后暴力枚举每种情况，O(1)判两个长度为n/K的数字大小即可

然后不能整除一定是有n-n/K*K个长度为n/K+1的数字

先二分答案排名，枚举起始位置，考虑一个贪心的做法，假如当前能够放一个长度为n/K+1的就立刻放，这样能够保证最优

why？放n/K+1和n/K的区别相当于多往前走了1步，但是多用了一组n/K+1，假如走n/K在未来某一决策能够走n/K+1，那么当前走n/K+1只需付出少用一组就可以拉回同一起跑线。其实也就是具有决策包容性

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
const int _=1e2;
const int maxn=*2e5+_;
int n,K,U;char ss[maxn];

int sa[maxn],Rank[maxn];
namespace SA
{
    int Rsort[maxn],tt[maxn];
    void getsa(int n,int m)
    {
        for(int i=;i<=n;i++)Rank[i]=ss[i]-'';

        memset(Rsort,,sizeof(Rsort));
        for(int i=;i<=n;i++)Rsort[Rank[i]]++;
        for(int i=;i<=m;i++)Rsort[i]+=Rsort[i-];
        for(int i=n;i>=;i--)sa[Rsort[Rank[i]]--]=i;

        int ln=,p=;
        while(p<n)
        {
            int k=;for(int i=n-ln+;i<=n;i++)tt[++k]=i;
            for(int i=;i<=n;i++)
                if(sa[i]-ln>)tt[++k]=sa[i]-ln;

            memset(Rsort,,sizeof(Rsort));
            for(int i=;i<=n;i++)Rsort[Rank[tt[i]]]++;
            for(int i=;i<=m;i++)Rsort[i]+=Rsort[i-];
            for(int i=n;i>=;i--)sa[Rsort[Rank[tt[i]]]--]=tt[i];

            for(int i=;i<=n;i++)tt[i]=Rank[i];

            p=;Rank[sa[]]=;
            for(int i=;i<=n;i++)
            {
                if(tt[sa[i]]!=tt[sa[i-]]||tt[sa[i]+ln]!=tt[sa[i-]+ln])p++;
                Rank[sa[i]]=p;
            }

            m=p;ln*=;
        }
    }
    void GetTrueRank(int n)
    {
        int p=;
        for(int i=;i<=*n-;i++)
            if(sa[i]<=n)Rank[sa[i]]=Rank[sa[i]+n]=++p;

        for(int i=;i<=n;i++)sa[Rank[i]]=i;
    }
    void main(int n){getsa(*n-,);GetTrueRank(n);}
}

namespace PRINT
{
    void main(int ans,int len)
    {
        int st=sa[ans];
        for(int i=;i<=len;i++)putchar(ss[i+st-]);
        puts("");
    }
}
namespace SOL1
{
    void main()
    {
        int ans=(<<);
        for(int i=;i<=U;i++)
        {
            int num=,tt=;
            for(int j=;j<=K;j++)
                num=max(num,Rank[i+tt]),tt+=U;
            ans=min(ans,num);
        }
        PRINT::main(ans,U);
    }
}
namespace SOL2
{
    int calc(int em)
    {
        int mmax=;
        for(int i=;i<=U+;i++)
        {
            int p=i,num=;
            for(int j=;j<=K;j++)
            {
                if(Rank[p]<=em)p+=U+,num++;
                else p+=U;
                if(p>=i+n-)break;
            }
            mmax=max(mmax,num);
        }
        return mmax;
    }
    void main()
    {
        int res=n-U*K;
        int el=,er=n,em,ans;
        while(el<=er)
        {
            em=(el+er)/;
            if(calc(em)>=res)
            {
                er=em-;
                ans=em;
            }
            else el=em+;
        }
        PRINT::main(ans,U+);
    }
}

int main()
{
    scanf("%d%d%s",&n,&K,ss+); U=n/K;
    for(int i=;i<n;i++)ss[i+n]=ss[i];
    SA::main(n);

    if(n%K==)SOL1::main();
    else SOL2::main();

    return ;
}