Uva 11992 Fast Matrix Operations 【线段树】【双标记pushdown】

博客围绕Uva 11992的矩阵快速操作展开,涉及线段树及双标记pushdown技术,这些信息技术可用于高效处理矩阵操作问题。
#include<bits/stdc++.h>
using namespace std;
#define lson rt << 1, l, mid
#define rson rt << 1 | 1, mid + 1, r
const int MAX = 1e5 + 7;
const int INF = 0x3f3f3f3f;
struct node{
    int sum, mn, mx, add, sett;
} a[25][MAX << 2];
void pushdown(int rt, int l, int r, int num){
    int mid = (l + r) >> 1;
    if(a[num][rt].sett >= 0){
        a[num][rt << 1].sett = a[num][rt << 1 | 1].sett = a[num][rt].sett;
        a[num][rt << 1].mn = a[num][rt << 1 | 1].mn = a[num][rt].sett;
        a[num][rt << 1].mx = a[num][rt << 1 | 1].mx = a[num][rt].sett;
        a[num][rt << 1].add = a[num][rt << 1 | 1].add = 0;
        a[num][rt << 1].sum = (mid - l + 1) * a[num][rt].sett;
        a[num][rt << 1 | 1].sum = (r - mid) * a[num][rt].sett;
        a[num][rt].sett = -1;
    }
    if(a[num][rt].add > 0){
        a[num][rt << 1].add += a[num][rt].add;
        a[num][rt << 1 | 1].add += a[num][rt].add;
        a[num][rt << 1 | 1].mn += a[num][rt].add;
        a[num][rt << 1 | 1].mx += a[num][rt].add;
        a[num][rt << 1].mx += a[num][rt].add;
        a[num][rt << 1].mn += a[num][rt].add;
        a[num][rt << 1].sum += (mid - l + 1) * a[num][rt].add;
        a[num][rt << 1 | 1].sum += (r - mid) * a[num][rt].add;
        a[num][rt].add = 0;
    }
}
void pushup(int rt, int num){
    a[num][rt].sum = a[num][rt << 1].sum +  a[num][rt << 1 | 1].sum;
    a[num][rt].mx = max(a[num][rt << 1].mx, a[num][rt << 1 | 1].mx);
    a[num][rt].mn = min(a[num][rt << 1].mn, a[num][rt << 1 | 1].mn);
}
void build(int rt, int l, int r, int num){
    a[num][rt].sett = -1, a[num][rt].sum = a[num][rt].mx = a[num][rt].mn = a[num][rt].add = 0;
    if(l == r) return;
    int mid = (l + r) >> 1;
    build(lson, num), build(rson, num);
}
void update_fuck(int rt, int l, int r, int x, int y, int v, int num){
    if(x <= l && r <= y){
        a[num][rt].sum += (r - l + 1) * v;
        a[num][rt].mn += v;
        a[num][rt].mx += v;
        a[num][rt].add += v;
        return;
    }
    pushdown(rt, l, r, num);
    int mid = (l + r) >> 1;
    if(x <= mid) update_fuck(lson, x, y, v, num);
    if(mid < y) update_fuck(rson, x, y, v, num);
    pushup(rt, num);
}
void update_shit(int rt, int l, int r, int x, int y, int v, int num){
    if(x <= l && r <= y){
        a[num][rt].sum = (r - l + 1) * v;
        a[num][rt].mn = v;
        a[num][rt].mx = v;
        a[num][rt].add = 0;
        a[num][rt].sett = v;
        return;
    }
    pushdown(rt, l, r, num);
    int mid = (l + r) >> 1;
    if(x <= mid) update_shit(lson, x, y, v, num);
    if(mid < y) update_shit(rson, x, y, v, num);
    pushup(rt, num);
}
int querySum(int rt, int l, int r, int x, int y, int num){
    if(x <= l && r <= y) return a[num][rt].sum;
    int mid = (l + r) >> 1;
    pushdown(rt, l, r, num);
    int ans = 0;
    if(x <= mid) ans += querySum(lson, x, y, num);
    if(mid < y) ans += querySum(rson, x, y, num);
    return ans;
}
int queryMin(int rt, int l, int r, int x, int y, int num){
    if(x <= l && r <= y) return a[num][rt].mn;
    pushdown(rt, l, r, num);
    int ans = INF;
    int mid = (l + r) >> 1;
    if(x <= mid) ans = min(ans, queryMin(lson, x, y, num));
    if(mid < y) ans = min(ans, queryMin(rson, x, y, num));
    return ans;
}
int queryMax(int rt, int l, int r, int x, int y, int num){
    if(x <= l && r <= y) return a[num][rt].mx;
    pushdown(rt, l, r, num);
    int ans = -INF;
    int mid = (l + r) >> 1;
    if(x <= mid) ans = max(ans, queryMax(lson, x, y, num));
    if(mid < y) ans = max(ans, queryMax(rson, x, y, num));
    return ans;
}
int main(){
    int n, m, q;
    while(scanf("%d%d%d", &n, &m, &q) != EOF){
        for(int i = 1; i <= n; i++) build(1, 1, m, i);
        for(int i = 0, x1, y1, x2, y2, hhd, op; i < q; i++){
            scanf("%d%d%d%d%d", &op, &x1, &y1, &x2, &y2);
            if(op != 3) scanf("%d", &hhd);
            if(op == 1) for(int i = x1; i <= x2; i++) update_fuck(1, 1, m, y1, y2, hhd, i);
            else if(op == 2) for(int i = x1; i <= x2; i++) update_shit(1, 1, m, y1, y2, hhd, i);
            else{
                int sum = 0, mx = -INF, mn = INF;
                for(int i = x1; i <= x2; i++){
                    sum += querySum(1, 1, m, y1, y2, i);
                    mn = min(mn, queryMin(1, 1, m, y1, y2, i));
                    mx = max(mx, queryMax(1, 1, m, y1, y2, i));
                }
                printf("%d %d %d\n", sum, mn, mx);
            }
        }
    }
    return 0;
}

 

内容概要:本文提出了一种基于非合作博弈理论的居民负荷分层调度模型,并结合双层鲸鱼优化算法(Two-level Whale Optimization Algorithm)进行高效求解,模型与算法均通过Matlab代码实现。研究针对电力系统中居民侧用电负荷的复杂调度问题,引入非合作博弈机制刻画各用户之间的利益竞争关系,实现负荷的分层优化分配;同时设计双层优化架构,上层优化资源配置,下层模拟用户自主决策行为,提升了模型的实用性与合理性。通过智能优化算法求解多层级、非凸非线性的博弈模型,有效提高了调度方案的收敛性与全局寻优能力,适用于现代智能电网中的需求侧管理与能源优化场景。; 适合人群:具备电力系统基础理论知识和Matlab编程能力,从事智能电网、能源优化调度、需求侧管理、博弈论应用等方向的科研人员、高校研究生及工程技术人员。; 使用场景及目标:①应用于居民区电力负荷的分层优化调度系统设计与仿真分析;②为非合作博弈在多主体能源系统建模中的应用提供方法论支持;③利用双层鲸鱼算法解决具有嵌套结构的复杂双层优化问题,提升求解效率与调度方案的可行性。; 阅读建议:建议读者结合提供的Matlab代码深入理解模型构建逻辑与算法实现流程,重点关注博弈模型的效用函数设计、纳什均衡求解思路以及双层优化结构的迭代机制,宜配合实际用电数据开展复现实验以验证模型有效性与鲁棒性。
内容概要:本文围绕基于自适应神经模糊推理系统(ANFIS)智能控制器的可再生能源微电网功率管理系统展开研究,结合Simulink仿真实现,深入探讨了微电网中功率的智能调控与经济机组组合调度问题。通过引入ANFIS控制器,有效应对风能、光伏等可再生能源出力的波动性与不确定性,提升系统运行的稳定性与电能质量。研究内容涵盖微电网多源协调控制策略、功率平衡管理、优化调度模型构建及仿真验证,实现了对分布式电源、储能系统和负荷的协同优化,兼顾经济性与可靠性目标,并通过仿真平台验证了所提方法的有效性与优越性。; 适合人群:具备电力系统、自动化或新能源相关专业背景,熟悉Matlab/Simulink仿真环境,从事微电网能量管理、智能控制、能源优化等领域研究的研究生、科研人员及工程技术人员。; 使用场景及目标:①用于高比例可再生能源接入场景下的微电网能量管理系统研发与教学实践;②为实现微电网功率稳定控制与经济高效运行提供先进的智能控制解决方案;③支撑高水平学术论文复现、科研课题攻关及实际工程项目的仿真验证与方案优化。; 阅读建议:建议结合提供的Simulink模型与相关代码进行动手实践,重点关注ANFIS控制器的设计流程、规则库构建与参数调优方法,并通过与传统PID或MPC控制策略的对比实验,深入理解其在动态响应与鲁棒性方面的优势。同时可进一步拓展文中提出的优化调度逻辑,应用于多目标、多约束的复杂实际应用场景中。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值