模板
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;
}