Hdu 5808 Price List Strike Back(分治,BestCoder Round #86 E)

传送门:Hdu 5808 Price List Strike Back

题意:一共有n家商店，编号依次为1到n。每家商店只会卖一种物品，其中第i家商店的物品单价为v,且它到Byteasar的家的距离为di
有m组询问(1≤n≤20000,1≤m≤100000,1≤vi<=100,1<=di<=10^9),第i天他会选择一个区间[li,ri]并给自己设定一个距离上限ci,
然后他会在编号在该区间内每家到自己家的距离不超过ci的商店购买最多一件物品，当然他也可以选择什么都不买。回家之后，Byteasar会把今天购物所花的钱的总数sum记录在账本上。
判断每一条记录是否会有错

思路:摘自bc题解
考虑对序列进行分治，设当前分治区间为[l,r],取mid=(l+r)/2,那么所有在[l,mid)的询问和(mid,r]的询问可以递归分治求解，故只需考虑必然经过midmid的询问。

设f[i][j]表示考虑了[i,mid],目前选出的物品和为j时,所选商店到家的距离的最大值最小是多少;g[i][j]表示考虑了(mid,i],目前选出的物品和为j时,所选商店到家的距离的最大值最小是多少,f和g都能在O(100(r-l+1))的复杂度内求出。

那么对于一个询问l,r,s,c,只需要求出min(max(f[j],g[s-j]))然后和c比较一下大小就好了。

时间复杂度O(100(nlog n+m)+m log n)

思维拓展:还可以限制每个店的上限(单调队列维护多重背包)或者每个店容量无限(完全背包)或者依赖性背包

---

#include<bits/stdc++.h>
using namespace std;
const int INF=0x3f3f3f3f;
const int N=2e4+5;
const int M=1e5+5;
char ans[M];
int v[N],d[N];
int f[N][101],g[N][101];
struct query{
    int l,r,c,s,id; //c表示距离,s表示价值和
}Query[M];

void getFG(int L,int R,int mid){
    f[mid+1][0]=0;
    for(int i=1;i<=100;i++)
        f[mid+1][i]=INF;
    for(int i=mid;i>=L;i--)
        for(int j=100;j>=0;j--){
            f[i][j]=f[i+1][j];
            if(j>=v[i])
                f[i][j]=min(f[i][j],max(f[i+1][ j-v[i] ],d[i]));
        }
    g[mid][0]=0;
    for(int i=1;i<=100;i++)
        g[mid][i]=INF;
    for(int i=mid+1;i<=R;i++)
        for(int j=100;j>=0;j--){
            g[i][j]=g[i-1][j];
            if(j>=v[i])
                g[i][j]=min(g[i][j],max(g[i-1][ j-v[i] ],d[i]));
        }
}

void getans(query Q[],int m){
    for(int i=0;i<m;i++){
        int Min=INF;
        for(int j=0;j<=Q[i].s;j++)
            Min=min(Min,max(f[ Q[i].l ][j],g[ Q[i].r ][ Q[i].s-j ]));
        ans[ Q[i].id ]=(Min <= Q[i].c ) ? '0':'1';
    }
}

void solve(int L,int R,query Q[],int m){
    if(L==R){
        for(int i=0;i<m;i++)
            ans[Q[i].id]=( (Q[i].s == v[L] && Q[i].c>=d[L])   ? '0' : '1');
        return ;
    }
    int mid=L+R>>1,CutL=0;
    for(int i=0;i<m;i++) //全部在左边
        if(Q[i].r<=mid)
            swap(Q[i],Q[CutL++]);
    int CutR=CutL;
    for(int i=CutL;i<m;i++)
        if(Q[i].l>mid)
            swap(Q[i],Q[CutR++]);
    if(CutR<m){
        getFG(L,R,mid);
        getans(Q+CutR,m-CutR);
    }
    if(CutL>0)
        solve(L,mid,Q,CutL);
    if(CutR>CutL)
        solve(mid+1,R,Q+CutL,CutR-CutL);
}

int main(){
    int _,n,m;
    scanf("%d",&_);
    while(_--){
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)   scanf("%d",&v[i]);
        for(int i=1;i<=n;i++)   scanf("%d",&d[i]);
        for(int i=0;i<m;i++)
            scanf("%d%d%d%d",&Query[i].l,&Query[i].r,&Query[i].c,&Query[i].s),Query[i].id=i;
        solve(1,n,Query,m);
        ans[m]='\0';
        puts(ans);
    }
    return 0;
}