李超树

模板

const int MAX_N=1010000;
struct line{
 long long k,b;
    int l,r;
};
struct node{
    line t;
    int flag;
}a[MAX_N<<1];
line ca;
int get_id(int l,int r){
    return (l+r)|(l!=r);
}
long long f(line t,long long pos){
    return t.k*pos+t.b;
}
void build(int l,int r){
    int now=get_id(l,r);
    a[now].flag=0;a[now].t=ca;
    if(l==r) return;
    int mid=(l+r)>>1;
    build(l,mid);
    build(mid+1,r);
}
void change(int l,int r,line k){
    int now=get_id(l,r);
    if(k.l<=l&&r<=k.r){
        if(!a[now].flag){
            a[now].flag=1;
            a[now].t=k;
        }
        else if(f(k,l)>=f(a[now].t,l)&&f(k,r)>=f(a[now].t,r))
        a[now].t=k;
        else if(f(k,l)>f(a[now].t,l)||f(k,r)>f(a[now].t,r)){
            int mid=(l+r)>>1;
            if(f(k,mid)>f(a[now].t,mid))
            swap(a[now].t,k);
            if(f(k,l)>f(a[now].t,l))
            change(l,mid,k);
            else
            change(mid+1,r,k);
        }
    }
    else{
        int mid=(l+r)>>1;
        if(k.l<=mid)
        change(l,mid,k);
        if(mid<k.r)
        change(mid+1,r,k);
    }
}
long long query(int l,int r,int pos){
    int now=get_id(l,r);
    if(l==r)
    return f(a[now].t,pos);
    int mid=(l+r)>>1;
    long long ans=f(a[now].t,pos);
    if(pos<=mid)
    return max(ans,query(l,mid,pos));
    else
    return max(ans,query(mid+1,r,pos));
}

P4254 [JSOI2008]Blue Mary开公司
MAX_N代表总共的边的个数,MAX_T代表x轴的长度。
最弱智的李超树,首先是直线,而且k是大于0的,所以比较好写,但是功能比较low。

#include<iostream>
#include<cstdio>
#include<string>
using namespace std;
const int MAX_N=100010;
const int MAX_T=50100;
struct line{
 double k,b;
}a[MAX_N];
int tr[MAX_T<<2];
int cnt;
double f(int id,int x){
 return a[id].k*(x-1)+a[id].b;
}
void change(int k,int l,int r,int id){
 if(l==r){
  if(f(id,l)>f(tr[k],l))
  tr[k]=id;
  return;
 }
 int mid=(l+r)>>1;
 if(a[tr[k]].k<a[id].k){
  if(f(id,mid)>f(tr[k],mid)){
   change(k<<1,l,mid,tr[k]);
   tr[k]=id;
  }
  else
  change(k<<1|1,mid+1,r,id);
 }
 if(a[tr[k]].k>a[id].k){
  if(f(id,mid)>f(tr[k],mid)){
   change(k<<1|1,mid+1,r,tr[k]);
   tr[k]=id;
  }
  else
  change(k<<1,l,mid,id);
 }
}
double query(int k,int l,int r,int x){
    if(l==r)
        return f(tr[k],x);
    int mid=(l+r)>>1;
    if(x<mid)
        return max(f(tr[k],x),query(k<<1,l,mid,x));
    else
        return max(f(tr[k],x),query(k<<1|1,mid+1,r,x));
}
int main(void){
    int i,n,x;
    int T=50000;
    scanf("%d",&n);
    for(i=0;i<n;i++){
        string s;
        cin>>s;
        if(s[0]=='P'){
            cnt++;
            scanf("%lf%lf",&a[cnt].b,&a[cnt].k);
            change(1,1,T,cnt);
        }
        else{
            scanf("%d",&x);
            printf("%d\n",(int)query(1,1,T,x)/100);
        }
    }
    return 0;
}

【2019南京icpc网络赛 I】Washing clothes
首先如果一个人用洗衣机洗衣服的话那么他后面的人都用洗衣机。因为如果这个人用洗衣机而他后面的人用手洗的话,那还不如让后面的人用洗衣机,他手洗。
对于第i个人y=t[i]+b和y=t[i]+(n-i+1)*x取小的两段。
从1到y询问当前x轴值的最大值。

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int MAX_N=1010000;
struct line{
 long long k,b;
 int l,r;
};
struct node{
 line t;
 int flag;
}a[MAX_N<<1];
line ca;
int get_id(int l,int r){
 return (l+r)|(l!=r);
}
long long f(line t,long long pos){
 return t.k*pos+t.b;
}
void build(int l,int r){
 int now=get_id(l,r);
 a[now].flag=0;a[now].t=ca;
 if(l==r) return;
 int mid=(l+r)>>1;
 build(l,mid);
 build(mid+1,r);
}
void change(int l,int r,line k){
 int now=get_id(l,r);
 if(k.l<=l&&r<=k.r){
  if(!a[now].flag){
   a[now].flag=1;
   a[now].t=k;
  }
  else if(f(k,l)>=f(a[now].t,l)&&f(k,r)>=f(a[now].t,r))
  a[now].t=k;
  else if(f(k,l)>f(a[now].t,l)||f(k,r)>f(a[now].t,r)){
   int mid=(l+r)>>1;
   if(f(k,mid)>f(a[now].t,mid))
   swap(a[now].t,k);
   if(f(k,l)>f(a[now].t,l))
   change(l,mid,k);
   else
   change(mid+1,r,k);
  }
 }
 else{
  int mid=(l+r)>>1;
  if(k.l<=mid)
  change(l,mid,k);
  if(mid<k.r)
  change(mid+1,r,k);
 }
}
long long query(int l,int r,int pos){
 int now=get_id(l,r);
 if(l==r)
 return f(a[now].t,pos);
 int mid=(l+r)>>1;
 long long ans=f(a[now].t,pos);
 if(pos<=mid)
 return max(ans,query(l,mid,pos));
 else
 return max(ans,query(mid+1,r,pos));
}
long long t[MAX_N];
int main(void){
 int n,i;
 long long y,b;
 ca.k=ca.b=0;ca.l=ca.r=0;
 while(scanf("%d%lld",&n,&y)!=EOF){
    build(1,y);
    for(i=1;i<=n;i++)
        scanf("%lld",&t[i]);
    sort(t+1,t+n+1);
 for(i=1;i<=n;i++){
        b=t[i];
  int x1,x2;
  if(y%(n-i+1)==0){
   x1=x2=y/(n-i+1);
  }
  else{
   x1=y/(n-i+1);
   x2=x1+1;
  }
  line s1,s2;
  s1.b=b;s1.k=(n-i+1);
  s1.l=1;s1.r=x1;
  s2.b=b+y;s2.k=0;
  s2.l=x2;s2.r=y;
        if(x1>=1)
  change(1,y,s1);
  if(x2<=y)
  change(1,y,s2);
 }
 for(i=1;i<=y;i++){
  printf("%lld",query(1,y,i));
  if(i!=y)
  printf(" ");
  else
  printf("\n");
 }
 }
 return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值