QOJ 6504. CCPC Final 2022 D Flower's Land 2题解

QOJ 6504. CCPC Final 2022 D Flower's Land 2题解

题意简述

给你一个只含 \(0,1,2\) 的序列,相邻两个相同的数字可以直接消掉。

询问包含两种

  • 区间所有数 \(+1\) 并对 \(3\) 取模。

  • 求一段区间能否用上述消除方式消完。

    样例输入

    8 9
    01211012
    2 4 5
    2 3 6
    1 6 8
    1 6 8
    2 3 6
    2 1 8
    1 1 1
    1 7 7
    2 1 8
    

    样例输出 #1

    Yes
    No
    Yes
    No
    Yes
    

提示

在我们做相邻两个能被消掉,判断一段区间能否被消掉时,常常用矩阵来考虑。

把每一种颜色用一种矩阵来表示,若当前位是偶数就设为这个矩阵,若当前位是奇数就设为这个矩阵的逆。

求解就把所有的矩阵乘起来,看最后结果矩阵是不是 \(I\)

为什么矩阵是正确的呢?因为矩阵满足结合律但不满足交换律。

这样就可以保证 \(1,2,3,1,2,3\) 会判断为错。

如果还没理解,下面再解释详细一点:

这是一段序列 \(0122221000\) 显然他是合法的。

在矩阵中,因为满足结合律,你先算中间那段 \(2222\) ,因为奇数和偶数个数相同,一定为 \(I\) ,相当于没有了,变成了 \(0110000\) 。一直向下就可以得到 \(I\)

题解

我们用线段树来维护矩阵乘法,这很容易,具体就是加了以后如何在矩阵中体现出来。

因为只有 \(0,1,2\) ,我们把当前,\(+1\) 后, \(+2\) 后的矩阵都记录下来。这样就可以了。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 5e5 + 10, mod = 998244353;
int n,q;
char x;
int a[N], opt, l, r;
inline ll mpow(ll x,int k){
	ll ans = 1;
	while(k){
		if(k & 1) ans = ans * x % mod;
		x = x * x % mod;
		k >>= 1;
	}
	return ans;
}
struct Mar{
	ll a[3][3];
	inline Mar operator *(const Mar b)const{
		Mar c;
		for(int i = 1; i <= 2; ++i){
			for(int j = 1; j <= 2; ++j){
				c.a[i][j] = 0;
				for(int k = 1;k <= 2; ++k){
					c.a[i][j] = (c.a[i][j] + a[i][k] * b.a[k][j] % mod) % mod;
				}	
			}
		}
		return c;
	}
	inline bool check(){
		if(a[1][1] != 1) return 0;
		if(a[1][2] != 0) return 0;
		if(a[2][1] != 0) return 0;
		if(a[2][2] != 1) return 0;
		return 1;
	}
	inline Mar inv()const{
		Mar c, b; 
		c.a[1][1] = 1;
		c.a[1][2] = 0;
		c.a[2][1] = 0;
		c.a[2][2] = 1;
		for(int i = 1; i <= 2; ++i)for(int j = 1; j <= 2; ++j) b.a[i][j] = a[i][j];
		for(int i = 1; i <= 2; ++i){
			for(int j = 1; j <= 2; ++j){
				if(i == j) continue;
				ll w = b.a[j][i] * mpow(b.a[i][i],mod - 2) % mod;
				for(int k = 1; k <= 2; ++k){
					b.a[j][k] = (b.a[j][k] - b.a[i][k] * w % mod + mod) % mod;
				}
				for(int k = 1; k <= 2; ++k){
					c.a[j][k] = (c.a[j][k] - c.a[i][k] * w % mod + mod) % mod;
				}
			}
		}
		for(int i = 1; i <= 2; ++i){
			for(int j = 1; j <= 2; ++j){
				c.a[i][j] = c.a[i][j] * mpow(b.a[i][i],mod - 2) % mod;
			}
		}
		return c;
	}
	inline void print(){
		for(int i = 1; i <= 2; ++i){
			for(int j = 1; j <= 2; ++j){
				cout<<a[i][j]<<' ';
			}
			cout<<'\n';
		}
	}
}I;
struct node{
	Mar now,nxt,nnt;
	int tag;
}tr[N << 2];
Mar m[3], m_[3];
inline void pre(){
	I.a[1][1] = 1,
	I.a[1][2] = 0;
	I.a[2][1] = 0;
	I.a[2][2] = 1;

	m[0].a[1][1] = 2,m[0].a[1][2] = 3;
	m[0].a[2][1] = 5,m[0].a[2][2] = 7;

	m[1].a[1][1] = 11,m[1].a[1][2] = 13;
	m[1].a[2][1] = 17,m[1].a[2][2] = 19;

	m[2].a[1][1] = 23,m[2].a[1][2] = 29;
	m[2].a[2][1] = 31,m[2].a[2][2] = 37;

	m_[0] = m[0].inv();
	m_[1] = m[1].inv();
	m_[2] = m[2].inv();
}
inline void input(){
	cin>> n >> q;
	for(int i = 1; i <= n; ++i){
		cin>>x;
		a[i] = x - '0';
	}
}
inline void pd(int x){
	cout<<"now:"<<'\n';
	tr[x].now.print();
	cout<<"nxt:"<<'\n';
	tr[x].nxt.print();
	cout<<"nnt:"<<'\n';
	tr[x].nnt.print();
	cout<<"tag:"<<'\n'<<tr[x].tag<<'\n';
} 
inline void downdate(int x){
	tr[x << 1].tag = (tr[x << 1].tag + tr[x].tag) % 3;
	tr[x << 1 | 1].tag = (tr[x << 1 | 1].tag + tr[x].tag) % 3;
	while(tr[x].tag > 0){
		swap(tr[x << 1].now, tr[x << 1].nxt);
		swap(tr[x << 1].nnt, tr[x << 1].nxt);
		swap(tr[x << 1 | 1].now, tr[x << 1 | 1].nxt);
		swap(tr[x << 1 | 1].nnt, tr[x << 1 | 1].nxt);
		--tr[x].tag;
	}
}
inline void pushup(int x){
	tr[x].now = tr[x << 1].now * tr[x << 1 | 1].now;
	tr[x].nxt = tr[x << 1].nxt * tr[x << 1 | 1].nxt;
	tr[x].nnt = tr[x << 1].nnt * tr[x << 1 | 1].nnt;
}
inline void build(int x, int l, int r){
	if(l == r){
		if(l % 2){
			tr[x].now = m_[a[l]];
			tr[x].nxt = m_[(a[l] + 1) % 3];
			tr[x].nnt = m_[(a[l] + 2) % 3];
		}else{ 
			tr[x].now = m[a[l]];
			tr[x].nxt = m[(a[l] + 1) % 3];
			tr[x].nnt = m[(a[l] + 2) % 3];
		}
		return ;
	}
	int mid = (l + r) >> 1;
	build(x << 1, l, mid);
	build(x << 1 | 1, mid + 1, r);
	pushup(x);
}
inline void adtr(int x){
	tr[x].tag = (tr[x].tag + 1) % 3;
	swap(tr[x].now, tr[x].nxt);
	swap(tr[x].nnt, tr[x].nxt);
}
inline void add(int x, int l, int r, int L, int R){
	if(L <= l && r <= R){
		adtr(x);
		return ;
	}
	downdate(x);
	int mid = (l + r) >> 1;
	if(L <= mid) add(x << 1, l, mid, L, R);
	if(R > mid) add(x << 1 | 1, mid + 1, r, L, R);
	pushup(x);
}
inline Mar query(int x, int l, int r, int L, int R){
	if(L <= l && r <= R){
		// cout<<x<<'\n';
		// tr[x].now.print();
		return tr[x].now;
	}
	downdate(x);
	int mid = (l + r) >> 1;
	Mar ans = I;
	if(L <= mid) ans = ans * query(x << 1, l, mid, L, R);
	if(R > mid) ans = ans * query(x << 1 | 1,mid + 1, r, L, R);
	return ans; 
}
inline void op(){
	build(1,1,n);
	for(int i = 1; i <= q; ++i){
		cin>> opt >> l >> r;
		if(opt == 1){
			add(1, 1, n, l, r);
		}else if(opt == 2){
			if(query(1, 1, n, l, r).check()){
				cout<<"Yes"<<'\n';
			}else{
				cout<<"No"<<'\n';
			}
		}
	}
}

int main(){
	cin.tie(0)->sync_with_stdio(false);
	pre();
	input();
	op();
	return 0;
}

热门相关:恭喜你被逮捕了   戏精老公今天作死没   锦庭娇   锦庭娇   万古至尊