[曼哈顿距离,线段树]Educational Codeforces Round 56 G Multidimensional Queries

本文介绍了一种使用线段树解决多维空间中曼哈顿距离最远点对查询的方法,通过状态压缩和维护多个线段树来处理动态更新和查询操作。文章详细解释了算法原理,并提供了完整的代码实现。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目:[曼哈顿距离,线段树]Educational Codeforces Round 56G. Multidimensional Queries

原题链接:

http://codeforces.com/contest/1093/problem/G

题意:

给n个k维的点,q个询问,求曼哈顿距离最远的两个点的距离,操作是改点的坐标。

题解:

和poj2926是一样的,压缩下状态就好。

poj2926看这位大神的博客https://www.cnblogs.com/tanhehe/archive/2013/05/25/3099400.html

然后在线操作的话就维护下(1<<k)个线段树就好了。

总结:

查询的时候分别查一次最大值和最小值,不要直接查询最大值和最小值的差,因为线段树维护的不是这个差,而是极值。

代码:
#define push_back pb
#define make_pair mk
#define rd read()
#define mem(a,b) memset(a,b,sizeof(a))
#define bug printf("*********\n");
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define ls tr[u].ch[0]
#define rs tr[u].ch[1]
#define FIN freopen(D://code//in.txt", "r", stdin);
#define debug(x) cout<<"["<<x<<"]" <<endl;
#define IO ios::sync_with_stdio(flse),cin.tie(0);
//#pragma comment(linker,"/STACk:1024000000,1024000000")
//#include<bits/stdc++.h>
#include<time.h>
#include<iostream>
#include<stdlib.h>
#include<stdio.h>
#include<cmath>
#include<map>
#include<algorithm>
#include<string>
#include<string.h>
#include<set>
#include<queue>
#include<stack>
#include<functional>
using std::pair;


typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
//const double PI=acos(-1);


const int maxn = 2e5 + 10;
const int maxm = 1e6 + 10;
const int mod = 998244353 ;
const int inf = 0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3f;
const double dinf=1e20;
const double eps=1e-8;
using namespace std;

ll read() {
    ll X = 0, p = 1; char c = getchar();
    for(; c > '9' || c < '0'; c = getchar()) if(c == '-') p = -1;
    for(; c >= '0' && c <= '9'; c = getchar()) X = X * 10 + c - '0';
    return X * p;
}
//*********************************************************************

int n,dem;
int a[maxn][6];
int ma[maxn<<2][33],mi[maxn<<2][33];
void pushup(int rt)
{
    for(int i=0;i<(1<<dem);i++)
    {
        ma[rt][i]=max(ma[rt<<1][i],ma[rt<<1|1][i]);
        mi[rt][i]=min(mi[rt<<1][i],mi[rt<<1|1][i]);
    }
}
void build(int l,int r,int rt)
{
    if(l==r)
    {
        int sum=0;
        for(int i=0;i<(1<<dem);i++)
        {
            ma[rt][i]=-inf;
            mi[rt][i]=inf;
        }
        for(int i=0;i<(1<<dem);i++)
        {
            int sum=0;
            for(int k=0;k<dem;k++)
            {
                if((i>>k)&1) sum+=a[l][k];
                else sum-=a[l][k];
            }
            ma[rt][i]=max(ma[rt][i],sum);
            mi[rt][i]=min(mi[rt][i],sum);
        }
        return ;
    }
    int mid=(l+r)>>1;
    build(lson);
    build(rson);
    pushup(rt);
}
void update(int p,int l,int r,int rt)
{
    if(l==p&&r==p)
    {
        int sum=0;
        for(int i=0;i<(1<<dem);i++)
        {
            ma[rt][i]=-inf;
            mi[rt][i]=inf;
        }
        for(int i=0;i<(1<<dem);i++)
        {
            int sum=0;
            for(int k=0;k<dem;k++)
            {
                if((i>>k)&1) sum+=a[l][k];
                else sum-=a[l][k];
            }
            ma[rt][i]=max(ma[rt][i],sum);
            mi[rt][i]=min(mi[rt][i],sum);
        }
        return ;
    }
    int mid=(l+r)>>1;
    if(p<=mid) update(p,lson);
    if(mid<p) update(p,rson);
    pushup(rt);
}

int query(int L,int R,int l,int r,int rt,int i,int f)
{
    if(L<=l&&r<=R)
    {
        if(f)return ma[rt][i];
        else return mi[rt][i];
    }
    int mid=(l+r)>>1;
    int ret;
    if(f)
    {
        ret=-inf;
        if(L<=mid) ret=max(ret,query(L,R,lson,i,f));
        if(R>mid) ret=max(ret,query(L,R,rson,i,f));
    }
    else
    {
        ret=inf;
        if(L<=mid) ret=min(ret,query(L,R,lson,i,f));
        if(R>mid) ret=min(ret,query(L,R,rson,i,f));
    }
    return ret;
}
/*void dfs(int l,int r,int rt,int i)
{
    if(l>r) return ;
    if(l==r)
    {
        cout<<l<<" "<<r<<" "<<ma[rt][i]<<" "<<mi[rt][i]<<endl;
        return;
    }
    int mid=(l+r)>>1;
    dfs(lson,i);
    cout<<l<<" "<<r<<" "<<ma[rt][i]<<" "<<mi[rt][i]<<endl;
    dfs(rson,i);
}*/
int main()
{
    cin>>n>>dem;
    for(int i=1;i<=n;i++)
    {
        for(int j=0;j<dem;j++)
            scanf("%d",&a[i][j]);
    }
    build(1,n,1);
    int q;
    cin>>q;
    int op,l,r;
    while(q--)
    {
        scanf("%d%d",&op,&l);
        if(op==1)
        {
            for(int i=0;i<dem;i++) scanf("%d",&a[l][i]);
            update(l,1,n,1);
        }
        else
        {
            scanf("%d",&r);
            int ans=-inf;
            for(int i=0;i<(1<<dem);i++)
            {
                ans=max(ans,query(l,r,1,n,1,i,1)-query(l,r,1,n,1,i,0));
            }
            printf("%d\n",ans);
        }
    }
}
### Codeforces Educational Round 26 比赛详情 Codeforces是一个面向全球程序员的比赛平台,其中Educational Rounds旨在帮助参与者提高算法技能并学习新技巧。对于具体的Educational Round 26而言,这类比赛通常具有如下特点: - **时间限制**:每道题目的解答需在规定时间内完成,一般为1秒。 - **内存限制**:程序运行所占用的最大内存量被限定,通常是256兆字节。 - 输入输出方式标准化,即通过标准输入读取数据并通过标准输出打印结果。 然而,关于Educational Round 26的具体题目细节并未直接提及于提供的参考资料中。为了提供更精确的信息,下面基于以往的教育轮次给出一些常见的题目类型及其解决方案思路[^1]。 ### 题目示例与解析 虽然无法确切描述Educational Round 26中的具体问题,但可以根据过往的经验推测可能涉及的问题类别以及解决这些问题的一般方法论。 #### 类型一:贪心策略的应用 考虑一个问题场景,在该场景下需要照亮一系列连续排列的对象。假设存在若干光源能够覆盖一定范围内的对象,则可以通过遍历整个序列,并利用贪心的思想决定何时放置新的光源以确保所有目标都被有效照射到。这种情况下,重要的是保持追踪当前最远可到达位置,并据此做出决策。 ```cpp #include <bits/stdc++.h> using namespace std; bool solve(vector<int>& a) { int maxReach = 0; for (size_t i = 0; i < a.size(); ++i) { if (maxReach < i && !a[i]) return false; if (a[i]) maxReach = max(maxReach, static_cast<int>(i) + a[i]); } return true; } ``` #### 类型二:栈结构处理匹配关系 另一个常见问题是涉及到成对出现元素之间的关联性判断,比如括号表达式的合法性验证。这里可以采用`<int>`类型的栈来记录左括号的位置索引;每当遇到右括号时就弹出最近一次压入栈底的那个数值作为配对依据,进而计算两者间的跨度长度累加至总数之中[^2]。 ```cpp #include <stack> long long calculateParens(const string& s) { stack<long long> positions; long long num = 0; for(long long i = 0 ; i<s.length() ;++i){ char c=s[i]; if(c=='('){ positions.push(i); }else{ if(!positions.empty()){ auto pos=positions.top(); positions.pop(); num+=i-pos; } } } return num; } ``` #### 类型三:特定模式下的枚举法 针对某些特殊条件约束下的计数类问题,如寻找符合条件的三位整数的数量。此时可通过列举所有可能性的方式逐一检验是否符合给定规则,从而统计满足要求的结果数目。例如求解形如\(abc\)形式且不含重复数字的正整数集合大小[^3]。 ```cpp vector<int> generateSpecialNumbers(int n) { vector<int> result; for (int i = 1; i <= min(n / 100, 9); ++i) for (int j = 0; j <= min((n - 100 * i) / 10, 9); ++j) for (int k = 0; k <= min(n % 10, 9); ++k) if ((100*i + 10*j + k)<=n&&!(i==0||j==0)) result.emplace_back(100*i+10*j+k); sort(begin(result), end(result)); return result; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值