cf1000F. One Occurrence

本文介绍了一种使用线段树数据结构解决特定数列查询问题的方法。通过预处理数列中每个元素的下一个相同元素的位置,利用线段树进行区间查询,实现了高效的查询算法。文章详细阐述了算法流程,并提供了完整的代码实现。

链接

点击跳转

题解

对于每个位置,算出其后面和它相同的数字第一次出现在哪里,不妨称为 n e x t i next_i nexti

然后离线,倒着枚举 l l l,用线段树维护 n e x t i next_i nexti的最大值

同一时刻每种权值只保留最左边的,如果查询得到的最大值大于 r r r,那就找到了一个数字,否则就无解

代码

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#define iinf 0x3f3f3f3f
#define linf (1ll<<60)
#define eps 1e-8
#define maxn 1000010
#define maxe 1000010
#define cl(x) memset(x,0,sizeof(x))
#define rep(_,__) for(_=1;_<=(__);_++)
#define em(x) emplace(x)
#define emb(x) emplace_back(x)
#define emf(x) emplace_front(x)
#define fi first
#define se second
using namespace std;
using namespace __gnu_pbds;
typedef long long ll;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
ll read(ll x=0)
{
    ll c, f(1);
    for(c=getchar();!isdigit(c);c=getchar())if(c=='-')f=-f;
    for(;isdigit(c);c=getchar())x=x*10+c-0x30;
    return f*x;
}
struct SegmentTree
{
    #define inf 0x3f3f3f3f
    int mn[maxn<<2], mx[maxn<<2], sum[maxn<<2], add[maxn<<2], set[maxn<<2], L[maxn<<2], R[maxn<<2];
    void maketag_set(int o, int v)
    {
        add[o]=0;
        set[o]=v;
        mx[o]=mn[o]=v;
        sum[o]=(R[o]-L[o]+1)*v;
    }
    void maketag_add(int o, int v)
    {
        add[o]+=v;
        mx[o]+=v, mn[o]+=v;
        sum[o]+=(R[o]-L[o]+1)*v;
    }
    void pushdown(int o)
    {
        if(L[o]==R[o])return;
        if(~set[o])
        {
            maketag_set(o<<1,set[o]);
            maketag_set(o<<1|1,set[o]);
            set[o]=-1;
        }
        if(add[o])
        {
            maketag_add(o<<1,add[o]);
            maketag_add(o<<1|1,add[o]);
            add[o]=0;
        }
    }
    void pushup(int o)
    {
        mx[o]=max(mx[o<<1],mx[o<<1|1]);
        mn[o]=min(mn[o<<1],mn[o<<1|1]);
        sum[o]=sum[o<<1]+sum[o<<1|1];
    }
    void build(int o, int l, int r, int* array)
    {
        int mid(l+r>>1);
        L[o]=l, R[o]=r;
        add[o]=0;
        set[o]=-1;
        if(l==r)
        {
            mn[o]=mx[o]=sum[o]=array[l];
            return;
        }
        build(o<<1,l,mid,array);
        build(o<<1|1,mid+1,r,array);
        pushup(o);
    }
    void segset(int o, int l, int r, int v)
    {
        int mid(L[o]+R[o]>>1);
        if(l<=L[o] and r>=R[o]){maketag_set(o,v);return;}
        pushdown(o);
        if(l<=mid)segset(o<<1,l,r,v);
        if(r>mid)segset(o<<1|1,l,r,v);
        pushup(o);
    }
    void segadd(int o, int l, int r, int v)
    {
        int mid(L[o]+R[o]>>1);
        if(l<=L[o] and r>=R[o]){maketag_add(o,v);return;}
        pushdown(o);
        if(l<=mid)segadd(o<<1,l,r,v);
        if(r>mid)segadd(o<<1|1,l,r,v);
        pushup(o);
    }
    int segsum(int o, int l, int r)
    {
        pushdown(o);
        int mid(L[o]+R[o]>>1), ans(0);
        if(l<=L[o] and r>=R[o])return sum[o];
        if(l<=mid)ans+=segsum(o<<1,l,r);
        if(r>mid)ans+=segsum(o<<1|1,l,r);
        return ans;
    }
    int segmin(int o, int l, int r)
    {
        int mid(L[o]+R[o]>>1), ans(inf);
        if(l<=L[o] and r>=R[o])return mn[o];
        pushdown(o);
        if(l<=mid)ans=min(ans,segmin(o<<1,l,r));
        if(r>mid)ans=min(ans,segmin(o<<1|1,l,r));
        return ans;
    }
    int segmax(int o, int l, int r)
    {
        int mid(L[o]+R[o]>>1), ans(-inf);
        if(l<=L[o] and r>=R[o])return mx[o];
        pushdown(o);
        if(l<=mid)ans=max(ans,segmax(o<<1,l,r));
        if(r>mid)ans=max(ans,segmax(o<<1|1,l,r));
        return ans;
    }
    #undef inf
}segtree;
int n, a[maxn], pos[maxn], nex[maxn], ans[maxn], l[maxn], r[maxn], b[maxn];
vector<int> qlis[maxn];
int main()
{
    int i, q;
    n=read();
    rep(i,n)a[i]=read();
    rep(i,500000)pos[i]=n+i;
    rep(i,500000)a[n+i]=i;
    for(i=n;i;i--)
    {
        nex[i]=pos[a[i]];
        pos[a[i]]=i;
    }
    rep(i,n)b[i]=-iinf;
    segtree.build(1,1,n,b);
    q=read();
    rep(i,q)l[i]=read(), r[i]=read(), qlis[l[i]].emb(i);
    for(i=n;i;i--)
    {
        if(nex[i]<=n)segtree.segset(1,nex[i],nex[i],-iinf);
        segtree.segset(1,i,i,nex[i]);
        for(auto q:qlis[i])
        {
            int p=segtree.segmax(1,i,r[q]);
            if(p>r[q])ans[q]=a[p];
        }
    }
    rep(i,q)printf("%d\n",ans[i]);
    return 0;
}
调一下代码(CF1000F One Occurrence,WA了): #include<bits/stdc++.h> using namespace std; int n ,q, a[500010], bl[500010], kuai_len; struct node{ int l, r, tim; friend bool operator <(node t1, node t2){ return bl[t1.l] == bl[t2.l] ? (t1.r == t2.r ? 0 : (bl[t1.l] & 1) ^ (t1.r > t2.r)) : t1.l < t2.l; } }b[500010]; int cnt, sum[500010], s[1010], ans[500010]; int js(){ if(cnt == 0) return 0; for(int i = 0; i <= 800; ++ i){ if(s[i]){ for(int j = i * kuai_len; j <= (i + 1) * kuai_len - 1; ++ j){ if(sum[j] == 1) return j; } } } } int main(){ scanf("%d", &n); kuai_len = pow(n, 0.5); for(int i = 1; i <= n; ++ i){ scanf("%d", &a[i]); bl[i] = (i - 1) / kuai_len + 1; } scanf("%d", &q); for(int i = 1; i <= q; ++ i){ scanf("%d %d", &b[i].l, &b[i].r); b[i].tim = i; } sort(b + 1, b + q + 1); int pl = 1, pr = 0; for(int i = 1; i <= q; ++ i){ while(pl > b[i].l){ -- pl; if(sum[a[pl]] == 1) s[a[pl] / kuai_len] --, -- cnt; else if(sum[a[pl]] == 0) s[a[pl] / kuai_len] ++, ++ cnt; ++ sum[a[pl]]; } while(pr < b[i].r){ ++ pr; if(sum[a[pr]] == 1) s[a[pr] / kuai_len] --, -- cnt; else if(sum[a[pr]] == 0) s[a[pr] / kuai_len] ++, ++ cnt; ++ sum[a[pr]]; } while(pl < b[i].l){ -- sum[a[pl]]; if(sum[a[pl]] == 1) s[a[pl] / kuai_len] ++, ++ cnt; else if(sum[a[pl]] == 0) s[a[pl] / kuai_len] --, -- cnt; ++ pl; } while(pr > b[i].r){ -- sum[a[pr]]; if(sum[a[pr]] == 1) s[a[pr] / kuai_len] ++, ++ cnt; else if(sum[a[pr]] == 0) s[a[pr] / kuai_len] --, -- cnt; -- pr; } ans[b[i].tim] = js(); } for(int i = 1; i <= q; ++ i){ printf("%d\n", ans[i]); } return 0; }
最新发布
03-16
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值