bzoj 3341: [Ceoi2013]adriatic 动态规划

题意

琪琪在大海上搜索涛涛，琪琪把大海分为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;
}