bzoj2149拆迁队 斜率优化dp+分治

2149: 拆迁队

Time Limit: 10 Sec  Memory Limit: 259 MB
Submit: 397  Solved: 177
[Submit][Status][Discuss]

Description

lanxisi带领着他的拆迁队来整治一个街道。这个街道由N个旧房子组成，从左到右编号为1..N。每个旧房子i有一个正整数的美观度A_i。

lanxisi希望整个街道从左到右美观度严格递增，也就是保证A_i<A_j（i<j）。但是旧的街道明显不符合这个要求，于是lanxisi希望拆迁一些旧房子并在原地建立新房子来满足这一要求。但是与很多拆迁队一样，很多钉子户拒绝拆迁。所以lanxisi希望能保留最多的旧房子来满足某些钉子户，当然，保留一个旧房子需要给房子主人B_i的赔偿金。最后，总花费=整治好以后所有房子美观度之和+赔偿金之和。新的美观值也都必须是正整数。

现在，请你求出lanxisi最多能保留多少旧房子和整治这个街道所需要的最少总花费（当然是在保留旧房子最多这个前提下的）。

Input

第一行有1个整数N，表示旧房子个数。第二行有N个整数，第i个数Ai表示旧房子的美观度。第三行有N个整数，第i个数Bi表示保留旧房子的赔偿金。

Output

第一行输出2个整数，表示保留的旧房子数量、总花费的大小。【样例输入一】 4 2 1 7 4 1 5 4 3 【样例输入二】 6 1 6 3 4 7 7 1 2 1 1 1 9

Sample Input

【样例输出一】
2 25

【样例输出二】
4 29

【数据范围】
100%的数据中，1<=N<=100000，0<=Ai,Bi<=100000000。

Sample Output

 

HINT

 

Source

这道题2个问 第一问是减去下标求lis 第二问在第一问基础上维护转移

还是有点麻烦 不多哔哔了

这篇博客写得不错http://www.cnblogs.com/zj75211/p/8125193.html

贴个代码

 #include<cstdio>
 #include<algorithm>
 #define N 100005
 using namespace std;
 typedef long long ll;
 int n,m,i,j,k,a[N],b[N],d[N],e[N],bit[N],f[N],ans0,G[N],nxt[N];
 int cl,cr,L[N],R[N],top,q[N];
 ll c[N],g[N],ans1=1LL<<;
 struct E{int x,t;E(){}E(int _x,int _t){x=_x,t=_t;}}h[N];
 void read(int&a){char c;while(!(((c=getchar())>='')&&(c<='')));a=c-'';while(((c=getchar())>='')&&(c<=''))(a*=)+=c-'';}
 void umax(int&a,int b){if(a<b)a=b;}
 void umin(ll&a,ll b){if(a>b)a=b;}
 int lower(int x){
   int l=,r=n+,mid,t;
   while(l<=r)if(e[mid=(l+r)>>]<=x)l=(t=mid)+;else r=mid-;
   return t;
 }
 void ins(int x,int y){for(;x<=n+;x+=x&-x)umax(bit[x],y);}
 int ask(int x){int t=-N;for(;x;x-=x&-x)umax(t,bit[x]);return t;}
 void add(int x,int y){nxt[y]=G[x];G[x]=y;}
 bool cmp(int x,int y){return d[x]==d[y]?c[x]>c[y]:d[x]<d[y];}
 double pos(int x,int y){return 1.0*(c[y]-c[x])/(d[x]-d[y]);}
 void insert(int x){
   if(top&&d[x]==d[q[top]])top--;
   while(top>&&pos(q[top],x)>pos(q[top],q[top-]))top--;
   q[++top]=x;
 }
 void query(int x){
   if(!top)return;
   int l=,r=top-,mid,t=top;
   while(l<=r){
     mid=(l+r)>>;
     if(x>pos(q[mid],q[mid+]))r=(t=mid)-;else l=mid+;
   }
   umin(g[x],1LL*x*d[q[t]]+c[q[t]]);
 }
 void solve(int l,int r){
   if(l==r)return;
   int mid=(l+r)>>,i,j;
   solve(l,mid),solve(mid+,r);
   for(cl=,i=l;i<=mid;i++)if(!h[i].t)L[++cl]=h[i].x;
   for(cr=,i=r;i>mid;i--)if(h[i].t)R[++cr]=h[i].x;
   if(!cl||!cr)return;
   sort(L+,L+cl+,cmp),sort(R+,R+cr+,cmp);
   for(i=j=,top=;i<=cr;i++){
     while(j<=cl&&d[L[j]]<=d[R[i]])insert(L[j++]);
     query(R[i]);
   }
 }
 int main(){
   read(n);
   for(i=;i<=n;i++)read(a[i]);
   for(i=;i<=n;i++)read(b[i]),e[i+]=d[i]=a[i]-i;
   sort(e+,e+n+);
   for(i=;i<=n+;i++)bit[i]=-N;
   ins(lower(),);
   for(i=;i<=n;i++)g[i]=ans1,j=lower(d[i]),ins(j,f[i]=ask(j)+);
   ans0=ask(n+);
   for(i=;i<=n;i++)G[i]=-;
   for(i=n;~i;i--)if(f[i]>=)add(f[i],i);
   for(i=;i<ans0;i++){
     j=G[i],k=G[i+],m=;
     while(~j||~k){
       if(j<)h[++m]=E(k,),k=nxt[k];
       else if(k<||j<k)h[++m]=E(j,),j=nxt[j];
       else h[++m]=E(k,),k=nxt[k];
     }
     solve(,m);
     for(k=G[i+];~k;k=nxt[k]){
       g[k]+=b[k]+a[k]+1LL*k*(k-)/;
       c[k]=g[k]+1LL*k*(k+)/-1LL*a[k]*k-a[k];
     }
   }
   for(i=;i<=n;i++)if(f[i]==ans0)umin(ans1,g[i]+1LL*(n-i)*a[i]+1LL*(n-i)*(n-i+)/);
     printf("%d %lld",ans0,ans1);
   return ;
 }