Bzoj2850 巧克力王国

Time Limit: 60 Sec  Memory Limit: 512 MB
Submit: 505  Solved: 204

Description

巧克力王国里的巧克力都是由牛奶和可可做成的。但是并不是每一块巧克力都受王国人民的欢迎，因为大家都不喜

欢过于甜的巧克力。对于每一块巧克力，我们设x和y为其牛奶和可可的含量。由于每个人对于甜的程度都有自己的

评判标准，所以每个人都有两个参数a和b，分别为他自己为牛奶和可可定义的权重，因此牛奶和可可含量分别为x

和y的巧克力对于他的甜味程度即为ax + by。而每个人又有一个甜味限度c，所有甜味程度大于等于c的巧克力他都

无法接受。每块巧克力都有一个美味值h。现在我们想知道对于每个人，他所能接受的巧克力的美味值之和为多少

Input

第一行两个正整数n和m，分别表示巧克力个数和询问个数。接下来n行，每行三个整数x,y,h，含义如题目所示。再

接下来m行，每行三个整数a,b,c，含义如题目所示。

Output

输出m行，其中第i行表示第i个人所能接受的巧克力的美味值之和。

Sample Input

3 3
1 2 5
3 1 4
2 2 1
2 1 6
1 3 5
1 3 7

Sample Output

5
0
4

HINT

1 <= n, m <= 50000，1 <= 10^9，-10^9 <= a, b, x, y <= 10^9。

Source

Violet 0

K-D tree

以x和y为坐标，通过ax+by估价来优化查询

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<cstring>
 #include<queue>
 using namespace std;
 const int mxn=;
 int read(){
     int x=,f=;char ch=getchar();
     while(ch<'' || ch>''){if(ch=='-')f=-;ch=getchar();}
     while(ch>='' && ch<=''){x=x*-''+ch;ch=getchar();}
     return x*f;
 }
 struct node{
     int max[],min[];
     int l,r;
     int d[],v;
     long long sum;
 }t[mxn];
 int nowD=;
 int cmp(const node a,const node b){
     return (a.d[nowD]<b.d[nowD] || (a.d[nowD]==b.d[nowD] && a.d[nowD^]<b.d[nowD^]));
 }
 int root,nct=;
 int n,m;
 long long a,b,c,ans=;
 void pushup(int rt,int x){
     t[rt].max[]=max(t[rt].max[],t[x].max[]);
     t[rt].max[]=max(t[rt].max[],t[x].max[]);
     t[rt].min[]=min(t[rt].min[],t[x].min[]);
     t[rt].min[]=min(t[rt].min[],t[x].min[]);
     return;
 }
 int Build(int l,int r,int D){
     nowD=D;
     int mid=(l+r)>>;
     nth_element(t+l,t+mid,t+r+,cmp);
     t[mid].max[]=t[mid].min[]=t[mid].d[];
     t[mid].max[]=t[mid].min[]=t[mid].d[];
     if(l!=mid){t[mid].l=Build(l,mid-,D^);pushup(mid,t[mid].l);}
     if(r!=mid){t[mid].r=Build(mid+,r,D^);pushup(mid,t[mid].r);}
     t[mid].sum=t[t[mid].l].sum+t[t[mid].r].sum+t[mid].v;
     return mid;
 }
 inline bool pd(long long x,long long y){return (a*x+b*y<c);}
 int cnt(int rt){
     int res=;
     if(pd(t[rt].max[],t[rt].max[]))res++;
     if(pd(t[rt].min[],t[rt].max[]))res++;
     if(pd(t[rt].max[],t[rt].min[]))res++;
     if(pd(t[rt].min[],t[rt].min[]))res++;
     return res;
 }
 void query(int rt){
     if(pd(t[rt].d[],t[rt].d[]))ans+=t[rt].v;
     int L=,R=;
     if(t[rt].l)L=cnt(t[rt].l);
     if(t[rt].r)R=cnt(t[rt].r);
     if(L==){ans+=t[t[rt].l].sum;}//还要算R所以不能return
         else if(L)query(t[rt].l);
     if(R==){ans+=t[t[rt].r].sum;}
         else if(R)query(t[rt].r);
     return;
 }
 int main(){
     n=read();m=read();
     int i,j;
     for(i=;i<=n;i++){
         t[i].d[]=read();t[i].d[]=read();t[i].v=read();
     }
     root=Build(,n,);
     for(i=;i<=m;i++){
         a=read();b=read();c=read();
         ans=;
         query(root);
         printf("%lld\n",ans);
     }
     return ;
 }