(Relax 线段树1.3)POJ 3667 Hotel(互不相交的线段的更新与统计)

本文介绍了一种使用线段树进行区间更新和查询的高效算法,并利用懒惰标记减少不必要的操作。该算法适用于需要频繁进行区间修改和查询的问题。


#include <iostream>
#include <cstdio>

#define maxn 80010
#define Fup(i, s, t) for (int i = s; i <= t; i ++)
#define Fdn(i, s, t) for (int i = s; i >= t; i --)

using namespace std;

struct node {
	int ls;//左端空区间的长度
	int rs;//右端空区间的长度
	int  ms;//最长子区间的长度
	int pos;//开始位置
	int mark;//懒惰标记(0:未定 , 1:全空  2:全满)
} tree[4 * maxn];//线段树
int n, m;

void build_tree(int l, int r, int i) {//构建全空的线段树
	tree[i].ls = tree[i].rs = tree[i].ms = r - l + 1;
	tree[i].pos = l;
	if (l == r)
		return;
	int mid = (l + r) / 2;
	build_tree(l, mid, i + i);
	build_tree(mid + 1, r, i + i + 1);
}

bool all_space(int l, int r, int i) {//判断根为i的区间[l,r]是否全空...
	if (tree[i].ls == r - l + 1)
		return 1;
	return 0;
}

void update(int l, int r, int i) {//通过标记法维护线段树..
	if (!tree[i].mark)
		return;
	if (tree[i].mark == 1) {
		int len = r - l + 1;
		tree[i + i].ls = tree[i + i].rs = tree[i + i].ms = (len + 1) / 2;
		tree[i + i].pos = l;
		tree[i + i + 1].ls = tree[i + i + 1].rs = tree[i + i + 1].ms = len / 2;
		tree[i + i + 1].pos = (l + r) / 2 + 1;
		tree[i + i].mark = tree[i + i + 1].mark = 1;
	} else {
		tree[i + i].ls = tree[i + i].rs = tree[i + i].ms = 0;
		tree[i + i].pos = l;
		tree[i + i + 1].ls = tree[i + i + 1].rs = tree[i + i + 1].ms = 0;
		tree[i + i + 1].pos = (l + r) / 2 + 1;
		tree[i + i].mark = tree[i + i + 1].mark = 2;
	}
	tree[i].mark = 0;
}

/**
 * 查询根为i的区间[l,r]是否存在长度为d的区间..
 * 如果存在则返回返回其左指针..
 * 否则返回0
 */
int query(int d, int l, int r, int i) {
	update(l, r, i);
	if (tree[i].ms < d)
		return 0;
	if (tree[i].ms == d)
		return tree[i].pos;
	int mid = (l + r) / 2;
	if (tree[i + i].ms >= d)
		return query(d, l, mid, i + i);
	if (tree[i + i].rs + tree[i + i + 1].ls >= d)
		return mid - tree[i + i].rs + 1;
	return query(d, mid + 1, r, i + i + 1);
}

/**
 * 在根为i的区间[l,r]上插入或删除子区间[tl,tr],插删标记为flag
 */
void change(int tl, int tr, int l, int r, int i, bool flag) {
	if (tl > r || tr < l)
		return;
	if (tl <= l && r <= tr) {
		if (flag) {
			tree[i].ls = tree[i].rs = tree[i].ms = 0;
			tree[i].pos = l;
			tree[i].mark = 2;
		} else {
			tree[i].ls = tree[i].rs = tree[i].ms = r - l + 1;
			tree[i].pos = l;
			tree[i].mark = 1;
		}
		return;
	}
	update(l, r, i);
	int mid = (l + r) / 2;
	change(tl, tr, l, mid, i + i, flag);
	change(tl, tr, mid + 1, r, i + i + 1, flag);
	tree[i].ls = tree[i + i].ls;
	if (all_space(l, mid, i + i))
		tree[i].ls += tree[i + i + 1].ls;
	tree[i].rs = tree[i + i + 1].rs;
	if (all_space(mid + 1, r, i + i + 1))
		tree[i].rs += tree[i + i].rs;
	tree[i].ms = max(tree[i + i].rs + tree[i + i + 1].ls,
			max(tree[i + i].ms, tree[i + i + 1].ms));
	if (tree[i].ms == tree[i + i].ms)
		tree[i].pos = tree[i + i].pos;
	else if (tree[i].ms == tree[i + i].rs + tree[i + i + 1].ls)
		tree[i].pos = mid - tree[i + i].rs + 1;
	else
		tree[i].pos = tree[i + i + 1].pos;
}

int main() {

	while (scanf("%d%d", &n, &m) != EOF) {
		memset(tree, 0, sizeof(tree));
		build_tree(1, n, 1);

		int i;
		for (i = 1; i <= m; ++i) {
			int kind;
			scanf("%d", &kind);

			if (kind == 1) {
				int d;
				scanf("%d", &d);

				int ans = query(d, 1, n, 1);

				printf("%d\n", ans);
				if (ans) {
					change(ans, ans + d - 1, 1, n, 1, 1);
				}
			} else {
				int x, d;
				scanf("%d%d", &x, &d);

				change(x, x + d - 1, 1, n, 1, 0);
			}
		}
	}

	return 0;
}


内容概要:本文档介绍了基于3D FDTD(时域有限差分)方法在MATLAB平台上对微带线馈电的矩形天线进行仿真分析的技术方案,重点在于模拟超MATLAB基于3D FDTD的微带线馈矩形天线分析[用于模拟超宽带脉冲通过线馈矩形天线的传播,以计算微带结构的回波损耗参数]宽带脉冲信号通过天线结构的传播过程,并计算微带结构的回波损耗参数(S11),以评估天线的匹配性能和辐射特性。该方法通过建立三维电磁场模型,精确求解麦克斯韦方程组,适用于高频电磁仿真,能够有效分析天线在宽频带内的响应特性。文档还提及该资源属于一个涵盖多个科研方向的综合性MATLAB仿真资源包,涉及通信、信号处理、电力系统、机器学习等多个领域。; 适合人群:具备电磁场微波技术基础知识,熟悉MATLAB编程及数值仿真的高校研究生、科研人员及通信工程领域技术人员。; 使用场景及目标:① 掌握3D FDTD方法在天线仿真中的具体实现流程;② 分析微带天线的回波损耗特性,优化天线设计参数以提升宽带匹配性能;③ 学习复杂电磁问题的数值建模仿真技巧,拓展在射频无线通信领域的研究能力。; 阅读建议:建议读者结合电磁理论基础,仔细理解FDTD算法的离散化过程和边界条件设置,运行并调试提供的MATLAB代码,通过调整天线几何尺寸和材料参数观察回波损耗曲线的变化,从而深入掌握仿真原理工程应用方法。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

帅气的东哥

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值