[GDOI模拟2015.12.19][HEOI2013]ALO

这是一篇关于OI题目的解析,主要讨论如何在给定数列a1..n中找到一段区间[l,r],使得区间次大值与区间内元素异或操作的结果最大。文章提出了通过枚举次大值、寻找合法区间以及使用可持续化Trie求解的方法,并提供了代码实现。" 83150419,7962148,插入整数到递增链表,"['数据结构', '链表操作', '排序算法']

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目大意

对于数列a1..n(ij,aiaj),找到一段区间[l,r](l<r),设区间次大值为mx2,最大化max{m2ai|i[l,r]}
1n50000,0ai109


题目分析

对于这类题目,一个很通用的思路就是枚举次大值是哪个数,然后求出合法的最大区间,再进行求解。
首先我们讨论如何找出合法区间。
显然要使ai为区间次大值,那区间内只能有一个比它大的。我们求出ai左边比ai大的最右边的数l1,次右边的数l2;同理,在右边找出比ai大的最左边的数r1,次左边的数r2。显然答案区间为(l2,r1)(l1,r2)。找出这些位置可以用二分加RMQ解决,这里不再累赘。
然后问题转化为给定一个区间,求区间内的数与给定值异或的最大值,显然,如果我们将区间所有数二进制存入Trie中,每次贪心地选择与定值对应位相反的方向走,问题就可以迎刃而解。
但是怎么求给定区间的Trie呢,我们发现这个满足区间减法,直接用可持续化Trie解决即可。


代码实现

#include <iostream>
#include <cstdio>
#include <cctype>
#include <cmath>

using namespace std;

int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch))
    {
        if (ch=='-')
            f=-1;
        ch=getchar();
    }
    while (isdigit(ch))
    {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}

const int D=30;
const int N=50000;
const int S=N*(D+1);
const int LGN=15;

bool buf[D];

struct TRIE
{
    int branch[S+1][2];
    int root[N+1];
    int sum[S+1];
    int tot;

    int newnode()
    {
        ++tot;
        branch[tot][0]=branch[tot][1]=sum[tot]=0;
        return tot;
    }

    void insert(int rt1,int &rt2,int x)
    {
        rt2=newnode();
        branch[rt2][0]=branch[rt1][0],branch[rt2][1]=branch[rt1][1],sum[rt2]=sum[rt1]+1;
        if (x==-1)
            return;
        insert(branch[rt1][buf[x]],branch[rt2][buf[x]],x-1);
    }

    int query(int rt1,int rt2,int x)
    {
        if (x==-1)
            if (sum[branch[rt2][1^buf[x]]]-sum[branch[rt1][1^buf[x]]])
                return 1;
            else
                return 0;
        if (sum[branch[rt2][1^buf[x]]]-sum[branch[rt1][1^buf[x]]])
            return (1<<x)+query(branch[rt1][1^buf[x]],branch[rt2][1^buf[x]],x-1);
        else
            return query(branch[rt1][buf[x]],branch[rt2][buf[x]],x-1);
    }
}t;

int a[N+1];
int n,ans;

struct RMQ
{
    int num[N+1][LGN+1];
    int lgn;

    void preparation()
    {
        lgn=(int)floor(log(n)/log(2));
        for (int i=1;i<=n;i++)
            num[i][0]=a[i];
        for (int j=1;j<=lgn;j++)
            for (int i=1;i<=n-(1<<j)+1;i++)
                if (num[i][j-1]>num[i+(1<<j-1)][j-1])
                    num[i][j]=num[i][j-1];
                else
                    num[i][j]=num[i+(1<<j-1)][j-1];
    }

    int query(int l,int r)
    {
        int lgrg=(int)floor(log(r-l+1)/log(2));
        if (num[l][lgrg]>num[r-(1<<lgrg)+1][lgrg])
            return num[l][lgrg];
        else
            return num[r-(1<<lgrg)+1][lgrg];
    }
}rmq;

int search_left(int ed,int x)
{
    int ret=-1,l=1,r=ed,mid;
    while (l<=r)
    {
        mid=l+r>>1;
        if (rmq.query(mid,ed)>x)
        {
            ret=mid;
            l=mid+1;
        }
        else
            r=mid-1;
    }
    return ret;
}

int search_right(int ed,int x)
{
    int ret=-1,l=ed,r=n,mid;
    while (l<=r)
    {
        mid=l+r>>1;
        if (rmq.query(ed,mid)>x)
        {
            ret=mid;
            r=mid-1;
        }
        else
            l=mid+1;
    }
    return ret;
}

int main()
{
    freopen("alo.in","r",stdin);
    freopen("alo.out","w",stdout);
    n=read();
    for (int i=1;i<=n;i++)
        a[i]=read();
    rmq.preparation();
    for (int i=1;i<=n;i++)
    {
        for (int j=0;j<D;j++)
            buf[j]=(a[i]>>j)&1;
        t.insert(t.root[i-1],t.root[i],D-1);
    }
    ans=0;
    for (int i=1;i<=n;i++)
    {
        int l1,l2,r1,r2;
        l1=search_left(i-1,a[i]);
        if (l1!=-1)
        {
            l2=search_left(l1-1,a[i]);
            if (l2==-1)
                l2=0;
        }
        else
            l1=0;
        r1=search_right(i+1,a[i]);
        if (r1!=-1)
        {
            r2=search_right(r1+1,a[i]);
            if (r2==-1)
                r2=n+1;
        }
        else
            r1=n+1;
        for (int j=0;j<D;j++)
            buf[j]=(a[i]>>j)&1;
        if (l1)
            ans=max(ans,t.query(t.root[l2],t.root[r1-1],D-1));
        if (r1!=n+1)
            ans=max(ans,t.query(t.root[l1],t.root[r2-1],D-1));
    }
    printf("%d\n",ans);
    fclose(stdin);
    fclose(stdout);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值