19/11/06 CSP-S模拟 T2 我的订书机之恋（随机化）（LCA）

今天考了三道题，第一道数论要用杜教筛等一大堆本蒟蒻不会的知识点，凉了；
第二道考试时题都没读懂，没从样例中理解题意，又凉了；
第三题花了1个小时打表，成功获得了10分，使我成功的避免了连续爆零的惨剧发生。
但本蒟蒻到现在为止也只会T2。

T2：

题意指求出一个$r$的合法的$l$的个数。

随机化解法：

我们可以将一个区间的左右两端点看作左右括号。
用随机化给每一对括号打一个不同的值$val$的标记。
预处理一个前缀异或和，如果左右括号都在其中，会异或成零，否则会异或成$val$。
于是一个合法的区间异或和为零，求一个前缀和，找$r_1$之前为$r_1$的个数。
什么？$r_2$?我们可以完全不用考虑，因为$d[r_1]==d[r_2]$，否则它们就没有共同的 $l$ 满足题意了。

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define pb push_back

int read()
{
	int ct = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) ct = ct*10 + (ch-'0'), ch = getchar();
	return ct * f;
}

const int N=2e6+5;
vector<int>v[N];
int n,m,l[N],r[N],cnt,sz,ans[N],bin[N];
bool isr[N];
ll d[N],b[N];

ll rnd() {return (rand()|((ll)rand()<<15))*(rand()|((ll)rand()<<15));}
//rand最大只有2^15,这四个rand之后最大能达到2^60，基本不会rand出相同的值
int main()
{
	srand(time(0));
	n=read();m=read();
	for(int i=1;i<=(n<<1);i++)
	{
		int x=read();
		if(x>i) {l[++cnt]=i;r[cnt]=x;isr[x]=true;}
	}
	for(int i=1;i<=cnt;i++) d[l[i]]=d[r[i]]=rnd();//随机化
	for(int i=1;i<=(n<<1);i++) d[i]^=d[i-1];
		for(int i=1;i<=(n<<1);i++) b[++sz]=d[i];b[++sz]=d[0];
		sort(b+1,b+sz+1);
		sz=unique(b+1,b+sz+1)-(b+1);
		for(int i=0;i<=(n<<1);i++) d[i]=lower_bound(b+1,b+sz+1,d[i])-b;//离散化
		for(int i=1;i<=m;i++)
		{
			int r1,r2;
			r1=read();r2=read();
			if(r1>r2) swap(r1,r2);
			if(!r1||!isr[r1]||!isr[r2]||d[r1]!=d[r2])
			{
				ans[i]=0;
				continue;
			}
			v[r1].pb(i);
		}
		bin[d[0]]++;
		for(int i=1;i<=(n<<1);i++)
		{
			for(int j=0;j<v[i].size();j++)
			{
				int id=v[i][j];
				ans[id]=bin[d[i]];
			}
			bin[d[i]]++;
		}
		for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}

正解：

正解用的是LCA，我还不会，这里我就引用一下sslz_fsy的博客；
博客网址：https://blog.csdn.net/sslz_fsy/article/details/102937893

发现可以对每个点求出它作为右端点的极小答案区间 然后可以作为答案的区间就是这些极小区间并起来 有一个 显然 的策略是一个 rrr 跳到它的最小区间的左端点 lll，然后再往左跳一步 如果跳到左括号已经凉了不能继续跳因为跳了一定跨出了一个区间 如果是右括号就可以用上述策略继续跳 于是有一个建图策略是 rrr 向 lll 的前一个连边，显然是一棵树的的深度就是合法的左端点个数。

#include<bits/stdc++.h>
#define cs const
using namespace std;
typedef long long ll;
int read(){
	int cnt = 0, f = 1; char ch = 0;
	while(!isdigit(ch)){ ch = getchar(); if(ch == '-') f = -1; }
	while(isdigit(ch)) cnt = cnt*10 + (ch-'0'), ch = getchar();
	return cnt * f;
}
cs int N = 2e6 + 5;
int n, m, l[N], r[N], p[N];
int sta[N], top;
int fa[N], rt[N], f[N][20], dep[N];
vector<int> v[N];
void dfs(int u, int RT){
	rt[u] = RT;
	for(int i = 1; i <= 18; i++) f[u][i] = f[f[u][i-1]][i-1];
	for(int i = 0; i < v[u].size(); i++){ 
		int t = v[u][i]; 
		dep[t] = dep[u] + 1; f[t][0] = u; dfs(t, RT); 
	}
}
int lca(int x, int y){
	if(dep[x] < dep[y]) swap(x, y);
	for(int i = 18; i >= 0; i--) if(dep[f[x][i]] >= dep[y]) x = f[x][i];
	if(x == y) return x;
	for(int i = 18; i >= 0; i--) if(f[x][i] ^ f[y][i]) x = f[x][i], y = f[y][i];
	return f[x][0];
}
int main(){
	n = read(), m = read();
	for(int i = 1; i <= (n << 1); i++) p[i] = read();
	for(int i = 1; i <= (n << 1); i++){
		l[i] = min(i, p[i]);
		r[i] = max(i, p[i]);
		while(top && l[i] <= sta[top]){
			int x = sta[top--];
			l[i] = min(l[i], l[x]);
			r[i] = max(r[i], r[x]);
		} sta[++top] = i;
	}
	memset(fa, -1, sizeof(fa));
	memset(rt, -1, sizeof(rt));
	for(int i = 1; i <= (n << 1); i++){
		if(r[i] == i){
			fa[i] = l[i] - 1;
			v[fa[i]].push_back(i); 
		}
	} 
	for(int i = 0; i <= (n << 1); i++){
		if(v[i].size() && !(~rt[i])) dep[i] = 1, dfs(i, i);
	}
	for(int i = 1; i <= m; i++){
		int r1 = read(), r2 = read();
		if(r1 > r2) swap(r1, r2);
		if(!r1){ puts("0"); continue; }
		if(!(~rt[r1]) || !(~rt[r2])){ puts("0"); continue; }
		if(rt[r1] ^ rt[r2]){ puts("0"); continue; }
		cout << dep[lca(r1, r2)] - 1 << '\n';
	} 
	return 0;
}