BZOJ.3170.[TJOI2013]松鼠聚会(切比雪夫距离转曼哈顿距离)

题目链接

将原坐标系每个点的坐标\((x,y)\)变为\((x+y,x-y)\),则原坐标系中的曼哈顿距离等于新坐标系中的切比雪夫距离。
反过来,将原坐标系每个点的坐标\((x,y)\)变为\((\frac{x+y}{2},\frac{x-y}{2})\),则原坐标系中的切比雪夫距离等于新坐标系中的曼哈顿距离。
随便写两个点就可以验证这是对的。

将题目中每个点的坐标\((x,y)\)改为\((\frac{x+y}{2},\frac{x-y}{2})\),然后记\(dis(a,b)=\Delta X+\Delta Y\)表示两个点\(a,b\)间的曼哈顿距离。
枚举每一个点\(x\),则到\(x\)的距离之和为\(Sum=\sum_{i=1}^ndis(i,x)\)
\(dis(i,x)\)拆开,即\[\begin{aligned}Sum&=\sum_{i=1}^n\Delta X(i,x)+\Delta Y(i,x)\\&=\Delta X(1,x)+\Delta X(2,x)+...+\Delta Y(1,x)+\Delta Y(2,x)...\end{aligned}\]

把每个点按\(x\)\(y\)排序,枚举点时就可以\(O(1)\)计算\(\Delta X\)\(\Delta Y\)的变化量了。
当然没必要先把坐标\(x,y\)除以\(2\)。最后把答案除以\(2\)即可。

我否认在刷水题的事实。

//3072kb    264ms
#include <cstdio>
#include <cctype>
#include <algorithm>
//#define gc() getchar()
#define MAXIN 300000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
typedef long long LL;
const int N=1e5+5;

LL Ans[N];
char IN[MAXIN],*SS=IN,*TT=IN;
struct Point
{
    int x,y,id;
    Point(int tid=0,int ty=0,int tx=0) {id=tid, x=tx+ty, y=tx-ty;}
}p[N];

inline int read()
{
    int now=0,f=1;register char c=gc();
    for(;!isdigit(c);c=='-'&&(f=-1),c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now*f;
}
bool cmpx(Point a,Point b)
{
    return a.x<b.x;
}
bool cmpy(Point a,Point b)
{
    return a.y<b.y;
}

int main()
{
    int n=read();
    for(int i=1; i<=n; ++i) p[i]=Point(i,read(),read());

    std::sort(p+1,p+1+n,cmpx);
    LL sum=-1ll*(n-1)*p[1].x;
    for(int i=2; i<=n; ++i) sum+=p[i].x;
    for(int i=1; i<=n; ++i)
        Ans[p[i].id]+=sum, sum+=(2ll*i-n)*(p[i+1].x-p[i].x);//i*dx-(n-i)*dx

    std::sort(p+1,p+1+n,cmpy);
    sum=-1ll*(n-1)*p[1].y;
    for(int i=2; i<=n; ++i) sum+=p[i].y;
    for(int i=1; i<=n; ++i)
        Ans[p[i].id]+=sum, sum+=(2ll*i-n)*(p[i+1].y-p[i].y);

    LL ans=1e18;
    for(int i=1; i<=n; ++i) ans=std::min(ans,Ans[i]);
    printf("%lld\n",ans>>1);

    return 0;
}

猜你喜欢

转载自www.cnblogs.com/SovietPower/p/9822449.html