[Tjoi2013]松鼠聚會
阿新 • • 發佈:2018-11-22
putc ostream ans cpp getch 兩個 insert void class
Description
有N個小松鼠,它們的家用一個點x,y表示,兩個點的距離定義為:點(x,y)和它周圍的8個點即上下左右四個點和對角的四個點,距離為1。現在N個松鼠要走到一個松鼠家去,求走過的最短距離。
Input
第一行給出數字N,表示有多少只小松鼠。0<=N<=10^5
下面N行,每行給出x,y表示其家的坐標。
-10^9<=x,y<=10^9
Output
表示為了聚會走的路程和最小為多少。
Sample Input
6
-4 -1
-1 -2
2 -4
0 2
0 3
5 -2
Sample Output
20
如果是曼哈頓距離十分好求,我們可以分開考慮
將x排序,利用前綴後綴和,算出每個點在x軸方向到其他點的距離,將答案記錄下來,再按y排序,即可統計答案
但是這題並不是曼哈頓距離,而是切比雪夫距離,怎麽辦?
其實有個結論,將每個點的坐標改為\((\frac{x+y}{2},\frac{x-y}{2})\)後,兩點之間的曼哈頓距離等於切比雪夫距離
證明的話可以自己手推,記得考慮大小關系(其實是我懶了)
這種結論題。。。不知道結論根本不會寫吧。。。至少我是沒有當場推結論的水平。。。
/*program from Wolfycz*/ #include<cmath> #include<cstdio> #include<cstring> #include<iostream> #include<algorithm> #define inf 0x7f7f7f7f using namespace std; typedef long long ll; typedef unsigned int ui; typedef unsigned long long ull; inline char gc(){ static char buf[1000000],*p1=buf,*p2=buf; return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++; } inline int frd(){ int x=0,f=1; char ch=gc(); for (;ch<'0'||ch>'9';ch=gc()) if (ch=='-') f=-1; for (;ch>='0'&&ch<='9';ch=gc()) x=(x<<3)+(x<<1)+ch-'0'; return x*f; } inline int read(){ int x=0,f=1; char ch=getchar(); for (;ch<'0'||ch>'9';ch=getchar()) if (ch=='-') f=-1; for (;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0'; return x*f; } inline void print(int x){ if (x<0) putchar('-'),x=-x; if (x>9) print(x/10); putchar(x%10+'0'); } const int N=1e5; struct S1{ int x,y,ID; void insert(int _x,int _y,int _ID){x=_x,y=_y,ID=_ID;} }A[N+10]; bool cmpx(const S1 &x,const S1 &y){return x.x<y.x;} bool cmpy(const S1 &x,const S1 &y){return x.y<y.y;} ll pre[N+10],suf[N+10],v[N+10]; int main(){ int n=read(); ll Ans=1e18; for (int i=1;i<=n;i++){ int x=read(),y=read(); A[i].insert(x+y,x-y,i); } sort(A+1,A+1+n,cmpx); for (int i=1;i<=n;i++) pre[i]=pre[i-1]+A[i].x; for (int i=n;i>=1;i--) suf[i]=suf[i+1]+A[i].x; for (int i=1;i<=n;i++) v[A[i].ID]=(1ll*i*A[i].x-pre[i])+(suf[i]-1ll*(n-i+1)*A[i].x); sort(A+1,A+1+n,cmpy); for (int i=1;i<=n;i++) pre[i]=pre[i-1]+A[i].y; for (int i=n;i>=1;i--) suf[i]=suf[i+1]+A[i].y; for (int i=1;i<=n;i++) Ans=min(Ans,(1ll*i*A[i].y-pre[i])+(suf[i]-1ll*(n-i+1)*A[i].y)+v[A[i].ID]); printf("%lld\n",Ans>>1); return 0; }
[Tjoi2013]松鼠聚會