首页 > 代码库 > BZOJ4700: 适者
BZOJ4700: 适者
先排序,枚举删一个点,在前面找出最优的另一个点,容易推出斜率方程,平衡树维护凸包。
#include<bits/stdc++.h> using namespace std; typedef long long ll; const int N=3e5+5; typedef struct node*ptr; struct node{ ptr i,j,s,t; ll w,x,y; node(){w=rand();} }*b,e[N]; void zag(ptr&o){ptr s=o->j;o->j=s->i,s->i=o,o=s;} void zig(ptr&o){ptr s=o->i;o->i=s->j,s->j=o,o=s;} void ins(ptr j,ptr&o=b){ if(!o)o=j,o->s?o->s->t=o:0,o->t?o->t->s=o:0; else if(j->x>o->x) {ins(j,o->j);if(o->j->w>o->w)zag(o);} else if(o->x>j->x) {ins(j,o->i);if(o->i->w>o->w)zig(o);} } void del(ptr j,ptr&o=b){ if(j->x>o->x)del(j,o->j); else if(o->x>j->x)del(j,o->i); else if(!o->i)o=o->j; else if(!o->j)o=o->i; else if(o->i->w>o->j->w)zig(o),del(j,o->j); else zag(o),del(j,o->i); } ptr pre(ll x,ptr o=b){ ptr s=0; while(o)x>=o->x?s=o,o=o->j:o=o->i; return s; } ptr suc(ll x,ptr o=b){ ptr s=0; while(o)o->x>=x?s=o,o=o->i:o=o->j; return s; } ll cal(ptr o,ptr s,ptr t){ ll x1=o->x-s->x,y1=o->y-s->y; ll x2=t->x-s->x,y2=t->y-s->y; return x1*y2-x2*y1; } void upd(ptr o){ if(ptr&s=o->s=pre(o->x)) while(s->s&&cal(o,s->s,s)<=0) del(s),s=pre(o->x); if(ptr&t=o->t=suc(o->x)) while(t->t&&cal(o,t,t->t)<=0) del(t),t=suc(o->x); if(!o->s||!o->t)ins(o); else if(cal(o->s,o,o->t)>0)ins(o); } ll slo(ptr s,ptr t){return(s->y-t->y)/(s->x-t->x);} ptr sol(ll k,ptr o=b){ if(o->s&&k>slo(o->s,o))return sol(k,o->i); if(o->t&&slo(o,o->t)>k)return sol(k,o->j); return o; } struct foo{int s,t;}c[N]; bool operator<(foo s,foo t){return s.t*t.s<t.t*s.s;} int n,m,d; ll k,q,s[N],r[N]; int main(){ scanf("%d%d",&n,&m); for(int i=1;i<=n;++i){ scanf("%d%d",&c[i].s,&d); c[i].t=(d+m-1)/m; } sort(c+1,c+n+1); for(int i=1;i<=n;++i) s[i]=s[i-1]+c[i].t; for(int i=n;i;--i){ k+=(s[i]-1)*c[i].s; r[i-1]=r[i]+c[i].s; e[i].x=c[i].t; e[i].y=(s[i]-1)*c[i].s+c[i].t*r[i]; } upd(e+1); for(int i=2;i<=n;++i){ ptr j=sol(c[i].s); q=max(q,e[i].y+j->y-j->x*c[i].s); upd(e+i); } printf("%lld\n",k-q); }
BZOJ4700: 适者
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。