首页 > 其他 > 详细

BZOJ2924 : [Poi1998]Flat broken lines

时间:2015-08-07 23:50:43      阅读:336      评论:0      收藏:0      [点我收藏+]

首先旋转坐标系

$x‘=x-y$

$y‘=-x-y$

则对于一个点,它下一步可以往它左上角任意一个点连线。

根据Dilworth定理,答案=这个偏序集最长反链的长度。

设f[i]为到i点为止的最长反链长度,则

f[i]=max(f[j])+1,j在i的左下角

按x坐标排序后用树状数组优化DP即可,时间复杂度$O(n\log n)$。

 

#include<cstdio>
#include<algorithm>
#define N 30010
int n,i,j,x,y,b[N],bit[N],f[N],ans;
struct P{int x,y;}a[N];
inline bool cmp(P a,P b){return a.x<b.x;}
inline void up(int&a,int b){if(a<b)a=b;}
inline void ins(int x,int y){for(;x<=n;x+=x&-x)up(bit[x],y);}
inline int ask(int x){int t=0;for(;x;x-=x&-x)up(t,bit[x]);return t;}
inline int lower(int x){
  int l=1,r=n,mid,t;
  while(l<=r)if(b[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
  return t;
}
inline void read(int&a){char c;while(!(((c=getchar())>=‘0‘)&&(c<=‘9‘)));a=c-‘0‘;while(((c=getchar())>=‘0‘)&&(c<=‘9‘))(a*=10)+=c-‘0‘;}
int main(){
  for(read(n),i=1;i<=n;i++)read(x),read(y),a[i].x=x-y,a[i].y=b[i]=-x-y;
  for(std::sort(a+1,a+n+1,cmp),std::sort(b+1,b+n+1),i=1;i<=n;i=j){
    for(j=i;j<=n&&a[j].x==a[i].x;j++)up(ans,f[j]=ask((a[j].y=lower(a[j].y))-1)+1);
    for(j=i;j<=n&&a[j].x==a[i].x;j++)ins(a[j].y,f[j]);
  }
  return printf("%d",ans),0;
}

  

BZOJ2924 : [Poi1998]Flat broken lines

原文:http://www.cnblogs.com/clrs97/p/4712249.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!