最远曼哈顿距离小结 poj 2926 Requirements&hdu 4666 Hyperspace

48 篇文章 0 订阅
30 篇文章 0 订阅

首先对于两点(x1,y1),(x2,y2)二维的曼哈顿距离为|x1-x2|+|y1-y2|
将绝对值去掉,再将同一点的坐标归于一处可得四种情况:
1.x1-x2+y1-y2->(x1+y1)-(x2+y2)
2.x2-x1+y1-y2->(-x1+y1)-(-x2+y2)
3.x1-x2+y2-y1->(x1-y1)-(x2-y2)
4.x2-x1+y2-y1->(-x1-y1)-(-x2-y2)
在右边式子中可以发现左半部分和右半部分的符号完全相同.
我们假设0代表’+’,1代表’-‘,则有
00 01 10 11.
这里和容斥原理中求区间和给点数互质个数时的二进制枚举没有任何差别.
我们可以处理处一个点所有的axi+byi的情况,若维数为k,易知有2^k种.
而最大距离必定在一个点和另一个点的2^k种情况下差值的最大值中形成.
静态:
v[i][j]: 第i个数的第j中状态
Double get(int n,int K)
{
for(int i=0;i<n;i++)
       {
           for(int j=0;j<(1<<K);j++)
           {
                double sum=0;
                for(int p=0;p<5;p++)
                {
                    if(j&(1<<p))
                    sum+=a[p];
                    else sum-=a[p];
                }
                v[j][i]=sum;
           }
       }
       double ans=-1;
       for(int i=0;i<(1<<K);i++)
       {
           double a=-inf,b=inf;
           for(int j=0;j<n;j++)
           {
                if(v[i][j]>a) a=v[i][j];
                if(v[i][j]<b) b=v[i][j];
           }
           if(a-b>ans) ans=a-b;
       }
}
 
 
动态:
struct st
{
   int a[6],b[35];
}s[maxn];
结构体保存每个点的维数距离(a),处理出来的状态(b)
multiset<int> data[35];
data动态维护,对每一组状态保存一个域,里面存所有点
work每次处理处所有的状态,add为增加一个点,del为删除一个点
void work(int id)
{
   for(int i=0;i<(1<<k);i++)
    {
       int sum=0;
       for(int j=0;j<k;j++)
       {
           if(i&(1<<j))
           {
                sum+=s[id].a[j];
           }
           else
           {
                sum-=s[id].a[j];
           }
       }
       s[id].b[i]=sum;
    }
}
void add(int id)
{
   work(id);
   for(int i=0;i<(1<<k);i++)
    {
       data[i].insert(s[id].b[i]);
    }
}
void del(int id)
{
   for(int i=0;i<(1<<k);i++)
    {
       data[i].erase(data[i].find(s[id].b[i]));
    }
}
 
int get()
{
   for(int i=0; i<(1<<k); i++)
    {
       int a=*(data[i].rbegin());
       int b=*(data[i].begin());
       if(a-b>ans) ans=a-b;
    }
}

poj2926
#include <stdio.h>
#include <cstring>
#define inf 999999999999999
using namespace std;
double v[35][100010];
double a[6];
int n;
//const double inf=1000000000;
int main()
{
    double sum;
    while(~scanf("%d",&n))
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<5;j++) scanf("%lf",&a[j]);
            for(int j=0;j<32;j++)
            {
                sum=0.0;
                for(int k=0;k<5;k++)
                {
                    if(j&(1<<k))
                    sum+=a[k];
                    else sum-=a[k];
                }
                v[j][i]=sum;
            }
        }
        double ans=-1;
        for(int i=0;i<32;i++)
        {
            double a=-inf,b=inf;
            for(int j=0;j<n;j++)
            {
                if(v[i][j]>a) a=v[i][j];
                if(v[i][j]<b) b=v[i][j];
            }
            if(a-b>ans) ans=a-b;
        }
        printf("%.2lf\n",ans);
    }
    return 0;
}

hdu 4666
#include <stdio.h>
#include <cstring>
#include <set>
#include <algorithm>
using namespace std;
const int maxn=60010;
struct st
{
    int a[6],b[35];
}s[maxn];
multiset<int> data[35];
int n,k;
const int inf=1000000000;
void work(int id)
{
    for(int i=0;i<(1<<k);i++)
    {
        int sum=0;
        for(int j=0;j<k;j++)
        {
            if(i&(1<<j))
            {
                sum+=s[id].a[j];
            }
            else
            {
                sum-=s[id].a[j];
            }
        }
        s[id].b[i]=sum;
    }
}
void add(int id)
{
    work(id);
    for(int i=0;i<(1<<k);i++)
    {
        data[i].insert(s[id].b[i]);
    }
}
void del(int id)
{
    for(int i=0;i<(1<<k);i++)
    {
        data[i].erase(data[i].find(s[id].b[i]));
    }
}
int main()
{
   // freopen("in.txt","r",stdin);
    int op,a,b;
    while(scanf("%d%d",&n,&k)==2)
    {
        for(int i=0;i<(1<<k);i++) data[i].clear();
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&op);
            if(op==0)
            {
                for(int j=0;j<k;j++)
                    scanf("%d",&s[i].a[j]);
                add(i);
            }
            else if(op==1)
            {
                scanf("%d",&a);
                del(a);
            }
            int ans=-inf;
           // printf("sz=%d\n",data[0].size());
            if(data[0].size()<2)
            {
                printf("0\n");
                continue;
            }
            for(int i=0;i<(1<<k);i++)
            {
                int a=*(data[i].rbegin());
                int b=*(data[i].begin());
                if(a-b>ans) ans=a-b;
            }
            printf("%d\n",ans);
        }
    }
    return 0;
}


  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值