HDU 3397 Sequence operation

本文深入探讨了一种复杂的线段树实现,包括更新、查询和下推操作,旨在提高数据结构的理解和应用能力。

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

很麻烦的一道线段树。。。


#include <algorithm>
#include <cstring>
#include <cstdio>
#include <vector>
#include <queue>
#include <cmath>
using namespace std;
struct SEG
{
    int l, r, va;
    int le, re, ml, sum1;
    int lz, rz, mz, sum2;
    int rev;
} seg[100005<<3];
void update(int rt)
{
    int le, re, ml, ll, rr;
    if(seg[rt].l==seg[rt].r) return;
    le=seg[rt<<1].le, re=seg[rt<<1|1].re;
    ml=max(seg[rt<<1].ml,seg[rt<<1|1].ml);
    ml=max(ml,seg[rt<<1].re+seg[rt<<1|1].le);
    seg[rt].le=seg[rt<<1].le;
    seg[rt].re=seg[rt<<1|1].re;
    ll=seg[rt<<1].r-seg[rt<<1].l+1, rr=seg[rt<<1|1].r-seg[rt<<1|1].l+1;
    if(le==ll)
        seg[rt].le+=seg[rt<<1|1].le;
    if(re==rr)
        seg[rt].re+=seg[rt<<1].re;
    seg[rt].ml=ml;
    //处理0
    int lz, rz, mz;
    lz=seg[rt<<1].lz, rz=seg[rt<<1|1].rz;
    mz=max(seg[rt<<1].mz,seg[rt<<1|1].mz);
    mz=max(mz,seg[rt<<1].rz+seg[rt<<1|1].lz);
    seg[rt].lz=seg[rt<<1].lz;
    seg[rt].rz=seg[rt<<1|1].rz;
    ll=seg[rt<<1].r-seg[rt<<1].l+1, rr=seg[rt<<1|1].r-seg[rt<<1|1].l+1;
    if(lz==ll)
        seg[rt].lz+=seg[rt<<1|1].lz;
    if(rz==rr)
        seg[rt].rz+=seg[rt<<1].rz;
    seg[rt].mz=mz;

    seg[rt].sum1=seg[rt<<1].sum1+seg[rt<<1|1].sum1;
    seg[rt].sum2=seg[rt].r-seg[rt].l+1-seg[rt].sum1;
}
void push_down(int rt)
{
    int tl, tr, m;
    tl=seg[rt].l, tr=seg[rt].r;
    m=(tl+tr)>>1;
    if(tl==tr) return;
    if(seg[rt].va!=-1)
    {
        seg[rt<<1].va=seg[rt<<1|1].va=seg[rt].va;
        seg[rt<<1].rev=seg[rt<<1|1].rev=0;
        if(seg[rt].va==1)
        {
            seg[rt<<1].sum1=seg[rt<<1].le=seg[rt<<1].re=seg[rt<<1].ml=m-tl+1;
            seg[rt<<1|1].sum1=seg[rt<<1|1].le=seg[rt<<1|1].re=seg[rt<<1|1].ml=tr-m;
            seg[rt<<1].sum2=seg[rt<<1].lz=seg[rt<<1].rz=seg[rt<<1].mz=0;
            seg[rt<<1|1].sum2=seg[rt<<1|1].lz=seg[rt<<1|1].rz=seg[rt<<1|1].mz=0;
        }
        else if(seg[rt].va==0)
        {
            seg[rt<<1].sum1=seg[rt<<1].le=seg[rt<<1].re=seg[rt<<1].ml=0;
            seg[rt<<1|1].sum1=seg[rt<<1|1].le=seg[rt<<1|1].re=seg[rt<<1|1].ml=0;
            seg[rt<<1].sum2=seg[rt<<1].lz=seg[rt<<1].rz=seg[rt<<1].mz=m-tl+1;
            seg[rt<<1|1].sum2=seg[rt<<1|1].lz=seg[rt<<1|1].rz=seg[rt<<1|1].mz=tr-m;
        }
        seg[rt].va=-1;
    }
    if(seg[rt].rev)
    {
        seg[rt].rev=0;
        seg[rt<<1].rev^=1, seg[rt<<1|1].rev^=1;
        swap(seg[rt<<1].le,seg[rt<<1].lz);
        swap(seg[rt<<1].re,seg[rt<<1].rz);
        swap(seg[rt<<1].ml,seg[rt<<1].mz);
        swap(seg[rt<<1].sum1,seg[rt<<1].sum2);
        swap(seg[rt<<1|1].le,seg[rt<<1|1].lz);
        swap(seg[rt<<1|1].re,seg[rt<<1|1].rz);
        swap(seg[rt<<1|1].ml,seg[rt<<1|1].mz);
        swap(seg[rt<<1|1].sum1,seg[rt<<1|1].sum2);
    }
}
void build(int l,int r,int rt)
{
    seg[rt].l=l, seg[rt].r=r, seg[rt].va=-1, seg[rt].rev=0;
    if(l==r)
    {
        scanf("%d",&seg[rt].va);
        seg[rt].sum1=seg[rt].va;
        seg[rt].sum2=1-seg[rt].va;
        seg[rt].ml=seg[rt].re=seg[rt].le=seg[rt].va;
        seg[rt].mz=seg[rt].rz=seg[rt].lz=1-seg[rt].va;
        return;
    }
    int m=(l+r)>>1;
    build(l,m,rt<<1);
    build(m+1,r,rt<<1|1);
    update(rt);
}
void change(int l,int r,int va,int rt)
{
    push_down(rt);
    int tl=seg[rt].l, tr=seg[rt].r;
    int m=(tl+tr)>>1;
    if(l<=tl&&tr<=r)
    {
        if(va==1)
        {
            seg[rt].va=va;
            seg[rt].sum1=seg[rt].le=seg[rt].re=seg[rt].ml=tr-tl+1;
            seg[rt].sum2=seg[rt].lz=seg[rt].rz=seg[rt].mz=0;
        }
        else if(va==0)
        {
            seg[rt].va=va;
            seg[rt].sum1=seg[rt].le=seg[rt].re=seg[rt].ml=0;
            seg[rt].sum2=seg[rt].lz=seg[rt].rz=seg[rt].mz=tr-tl+1;
        }
        else
        {
            seg[rt].rev=1;
            swap(seg[rt].le,seg[rt].lz);
            swap(seg[rt].re,seg[rt].rz);
            swap(seg[rt].ml,seg[rt].mz);
            swap(seg[rt].sum1,seg[rt].sum2);
        }
        return;
    }
    if(l<=m)
        change(l,r,va,rt<<1);
    if(r>m)
        change(l,r,va,rt<<1|1);
    update(rt);
}
int query(int l,int r,int rt)
{
    push_down(rt);
    int tl=seg[rt].l, tr=seg[rt].r;
    if(l<=tl&&tr<=r)
        return seg[rt].ml;
    int m=(tl+tr)>>1;
    int ans;
    if(r<=m)
        ans=query(l,r,rt<<1);
    else if(l>m)
        ans=query(l,r,rt<<1|1);
    else
        ans=max(min(m-l+1,seg[rt<<1].re)+min(r-m,seg[rt<<1|1].le),max(query(l,r,rt<<1),query(l,r,rt<<1|1)));
    return ans;
}
int query2(int l,int r,int rt)
{
    push_down(rt);
    int tl=seg[rt].l, tr=seg[rt].r;
    if(l<=tl&&tr<=r)
        return seg[rt].sum1;
    int m=(tl+tr)>>1;
    int ans=0;
    if(l<=m)
        ans+=query2(l,r,rt<<1);
    if(r>m)
        ans+=query2(l,r,rt<<1|1);
    return ans;
}
int main()
{
    int i, j, n, m, op, a, b, T;
    while(~scanf("%d",&T))
    {
        while(T--)
        {
            scanf("%d%d",&n,&m);
            build(1,n,1);
            while(m--)
            {
                scanf("%d%d%d",&op,&a,&b);
                if(op==0)
                    change(a+1,b+1,0,1);
                else if(op==1)
                    change(a+1,b+1,1,1);
                else if(op==2)
                    change(a+1,b+1,2,1);
                else if(op==3)
                    printf("%d\n",query2(a+1,b+1,1));
                else
                    printf("%d\n",query(a+1,b+1,1));
            }
        }
    }
    return 0;
}


内容概要:本文探讨了在MATLAB/SimuLink环境中进行三相STATCOM(静态同步补偿器)无功补偿的技术方法及其仿真过程。首先介绍了STATCOM作为无功功率补偿装置的工作原理,即通过调节交流电压的幅值和相位来实现对无功功率的有效管理。接着详细描述了在MATLAB/SimuLink平台下构建三相STATCOM仿真模型的具体步骤,包括创建新模型、添加电源和负载、搭建主电路、加入控制模块以及完成整个电路的连接。然后阐述了如何通过对STATCOM输出电压和电流的精确调控达到无功补偿的目的,并展示了具体的仿真结果分析方法,如读取仿真数据、提取关键参数、绘制无功功率变化曲线等。最后指出,这种技术可以显著提升电力系统的稳定性与电能质量,展望了STATCOM在未来的发展潜力。 适合人群:电气工程专业学生、从事电力系统相关工作的技术人员、希望深入了解无功补偿技术的研究人员。 使用场景及目标:适用于想要掌握MATLAB/SimuLink软件操作技能的人群,特别是那些专注于电力电子领域的从业者;旨在帮助他们学会建立复杂的电力系统仿真模型,以便更好地理解STATCOM的工作机制,进而优化实际项目中的无功补偿方案。 其他说明:文中提供的实例代码可以帮助读者直观地了解如何从零开始构建一个完整的三相STATCOM仿真环境,并通过图形化的方式展示无功补偿的效果,便于进一步的学习与研究。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值