【bzoj1568】【JSOI2008】【Blue Mary开公司】【线段树】

Description

Input

第一行 ：一个整数N ，表示方案和询问的总数。 接下来N行，每行开头一个单词“Query”或“Project”。 若单词为Query，则后接一个整数T，表示Blue Mary询问第T天的最大收益。 若单词为Project，则后接两个实数S，P，表示该种设计方案第一天的收益S，以及以后每天比上一天多出的收益P。

Output

对于每一个Query，输出一个整数，表示询问的答案，并精确到整百元（以百元为单位，例如：该天最大收益为210或290时，均应该输出2）。

Sample Input

10
Project 5.10200 0.65000
Project 2.76200 1.43000
Query 4
Query 2
Project 3.80200 1.17000
Query 2
Query 3
Query 1
Project 4.58200 0.91000
Project 5.36200 0.39000

Sample Output

0
0
0
0
0

HINT

约定： 1 <= N <= 100000 1 <= T <=50000 0 < P < 100，| S | <= 10^6 提示：本题读写数据量可能相当巨大，请选手注意选择高效的文件读写方式。

题解：同bzoj3165

#include<iostream>#include<cstring>#include<cmath>#include<cstdio>#define N 500000#define P 500010#define MOD 1000000000#define eps 1e-10using namespace std;int t[P<<2],ansa,n,xx,kind,top,a[P],lastans;double b[P],ansb,x,y,x2,y2,p,q;char ch[10];struct use{   double k,b;   int l,r;   double f(int x){return k*x+b;}}line[N];int cross(use a,use b){  return floor((b.b-a.b)/(a.k-b.k));}use cal(double x,double y,double x2,double y2){   use ans;  ans.r=max(x,x2);ans.l=min(x,x2);  if (x2!=x) ans.k=(y2-y)/(x2-x),ans.b=y-ans.k*x;  else ans.k=0.0,ans.b=max(y,y2);  return ans;}int check(double x){  return (x>-eps)-(x<eps); }void up(int k,int x){   double t=line[k].f(x);   int ff=check(t-b[x]);   if (!a[x]||(ff>0||(ff==0&&k<a[x]))){     a[x]=k;b[x]=t;   }    }void insert(int k,int l,int r,int x){  if (line[x].l<=l&&r<=line[x].r){    if (!t[k]){t[k]=x;return;}    else{      int fl=check(line[x].f(l)-line[t[k]].f(l))>0;      int fr=check(line[x].f(r)-line[t[k]].f(r))>0;       if (fl&&fr) t[k]=x;      else       if(fl||fr){         int mid=(l+r)>>1;         int p=cross(line[x],line[t[k]]);         if (p<=mid&&fl) insert(k<<1,l,mid,x);         if (p<=mid&&fr) insert(k<<1,l,mid,t[k]),t[k]=x;         if (p>mid&&fl) insert(k<<1|1,mid+1,r,t[k]),t[k]=x;         if (p>mid&&fr) insert(k<<1|1,mid+1,r,x);        }      else up(x,l),up(x,r);    }   return;  }  int mid=(l+r)>>1;  if (line[x].l<=mid) insert(k<<1,l,mid,x);  if (line[x].r>mid) insert(k<<1|1,mid+1,r,x);}void query(int k,int l,int r,int x){   if (t[k]){     double tt=line[t[k]].f(x);     int ff=check(tt-ansb);     if ((ff>0||(ff==0&&ansa<t[k]))){        ansa=t[k];ansb=tt;     }   }   if (l==r) return;   int mid=(l+r)>>1;   if (x<=mid) query(k<<1,l,mid,x);   else query(k<<1|1,mid+1,r,x); }int main(){ scanf("%d",&n); for (int i=1;i<=n;i++){   scanf("%s",ch);   if (ch[0]=='P'){     scanf("%lf%lf",&p,&q);     x=1;y=p;x2=P;y2=y+(P-1)*q;     line[++top]=cal(x,y,x2,y2);     insert(1,1,P,top);   }   else{     scanf("%d",&xx);ansa=0;ansb=-1.0;     query(1,1,P,xx);     int ff=check(b[xx]-ansb);     if ((ff>0||(ff==0&&a[xx]<ansa))){        ansa=a[xx];     }     p=line[ansa].f(xx);p=p/100.0;     printf("%d\n",(int)p);   }  }}