3341: [Ceoi2013]adriatic

题目链接

题目大意：平面上有一些点，只能想正左上方或正右下方走，问每个点到所有其他点的距离和

题解：比较玄妙的分治思想……题解

我的收获：思想

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
using namespace std;

const int M=2505;

int n;
bool on[M][M];//这个点是否有船 
int maxx[M][M],maxy[M][M];//(i,j)~（n,n）中最远的x坐标，和最远的y坐标。
int minx[M][M],miny[M][M];//(1,1)~（i,j）中最小的x坐标，和最小的y坐标。
int cnt[M][M];//前缀和 
int mur[M][M],mdl[M][M];//记忆化

struct point{int x,y;}p[M*100];

int ur(int x,int y)//计算访问(x,y) 右上方的额外代价 
{
    int &ans=mur[x][y];
    if(ans!=-1) return ans;
    if(cnt[x][2500]-cnt[x][y-1]==0) return ans=0;//没有点要访问
    return ans=cnt[x][2500]-cnt[x][y-1]+ur(min(x,minx[x][y-1]),max(y,maxy[x+1][y]));//还要加上当前一次拓展都到不了的额外代价 
}

int dl(int x,int y)//计算访问(x,y) 左下方的额外代价，大致上同上。
{
    int &ans=mdl[x][y];
    if(ans!=-1) return ans;
    if(cnt[2500][y]-cnt[x-1][y]==0) return ans=0;
    return ans=cnt[2500][y]-cnt[x-1][y]+dl(max(x,maxx[x][y+1]),min(y,miny[x-1][y]));
} 

void pre()
{
    memset(minx,0x3f,sizeof(minx));
    memset(miny,0x3f,sizeof(miny));
    for(int i=1;i<=2500;i++)
        for (int j=1;j<=2500;j++)
        {
            cnt[i][j]=on[i][j]+cnt[i-1][j]+cnt[i][j-1]-cnt[i-1][j-1];
            minx[i][j]=on[i][j]?min(i,minx[i-1][j]):min(minx[i-1][j],minx[i][j-1]);
            miny[i][j]=on[i][j]?min(j,miny[i][j-1]):min(miny[i-1][j],miny[i][j-1]);
        }
    for(int i=2500;i>=1;i--)
        for(int j=2500;j>=1;j--)
        {
            maxx[i][j]=on[i][j]?max(i,maxx[i+1][j]):max(maxx[i+1][j],maxx[i][j+1]);
            maxy[i][j]=on[i][j]?max(j,maxy[i][j+1]):max(maxy[i+1][j],maxy[i][j+1]);
        }
}

void work()
{
    pre();
    memset(mur,-1,sizeof(mur));
    memset(mdl,-1,sizeof(mdl));
    for(int i=1;i<=n;i++)
        printf("%d\n",cnt[2500][2500]+ur(p[i].x,p[i].y)+dl(p[i].x,p[i].y)-3);//-3因为i这个点被多算了3次 
}

void init()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&p[i].x,&p[i].y);
        on[p[i].x][p[i].y]=true;
    }
}

int main()
{
    init(); 
    work();
    return 0;
}