J. 二进制与、平方和

最新推荐文章于 2024-09-13 16:27:15 发布

zzzyyzz_

最新推荐文章于 2024-09-13 16:27:15 发布

阅读量196

点赞数

分类专栏：数据结构文章标签：算法

本文链接：https://blog.csdn.net/zzzyyzz_/article/details/131689989

版权

数据结构专栏收录该内容

13 篇文章 1 订阅

订阅专栏

文章描述了一种解决Codeforces竞赛中关于二进制与操作和平方和计算的问题的方法。通过使用线段树和bitset数据结构，优化处理每一位的更新和查询操作，以提高效率。在处理二进制与操作时，忽略x为1的位，仅处理x为0的位，同时维护区间内的非零位信息。对于平方和，通过重新计算受影响部分来更新结果。

摘要由CSDN通过智能技术生成

Problem - J - Codeforces

思路：因为与操作是按位进行的，我们可以考虑拆开每一位，并且我们发现了一个性质，如果x的某一位是1，那么x对这一位是无影响的，因为任何数与1都是它本身，所以对于x的二进制中为1的位不用管，如果某一位为零，那么我们就把l,r的所有数的该位设置为0，我们可以用一个线段树来维护，二进制的位则用一个bitset来维护，同时我们用一个bitset来维护当前区间内第i位是否存在不为零的数，这样当整个区间都为零时我们可以不用进行操作了，再就是维护平方和时，我们可以直接在修改时重新计算一边即可，因为每个数最多会被与24次，每次枚举24位，所以重新只需要24*24*n次，是可以完成的

// Problem: J. 二进制与、平方和
// Contest: Codeforces - 2020 CCPC Henan Provincial Collegiate Programming Contest
// URL: https://codeforces.com/gym/104095/problem/J
// Memory Limit: 512 MB
// Time Limit: 2000 ms

#include<iostream>
#include<cstring>
#include<string>
#include<sstream>
#include<cmath>
#include<bitset>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<map>
#include<stack>
#include<vector> 
#include<set>
#include<unordered_map>
#include<ctime>
#include<cstdlib>
#define fi first
#define se second
#define i128 __int128
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<int,int> PII;
typedef pair<int,pair<int,int> > PIII;
const double eps=1e-7;
const int N=5e5+7 ,M=5e5+7, INF=0x3f3f3f3f,mod=1e9+7,mod1=998244353;
const long long int llINF=0x3f3f3f3f3f3f3f3f;
inline ll read() {ll x=0,f=1;char c=getchar();while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
while(c>='0'&&c<='9') {x=(ll)x*10+c-'0';c=getchar();} return x*f;}
inline void write(ll x) {if(x < 0) {putchar('-'); x = -x;}if(x >= 10) write(x / 10);putchar(x % 10 + '0');}
inline void write(ll x,char ch) {write(x);putchar(ch);}
void stin() {freopen("in_put.txt","r",stdin);freopen("my_out_put.txt","w",stdout);}
bool cmp0(int a,int b) {return a>b;}
template<typename T> T gcd(T a,T b) {return b==0?a:gcd(b,a%b);}
template<typename T> T lcm(T a,T b) {return a*b/gcd(a,b);}
void hack() {printf("\n----------------------------------\n");}

int T,hackT;
int n,m,k;
int w[N];
struct Node{
	int l,r;
	int sum;
	bitset<24> s;
	bitset<24> st;
};

Node tr[N*4];

void pushup(int u) {
	tr[u].sum=(tr[u<<1].sum+tr[u<<1|1].sum)%mod1;
	tr[u].st=tr[u<<1].st|tr[u<<1|1].st;
}

void build(int u,int l,int r) {
	if(l==r) {
		bitset<24> temp;
		for(int i=23;i>=0;i--) if((w[l]>>i)&1) temp[i]=1;
		tr[u]={l,r,(ll)w[l]*w[l]%mod1,temp,temp};
	}else {
		tr[u]={l,r};
		int mid=l+r>>1;
		build(u<<1,l,mid),build(u<<1|1,mid+1,r);
		
		pushup(u);
	}
	
}

void modify(int u,int l,int r,int x) {
	if(tr[u].st[x]==0) return ;
	if(tr[u].l==tr[u].r) {
		tr[u].s[x]=0;
		tr[u].st[x]=0;
		int res=0;
		for(int i=23;i>=0;i--) if(tr[u].s[i]==1) res+=(1<<i);
		tr[u].sum=(ll)res*res%mod1;
	}else {
		int mid=tr[u].l+tr[u].r>>1;
		if(l<=mid) modify(u<<1,l,r,x);
		if(r>mid) modify(u<<1|1,l,r,x);
		
		pushup(u);
	}
}

int query(int u,int l,int r) {
	if(tr[u].l>=l&&tr[u].r<=r) return tr[u].sum;
	else {
		int mid=tr[u].l+tr[u].r>>1;
		int res=0;
		if(l<=mid) res=(res+query(u<<1,l,r))%mod1;
		if(r>mid) res=(res+query(u<<1|1,l,r))%mod1;
		
		return res;
	}
}

void solve() {
	n=read();
	
	for(int i=1;i<=n;i++) w[i]=read();
	
	build(1,1,n);
	
	int q=read();
	while(q--) {
		int op=read();
		if(op==1) {
			int l=read(),r=read();
			int x=read();
			for(int i=23;i>=0;i--) {
				if(!((x>>i)&1)) modify(1,l,r,i);
			}
		}else if(op==2) {
			int l=read(),r=read();
			printf("%d\n",query(1,l,r));
		}
	}
}   

int main() {
    // init();
    // stin();

    // scanf("%d",&T);
    T=1; 
    while(T--) hackT++,solve();
    
    return 0;       
}