hdu 4666 Hyperspace（多维度最远曼哈顿距离）

献上博文一篇http://hi.baidu.com/byplane747/item/53ca46c159e654bc0d0a7b8d

设维度为k，维护（1<<k）个优先队列，用来保存0~(1<<k)-1 种状态（状态压缩），设状态1为“+”，状态0为“0”。

对于命令0，求出每个状态j的值，并与优先队列(s-j)的top()值相加，计算最大值。

对于命令1，标记消除的点，对每个队列pop()到存在的点。重新算一遍最大值。

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#define rep(i,a,b) for(int i=a;i<=b;i++)
#define clr(a,m) memset(a,m,sizeof(a))
using namespace std;

const int MAXN=666666;

struct Q{
    int t,id;
    Q(){}
    Q(int _t,int _id):t(_t),id(_id){}
    bool operator < (const Q &tmp) const{
        return t<tmp.t;
    }
};

priority_queue<Q>q[1<<5];
int vis[MAXN];
int a[5],b[1<<5];

int main()
{
    int n,cnt;
    while(~scanf("%d%d",&n,&cnt))
    {
        int op,s=(1<<cnt)-1;
        int ans=0;
        clr(vis,0);
        rep(i,0,s)
            while(!q[i].empty())
                q[i].pop();
        rep(i,1,n){
            scanf("%d",&op);

            if(!op){
                rep(j,0,cnt-1)
                    scanf("%d",&a[j]);
                rep(j,0,s){
                    int ad=0;
                    rep(k,0,cnt-1){
                        if(j&(1<<k))
                            ad+=a[k];
                        else
                            ad-=a[k];
                    }
                    if(!q[s-j].empty()){
                        ans=max(ans,ad+q[s-j].top().t);
                    }
                    b[j]=ad;
                }
                rep(j,0,s)
                    q[j].push(Q(b[j],i));
            }else {
                int pos;
                scanf("%d",&pos);
                vis[pos]=1;
                rep(j,0,s){
                    while(!q[j].empty()&&vis[q[j].top().id])
                        q[j].pop();
                }
                ans=0;
                rep(j,0,s)
                    if(!q[j].empty()&&!q[s-j].empty())
                        ans=max(ans,q[j].top().t+q[s-j].top().t);
            }
            printf("%d
",ans);
        }
    }
    return 0;
}