【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-10
using 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);
   } 
 }
}