【bzoj1269】【AHOI2006】文本编辑器editor【Splay】

本文讨论了使用平衡树解决字符串插入问题时遇到的挑战,并提供了有效的解决方案。通过详细解释平衡树的操作,如翻转标记、节点旋转等,作者成功地解决了字符串插入过程中的空字符跳过问题。此外,文章还介绍了如何使用特定函数实现字符串的建树和树根挂接,以及如何通过调整代码结构避免低级错误。通过实例代码演示,读者可以更好地理解并应用这些技巧。

被文艺平衡树折磨了一天以后发现这道题就很好做啦~~~

传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1269
但是c++喜闻乐见的gets十分不好使。

如果insert的字符串有空格,gets会跳过去不读空格。

好在题目给了ASCII码的范围。

这样就可做了。
写的时候不要把功能一次性全写完,先写个基本的,排除一下低级错误……我把Move和insert写完以后发现Splay又出了几个沙茶错误- -

还有居然把Next和Prev搞反了- -……

Splay中cur代表当前光标位置。

flip是标记是否反转的标记~

p表示当前节点是他爸的左孩子还是右孩子。

build是用来给字符串建树,建完之后把这棵树的树根直接挂到上面去~(预先把光标两侧的节点Splay到一块去)。

代码在下面。

#include<cstring>
#include<iostream>
#include<cstdio>
using namespace std;
char str[3000002],cmd[10];
int n,x;
struct Splay{
	static const int maxn=3000001;
	struct node{
		int ch[2],f,s;
		char v;
		bool p,flip;
	}t[maxn];
	int rt,tot,cur;
	inline void upd(int x){
		t[x].s=t[t[x].ch[0]].s+t[t[x].ch[1]].s+1;
	}
	inline void pushdown(int x){
		if(!t[x].flip) return;
		t[t[x].ch[0]].flip^=1;
		t[t[x].ch[0]].p^=1;
		t[t[x].ch[1]].flip^=1;
		t[t[x].ch[1]].p^=1;
		swap(t[x].ch[0],t[x].ch[1]);
		t[x].flip=0;
	}
	inline void sc(int a,int b,bool c){//把节点a的c孩子置为b(c=0:左孩子,c=1:右孩子)
		t[a].ch[c]=b;
		t[b].f=a;
		t[b].p=c;
	}
	inline void rot(int x){
		int y=t[x].f,p=t[x].p;
		if(y==rt) t[rt=x].f=0;
		else sc(t[y].f,x,t[y].p);
		sc(y,t[x].ch[!p],p);
		sc(x,y,!p);
		upd(y);
	}
	void splay(int x,int r){
		int f;
		while((f=t[x].f)!=r){
			if(t[f].f==r) rot(x);
			else if(t[f].p==t[x].p) rot(f),rot(x);
			else rot(x),rot(x);
//			cout<<x<<"\n";
		}
		upd(x);
	}
	int build(int l,int r,int f,int p){
		if(l>r) return 0;
		int now=++tot;
		t[now].f=f;
		t[now].p=p;
		t[now].flip=false;
		t[now].v=str[l+r>>1];
		t[now].ch[0]=build(l,(l+r>>1)-1,now,0);
		t[now].ch[1]=build((l+r>>1)+1,r,now,1);
		upd(now);
		return now;
	}
	int select(int x,int k){
		while(1){
			pushdown(x);
			if(k==t[t[x].ch[0]].s+1)return x;
			else if(k<t[t[x].ch[0]].s+1) x=t[x].ch[0];
			else k-=t[t[x].ch[0]].s+1,x=t[x].ch[1];
		}
	}
	inline int insert(int subtree){
		splay(select(rt,cur+1),0);
		splay(select(rt,cur+2),rt);
		sc(t[rt].ch[1],subtree,0);
		upd(t[rt].ch[1]);
		upd(rt);
	}
	inline int Flip(int length){
		splay(select(rt,cur+1),0);
		splay(select(rt,cur+2+length),rt);
		t[t[t[rt].ch[1]].ch[0]].flip^=1;
	}
	inline void del(int length){
		splay(select(rt,cur+1),0);
		splay(select(rt,cur+2+length),rt);
		sc(t[rt].ch[1],0,0);
		upd(t[rt].ch[1]);
		upd(rt);
	}
	inline char get(){
		return t[select(rt,cur+2)].v;
	}
	Splay():rt(0),tot(0),cur(0){}
}editor;
int main(){
	scanf("%d\n",&n);
	editor.rt=editor.build(0,1,0,0);
	while(n--){
		scanf("%s",cmd);
		switch(cmd[0]){
			case 'M':
				scanf("%d",&x);
				editor.cur=x;
				break;
			case 'I':
				scanf("%d",&x);
				do{str[0]=getchar();}while(str[0]<32||str[0]>126);
				for(int i=1;i<x;++i) str[i]=getchar();
				getchar();
				editor.insert(editor.build(0,x-1,0,0));
				break;
			case 'D':
				scanf("%d",&x);
				editor.del(x);
				break;
			case 'R':
				scanf("%d",&x);
				editor.Flip(x);
				break;
			case 'G':
				putchar(editor.get());
				putchar('\n');
				break;
			case 'P':
				editor.cur--;
				break;
			case 'N':
				editor.cur++;
				break;
		}
	}
	return 0;
}

### 题目内容 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在被困在了这个n维球体中,仅知道球面上n + 1个点的坐标,需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器[^1][^2]。 ### 解题思路 - **原理依据**:根据球心的定义,球心到球面上任意一点距离都相等。设球心坐标为\((x_1,x_2,\cdots,x_n)\),球面上一点坐标为\((a_1,a_2,\cdots,a_n)\),两点距离公式为\(dist = \sqrt{(x_1 - a_1)^2+(x_2 - a_2)^2+\cdots+(x_n - a_n)^2}\)。 - **构建方程**:设球面上\(n + 1\)个点的坐标分别为\((a_{i1},a_{i2},\cdots,a_{in})\),\(i = 1,2,\cdots,n + 1\)。以第一个点和其他点为例,根据球心到各点距离相等,可得: \((x_1 - a_{11})^2+(x_2 - a_{12})^2+\cdots+(x_n - a_{1n})^2=(x_1 - a_{21})^2+(x_2 - a_{22})^2+\cdots+(x_n - a_{2n})^2\) 展开式子: \(x_1^2 - 2a_{11}x_1+a_{11}^2+x_2^2 - 2a_{12}x_2+a_{12}^2+\cdots+x_n^2 - 2a_{1n}x_n+a_{1n}^2=x_1^2 - 2a_{21}x_1+a_{21}^2+x_2^2 - 2a_{22}x_2+a_{22}^2+\cdots+x_n^2 - 2a_{2n}x_n+a_{2n}^2\) 消去\(x_1^2,x_2^2,\cdots,x_n^2\)后可得: \(2(a_{21}-a_{11})x_1 + 2(a_{22}-a_{12})x_2+\cdots+2(a_{2n}-a_{1n})x_n=a_{21}^2 - a_{11}^2+a_{22}^2 - a_{12}^2+\cdots+a_{2n}^2 - a_{1n}^2\) 同理,用第一个点和第\(i\)个点\((i = 3,\cdots,n + 1)\)可得到\(n\)个线性方程,构成一个\(n\)元一次方程组。 - **求解方程组**:使用高斯消元法求解这个\(n\)元一次方程组,得到的解就是球心的\(n\)维坐标。 ### 代码实现 ```python n = int(input()) points = [] for _ in range(n + 1): points.append(list(map(float, input().split()))) # 构建方程组的系数矩阵和常数项 a = [[0] * (n + 1) for _ in range(n)] b = [0] * n for i in range(n): for j in range(n): a[i][j] = 2 * (points[i + 1][j] - points[0][j]) b[i] += points[i + 1][j] ** 2 - points[0][j] ** 2 # 高斯消元 for i in range(n): # 选主元 max_row = i for j in range(i + 1, n): if abs(a[j][i]) > abs(a[max_row][i]): max_row = j a[i], a[max_row] = a[max_row], a[i] b[i], b[max_row] = b[max_row], b[i] # 消元 for j in range(i + 1, n): factor = a[j][i] / a[i][i] for k in range(i, n): a[j][k] -= factor * a[i][k] b[j] -= factor * b[i] # 回代求解 x = [0] * n for i in range(n - 1, -1, -1): s = 0 for j in range(i + 1, n): s += a[i][j] * x[j] x[i] = (b[i] - s) / a[i][i] # 输出结果 print(" ".join(map(lambda num: "{:.3f}".format(num), x))) ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值