HDU 5316 Magician

本文介绍了一种使用线段树解决特定问题的方法,包括如何在节点中存储关键信息并进行更新和查询操作。
题目不多说,就是更新和查询,然后用线段树解决问题,每个节点储存四个信息,分别是以奇数位开头奇数位结尾,奇数位开头偶数位结尾,偶数位开头奇数位结尾,偶数位开头偶数位结尾的最大和,然后更新就行没具体看代码
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <functional>
#include <cstdio>
#include <queue>
#include <map>
#include <algorithm>
#include <stack>
#include <utility>
typedef long long ll;
using namespace std;
const int mx = 100009;
const ll inf = 0x3f3f3f3f3f3f3f3f;
ll M[mx << 2][2][2],num[mx],k[2][2],cc[2][2];
int n,ql,qr,position,date;
void pushup(int o)
{
    M[o][1][1] = max(M[o << 1][1][1] + M[o << 1 | 1][0][1],M[o << 1][1][0] + M[o << 1 | 1][1][1]);
    M[o][1][0] = max(M[o << 1][1][0] + M[o << 1 | 1][1][0],M[o << 1][1][1] + M[o << 1 | 1][0][0]);
    M[o][0][0] = max(M[o << 1][0][0] + M[o << 1 | 1][1][0],M[o << 1][0][1] + M[o << 1 | 1][0][0]);
    M[o][0][1] = max(M[o << 1][0][1] + M[o << 1 | 1][0][1],M[o << 1][0][0] + M[o << 1 | 1][1][1]);
    for(int i = 0; i < 2; i++)
        for(int j = 0; j < 2; j++)
            M[o][i][j] = max(M[o][i][j],max(M[o << 1][i][j],M[o << 1 | 1][i][j]));
}
void add(int o)
{
    cc[1][1] = max(k[1][1],max(k[1][1] + M[o][0][1],k[1][0] + M[o][1][1]));
    cc[1][0] = max(k[1][0],max(k[1][0] + M[o][1][0],k[1][1] + M[o][0][0]));
    cc[0][0] = max(k[0][0],max(k[0][0] + M[o][1][0],k[0][1] + M[o][0][0]));
    cc[0][1] = max(k[0][1],max(k[0][1] + M[o][0][1],k[0][0] + M[o][1][1]));
    for(int i = 0; i < 2; i++)
        for(int j = 0; j < 2; j++)
            k[i][j] = max(cc[i][j],M[o][i][j]);

}
void build(int o = 1,int L = 1,int R = n)
{
    if(L == R)
    {
        for(int i = 0; i < 2; i++)
            for(int j = 0; j < 2; j++)
                M[o][i][j] = -inf;
        M[o][1&L][1&L] = num[L];
        return;
    }
    int mid = (L + R) >> 1;
    build(o << 1,L,mid);
    build(o << 1 | 1,mid + 1,R);
    pushup(o);
}

void update(int o = 1,int L = 1,int R = n)
{
    if(L == R)
    {
        num[position] = date;
        M[o][1&L][1&L] = num[L];
        return;
    }
    int mid = (L + R) >> 1;
    if(position <= mid) update(o << 1,L,mid);
    else update(o << 1 | 1,mid + 1,R);
    pushup(o);
}
void que(int o = 1,int L = 1,int R = n)
{
    if(ql <= L && qr >= R)
    {
        if(L == ql)
            for(int i = 0; i < 2; i++)
                for(int j = 0; j < 2; j++)
                    k[i][j] = M[o][i][j];
        else
            add(o);
        return;
    }
    int mid = (L + R) >> 1;
    if(ql <= mid)  que(o << 1,L,mid);
    if(qr > mid)  que(o << 1 | 1,mid + 1,R);
}
int main ()
{
    int T,i,j,ki,m,x,y;
    ll ans;
    scanf("%d",&T);
    while (T--)
    {
        scanf("%d%d",&n,&m);
        for(i = 1; i <= n; i++)
            scanf("%I64d",&num[i]);
        build();
        while(m--)
        {
            scanf("%d",&ki);
            if(ki)
            {
                scanf("%d%d",&position,&date);
                update();
            }
            else
            {
                ans = -inf;
                scanf("%d%d",&ql,&qr);
                que();
                for(i = 0; i < 2; i++)
                    for(j = 0; j < 2; j++)
                        ans = max(ans,k[i][j]);
                printf("%I64d\n",ans);
            }
        }
    }
    return 0;
}

内容概要:本文围绕“基于交流潮流的电力系统多元件N-k故障模型研究”展开,深入探讨了利用Matlab代码实现电力系统在发生多个关键元件同时故障(即N-k故障)情况下的交流潮流计算与故障分析方法。该模型不仅考虑了传统潮流方程的非线性特性,还引入了故障约束条件,能够精确模拟复杂多样的故障场景,如短路、断线等,进而评估电网在极端运行条件下的稳态与动态行为。研究通过构建典型电力系统算例,验证了所提模型在故障筛选、脆弱性识别及系统恢复策略制定方面的有效性,为电力系统安全评估、风险预警和防御体系构建提供了坚实的理论依据和技术支撑。此外,模型具备良好的扩展性,可进一步应用于连锁故障传播分析、恶意攻击模拟等高级安全分析领域。; 适合人群:具备电力系统分析基础理论知识和Matlab编程能力的高校研究生、科研院所研究人员以及电力公司从事电网规划、运行与安全管理的技术人员,特别适用于开展电力系统安全稳定、可靠性评估与应急响应机制研究的专业人士。; 使用场景及目标:①开展电力系统在多重故障条件下的交流潮流仿真,评估系统电压稳定性、线路过载风险及负荷损失程度;②识别电网中的关键薄弱环节与脆弱元件,支撑电网加固改造与防御资源配置;③用于科研项目中的故障场景建模与算法验证,或作为教学案例帮助学生理解复杂故障下的系统响应机制。; 阅读建议:此资源以Matlab代码为核心实现手段,建议读者结合理论推导与代码实现进行对照学习,重点关注故障建模过程中雅可比矩阵的修正方法、故障注入方式及收敛性处理策略,建议在仿真中逐步增加故障数量与复杂度,深入理解N-k故障对系统潮流分布的影响规律,并尝试将其拓展至含新能源接入的现代电力系统场景中进行验证与优化。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值