Luogu P4402 [Cerc2007]robotic sort 机械排序

本文详细介绍了Splay树的区间翻转实现方法,包括关键的数据结构定义、旋转操作、下推操作等核心步骤,并提供了完整的代码示例。

题目描述

思路:

我们要知道 splay 区间翻转的本质是什么,
本质就是 每个节点编号可以随便改变, 想怎么变怎么变,变完之后会落在一个下标上,
然而splay内在的数列下标是满足splay 的性质,

当前节点的左子树的下标小于当前节点,
当前节点的右子树的下标大于当前节点.

#include<bits/stdc++.h>
using namespace std;
const int N = 1e5 + 100;
int sz[N], cnt[N], fa[N], ch[N][2], rt, n, m, lazy[N], ans[N], tot;
stack<int>S;
struct node {
    int a, id;
    bool operator < (node A)const {
        if (a != A.a) return a < A.a;
        return id < A.id;
    }
} f[N];
void pushup(int x) {
    sz[x] = sz[ch[x][0]] + sz[ch[x][1]] + 1;
}
bool get(int x) {
    return ch[fa[x]][1] == x;
}
void pushdown(int x) {
    if (lazy[x]) {
        lazy[ch[x][0]] ^= 1;
        lazy[ch[x][1]] ^= 1;
        swap(ch[x][0], ch[x][1]);
        lazy[x] = 0;
    }
}
void rotate(int x) {
    int y = fa[x], z = fa[y], chk = get(x);
    pushdown(y);
    pushdown(x);
    ch[y][chk] = ch[x][chk ^ 1];
    fa[ch[x][chk ^ 1]] = y;
    ch[x][chk ^ 1] = y;
    fa[y] = x;
    fa[x] = z;
    if (z) ch[z][y == ch[z][1]] = x;
    pushup(x);
    pushup(y);
}
void splay(int x, int goal) {
    while(fa[x] != goal) {
        int y = fa[x], z = fa[y];
        if (z != goal)
            (ch[z][1] == y) ^ (ch[y][1] == x) ? rotate(x) : rotate(y);
        rotate(x);
    }
    if (goal == 0) rt = x;
}
int find(int k) {
    int now = rt;
    while(1) {
        pushdown(now);
        if (ch[now][0] && k <= sz[ch[now][0]]) {
            now = ch[now][0];
        } else {
            if (ch[now][0]) k -= sz[ch[now][0]];
            k --;
            if (k <= 0) return now;
            now = ch[now][1];
        }
    }
}
int build(int l, int r, int root) {
    int now = (l + r) >> 1;
    fa[now] = root;
    if (l < now) ch[now][0] = build(l, now - 1, now);
    if (r > now) ch[now][1] = build(now + 1, r, now);
    pushup(now);
    return now;
}
void solve(int x) {
    int y = x;
    while(!S.empty()) S.pop();
    while(fa[y] != 0) {
        S.push(fa[y]);
        y = fa[y];
    }
    while(!S.empty()) {
        y = S.top();
        S.pop();
        pushdown(y);
    }
    splay(x, 0);
}

int main() {
    scanf("%d", &n);
    for (int i = 1; i <= n; ++i) {
        scanf("%d", &f[i].a);
        f[i].id = i;
    }
    int x, y;
    sort(f + 1, f + n + 1); //每个节点的节点编号就是 id + 1, 因为build 的时候比较巧妙.
    f[0].id = 0, f[n + 1].id = n + 1;
    build(1, n + 2, 0);
    for (int i = 1; i <= n; ++i) {
        solve(f[i].id + 1);  // 1 号节点,n + 1 号节点 是两边节点,所以要加一
        ans[i] = sz[ch[rt][0]];
        x = find(i);  // 区间是   [i+ 1, ans[i] + 1] , 要左边减一,右边加一 .
        y = find(ans[i] + 2);
        splay(x, 0); 
        splay(y, x);
        lazy[ch[ch[rt][1]][0]] ^= 1;
    }
    for (int i = 1; i <= n; ++i)
        printf("%d ", ans[i]);
    printf("\n");
    return 0;
}
/*
6
3 4 5 1 6 2

*/
#include <iostream>
#include <cstdio>
#include <algorithm>
using namespace std;
const int N=100005;
int n,fa[N],ch[N][2],val[N],siz[N],tot,pos[N],rt;
struct Node{
    int id,v;
}a[N];
bool cmp1(Node x,Node y) {return x.v==y.v?x.id<y.id:x.v<y.v;}
bool cmp2(Node x,Node y) {return x.id<y.id;}
bool rev[N];
inline void pushup(int x) {siz[x]=siz[ch[x][0]]+siz[ch[x][1]]+1;}
void pushdown(int x) {
    if(rev[x]) {
        swap(ch[x][0],ch[x][1]);
        rev[ch[x][0]]^=1,rev[ch[x][1]]^=1;
        rev[x]=0;
    }
}
int build(int f,int l,int r) {
    if(l>r) return 0;
    int now=++tot,mid=l+r>>1;
    val[now]=a[mid].v;fa[now]=f;
    pos[a[mid].v]=now;
    ch[now][0]=build(now,l,mid-1);
    ch[now][1]=build(now,mid+1,r);
    pushup(now);
    return now;
}
inline void rotate(int x) {
    int f=fa[x],ff=fa[f];
    pushdown(f);pushdown(x);
    bool tag=ch[fa[x]][1]==x;
    ch[f][tag]=ch[x][tag^1];
    fa[ch[f][tag]]=f;
    fa[f]=x;
    ch[x][tag^1]=f;
    fa[x]=ff;
    if(ff) ch[ff][f==ch[ff][1]]=x;
    pushup(f);pushup(x);
}
inline void splay(int x,int tar){
    for(int f;(f=fa[x])!=tar;rotate(x)) if(fa[f]!=tar)rotate((x==ch[fa[x]][0])==(f==ch[fa[f]][0])?f:x);
    if(!tar) rt=x;
}
int nxt(){
    pushdown(rt);
    int x=ch[rt][1];
    while(pushdown(x),ch[x][0]) x=ch[x][0];
    return x;
}
int main() {
    scanf("%d",&n);
    for(int i=1;i<=n;i++) scanf("%d",&a[i+1].v),a[i+1].id=i+1;
    a[1].v=0;a[n+2].v=n+1;
    sort(a+2,a+2+n,cmp1);
    for(int i=2;i<=n+1;i++) a[i].v=i-1;
    sort(a+2,a+2+n,cmp2);
    rt=build(0,1,n+2);
    for(int i=1;i<=n;i++) {
        int x=pos[i];splay(x,0);
        printf("%d ",siz[ch[x][0]]);        
        x=nxt();
        int y=pos[i-1];
        splay(y,0);
        splay(x,y);
        rev[ch[x][0]]^=1;
    }
}
内容概要:本文围绕“光伏-储能-数据中心”系统的容量优化配置展开研究,旨在通过构建数学模型并结合Matlab编程实现,对园区内光伏发电、储能设备与数据中心算力负荷之间的多能互补关系进行协同优化。研究综合考虑可再生能源出力的波动性、储能系统的充放电特性以及数据中心的动态用电需求与算力调度特性,建立以最小化综合成本、降低碳排放强度、提升能源自给率等为目标的多目标优化模型,并采用先进的智能优化算法(如粒子群、灰狼、鲸鱼等)进行求解,从而确定光伏装机容量与储能系统配置的最优方案。文中提供了完整的Matlab代码实现流程,涵盖数据预处理、模型构建、算法求解与结果可视化全过程,属于综合能源系统与信息基础设施深度融合的前沿交叉课题。; 适合人群:具备电力系统、能源工程、自动化或计算机等相关专业背景,熟悉Matlab编程与基本优化算法原理,从事新能源利用、绿色数据中心、综合能源系统规划与低碳调度等方向的研究生、科研人员及工程技术人员。; 使用场景及目标:① 掌握光伏-储能-数据中心耦合系统的建模思路与多能流协同机制;② 学习基于Matlab的多目标容量优化配置方法与智能算法应用技巧;③ 为工业园区、数字经济园区及绿色数据中心的能源系统规划、节能减排与可持续运行提供科学决策依据和技术解决方案。; 阅读建议:建议结合所提供的Matlab代码,深入理解目标函数设计、约束条件设定及算法实现细节,可通过调整负荷参数、改变气候数据、引入新的优化目标或约束条件等方式进行拓展性实验,以深化对系统特性的认知并提升实际科研与工程应用能力。
内容概要:本文档聚焦于“考虑源网荷储协调的主动配电网优化调度方法”的Matlab代码复现研究,属于电力系统与综合能源系统领域的高水平科研资源。文档系统阐述了融合电源侧、电网侧、负荷侧及储能系统(源-网-荷-储)协同优化的主动配电网调度模型,旨在通过先进的数学优化方法提升电力系统的经济性、可靠性与低碳化水平。研究内容涵盖多时间尺度调度架构、不确定性建模(如新能源出力波动)、鲁棒优化或分布鲁棒优化(DRO)等核心技术,并依托Matlab平台完成模型构建与仿真验证。同时,文档列举了多个前沿研究方向,如光储充一体化、主从博弈、微电网双层优化、电氢耦合系统、电动汽车集群调度等,体现出其在现代智能电网与综合能源系统优化调度中的代表性与前瞻性。; 适合人群:具备电力系统基础理论知识和Matlab编程能力,从事能源、电气工程、自动化、可再生能源等方向研究的研究生、科研人员及工程技术人员,特别适合致力于优化建模、智能算法应用与综合能源系统仿真的研究人员。; 使用场景及目标:① 复现高水平期刊论文中的主动配电网优化调度模型;② 掌握源网荷储协同优化的建模方法与Matlab实现技术;③ 借助所提供的完整代码资源开展二次开发、算法改进或新场景拓展,服务于科研创新、项目申报与学术论文撰写。; 阅读建议:建议结合文档中提及的相关文献与案例,按照目录结构循序渐进地学习,重点理解模型构建的逻辑框架与代码实现的关键细节,并利用附带的百度网盘资源获取全部代码与数据文件,通过动手实践加深对优化算法与工程应用的理解。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值