洛谷 P5012 水の数列

洛谷 P5012 水の数列

Problem

给你一个长度为\(n(n\le10^6)\)的数列，有\(T(T\le 10^3)\)组询问，每一组询问查询区间\([l,r]\)，请选择一个\(x\)，将所有\(\le x\)的数标记后得到的段数\(\in[l,r]\)，得分为每一段连续标记的区间的长度的平方和除以\(x\)，最大化得分。

\(0\le a_i\le 10^6\).

强制在线。

Solution

注意到没有修改，那么就可以愉快的预处理。

首先从小到大枚举\(x(1\le x\le 10^6)\)同时用并查集维护连通性，计算答案为\(x\)时有多少段，得分是多少。

把这个\((得分,x)\)二元组放到一个数组中(下标为段数)，用线段树维护即可。

结果大概可能是像我一样辛辛苦苦打完线段树，发现线段树无论如何都会\(MLE\)……线段树的空间复杂度是\(4\times n\)左右，而分块的空间复杂度只有\(n+\sqrt{n}\)，于是换成分块就可以愉快的AC啦！这是一种时间换空间的方法。

Code

代码实属远古且肉麻，见谅。

/**************************************************************
 * Problem: 5012
 * Author: Vanilla_chan
 * Date: 20210331
 * E-Mail: Vanilla_chan@outlook.com
**************************************************************/
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define IL inline
#define re register
#define LL long long
#define ULL unsigned long long
#ifdef TH
#define debug printf("Now is %d\n",__LINE__);
#else
#define debug
#endif
#ifdef ONLINE_JUDGE
char buf[1<<23],* p1=buf,* p2=buf,obuf[1<<23],* O=obuf;
#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
#endif
using namespace std;

namespace oi
{
    inline bool isdigit(const char& ch)
    {
        return ch<='9'&&ch>='0';
    }
    inline bool isalnum(const char& ch)
    {
        return (ch>='a'&&ch<='z')||(ch>='A'&&ch<='Z')||(ch>='0'&&ch<='9');
    }
    struct istream
    {
        char ch;
        bool fu;
        templateinline istream& operator>>(T& d)
        {
            fu=d=0;
            while(!isdigit(ch)&&ch!='-') ch=getchar();
            if(ch=='-') fu=1,ch=getchar();
            d=ch-'0';ch=getchar();
            while(isdigit(ch))
                d=(d<<3)+(d<<1)+(ch^'0'),ch=getchar();
            if(fu) d=-d;
            return *this;
        }
        inline istream& operator>>(char& ch)
        {
            ch=getchar();
            for(;!isalnum(ch);ch=getchar());
            return *this;
        }
        inline istream& operator>>(string& str)
        {
            str.clear();
            for(;!isalnum(ch);ch=getchar());
            while(isalnum(ch))
                str+=ch,ch=getchar();
            return *this;
        }
    }cin;
    inline int read()
    {
        int x=0,fu=1;
        char ch=getchar();
        while(!isdigit(ch)&&ch!='-') ch=getchar();
        if(ch=='-') fu=-1,ch=getchar();
        x=ch-'0';ch=getchar();
        while(isdigit(ch)) { x=x*10+ch-'0';ch=getchar(); }
        return x*fu;
    }
    int G[55];
    templateinline void write(T x)
    {
        int g=0;
        if(x<0) x=-x,putchar('-');
        do { G[++g]=x%10;x/=10; } while(x);
        for(int i=g;i>=1;--i)putchar('0'+G[i]);putchar('\n');
    }
};
using namespace oi;
#define N 1000010

int n,m;
int maxnum;
int a[N],f[N];
LL sze[N];
int getf(int x)
{
    if(f[x]==x) return x;
    return f[x]=getf(f[x]);
}
LL now;
int duan;
void merge(int x,int y)
{
    x=getf(x);
    y=getf(y);
    if(x==y) return;
    duan--;
    now-=sze[x]*sze[x];
    now-=sze[y]*sze[y];
    f[y]=x;
    sze[x]+=sze[y];
    sze[y]=0;
    now+=sze[x]*sze[x];
}
vectorpos[N];
struct node
{
    LL ans;
    LL x;
    node(LL defen=0,int xx=0)
    {
        ans=defen;
        x=xx;
    }
    IL bool operator>(const node &z)const
    {
        if(ans==0) return 0;
        if(z.ans==0) return 1;
        if(ans*z.x==z.ans*x) return x>z.x;
        return ans*z.x>z.ans*x;
    }
};
IL node max(const node &x,const node &y)
{
    return x>y?x:y;
}
/*
namespace Tree
{
    ::node a[N];
    struct node
    {
        ::node mx;
        #define mx(x) b[x].mx
    }b[N<<2];
    void upd(int p)
    {
        mx(p)=::max(mx(p<<1),mx(p<<1|1));
    }
    void build(int p,int l,int r)
    {
        if(l==r)
        {
            mx(p)=a[l];
            return;
        }
        int mid=l+r>>1;
        build(p<<1,l,mid);
        build(p<<1|1,mid+1,r);
        upd(p);
    }
    ::node ask(int p,int L,int R,int l,int r)
    {
        if(l<=L&&R<=r)
        {
            return mx(p);
        }
        int mid=L+R>>1;
        if(l<=mid&&r>mid) return max(ask(p<<1,L,mid,l,r),ask(p<<1|1,mid+1,R,l,r));
        if(l<=mid) return ask(p<<1,L,mid,l,r);
        else return ask(p<<1|1,mid+1,R,l,r);
    }
};
*/
namespace Block
{
    ::node a[N],mx[2010];
    int L[2010],R[2010],belong[N];
    int s,block;
    void build()
    {
        s=sqrt(n);
        block=n/s;
        for(int i=1;i<=s;i++)
        {
            L[i]=block*(i-1)+1;
            R[i]=block*i;
        }
        R[s]=n;
        for(int i=1;i<=s;i++)
        {
            for(int j=L[i];j<=R[i];j++)
            {
                belong[j]=i;
                mx[i]=max(mx[i],a[j]);
            }
        }
    }
    ::node ask(int l,int r)
    {
        ::node ans;
        if(belong[r]-belong[l]<=3)
        {
            for(int i=l;i<=r;i++) ans=::max(ans,a[i]);
        }
        else
        {
            for(int i=belong[l]+1;i<=belong[r]-1;i++) ans=::max(ans,mx[i]);
            for(int i=l;i<=R[belong[l]];i++) ans=::max(ans,a[i]);
            for(int i=L[belong[r]];i<=r;i++) ans=::max(ans,a[i]);
        }
        return ans;
    }
}
void pre()
{
    for(int i=1;i<=n;i++)
    {
        f[i]=0;
        sze[i]=0;
    }
    for(int i=1;i<=n;i++) pos[a[i]].push_back(i);
    for(int i=1;i<=maxnum;i++)
    {
        for(unsigned int j=0;jr) swap(l,r);
        ans=Block::ask(l,r);
        if(ans.ans)
        {
            cout<

Other

如果您最后一排\(RE\)/全\(WA\)，除了在计算数组大小，检查算法的正确性，别忘记看看加密方式QAQ

还有，记得多%。

时隔几个月，在大家的帮助下（见讨论区）终于完成了历史性的突破，\(95pts\to AC\)，通过了这道题。

无可奈何花落去，似曾相识燕归来。小园香径独徘徊。

分块2021准备省选

相关

题解 Count on a tree II/【模板】树分块

LOJ 数列分块入门 8

分块矩阵进行分块初等变换后秩不变的证明

【CodeForces】700 D. Huffman Coding on Segment 哈夫曼树+莫队+分块

「 Luogu P2801 」 教主的魔法——分块

#分块，可撤销并查集#洛谷 3247 [HNOI2016]最小公倍数

#分块，可撤销并查集#洛谷 3247 [HNOI2016]最小公倍数

#分块，可撤销并查集#洛谷 3247 [HNOI2016]最小公倍数

#分块，可撤销并查集#洛谷 3247 [HNOI2016]最小公倍数

牛客练习赛95Dgcd 题解(整除分块)

【题解】LibreOJ #6281. 数列分块入门 5

LOJ6277 数列分块入门 1

标签