洛谷 P2178 品酒大会(bzoj P4199 [Noi2015]品酒大会/uoj P131【NOI2015】品酒大会)

这道题一看就是一种非常恶心的感觉，事实的确如此恶心到我饭都吃不下去。。

我恐怕真的不适合这个竞赛。。

这道题我想的是sa+单调栈然后全世界都比我快，全世界都用的是并查集。。

好吧说下zz的做法：首先我们想想求lcp的过程两个后缀的lcp=rmq min height(i,j)，那么我们考虑height[i]并维护一个单调递增栈，如果height[i]<=height[stack[top]],那么height[stack[top]]对后面的height[j](j>i)没有贡献那么就弹出，只要知道它待在栈内的时间新加入了多少个后缀，就能知道它会被累加多少遍。很显然打一个时间戳就能记录height[stack[top]]对多少元素累加的次数。设L[i]表示i左边有多少元素height比它大对于比他的lcm肯定等于自身的height，所以两者相乘就是对于第i个height相似度为height[i]的个数，把所有累加起来就是相似度为height[i]的个数，至于最大值呢呢就用st表维护就好了。说的简单（思考的真是无比困难。。程序也是无比难打。。菜飞了我。。

PS：介于楼主表达能力有限所以有不懂的地方可以评论哟！！

注意：1.long long 2.s[i]可能为负，对于负的情况要维护最小值，所以不一定维护最大值即最大，两个负数乘起来有可能比两个正数乘起来大，神坑。。

代码：

#include<iostream>
#include<stdio.h>
#define ll long long
using namespace std;
const int Maxn=300005;
const long long oo=2e18;
char s[Maxn];
int tax[Maxn],sa[Maxn],ranks[Maxn],height[Maxn],tp[Maxn];
int log2[Maxn],st[Maxn];
int Min[Maxn][20],Max[Maxn][20],a[Maxn],L[Maxn],times[Maxn];
long long sum[Maxn],Maxa[Maxn],Mina[Maxn],ans[Maxn];
int m=300,n,now=1,top;
inline long long max(long long x,long long y)
{
    return x>=y?x:y;
}
inline long long min(long long x,long long y)
{
    return x>=y?y:x;
}
inline void Rsort()
{
    for(int i=1;i<=m;i++)
        tax[i]=0;
    for(int i=1;i<=n;i++)
        tax[ranks[tp[i]]]++;
    for(int i=1;i<=m;i++)
        tax[i]+=tax[i-1];
    for(int i=n;i>=1;i--)
        sa[tax[ranks[tp[i]]]--]=tp[i];
}
inline int cmp(int *f,int x,int y,int w)
{
    return ((f[x]==f[y])&&(f[x+w]==f[y+w]));
}
inline void suffix()
{
    for(int i=1;i<=n;i++)
        ranks[i]=s[i],tp[i]=i;
    Rsort();
    for(int p=1,w=1,i;p<n;w+=w,m=p)
    {
        for(p=0,i=n-w+1;i<=n;i++)
            tp[++p]=i;
        for(i=1;i<=n;i++)
            if(sa[i]>w)
                tp[++p]=sa[i]-w;
        Rsort();
        swap(ranks,tp);
        ranks[sa[1]]=p=1;
        for(i=2;i<=n;i++)
            ranks[sa[i]]=cmp(tp,sa[i],sa[i-1],w)?p:++p;
    }
    for(int i=1,p,j;i<=n;height[ranks[i]]=j,i++)
        for(j=max(height[ranks[i-1]]-1,0),p=sa[ranks[i]-1];s[p+j]==s[i+j];j++);
}
inline long long rmqmin(int q1,int q2)
{
    return q1>q2?0x3f3f3f3f:min(Min[q1][log2[q2-q1+1]],Min[q2-(1<<log2[q2-q1+1])+1][log2[q2-q1+1]]);
}
inline long long rmqmax(int q1,int q2)
{
    return q1>q2?-0x3f3f3f3f:max(Max[q1][log2[q2-q1+1]],Max[q2-(1<<log2[q2-q1+1])+1][log2[q2-q1+1]]);
}
int main()
{
    scanf("%d",&n);
    scanf("%s",s+1);
    suffix();
    for(int i=0;i<=n;i++)
    {
        ans[i]=-oo;
        Maxa[i]=-oo;
        Mina[i]=oo;
    }
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    log2[0]=-1;
    for(int i=1;i<=n;i++)
        log2[i]=log2[i>>1]+1;
    for(int i=1;i<=n;i++)
    {
        Min[i][0]=a[sa[i]];
        Max[i][0]=a[sa[i]];
    }
    for(int j=1;j<=log2[n];j++)
        for(int i=1;i+(1<<j)-1<=n;i++)
        {
            Min[i][j]=min(Min[i][j-1],Min[i+(1<<j-1)][j-1]);
            Max[i][j]=max(Max[i][j-1],Max[i+(1<<j-1)][j-1]);
        }
    for(int i=2;i<=n;i++)
    {
        L[i]=i-1;
        Mina[i]=min(Mina[i],a[sa[i-1]]);
        Maxa[i]=max(Maxa[i],a[sa[i-1]]);
        while(top!=0&&height[st[top]]>=height[i])
        {
            int len=st[top]-L[st[top]],time_now=now-times[st[top]];
            L[i]=L[st[top]];
            sum[height[st[top]]]+=(long long)len*time_now;
            long long Ma=rmqmax(times[st[top]]+1,now),Mi=rmqmin(times[st[top]]+1,now);
            ans[height[st[top]]]=max(ans[height[st[top]]],Ma*Maxa[st[top]]);
            ans[height[st[top]]]=max(ans[height[st[top]]],Mi*Mina[st[top]]);
            Mina[i]=min(Mina[st[top]],Mina[i]);
               Maxa[i]=max(Mina[st[top]],Maxa[i]);
            Mina[i]=min(Maxa[st[top]],Mina[i]);
            Maxa[i]=max(Maxa[st[top]],Maxa[i]);
            top--;
        }
        st[++top]=i;
        times[i]=now;
        now++;
    }
    while(top!=0)
    {
        int len=st[top]-L[st[top]],time_now=now-times[st[top]];
        sum[height[st[top]]]+=(long long)len*time_now;
        long long Ma=rmqmax(times[st[top]]+1,now),Mi=rmqmin(times[st[top]]+1,now);
        ans[height[st[top]]]=max(ans[height[st[top]]],Ma*Maxa[st[top]]);
        ans[height[st[top]]]=max(ans[height[st[top]]],Mi*Mina[st[top]]);
        top--;
    }
    for(int i=n-1;i>=0;i--)
    {
        ans[i]=max(ans[i+1],ans[i]);
        sum[i]+=sum[i+1];
    }
    for(int i=0;i<n;i++)
        if(ans[i]==-oo)
        {
            ans[i]=0;
        }
    for(int i=0;i<n;i++)
        printf("%lld %lld\n",sum[i],ans[i]);
    return 0;
}