bzoj 3674

在这里插入图片描述
主席树加并查集

#include <cstdio>
#include <iostream>
#include <iomanip>
#include <string>
#include <cstdlib>
#include <cstring>
#include <queue>
#include <set>
#include <vector>
#include <map>
#include <algorithm>
#include <math.h>
#include <stack>

#define INF 0x3f3f3f3f
#define IMAX 2147483646
#define LINF 0x3f3f3f3f3f3f3f3f
#define ll long long
#define ull unsigned long long
#define uint unsigned int

using namespace std;

int n, m, tot = 0, root[200006];
struct T {
	int l, r, f, d;
}t[200006 << 5];

int build(int l, int r) {
	int p = ++tot;
	if (l == r) {
		t[p].f = l;
		t[p].d = 1;
		return p;
	}
	int mid = (l + r) >> 1;
	t[p].l = build(l, mid);
	t[p].r = build(mid + 1, r);
	return p;
}

int ask(int p, int l, int r, int x) {
	while (l < r) {
		int mid = (l + r) >> 1;
		if (x <= mid) {
			r = mid;
			p = t[p].l;
		}
		else {
			l = mid + 1;
			p = t[p].r;
		}
	}
	return p;
}

int add(int pre, int l, int r, int x) {
	int p = ++tot;
	if(l == r){
		t[p].f = x;
		t[p].d = t[pre].d + 1;
		return p;
	}
	t[p].l = t[pre].l;
	t[p].r = t[pre].r;
	int mid = (l + r) >> 1;
	if (x <= mid)t[p].l = add(t[pre].l, l, mid, x);
	else t[p].r = add(t[pre].r, mid + 1, r, x);
	return p;
}

int add(int pre, int l, int r, int s, int f) {
	int p = ++tot;
	if (l == r) {
		t[p].f = f;
		t[p].d = t[pre].d;
		return p;
	}
	t[p].l = t[pre].l;
	t[p].r = t[pre].r;
	int mid = (l + r) >> 1;
	if (s <= mid) t[p].l = add(t[pre].l, l, mid, s, f);
	else t[p].r = add(t[pre].r, mid + 1, r, s, f);
	return p;
}

int get(int p, int x) {
	int k = ask(p, 1, n, x);
	if (t[k].f == x) return k;
	return get(p, t[k].f);
}

int main() {
	cin >> n >> m;
	root[0] = build(1, n);
	int ans = 0, o, x, y;
	for (int i = 1; i <= m; i++) {
		root[i] = root[i - 1];
		scanf("%d", &o);
		if (o == 1) {
			scanf("%d %d", &x, &y);
			int rx = get(root[i], x ^= ans);
			int ry = get(root[i], y ^= ans);
			if (t[rx].f == t[ry].f) continue;
			if (t[rx].d > t[ry].d) swap(rx, ry);
			root[i] = add(root[i - 1], 1, n, t[rx].f, t[ry].f);
			if (t[rx].d == t[ry].d)
				root[i] = add(root[i], 1, n, t[ry].f);
		}
		else if (o == 2) {
			scanf("%d", &x);
			root[i] = root[x ^= ans];
		}
		else {
			scanf("%d %d", &x, &y);
			int rx = get(root[i], x ^= ans);
			int ry = get(root[i], y ^= ans);
			printf("%d\n", ans = rx == ry);
		}
	}
	return 0;
}
/*
5 6
1 1 2
3 1 2
2 1
3 0 3
2 1
3 1 2
*/
【电动汽车充电站有序充电调度的分散式优化】基于蒙特卡诺和拉格朗日的电动汽车优化调度(分时电价调度)(Matlab代码实现)内容概要:本文介绍了基于蒙特卡洛和拉格朗日方法的电动汽车充电站有序充电调度优化方案,重点在于采用分散式优化策略应对分时电价机制下的充电需求管理。通过构建数学模型,结合不确定性因素如用户充电行为和电网负荷波动,利用蒙特卡洛模拟生成大量场景,并运用拉格朗日松弛法对复杂问题进行分解求解,从而实现全局最优或近似最优的充电调度计划。该方法有效降低了电网峰值负荷压力,提升了充电站运营效率与经济效益,同时兼顾用户充电便利性。 适合人群:具备一定电力系统、优化算法和Matlab编程基础的高校研究生、科研人员及从事智能电网、电动汽车相关领域的工程技术人员。 使用场景及目标:①应用于电动汽车充电站的日常运营管理,优化充电负荷分布;②服务于城市智能交通系统规划,提升电网与交通系统的协同水平;③作为学术研究案例,用于验证分散式优化算法在复杂能源系统中的有效性。 阅读建议:建议读者结合Matlab代码实现部分,深入理解蒙特卡洛模拟与拉格朗日松弛法的具体实施步骤,重点关注场景生成、约束处理与迭代收敛过程,以便在实际项目中灵活应用与改进。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值