bzoj 3341: [Ceoi2013]adriatic 动态规划
2017-12-16 15:46
337 查看
题意
琪琪在大海上搜索涛涛,琪琪把大海分为2500*2500的网格,现在琪琪拥有N条小帆船,她可以在船之间移动,由于风向的原因,她只能从一条船移动到位于这条船严格左上方或严格右下方的船。对于每条帆船,琪琪想要知道从这条船跳到其他各条船的最小跳转次数和。3 ≤ N ≤ 250 000
分析
咋一看还以为是KDtree。。。一开始很直观的想法就是先把左上和右下两块搞掉,然后对于左下的一块,找到最左边和左下边的一个点来跳,这样就又把矩形变成了一个更小的矩形。但是对每个点这样切的话会T,那么只要记忆化一下即可。
代码
#include<iostream> #include<cstdio> #include<cstdlib> #include<cstring> #include<algorithm> using namespace std; const int N=2505; const int inf=1000000000; const int M=2500; int n,s[N][N],ur[N][N],dl[N][N],ul[N][N],dr[N][N],drd[N][N],drr[N][N],ulu[N][N],ull[N][N]; struct data{int x,y;}a[250005]; int read() { int x=0,f=1;char ch=getchar(); while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } int get_dl(int x,int y) { if (x<1||x>M||y<1||y>M) return 0; if (dl[x][y]>-1) return dl[x][y]; dl[x][y]=0; int d=min(drd[x][y+1],x),l=min(ull[x+1][y],y); return dl[x][y]=get_dl(d,l)+dr[x][1]-dr[x][y+1]; } int get_ur(int x,int y) { if (x<1||x>M||y<1||y>M) return 0; if (ur[x][y]>-1) return ur[x][y]; ur[x][y]=0; int u=max(ulu[x][y-1],x),r=max(drr[x-1][y],y); return ur[x][y]=get_ur(u,r)+ul[x][M]-ul[x][y-1]; } int main() { n=read(); for (int i=1;i<=n;i++) a[i].x=read(),a[i].y=read(),a[i].y=M-a[i].y+1,s[a[i].x][a[i].y]++; for (int i=0;i<=M+1;i++) for (int j=0;j<=M+1;j++) { ull[i][j]=drd[i][j]=inf; ulu[i][j]=drr[i][j]=-1; dl[i][j]=ur[i][j]=-1; } for (int i=M;i>=1;i--) for (int j=1;j<=M;j++) { ull[i][j]=min(ull[i][j-1],ull[i+1][j]); ulu[i][j]=max(ulu[i][j-1],ulu[i+1][j]); if (s[i][j]) ull[i][j]=min(ull[i][j],j),ulu[i][j]=max(ulu[i][j],i); ul[i][j]=ul[i][j-1]+ul[i+1][j]-ul[i+1][j-1]+s[i][j]; } for (int i=1;i<=M;i++) for (int j=M;j>=1;j--) { drr[i][j]=max(drr[i][j+1],drr[i-1][j]); drd[i][j]=min(drd[i][j+1],drd[i-1][j]); if (s[i][j]) drr[i][j]=max(drr[i][j],j),drd[i][j]=min(drd[i][j],i); dr[i][j]=dr[i][j+1]+dr[i-1][j]-dr[i-1][j+1]+s[i][j]; } for (int i=1;i<=n;i++) { int x=a[i].x,y=a[i].y; int ans=get_dl(x,y)+get_ur(x,y)+n-3; printf("%d\n",ans); } return 0; }
相关文章推荐
- bzoj 3341: [Ceoi2013]adriatic
- 3341: [Ceoi2013]adriatic
- BZOJ 1948 CEOI 2006 Connect 状态压缩动态规划
- jzoj3578 【CEOI2013】Adriatic 巧妙dp
- BZOJ2013 [Ceoi2010]A huge tower 乱搞
- BZOJ_P4416 [SHOI2013]阶乘字符串(状态压缩动态规划)
- [BZOJ 3111] ZJOI 2013 蚂蚁寻路 · 动态规划
- 【CEOI2010】【BZOJ2013】A huge tower
- BZOJ3107 CQOI2013二进制A+B(动态规划)
- bzoj 3233: [Ahoi2013]找硬币 动态规划
- bzoj 3111: [Zjoi2013]蚂蚁寻路 动态规划
- [bzoj3126][USACO2013]Photo_动态规划_单调队列
- BZOJ 3233 Ahoi2013 找硬币 动态规划
- bzoj2013 [Ceoi2010]A huge tower
- bzoj3341: [Ceoi2013]adriatic 记忆化搜索
- 【CEOI2013】Adriatic
- bzoj2013[CEOI2010] A huge tower
- [BZOJ3242][NOI2013]快餐店-基环树-动态规划
- BZOJ 3233: [Ahoi2013]找硬币|动态规划
- [bzoj3165][Heoi2013]Segment——李超线段树