码迷,mamicode.com
首页 > 其他好文 > 详细

bzoj2716: [Violet 3]天使玩偶

时间:2016-05-20 13:25:17      阅读:259      评论:0      收藏:0      [点我收藏+]

标签:

题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2716

由于距离定义为曼哈顿距离,我们先将绝对值符号去掉。这时可以分为四种情况,但是我们可以对点的横纵坐标进行一些处理(处理方法可以看程序),可以使四种情况都变为一种:Ax>Bx,Ay>By。这时距离表示为dist(A,B)=(Ax+Ay)-(Bx+By)。对于一个点,我们只需要去找位于它左下方的x+y最大的点就行了。那么问题就变成了一道简单的CDQ分治题。

不过像我这种渣渣,TLE了好几遍,把读入优化和一些小优化加进去才过,哎,太弱了。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #include<algorithm>
 5 #define inf 4000000
 6 #define maxn 1000010
 7 using namespace std;
 8 int n,m,cnt,tot,pos[maxn],ans[maxn],t[maxn],maxx,maxy;
 9 struct fuck{int x,y,id,op;}e[maxn],f[maxn];
10 int read(){
11     int x=0,f=1; char ch;
12     for(ch=getchar();ch<0||ch>9;ch=getchar()) if(ch==-) f=-1;
13     for(;ch>=0&&ch<=9;ch=getchar()) x=x*10+ch-0;
14     return x*f;
15 }
16 void change(int x,int y){for(int i=x;i<=maxy;i+=i&-i) t[i]=max(t[i],y);}
17 int query(int x){int y=0; for(int i=x;i;i-=i&-i) y=max(y,t[i]); return y==0?-inf:y;}
18 bool comp(fuck a,fuck b){return a.x<b.x;}
19 void solve(int l,int r){
20     if(l==r) return;
21     int mid=(l+r)>>1;
22     solve(l,mid); solve(mid+1,r);
23     sort(f+l,f+mid+1,comp); sort(f+mid+1,f+r+1,comp);
24     int i=l,j=mid+1,last=0;
25     while(j<=r){
26         while(i<=mid&&f[i].op==2) i++;
27         while(j<=r&&f[j].op==1) j++;
28         if(i<=mid&&f[i].x<=f[j].x) change(f[i].y,f[i].x+f[i].y),last=i++;
29         else if(j<=r) ans[f[j].id]=min(ans[f[j].id],f[j].x+f[j].y-query(f[j].y)),j++;
30     }
31     for(int i=l;i<=last;i++){
32         if(f[i].op==1)
33             for(int j=f[i].y;j<=maxy;j+=j&-j) t[j]=0;
34     }
35 }
36 int main(){
37     n=read(); m=read();
38     for(int i=1;i<=n+m;i++) ans[i]=inf;
39     for(int i=1;i<=n;i++){
40         e[i].op=1; e[i].id=i; e[i].x=read()+1; e[i].y=read()+1;
41         maxx=max(maxx,e[i].x); maxy=max(maxy,e[i].y);
42     }
43     int type,x,y;
44     for(int i=n+1;i<=n+m;i++){
45         e[i].op=read(); e[i].id=i; e[i].x=read()+1; e[i].y=read()+1;
46         maxx=max(maxx,e[i].x); maxy=max(maxy,e[i].y);
47         if(e[i].op==2) pos[++cnt]=i;
48     }
49     maxx++; maxy++; tot=n+m;
50     for(int i=1;i<=tot;i++) f[i]=e[i]; solve(1,tot);
51     for(int i=1;i<=tot;i++) f[i]=e[i],f[i].x=maxx-f[i].x; solve(1,tot);
52     for(int i=1;i<=tot;i++) f[i]=e[i],f[i].y=maxy-f[i].y; solve(1,tot);
53     for(int i=1;i<=tot;i++) f[i]=e[i],f[i].x=maxx-f[i].x,f[i].y=maxy-f[i].y; solve(1,tot);
54     for(int i=1;i<=cnt;i++) printf("%d\n",ans[pos[i]]);
55     return 0;
56 }

 

bzoj2716: [Violet 3]天使玩偶

标签:

原文地址:http://www.cnblogs.com/longshengblog/p/5511593.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!