HDU 5808 BestCoder Round #86 Price List Strike Back （树状数组）

Price List Strike Back  

题目链接：点我打开链接

Source

BestCoder Round #86

题意：

 在Byteland一共有n家商店，编号依次为1到n。每家商店只会卖一种物品，其中第i家商店的物品单价为v_i，且它到Byteasar的家的距离为d_i。
Byteasar每天都会进行一次购物，第i天他会选择一个区间[l_i,r_i]，并给自己设定一个距离上限c_i，然后他会在编号在该区间内每家到自己家的距离不超过c_i的商店购买最多一件物品，当然他也可以选择什么都不买。回家之后，Byteasar会把今天购物所花的钱的总数sum_i记录在账本上。
Byteasar的数学不好，他可能会把花的钱记错。
请写一个程序，帮助Byteasar判断每条记录是否一定是错的。
注意：记多或者记少都算记错。

输入的第一行包含一个正整数T，表示测试数据的组数。
对于每组数据，第一行包含两个正整数n,m，表示商店的个数和记录的个数。
第二行包含n个正整数v_i，依次表示每家商店的物品的单价。
第三行包含n个正整数d_i，依次表示每家商店到Byteasar的家的距离。
接下来m行，每行包含四个整数l_i,r_i,c_i,sum_i，表示一条记录。

对于每组数据，输出一行m个字符，依次回答每个询问。如果一定记错了，请输出'1'，否则输出'0'。

题解：这道题涉及到一个距离的影响。 所以要把距离的影响去掉，可以考虑对"所有商店距离家的距离" 和"所有询问距离家的距离" 分别进行升序排序，这就可以使得，在处理一个询问的时候，我们可以只考虑了所有距离范围内的商店。 而在过程中需要用树状数组维护前缀和，然后求出区间范围内某种价值的物品有多少件。（这题比赛时没有做出来。。。我还是鶸）

官方题解：

那么所有在[ l , mid )的询问和( mid , r ]的询问可以递归分治求解，故只需考虑必然经过mid的询问。 设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( fi , gs − i ))，然后和c比较一下大小就好了。 时间复杂O(100(nlogn+m)+mlogn)。

AC代码：

//#include<bits/stdc++.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
#include<map>
#include<queue>
#include<set>
#include<stack>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int N = 20000+10;
const int mod=998244353;
const int M=1e5+10;
const int INF= 1e9;
struct A
{
	int v,d,ord;
}a[N];

struct Q
{
	int l,r,d,sum,ord;
}q[M];

int n,m,C[N][110],ans[M];

int cmp(A a, A b)
{
	return a.d < b.d;
}

int cmp2(Q a, Q b)
{
	return a.d < b.d;
}

int lowbit(int x)
{
	return x & (-x);
}

void add(int x, int y)
{
	for(int i = x ;i <= n; i += lowbit(i)){
		C[i][y]++;
	}
}

int query(int x,int y)
{
	int ans = 0;
	for(int i = x; i >= 1; i -= lowbit(i)){
		ans+=C[i][y];
	}
	return ans;
}

bool e[110];
int b[3010];
int f[110];
int main()
{
	//freopen("in.txt","r",stdin);
	// freopen("out.txt","w",stdout);
	int t;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++)scanf("%d",&a[i].v);
		for(int i=1;i<=n;i++){
			scanf("%d",&a[i].d), a[i].ord=i;
		}
		sort(a+1,a+n+1,cmp);
		memset(C,0,sizeof(C));
		for(int i=1; i <= m; i++){
			scanf("%d%d%d%d",&q[i].l, &q[i].r, &q[i].d, &q[i].sum);
			q[i].ord = i;
			
		}
		sort( q+1 , q+m+1 , cmp2);
		a[n+1].d = INF+10;
		int l = 1;
		for(int i=1; i<=m; i++)
		{
			while(a[l].d <= q[i].d && l <= n)
			{
				add(a[l].ord,a[l].v);
				l++;
			}
			int Ans = query(q[i].r, q[i].sum)- query(q[i].l-1 , q[i].sum);
			if(Ans) ans[q[i].ord] = 0;
			//一定买不了 
			else{
				int cnt = 0;
				for(int j = 1; j < q[i].sum; j++){
					int y = query(q[i].r , j) - query(q[i].l - 1, j);
					int p = 1;
					while(p < y){
						b[++cnt] = j * p;
						y -= p;
						p <<= 1;
					}
					if(y) b[++cnt] = y*j;
				}
				memset(f,0,sizeof(f));
				f[0] = 1;
				for(int j = 1; j <= cnt; j++){
					for(int v = q[i].sum; v >= b[j]; --v){
						f[v]= max(f[v],f[v-b[j]]);
					}
				}
				if(f[q[i].sum]) ans[q[i].ord] = 0;
				else ans[q[i].ord] = 1;
			}
		}
		for(int i=1; i<=m ; i++) printf("%d",ans[i]);
		printf("\n");
	}
	return 0;
}