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

分析：
超哥线段树，标记永久化处理
一开始按照自己的想法写的，但是发现有两个小地方理解的不好

if (x<=mid)         //交点在左区间 
{
    if (z4<z2) 
        add(bh<<1,l,mid,a,b);
    else 
    {
        add(bh<<1,l,mid,T[bh].a,T[bh].b);
        T[bh].a=a; T[bh].b=b;
    }
} 
else               //交点在右区间 
{
    if (z2<z4) 
        add((bh<<1)+1,mid+1,r,a,b);
    else 
    {
        add((bh<<1)+1,mid+1,r,T[bh].a,T[bh].b);
        T[bh].a=a; T[bh].b=b;
    }   
}

为了加深印象，就详细分析一下把：
在区间添加一条线段，我们要分情况讨论

• 新线段完全在旧线段上，或者区间没有线段覆盖，直接修改成新线段
• 新线段完全在旧线段下，不用修改
• 新线段与旧线段有交点
  
  • 交点在左儿子
    
    
    • 如左图，右区间的新线段是肯定影响不到的，我们只用递归处理左区间，最左端在上边的新线段，所以我们把新线段传下去：
      add(bh<<1,l,mid,a,b);
    • 如右图，右区间是可以完全覆盖的，但是因为我们贯彻了标记永久化的思想，我们只用在当前结点上打一个标记修改，但是左区间并不是完全覆盖，所以我们还是需要递归处理左区间，最左端在上边的旧线段，所以我们把旧线段传下去：
      add(bh<<1,l,mid,T[bh].a,T[bh].b);
T[bh].a=a; T[bh].b=b;
  • 交点在右儿子
    
    
    • 如左图，新线段在左区间是没有影响的，最右端在上面的新线段，所以我们把新线段传下去：
      add((bh<<1)+1,mid+1,r,a,b);
    • 如右图，左区间是可以完全覆盖的，但是因为我们贯彻了标记永久化的思想，我们只用在当前结点上打一个标记修改，但是右区间并不是完全覆盖，所以我们还是需要递归处理右区间，最右端在上边的旧线段，所以我们把旧线段传下去：
      add((bh<<1)+1,mid+1,r,T[bh].a,T[bh].b);
T[bh].a=a; T[bh].b=b;

tip

一开始把K和B搞错了
卡了一个小时。。。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define ll long long
#define LD long double

using namespace std;

const int N=50000;
int n,pos,tot=0;
double x,y,ans;
struct node{
    double a,b;
};
node T[(N+10)<<2];
bool flag[(N+10)<<2];

void add(int bh,int l,int r,double a,double b)
{
    if (!flag[bh])
    {
        T[bh].a=a; T[bh].b=b;
        flag[bh]=1;
        return;
    }

    int mid=(l+r)>>1;
    double z1=(double)l*T[bh].a+T[bh].b;    //原线段 
    double z2=(double)r*T[bh].a+T[bh].b;
    double z3=(double)l*a+b;                //新线段 
    double z4=(double)r*a+b;

    if (z1>=z3&&z2>=z4) return;               //不用修改 
    if (z1<z3&&z2<z4)                         //直接覆盖 
    {
        T[bh].a=a; T[bh].b=b;
        return;
    }
    double x=(double)(T[bh].b-b)/(a-T[bh].a);     //交点 
    if (x<=mid)         //交点在左区间 
    {
        if (z4<z2) 
            add(bh<<1,l,mid,a,b);
        else 
        {
            add(bh<<1,l,mid,T[bh].a,T[bh].b);
            T[bh].a=a; T[bh].b=b;
        }
    } 
    else               //交点在右区间 
    {
        if (z2<z4) 
            add((bh<<1)+1,mid+1,r,a,b);
        else 
        {
            add((bh<<1)+1,mid+1,r,T[bh].a,T[bh].b);
            T[bh].a=a; T[bh].b=b;
        }   
    }
}

double get_y(node x,double p)
{
    return (double)p*x.a+x.b;
}

void ask(int bh,int l,int r,int pos)
{
    if (flag[bh]) 
        ans=max(ans,get_y(T[bh],pos));    //这个区间上有线段 
    if (l==r) return;
    int mid=(l+r)>>1;
    if (pos<=mid) ask(bh<<1,l,mid,pos);
    else ask((bh<<1)+1,mid+1,r,pos);
}

int main()
{
    scanf("%d",&n);
    char s[20];
    for (int i=1;i<=n;i++)
    {
        scanf("%s",s);
        if (s[0]=='P')
        {
            scanf("%lf%lf",&x,&y);
            add(1,1,N,y,x-y);
        }
        else
        {
            scanf("%d",&pos);
            ans=0.0; ask(1,1,N,pos);
            printf("%lld\n",(ll)floor(ans/100.0));
        }
    }
    return 0;
}