BZOJ2794/POI2012 Cloakroom

Task
有n件物品，每件物品有三个属性a[i], b[i], c[i] .(a[i]＜b[i])
再给出q个询问，每个询问由非负整数m, k, s组成，问是否能够选出某些物品使得：
1. 对于每个选的物品i，满足a[i]<=m且b[i]>m+s。
2. 所有选出物品的c[i]的和正好是k。
n<=1,000, q<=1,000,000.
c[i]<=1,000, 1<=a[i]＜b[i]<=10^9
1<=m<=10^9, 1<=k<=100,000, 0<=s<=10^9

Solution
　离线.
　把询问和节点都根据左端点排序,保证询问[l,r]时,所有考虑的物品的左端点都满足条件.满足左端点的物品个数是不断增加的,一旦满足,接下来所有询问都会满足.那么现在只用考虑右端点了.
　直接确定所有物品再进行01背包肯定不行.那么考虑新思路.假如因为现在可选物品的范围一直增大,我们可以对当前物品的信息进行记录.由于k并不大,可以考虑把k存入下标:
　设dp[x]为当前价值构成x最远的b[i].
　对于询问可以直接判断dp[k]与m+s的关系,
　对于更新,用01背包进行更新即可.

题外话:
这道题的思路来源于BestCoder round#86 1005(hdu5808,题解戳这,题目链接戳这)的思路,都是询问一些符合范围的物品价值是否能够达到k.而且思路也是相近的,只不过BC的这道题信息多一维,用了分治来优化,思路更复杂一些.

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#define ll long long
#include<queue>
#include<set>
using namespace std;
inline void rd(int &res){
    res=0;char c;
    while(c=getchar(),c<48);
    do res=(res<<1)+(res<<3)+(c^48);
    while(c=getchar(),c>=48);
}
inline void print(ll x){
    if(!x)return ;
    print(x/10);
    putchar((x%10)^48);
}
inline void sc(ll x){
    if(x<0){x=-x;putchar('-');}
    print(x);
    if(!x)putchar('0');
    putchar('\n');
}
inline void Max(int &x,int y){if(x<y)x=y;}
inline void Min(int &x,int y){if(x>y)x=y;}
const int M=1005;
const int N=1e6+5;
const int P=100000;
int dp[P+5],n,m,ans[N];
char s[2][3];
struct node{
    int l,r,v,id;
    node(){id=-1;}
    bool operator<(const node &tmp)const{
        if(l!=tmp.l)return l<tmp.l;
        return id<tmp.id;
    }
}A[M+N];
void pt(int c){
    for(int i=0;i<3;i++)putchar(s[c][i]);
    putchar('\n');
}
int main(){
    int n,i,j,mx=0,sum=0,cas,a,b,k;
    memset(dp,-1,sizeof(dp));
    rd(n);
    s[0][0]='N',s[0][1]='I',s[0][2]='E';
    s[1][0]='T',s[1][1]='A',s[1][2]='K';
    for(i=1;i<=n;i++){
        rd(A[i].v);rd(A[i].l);rd(A[i].r);
        Max(mx,A[i].r);
    }
    rd(m);
    for(i=1;i<=m;i++){
        rd(A[i+n].l);rd(A[i+n].v);rd(A[i+n].r);
        A[i+n].r+=A[i+n].l;
        A[i+n].id=i;
    }
    sort(A+1,A+1+n+m);
    dp[0]=mx+1;
    for(i=1;i<=n+m;i++){
        int l=A[i].l,r=A[i].r,v=A[i].v,id=A[i].id;
        if(~id){//询问 
            if(~dp[v]&&dp[v]>r)ans[id]=1;
            else ans[id]=0;
            continue;
        }//update (j+v)<=P  j<=P-v 
        for(j=min(P-v,sum);j>=0;j--){
            if(~dp[j])Max(dp[v+j],min(dp[j],r));//更新dp值
        }
        sum+=v;
    }
    for(i=1;i<=m;i++)pt(ans[i]);
    return 0;
}