码迷,mamicode.com
首页 > 其他好文 > 详细

BZOJ2924 : [Poi1998]Flat broken lines

时间:2015-08-07 23:52:20      阅读:277      评论: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
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!