标签：int top st read 1ll HNOI2016 序列 s1

传送门

Description

有 \(q\) 个询问，每个询问给定两个数\(l\) 和\(r\)，求 \(a[l:r]\) 的不同子序列的最小值之和

Solution 

校内模拟赛用了这道题，但是莫队只能拿\(80\)分，正解是猫树

当然还是莫队啦

考虑一个数加入时的贡献，就是以它为端点的区间的贡献

发现可以将现有的区间分成一段一段的，每段的最小值不同

这个可以用单调栈\(+\)倍增解决

于是就有了\(O(n\sqrt n\log n)\)的做法

但是显然过不了原题

所以想办法把\(\log\)去掉

以右端点为例，令\(f[i]\)表示\(a_i\)对区间\([1,i-1]\)的贡献

发现答案可以看成是\(f[i]-f[k],(k\leq i)\)然后再加上\(k\)以后的贡献

但是，这里的\(k\)得满足它是端点之一，且它往左的贡献很好算

那是什么？就是区间最小值啦，它往左的贡献就是\(a[k]*(k-nowl+1)\)

所以可以先预处理所有的\(f[i]\)，用\(RMQ\)算法\(O(1)\)求出区间最小值就行了

复杂度\(O(n\sqrt n)\)

附上不能看的代码

Code 

#include<bits/stdc++.h>
#define ll long long
#define db double
#define dbg1(x) cerr<<#x<<"="<<(x)<<" "
#define dbg2(x) cerr<<#x<<"="<<(x)<<"\n"
#define dbg3(x) cerr<<#x<<"\n"
#define int ll
using namespace std;
inline int read()
{
    int f=1,x=0;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
    return x*f;
}
const int MN=1e5+5,Mod=1e9+7;
int A,B,C,P;
ll lastAns=0;
inline int rnd(){return A=(A*B+(C^(int)(lastAns&0x7fffffffll))%P)%P;}
int n,q,l,r,a[MN];
bool f1,f2;
namespace solve1
{
    int L[MN],R[MN],st[MN],top;
    void Main()
    {
        int i;
        st[top=0]=0;
        for(i=1;i<=n;++i)
        {
            while(top&&a[st[top]]>=a[i]) --top;
            L[i]=st[top]+1;
            st[++top]=i;
        }
        st[top=0]=n+1;
        for(i=n;i;--i)
        {
            while(top&&a[st[top]]>a[i]) --top;
            R[i]=st[top]-1;
            st[++top]=i;
        }
        while(q--)
        {
            l=read(),r=read();
            ll ans=0ll;
            for(i=l;i<=r;++i)
                ans+=1ll*a[i]*1ll*(min(r,R[i])-i+1)*1ll*(i-max(l,L[i])+1);
            printf("%lld\n",ans);
        }
    }
}
ll suf1[MN],suf2[MN];
void init()
{
    int i;suf1[n+1]=suf2[n+1]=0ll;
    for(i=n;i;--i)
        suf1[i]=suf1[i+1]+1ll*a[i]*(n-i+1),
        suf2[i]=suf2[i+1]+a[i];
}
ll cal(int x,int y){return suf1[x]-suf1[y+1]-1ll*(suf2[x]-suf2[y+1])*1ll*(n-y);}
struct Ques{int l,r,id;}_[MN];int tt;
namespace solve2
{
    int Bl,bel[MN],nowl,nowr;
    ll ans[MN],Sum,nowans;
    ll s1[20][MN],s2[20][MN],p1[20][MN],p2[20][MN],f[MN],g[MN];
    bool cmp(Ques x,Ques y)
    {
        return bel[x.l]!=bel[y.l]?bel[x.l]<bel[y.l]
            :(bel[x.l]&1?x.r<y.r:x.r>y.r);
    }
    int st[MN],top;
    ll cal(int x,bool ty)
    {
        int tmp=0;
        if(!ty)
        {
            for(int i=17;~i;--i)
                if(s1[i][x]==nowr) return tmp+s2[i][x];
                else if(s1[i][x]<nowr)tmp+=s2[i][x],x=s1[i][x]+1;
            return tmp+1ll*a[x]*(nowr-x+1);
        }
        else
        {
            for(int i=17;~i;--i)
                if(p1[i][x]==nowl) return tmp+p2[i][x];
                else if(p1[i][x]>nowl)tmp+=p2[i][x],x=p1[i][x]-1;
            return tmp+1ll*a[x]*(x-nowl+1);
        }
    }
    int Log[MN],mi[20][MN];
    int cmp_(int x,int y){return a[x]<a[y]?x:y;}
    int gmin(int x,int y){int og=Log[y-x+1];return cmp_(mi[og][x],mi[og][y-(1<<og)+1]);}
    void Init()
    {
        int i,j;
        for(i=2;i<=n;++i) Log[i]=Log[i>>1]+1;
        for(i=1;i<=n;++i) mi[0][i]=i;
        for(j=1;j<18;++j)for(i=1;i<=n;++i)if(i-1+(1<<j)<=n)mi[j][i]=cmp_(mi[j-1][i],mi[j-1][i+(1<<j-1)]);
        st[top=0]=0;
        for(i=1;i<=n;++i)
        {
            while(top&&a[st[top]]>=a[i]) --top;
            p1[0][i]=st[top]+1;
            p2[0][i]=1ll*(i-p1[0][i]+1)*a[i];
            st[++top]=i;
        }
        st[top=0]=n+1;
        for(i=n;i;--i)
        {
            while(top&&a[st[top]]>a[i]) --top;
            s1[0][i]=st[top]-1;
            s2[0][i]=1ll*(s1[0][i]-i+1)*a[i];
            st[++top]=i;
        }
        for(i=1;i<18;++i)
            for(j=1;j<=n;++j)
            {
                if(s1[i-1][j]>=n) s1[i][j]=n+1;
                else s1[i][j]=s1[i-1][s1[i-1][j]+1],
                     s2[i][j]=s2[i-1][j]+s2[i-1][s1[i-1][j]+1];
                if(p1[i-1][j]<=1) p1[i][j]=0;
                else p1[i][j]=p1[i-1][p1[i-1][j]-1],
                     p2[i][j]=p2[i-1][j]+p2[i-1][p1[i-1][j]-1];
            }
        nowl=1;nowr=n;
        for(i=1;i<=n;++i) f[i]=cal(i,1),g[i]=cal(i,0);
    }
    ll Cal(int x,int ty)
    {
        ll tmp=0ll;
        if(ty)
        {
            int pos=gmin(nowl,x);
            tmp=f[x]-f[gmin(nowl,x)];
            tmp+=1ll*a[pos]*(pos-nowl+1);
        }
        else
        {
            int pos=gmin(x,nowr);
            tmp=g[x]-g[gmin(x,nowr)];
            tmp+=1ll*a[pos]*(nowr-pos+1);
        }
        return tmp;
    }
    void Add(int x,bool ty){nowans+=Cal(x,ty);}
    void Del(int x,bool ty){nowans-=Cal(x,ty);}
    void Main()
    {
        Init();Bl=sqrt((db)(n));
        for(int i=1;i<=n;++i)bel[i]=(i+Bl-1)/Bl;
        std::sort(_+1,_+tt+1,cmp);
        nowl=1,nowr=0;Sum=0;nowans=0;
        for(int i=1;i<=tt;++i)
        {
            while(nowr<_[i].r) Add(++nowr,1);
            while(nowl>_[i].l) Add(--nowl,0);
            while(nowr>_[i].r) Del(nowr--,1);
            while(nowl<_[i].l) Del(nowl++,0);
            if(f2) (Sum+=nowans)%=Mod;
            else ans[_[i].id]=nowans;
        }
        if(f2) printf("%lld\n",Sum);
        else for(int i=1;i<=tt;++i) printf("%lld\n",ans[i]);
    }
}
signed main()
{   
    n=read();q=read();
    if(f2) A=read(),B=read(),C=read(),P=read();
    int i;bool c1=1,c2=1;a[0]=-Mod;
    for(i=1;i<=n;++i) a[i]=read(),c1&=(a[i]==a[1]),c2&=(a[i]>a[i-1]);
    if(c2) init();
    if(n*q<=Mod&&!f1&&!f2) {solve1::Main();return 0;}
    while(q--)
    {
        if(f2)
        {
            l=rnd()%n+1;r=rnd()%n+1;
            if(l>r)swap(l,r);
        }
        else if(f1)
        {
            l=read();r=read();
            l=((l+lastAns)%n+n)%n+1;
            r=((r+lastAns)%n+n)%n+1;
            if(l>r)swap(l,r);
        }
        else l=read(),r=read();
        if(c1) printf("%lld\n",lastAns=1ll*(1ll*(r-l+1)*1ll*(r-l+2)/2)*1ll*a[1]);
        else if(c2) printf("%lld\n",lastAns=cal(l,r));
        else _[++tt].l=l,_[tt].r=r,_[tt].id=tt;
    }
    if(!c1&&!c2){solve2::Main();return 0;}
    return 0;
}

Blog来自PaperCloud，未经允许，请勿转载，TKS！

标签：int,top,st,read,1ll,HNOI2016,序列,s1
来源： https://www.cnblogs.com/PaperCloud/p/11341650.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。