【洛谷P4169】天使玩偶/SJY摆棋子【CDQ分治】
生活随笔
收集整理的這篇文章主要介紹了
【洛谷P4169】天使玩偶/SJY摆棋子【CDQ分治】
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
傳送門
題意:動態加點,給定點詢問曼哈頓距離最近的點
N,M≤3e5,x,y≤1e6N,M \leq 3e5,x,y \leq 1e6N,M≤3e5,x,y≤1e6
經(kan)過(le)分(ti)析(jie),這是一道cdqcdqcdq分治
考慮當前區間左半邊修改對右半邊的詢問的影響
設左邊某個修改為(x1,y1)(x_1,y_1)(x1?,y1?),右邊的某個詢問為(x2,y2)(x_2,y_2)(x2?,y2?)
考慮x1≤x2,y1≤y2x_1 \leq x_2,y_1 \leq y_2x1?≤x2?,y1?≤y2?的情況,答案為x2+y2?x1?y1x_2+y_2-x_1-y_1x2?+y2??x1??y1?
因為坐標是線性級別,所以需要排序搞掉一維
左右分別以xxx坐標排序
這樣可以用雙指針搞掉xxx的限制
yyy坐標開個樹狀數組,記錄不超過yyy的修改中最大的x+yx+yx+y
然后很容易計算答案
其他三種情況類似
為了實現方便,直接旋轉坐標系,即用infinfinf減
樹狀數組已經有了logloglog,我們已經無所畏懼,所以可以直接sortsortsort
復雜度O(NlogN2)O(Nlog_N^2)O(NlogN2?)
#include <iostream> #include <cstdio> #include <cstring> #include <cctype> #include <algorithm> #define MAXN 600005 #define MAXM 2000005 #define MAX 1000000 using namespace std; inline int read() {int ans=0;char c=getchar();while (!isdigit(c)) c=getchar();while (isdigit(c)) ans=(ans<<3)+(ans<<1)+(c^48),c=getchar();return ans; } int n,m,siz; struct BIT {int s[MAXM];inline int lowbit(const int& x){return x&-x;}inline void modify(int x,const int& v){for (;x<=(MAX<<1);s[x]=max(s[x],v),x+=lowbit(x));}inline int query(int x){int ans=-0x3f3f3f3f;for (;x;ans=max(ans,s[x]),x-=lowbit(x));return ans;}inline void clear(int x){for (;x<=(MAX<<1);s[x]=-0x3f3f3f3f,x+=lowbit(x));} }bit; int ans[MAXN]; struct query{int type,x,y,pos;}q[MAXN]; inline bool operator <(const query& a,const query& b){if (a.x==b.x) return a.y<b.y;return a.x<b.x;} void calc(int l,int r) {int now=l-1,mid=(l+r)>>1;for (int i=mid+1;i<=r;i++){if (q[i].type==1) continue;while (now<mid&&q[now+1].x<=q[i].x) {++now;if (q[now].type==1) bit.modify(q[now].y,q[now].x+q[now].y);}ans[q[i].pos]=min(ans[q[i].pos],q[i].x+q[i].y-bit.query(q[i].y));}for (int i=l;i<=mid;i++) bit.clear(q[i].y); } void cdq(int l,int r) {if (l==r) return;int mid=(l+r)>>1;cdq(l,mid);cdq(mid+1,r);calc(l,r);for (int i=l;i<=r;i++) q[i].x=MAX-q[i].x;sort(q+l,q+mid+1),sort(q+mid+1,q+r+1);calc(l,r);for (int i=l;i<=r;i++) q[i].y=MAX-q[i].y;sort(q+l,q+mid+1),sort(q+mid+1,q+r+1);calc(l,r);for (int i=l;i<=r;i++) q[i].x=MAX-q[i].x;sort(q+l,q+mid+1),sort(q+mid+1,q+r+1);calc(l,r);for (int i=l;i<=r;i++) q[i].y=MAX-q[i].y;sort(q+l,q+r+1); } int main() {n=read(),m=read();for (int i=1;i<=n;i++) q[i].type=1,q[i].x=read()+1,q[i].y=read()+1;for (int i=1;i<=m;i++) q[n+i].type=read(),q[n+i].x=read()+1,q[n+i].y=read()+1,q[n+i].pos=i;memset(ans,0x3f,sizeof(ans));for (int i=0;i<MAXM;i++) bit.s[i]=-0x3f3f3f3f;cdq(1,n+m);for (int i=1;i<=m;i++)if (ans[i]<0x3f3f3f3f)printf("%d\n",ans[i]);return 0; }總結
以上是生活随笔為你收集整理的【洛谷P4169】天使玩偶/SJY摆棋子【CDQ分治】的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 圣剑契约电脑版攻略
- 下一篇: lol2018s9英雄联盟放逐之刃 锐雯