NOI2015品酒大会 后缀数组

自己尝试敲后缀数组，发现难看（tiao）的不行，于是抄了板子

考虑建出hei以后转化出的问题：

对于一个数组中权值大于等于k的连续部分，求取两个数的方案数和两数积的最大值

（好气啊，可以有负数）

把询问倒序以后相当于连续部分之间会动态加元素，使他们连起来

维护一段极大连续段的最大值、最小值、长度，保存在左右两端——wtf我也不知道这算什么做法

然后瞎**乱搞一下，就能每次合并O(1)维护出来了（似乎没比并查集好多少）

 #include <bits/stdc++.h>
 #define N 300001
 #define MIN -9223372036854775807
 #define calc(x) (x)*((x)-1)/2
 #define ll long long
 using namespace std;
 ll n,ans,ans2;
 ll c[N],mi[N],ma[N],len[N],w[N];
 ll sa[N],Rank[N],height[N];
 ll wa[N],wb[N],wv[N],wd[N];
 ll out1[N],out2[N];
 struct node
 {
     ll x,y;
     node(ll p,ll q)
     {
         x=p;y=q;
     }
     node()
     {
     }
 } t[N];
 bool cmp(ll *r,ll a,ll b,ll l)
 {
     return r[a]==r[b]&&r[a+l]==r[b+l];
 }
 bool com(node a,node b)
 {
     return a.x>b.x;
 }
 void da(ll *r,ll n,ll m)
 {
     ll *x=wa,*y=wb,*t;
     for (ll i=;i<m;i++) wd[i]=;
     for (ll i=;i<n;i++) wd[x[i]=r[i]]++;
     for (ll i=;i<m;i++) wd[i]+=wd[i-];
     for (ll i=n-;i>=;i--) sa[--wd[x[i]]]=i;
     for (ll j=,p=,i;p<n;j<<=,m=p)
     {
         for (p=,i=n-j;i<n;i++) y[p++]=i;
         for (i=;i<n;i++) if (sa[i]>=j) y[p++]=sa[i]-j;
         for (i=;i<n;i++) wv[i]=x[y[i]];
         for (i=;i<m;i++) wd[i]=;
         for (i=;i<n;i++) wd[wv[i]]++;
         for (i=;i<m;i++) wd[i]+=wd[i-];
         for (i=n-;i>=;i--) sa[--wd[wv[i]]]=y[i];
         for (t=x,x=y,y=t,p=,x[sa[]]=,i=;i<n;i++)
             x[sa[i]]=cmp(y,sa[i-],sa[i],j)?p-:p++;
     }
 }
 void get_height(ll *r,ll n)
 {
     for (ll i=;i<=n;i++) Rank[sa[i]]=i;
         for (ll i=,k=,j;i<n;height[Rank[i++]]=k)
             for (k?k--:,j=sa[Rank[i]-];r[i+k]==r[j+k];k++);
 }
 void merge(ll x,ll y)
 {
     if(!len[x] || !len[y]) return;
     ans2+=calc(len[x]+len[y]+)-calc(len[x]+)-calc(len[y]+);
     ans=max(ans,w[sa[x-len[x]]]*mi[y]);ans=max(ans,mi[x]*mi[y]);
     ans=max(ans,w[sa[x-len[x]]]*ma[y]);ans=max(ans,ma[x]*ma[y]);
     ll MI=min(mi[x],mi[y]),MA=max(ma[x],ma[y]);
     mi[x-len[x]+]=mi[y+len[y]-]=MI;
     ma[x-len[x]+]=ma[y+len[y]-]=MA;
     len[x-len[x]+]=len[y+len[y]-]=len[x]+len[y];
 }
 void add(ll now)
 {
     if(!now) return;
     len[now]=;
     mi[now]=ma[now]=w[sa[now]];
     ans=max(ans,w[sa[now]]*w[sa[now-]]);
     ++ans2;
     if(now)
         merge(now-,now);
     if(now<n-)
         merge(now,now+);
 }
 int main()
 {
     scanf("%lld",&n);
     for(ll i=;i<n;i++)
     {
         char ch;
         while(ch=getchar(),!isalpha(ch));
         c[i]=ch-'a'+;
     }
     for(ll i=;i<n;i++)
         scanf("%lld",&w[i]);
     c[n]=;
     da(c,n+,);
     get_height(c,n);
     for(ll i=;i<n;i++)
         sa[i]=sa[i+],height[i]=height[i+];
     if()
     for(ll i=;i<n;i++)
     {
         printf("%lld %lld\n",sa[i],height[i]);
     }
     for(ll i=;i<n;i++)
         t[i]=node(height[i],i);
     sort(t,t+n,com);
     ll j=;
     ans=MIN;ans2=;
     for(ll i=;i<n;i++)
         len[i]=,ma[i]=MIN,mi[i]=-MIN;
     for(ll i=n-;i>=;i--)
     {
         while(j<n && t[j].x==i) add(t[j].y),j++;
         out1[i]=ans2;out2[i]=(ans==MIN)?:ans;
     }
     for(ll i=;i<n;i++)
         printf("%lld %lld\n",out1[i],out2[i]);
     return ;
 }