Codeforces 1093G Multidimensional Queries(线段树+状态压缩)

题目链接：Multidimensional Queries

题意：一个k（1<=k<=5）维数组n（1≤n≤2⋅10^5）每维的元素a（-10^6<=a<=10^6），进行q（1≤q≤2⋅10^5）次操作，操作有以下两种类型

1，修改数组中第x个元素

2.求出[l,r]区间内曼哈顿距离最大的两点间的曼哈顿距离

题解：我们可以发现，最终答案一定是两个元素的每一维分别有一个取正，一个取负最为最终答案，即答案中i，j的各维正负情况正好相反，并且k只有5维，于是我们可以用状态压缩表示各维的取正负情况，就可以知道每个数选取当前状态时对答案的贡献，每个状态用一颗线段树维护一下最大区间最大贡献即可，查找答案时因为两个数i，j的正负情况正好相反，枚举所有情况取最大即可

不懂细节见代码注释：

#include<iostream>
#include<stack>
#include<list>
#include<set>
#include<vector>
#include<algorithm>
#include<math.h>
#include<numeric>
#include<map>
#include<cstring>
#include<queue>
#include<iomanip>
#include<cmath>
#include<queue>
#include <bitset>
#include<unordered_map>
	#ifndef local
	#define endl '\n'
#endif */
#define mkp make_pair
using namespace std;
using std::bitset;
typedef long long ll;
typedef long double ld;
const int inf=0x3f3f3f3f;
const ll MAXN=2e6+10;
const ll N=1e5+100;
const ll mod=1e9+7;
const ll hash_p1=1610612741;
const ll hash_p2=805306457;
const ll hash_p3=402653189;
//-----------------------------------------------------------------------------------------------------------------*/
// ll head[MAXN],net[MAXN],to[MAXN],edge[MAXN]/*流量*/,cost[MAXN]//费用;
/* 
void add(ll u,ll v,ll w,ll s){
	to[++cnt]=v;net[cnt]=head[u];edge[cnt]=w;cost[cnt]=s;head[u]=cnt;
	to[++cnt]=u;net[cnt]=head[v];edge[cnt]=0;cost[cnt]=-s;head[v]=cnt;
}
struct elemt{
	int p,v;
};
struct comp{
	public:
		bool operator()(elemt v1,elemt v2){
			return v1.v<v2.v;
		}
};
-----------------------------------
求[1,MAXN]组合式和逆元 
ll mi(ll a,ll b){
	ll res=1;
	while(b){
		if(b%2){
			res=res*a%mod;
		}	
		a=a*a%mod;
	}
	return res;
}
ll fac[MAXN],inv[MAXN]
fac[0]=1;inv[0]=1;
for(int i=1;i<=MAXN;i){
	fac[i]=(fac[i-1]*i)%mod;
	inv[i]=mi(fac[i],mod-2);
}
ll C(int m,int n){//组合式C(m,n); 
	if(!n){
		return 1;
	}
	return fac[m]*(inv[n]*inv[m*-n]%mod)%mod;
}
---------------------------------
 unordered_map<int,int>mp;
//优先队列默认小顶堆 ， greater<int> --小顶堆  less<int> --大顶堆  
priority_queue<elemt,vector<elemt>,comp>q;
	set<int>::iterator it=st.begin();
*/
// vector<vector<int>>edge; 二维虚拟储存坐标 
//-----------------------------------------------------------------------------------------------------------------*/
  //map<int,bool>mp[N]; 
int dis[8*N][40];//维护当前状态下区间对应最大贡献，第二维表示各维度取正负状态
int n,k;
int a[2*N][5];//存元素
void pushin(int rt,int p){
	dis[rt][p]=max(dis[2*rt][p],dis[2*rt+1][p]);
}
void build(int l,int r,int rt,int p){
	if(l==r){
		for(int i=0;i<k;i++){
			if(p&(1<<i)){//取负时贡献
				dis[rt][p]-=a[l][i];
			}
			else{//取正时贡献
				dis[rt][p]+=a[l][i];
			}
		}
		return ;
	}
	int mid=(l+r)/2;
	build(l,mid,2*rt,p);
	build(mid+1,r,2*rt+1,p);
	pushin(rt,p);
}
int b[5];//存修改的目标位置
void update(int l,int r,int rt,int x,int p){//在当前状态下进行修改
	if(l==r){//查找到目标位置并进行修改
		dis[rt][p]=0;
		for(int i=0;i<k;i++){
			if(p&(1<<i)){//取负时贡献
				dis[rt][p]-=b[i];
			}
			else{//取正时贡献
				dis[rt][p]+=b[i];
			}
		}
		return ;
	}
	int mid=(l+r)/2;
	if(x<=mid){
		update(l,mid,2*rt,x,p);
	}
	else{
		update(mid+1,r,2*rt+1,x,p);
	}
	pushin(rt,p);
}
int query(int l,int r,int rt,int le,int ri,int p){//查找当前状态下的最大贡献
	if(le<=l&&ri>=r){
		return dis[rt][p];
	}
	int mid=(l+r)/2;
	int res=-inf;
	if(le<=mid){
		res=max(res,query(l,mid,2*rt,le,ri,p));
	}
	if(ri>mid){
		res=max(res,query(mid+1,r,2*rt+1,le,ri,p));
	}
	return res;
}
int main(){
/*cout<<setiosflags(ios::fixed)<<setprecision(8)<<ans<<endl;//输出ans（float）格式控制为8位小数（不含整数部分）*/
/*cout<<setprecision(8)<<ans<<endl;//输出ans（float）格式控制为8位小数（含整数部分）*/
	ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);//同步流
	cin>>n>>k;
	for(int i=1;i<=n;i++){
		for(int j=0;j<k;j++){
			cin>>a[i][j];
		}
	}
	for(int i=0;i<(1<<k);i++){//枚举状态，0为取正，1为取负
		build(1,n,1,i);
	}
	int q;
	cin>>q;
	while(q--){
		int lei,x,l,r;
		cin>>lei;
		if(lei==1){
			cin>>x;
			for(int i=0;i<k;i++){//读入当前要修改的点的各维坐标
				cin>>b[i];
			}
			for(int i=0;i<(1<<k);i++){//枚举所有状态并更新
				update(1,n,1,x,i);
			}
		}
		else{
			cin>>l>>r;
			int res=0;
			for(int i=0;i<(1<<k);i++){//枚举答案可能组成情况取最大
				res=max(res,query(1,n,1,l,r,i)+query(1,n,1,l,r,i^((1<<k)-1)/*相当于取反操作*/));
			}
			cout<<res<<endl;
		}
	}
	return 0;
}