hdu 4902 线段树双标记类型题

本文提供了一种解决HDU 4902问题的有效方法,通过使用段树和懒惰传播技术来优化对区间进行的两种操作:更新为指定值和更新为最大公约数。提供了完整的C++代码实现。

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

比赛时没敲出来,似乎现在敲得还是不好

后来找了份题解,http://www.cnblogs.com/qq1012662902/p/3883614.html  600多msAC  很高效了


参考了之后写了代码  近期在重写一遍

</pre><pre name="code" class="cpp">#include <cstdio>
#include <cstring>
#include <algorithm>
#include <string>
#include <iostream>
using namespace std;


#define ls(rt) rt*2
#define rs(rt) rt*2+1
#define ll long long
#define rep(i,s,e) for(int i=s;i<e;i++)
#define repe(i,s,e) for(int i=s;i<=e;i++)
#define IN(s) freopen(s,"r",stdin)
const int MAXN = 100000+100;


int gcd(int a, int b)
{
    if(!a)return b;////
    if(!b)return a;
    return gcd(b,a%b);
}
int num[MAXN];
struct Node{
    int l,r;
    int mx,val;
    //mx记录最大值,用于剪枝
    //val为-1是没有lazy,否则存储需要更新的值v--操作1
    //注意因为更新的值是大于等于0的,所以可以这么用
}nodes[MAXN*4];




void pushup(int rt)
{
    nodes[rt].mx=max(nodes[ls(rt)].mx,nodes[rs(rt)].mx);
}


void build(int rt, int l, int r)
{
    nodes[rt].l=l;
    nodes[rt].r=r;
    nodes[rt].mx=0;
    nodes[rt].val=-1;
    if(l == r)
    {
        nodes[rt].mx=nodes[rt].val=num[l];
        return;
    }
    int mid=(l+r)/2;
    build(ls(rt),l,mid);
    build(rs(rt),mid+1,r);
    pushup(rt);
}


void modify(int rt,int v)//2操作
{
    if(nodes[rt].mx<=v)return;
    if(nodes[rt].val!=-1)//仍然更新到段--到仅仅有一种更新lazy标记
    {
        nodes[rt].mx=nodes[rt].val=gcd(nodes[rt].val,v);
        //nodes[rt].val=-1;
        return;
    }
    modify(ls(rt),v);
    modify(rs(rt),v);
    pushup(rt);
}


void update(int rt, int l, int r,int op,int v)
{
    //找到区间之后,根据操作做更新
    //一种操作的话,当找到对应的
    //区间(nodes[rt].l==l && nodes[rt].r==r)之后,
    //只标记而不往下更新,直到下次update到这次才将这次的pushdown
    //此题中,操作一是赋值,会覆盖操作二,所以直接操作
    //如果是操作二,查看之前是不是已经有lazy操作,有的话,先做之前的lazy,
    //直到达到一个结点,该节点没有做过lazy标记,那么标记这个节点,停止update‘


    if(nodes[rt].l==l && nodes[rt].r==r)
    {
        if(op==1)
        {
            nodes[rt].val=v;
            nodes[rt].mx=v;//该区间所有的值都变成了v,所以最大值为v
        }
        else
        {
            modify(rt,v);
        }
        return;
    }
    //pushdown
    if(nodes[rt].val!=-1)
    {
        //此处把mx也更新了,因为gcd(a,b)<=a&&gcd(a,b)<=b
        //又因为题目中说了,当要更新的val<num[i]时才更新num[i]为gcd();
        //
        nodes[ls(rt)].mx=nodes[rs(rt)].mx=nodes[rt].val;
        nodes[ls(rt)].val=nodes[rs(rt)].val=nodes[rt].val;
        nodes[rt].val=-1;
    }
    int mid=(nodes[rt].l+nodes[rt].r)/2;
    if(r<=mid)update(ls(rt),l,r,op,v);
    else
    {
        if(l>mid)
            update(rs(rt),l,r,op,v);
        else
        {
            update(ls(rt),l,mid,op,v);
            update(rs(rt),mid+1,r,op,v);
        }
    }
    pushup(rt);
}


void dfs(int rt)
{
    if(nodes[rt].l==nodes[rt].r)
    {
        num[nodes[rt].l]=nodes[rt].val;
        return;
    }
    if(nodes[rt].val!=-1)
    {
        nodes[ls(rt)].val=nodes[rs(rt)].val=nodes[rt].val;
        nodes[rt].val=-1;
    }
    dfs(ls(rt));
    dfs(rs(rt));
}


int main()
{
    //IN("hdu4902.txt");
    int ncase,n,q,l,v,r;
    int op;
    scanf("%d",&ncase);
    while(ncase--)
    {
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d",&num[i]);
        build(1,1,n);
        scanf("%d",&q);
        while(q--)
        {
            scanf("%d%d%d%d",&op,&l,&r,&v);
            update(1,l,r,op,v);
        }
        dfs(1);
        for(int i=1;i<=n;i++)
            printf("%d ",num[i]);
        putchar('\n');
    }
    return 0;
}



现在不会TLE了,但还是WA了: #include <bits/stdc++.h> using namespace std; #define int long long const int INF = 0x3f3f3f3f; const int MAX_N = 2e5 + 50; int T, n, m, a[MAX_N]; #define ls(cur) cur << 1 #define rs(cur) cur << 1 | 1 namespace SegmentTree1 { // the case of a[i] <= x long long sum[MAX_N << 2], vtag[MAX_N << 2]; int cnt[MAX_N << 2], stag[MAX_N << 2]; void pushup(int cur) { sum[cur] = sum[ls(cur)] + sum[rs(cur)]; cnt[cur] = cnt[ls(cur)] + cnt[rs(cur)]; } void mark(int cur, int sign, long long val) { sum[cur] *= sign; sum[cur] += cnt[cur] * val; stag[cur] *= sign; vtag[cur] += val; } void pushdown(int cur) { if (stag[cur] != 1) { mark(ls(cur), stag[cur], 0); mark(rs(cur), stag[cur], 0); stag[cur] = 1; } if (vtag[cur]) { mark(ls(cur), 1, vtag[cur]); mark(rs(cur), 1, vtag[cur]); vtag[cur] = 0; } } void build(int cur, int l, int r) { stag[cur] = 1; vtag[cur] = 0; if (l == r) { sum[cur] = cnt[cur] = 0; return ; } int mid = l + r >> 1; build(ls(cur), l, mid); build(rs(cur), mid + 1, r); pushup(cur); } void insert(int cur, int l, int r, int idx, int val) { if (l == r) { sum[cur] = val; cnt[cur] = 1; return ; } pushdown(cur); int mid = l + r >> 1; if (idx <= mid) insert(ls(cur), l, mid, idx, val); else insert(rs(cur), mid + 1, r, idx, val); pushup(cur); } void modify(int cur, int l, int r, int L, int R, int val) { if (L <= l && r <= R) { mark(cur, -1, val); return ; } pushdown(cur); int mid = l + r >> 1; if (L <= mid) modify(ls(cur), l, mid, L, R, val); if (mid + 1 <= R) modify(rs(cur), mid + 1, r, L, R, val); pushup(cur); } long long query(int cur, int l, int r, int L, int R) { if (L <= l && r <= R) return sum[cur]; pushdown(cur); int mid = l + r >> 1; long long res = 0; if (L <= mid) res += query(ls(cur), l, mid, L, R); if (mid + 1 <= R) res += query(rs(cur), mid + 1, r, L, R); return res; } }; namespace SegmentTree2 { // the case of a[i] > x long long sum[MAX_N << 2], tag[MAX_N << 2]; int cnt[MAX_N << 2], mn[MAX_N << 2]; void pushup(int cur) { sum[cur] = sum[ls(cur)] + sum[rs(cur)]; cnt[cur] = cnt[ls(cur)] + cnt[rs(cur)]; mn[cur] = min(mn[ls(cur)], mn[rs(cur)]); } void mark(int cur, long long val) { sum[cur] -= cnt[cur] * val; mn[cur] -= bool(cnt[cur]) * val; tag[cur] += val; } void pushdown(int cur) { if (tag[cur]) { mark(ls(cur), tag[cur]); mark(rs(cur), tag[cur]); tag[cur] = 0; } } void build(int cur, int l, int r, int val[]) { tag[cur] = 0; if (l == r) { sum[cur] = mn[cur] = val[l]; cnt[cur] = 1; return ; } int mid = l + r >> 1; build(ls(cur), l, mid, val); build(rs(cur), mid + 1, r, val); pushup(cur); } void modify(int cur, int l, int r, int L, int R, int val) { if (l == r && mn[cur] <= val) { SegmentTree1::insert(1, 1, n, l, mn[cur]); sum[cur] = cnt[cur] = 0, mn[cur] = INF; return ; } if (L <= l && r <= R && mn[cur] > val) { mark(cur, val); return ; } pushdown(cur); int mid = l + r >> 1; if (L <= mid) modify(ls(cur), l, mid, L, R, val); if (mid + 1 <= R) modify(rs(cur), mid + 1, r, L, R, val); pushup(cur); } long long query(int cur, int l, int r, int L, int R) { if (L <= l && r <= R) return sum[cur]; pushdown(cur); int mid = l + r >> 1; long long res = 0; if (L <= mid) res += query(ls(cur), l, mid, L, R); if (mid + 1 <= R) res += query(rs(cur), mid + 1, r, L, R); return res; } }; signed main() { ios::sync_with_stdio(false); cin.tie(0), cout.tie(0); cin >> T; while (T--) { cin >> n >> m; for (int i = 1; i <= n; i++) cin >> a[i]; SegmentTree1::build(1, 1, n); SegmentTree2::build(1, 1, n, a); for (int opt, l, r, x; m--; ) { cin >> opt >> l >> r; if (opt == 1) { cin >> x; SegmentTree2::modify(1, 1, n, l, r, x); SegmentTree1::modify(1, 1, n, l, r, x); } else { long long ans1 = SegmentTree1::query(1, 1, n, l, r); long long ans2 = SegmentTree2::query(1, 1, n, l, r); cout << ans1 + ans2 << '\n'; } } } return 0; }
最新发布
07-29
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值