这道好(du)题(liu)还是很不错的 挺锻炼代码能力和不断优化 卡常的能力的。
对于 每次询问 我都可以将其分出方向 然后 写 也就是针对于4个方向 左下 左上 右下 右上
这样的话 就成功转换了问题 求4次 三维偏序即可 水题啊。
然后 打完代码 就提交 T飞了
//#include<bits/stdc++.h> #include<iostream> #include<iomanip> #include<ctime> #include<cstdio> #include<cstring> #include<string> #include<cmath> #include<cstdlib> #include<cctype> #include<utility> #include<queue> #include<deque> #include<vector> #include<stack> #include<algorithm> #include<set> #include<bitset> #include<map> #define INF 2147483646 #define ll long long #define ldb long double #define x(i) t[i].x #define y(i) t[i].y #define k(i) t[i].k #define t(i) t[i].t #define id(i) t[i].id using namespace std; char buf[1<<15],*ft,*fs; inline char getc() { return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++; } inline int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } inline void put(int x) { x<0?x=-x,putchar('-'):0; int num=0;char ch[50]; while(x)ch[++num]=x%10+'0',x/=10; num==0?putchar('0'):0; while(num)putchar(ch[num--]); putchar(10);return; } const int MAXN=300002,maxn=1000001; int n,m,tot,num,cnt; struct wy { int x,y; int t,k; int id; }t[MAXN<<1],tmp[MAXN<<1]; int ans[MAXN],c[maxn+10];//树状数组维护前缀最大值！ //注意时间戳也得排序 我的思维漏洞！ int cmp(wy x,wy y)//左下 { if(x.x==y.x&&x.y==x.y)return x.t<y.t; return x.x==y.x?x.y<y.y:x.x<y.x; } int cmp1(wy x,wy y)//左上 { if(x.x==y.x&&x.y==x.y)return x.t<y.t; return x.x==y.x?x.y>y.y:x.x<y.x; } int cmp2(wy x,wy y)//右下 { if(x.x==y.x&&x.y==x.y)return x.t<y.t; return x.x==y.x?x.y<y.y:x.x>y.x; } int cmp3(wy x,wy y)//右上 { if(x.x==y.x&&x.y==x.y)return x.t<y.t; return x.x==y.x?x.y>y.y:x.x>y.x; } inline int max(int x,int y){return x>y?x:y;} void add(int x,int y) { if(y==INF){for(;x<=maxn;x+=x&(-x))c[x]=-INF;return;} for(;x<=maxn;x+=x&(-x))c[x]=max(c[x],y); return; } inline int ask(int x) { int sum=-INF; for(;x;x-=x&(-x))sum=max(sum,c[x]); return sum; } void CDQ(int l,int r,int p) { if(l==r)return; int mid=(l+r)>>1; CDQ(l,mid,p); CDQ(mid+1,r,p); int i=l,j=mid+1; if(p==0) { for(int k=l;k<=r;++k) { if(j>r||(i<=mid&&t(i)<=t(j))) { if(k(i)!=1)add(y(i),x(i)+y(i)); tmp[k]=t[i]; ++i; } else { if(k(j)==0){tmp[k]=t[j];++j;continue;} //if(id(j)==1)cout<<x(j)<<' '<<y(j)<<endl; int maxx=ask(y(j)); if(maxx==-INF){tmp[k]=t[j];++j;continue;} ans[id(j)]=min(ans[id(j)],y(j)+x(j)-maxx); tmp[k]=t[j]; ++j; } } for(int k=l;k<=r;k++) { if(k<=mid&&k(k)!=1)add(y(k),INF); t[k]=tmp[k]; } } else { for(int k=l;k<=r;++k) { if(j>r||(i<=mid&&t(i)<=t(j))) { if(k(i)!=1)add(maxn-y(i),x(i)-y(i)); tmp[k]=t[i]; ++i; } else { if(k(j)==0){tmp[k]=t[j];++j;continue;} int maxx=ask(maxn-y(j)); if(maxx==-INF){tmp[k]=t[j];++j;continue;} //cout<<maxx<<endl; ans[id(j)]=min(ans[id(j)],x(j)-y(j)-maxx); tmp[k]=t[j]; ++j; } } for(int k=l;k<=r;k++) { if(k<=mid&&k(k)!=1)add(maxn-y(k),INF); t[k]=tmp[k]; } } return; } void CDQ1(int l,int r,int p) { if(l==r)return; int mid=(l+r)>>1; CDQ1(l,mid,p); CDQ1(mid+1,r,p); int i=l,j=mid+1; if(p==0) { for(int k=l;k<=r;++k) { if(j>r||(i<=mid&&t(i)<=t(j))) { if(k(i)!=1)add(y(i),y(i)-x(i)); tmp[k]=t[i]; ++i; } else { if(k(j)==0){tmp[k]=t[j];++j;continue;} //if(id(j)==1)cout<<x(j)<<' '<<y(j)<<endl; int maxx=ask(y(j)); if(maxx==-INF){tmp[k]=t[j];++j;continue;} ans[id(j)]=min(ans[id(j)],y(j)-x(j)-maxx); tmp[k]=t[j]; ++j; } } for(int k=l;k<=r;k++) { if(k<=mid&&k(k)!=1)add(y(k),INF); t[k]=tmp[k]; } } else { for(int k=l;k<=r;++k) { if(j>r||(i<=mid&&t(i)<=t(j))) { if(k(i)!=1)add(maxn-y(i),-x(i)-y(i));//此处不要想当然的写 tmp[k]=t[i]; ++i; } else { if(k(j)==0){tmp[k]=t[j];++j;continue;} int maxx=ask(maxn-y(j)); if(maxx==-INF){tmp[k]=t[j];++j;continue;} ans[id(j)]=min(ans[id(j)],-maxx-x(j)-y(j));//变形要正确啊！ tmp[k]=t[j]; ++j; } } for(int k=l;k<=r;k++) { if(k<=mid&&k(k)!=1)add(maxn-y(k),INF); t[k]=tmp[k]; } } return; } int main() { //freopen("1.in","r",stdin); //freopen("1.out","w",stdout); memset(ans,10,sizeof(ans)); n=read();m=read();num=n; for(int i=1;i<=n;i++)x(i)=read(),y(i)=read(); for(int i=1;i<=m;i++) { int p=read(); x(++num)=read();y(num)=read(); if(p==2)k(num)=1,id(num)=++cnt; t(num)=++tot; } //for(int i=1;i<=num;i++)cout<<x(i)<<' '<<y(i)<<' '<<t(i)<<' '<<k(i)<<endl; for(int i=1;i<=maxn;i++)c[i]=-INF; sort(t+1,t+1+num,cmp);//左下 CDQ(1,num,0);//左下 //for(int i=1;i<=cnt;i++)cout<<ans[i]<<endl; for(int i=1;i<=maxn;i++)c[i]=-INF; sort(t+1,t+1+num,cmp1);//左上 CDQ(1,num,1);//左上 //for(int i=1;i<=cnt;i++)cout<<ans[i]<<endl; for(int i=1;i<=maxn;i++)c[i]=-INF; sort(t+1,t+1+num,cmp2);//右下 CDQ1(1,num,0);//右下 //for(int i=1;i<=cnt;i++)cout<<ans[i]<<endl; for(int i=1;i<=maxn;i++)c[i]=-INF; sort(t+1,t+1+num,cmp3);//右上 CDQ1(1,num,1);//右上 //for(int i=1;i<=cnt;i++)cout<<ans[i]<<endl; for(int i=1;i<=cnt;i++)put(ans[i]); return 0; }
都别管我我要写4次CDQ 我要sort 4次然后不出意料只拿到50分。
原因是sort的次数多了 常数是别人的近乎4倍 （有的人一次sort都不做 按照时间直接排）
然后脑残的写了这么段代码 一直T 照着书上的 标程改了一下是一种比较简洁的做法。
但是其CDQ的每层都需要sort 一下比较浪费时间但是还是能拿到60分的 开O2 能A呢。
然后学了一下书上的做法 发现真心简单 代码复杂度 时间复杂度 空间复杂度包括常数都远远低于我
// luogu-judger-enable-o2 //#include<bits/stdc++.h> #include<iostream> #include<iomanip> #include<ctime> #include<cstdio> #include<cstring> #include<string> #include<cmath> #include<cstdlib> #include<cctype> #include<utility> #include<queue> #include<deque> #include<vector> #include<stack> #include<algorithm> #include<set> #include<bitset> #include<map> #define INF 2147483646 #define ll long long #define ldb long double #define x(i) t[i].x #define y(i) t[i].y #define t(i) t[i].t #define id(i) t[i].id using namespace std; char buf[1<<15],*ft,*fs; inline char getc() { return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++; } inline int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } inline void put(int x) { x<0?x=-x,putchar('-'):0; int num=0;char ch[50]; while(x)ch[++num]=x%10+'0',x/=10; num==0?putchar('0'):0; while(num)putchar(ch[num--]); putchar(10);return; } const int MAXN=600002,maxn=1000010; int n,m,cnt,num,tot,h; int ans[MAXN],c[maxn],maxx; struct wy { int x,y; int t; int id; friend int operator <(const wy &a,const wy &b) { return a.x==b.x?a.y<b.y:a.x<b.x;//此时排序不用时间戳(经过实验 //考虑一次CDQ 分治 时间复杂度 ≈(n+m)*log(n+m)^3左右 } }t[MAXN],tmp[MAXN]; inline int max(int x,int y){return x>y?x:y;} inline int min(int x,int y){return x>y?y:x;} inline void add(int x,int y){for(;x<=maxx;x+=x&(-x))c[x]=max(c[x],y);} inline int abs(int x){return x<0?-x:x;} inline int ask(int x) { int maxx1=-INF; for(;x;x-=x&(-x))maxx1=max(maxx1,c[x]); return maxx1; } inline void remove1(int x) { for(;x<=maxx;x+=x&(-x))c[x]=-INF; return; } void calculate(int st,int en,int w,int kx,int ky)//属性 { for(int i=st;i!=en;i+=w) { int yy=(ky==1?tmp[i].y:maxx-tmp[i].y); int sum=kx*tmp[i].x+ky*tmp[i].y; if(tmp[i].t==1)add(yy,sum); else { int an=ask(yy); if(an==-INF)continue; ans[tmp[i].id]=min(ans[tmp[i].id],sum-an); } //if(tmp[i].id==2)cout<<sum<<' '<<ask(yy)<<endl; } for(int i=st;i!=en;i+=w)if(tmp[i].t==1)remove1(ky==1?tmp[i].y:maxx-tmp[i].y); } void CDQ(int l,int r) { int mid=(l+r)>>1; if(l<mid)CDQ(l,mid); if(mid+1<r)CDQ(mid+1,r); h=0; for(int i=l;i<=r;i++)if((i<=mid&&t[i].t==1)||(i>mid&&t[i].t==2))tmp[++h]=t[i]; sort(tmp+1,tmp+1+h); calculate(1,h+1,1,1,1);//左下 calculate(1,h+1,1,1,-1);//左上 calculate(h,0,-1,-1,1);//右下 calculate(h,0,-1,-1,-1);//右上 } int main() { //freopen("1.in","r",stdin); //freopen("1.out","w",stdout); n=read();m=read();num=n; for(int i=1;i<=n;++i) { x(i)=read(); y(i)=read(); t(i)=1; maxx=max(maxx,y(i)); } for(int i=1;i<=m;++i) { t(++num)=read(); x(num)=read(); y(num)=read(); id(num)=t(num)==2?++cnt:0; maxx=max(maxx,y(num)); } //put(cnt); ++maxx;//巧妙一招 for(int i=1;i<=maxx;++i)c[i]=-INF; memset(ans,10,sizeof(ans)); CDQ(1,num); for(int i=1;i<=cnt;i++)put(ans[i]); return 0; }
比较难受 因为这浪费了我很多的时间 然后最后还是T 开了O2 死T一个点（原因是 我是个NC啊。。。
开O2 过了不算什么我要不开O2 过
点开题解发现我的第一个代码的思想及其接近其 题解中跑的最快的。
所以发现可以不用sort 于是顺利不开O2 第一个代码顺利得到80分。
发现可以加一个优化是 能不加到树状数组里就不加 到需要的时候再加 还有一些卡常的东西
最终不开O2 直接A （一直T的那个点是因为下标为0 树状数组原地GG了。。。）
//#include<bits/stdc++.h> #include<iostream> #include<iomanip> #include<ctime> #include<cstdio> #include<cstring> #include<string> #include<cmath> #include<cstdlib> #include<cctype> #include<utility> #include<queue> #include<deque> #include<vector> #include<stack> #include<algorithm> #include<set> #include<bitset> #include<map> #define INF 2147483646 #define ll long long #define ldb long double #define x(i) t[i].x #define y(i) t[i].y #define k(i) t[i].k #define t(i) t[i].t #define id(i) t[i].id #define R register using namespace std; //char buf[1<<15],*ft,*fs; //inline char getc(){return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;} inline int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } inline void put(int x) { x<0?x=-x,putchar('-'):0; int num=0;char ch[50]; while(x)ch[++num]=x%10+'0',x/=10; num==0?putchar('0'):0; while(num)putchar(ch[num--]); putchar(10);return; } const int maxn=1000010; int n,m,tot,num,cnt,maxx1,maxx2; struct wy { int x,y; int k,id; inline bool operator <(const wy &w)const { return x<w.x; } }t[maxn],tmp[maxn],a[maxn]; int ans[maxn],c[maxn+10]; inline int max(int x,int y){return x>y?x:y;} inline void add(int x,int y) { for(;x<=maxx1;x+=x&(-x)) { if(c[x]>y)break; else c[x]=y; } return; } inline void remove1(int x) { for(;x<=maxx1;x+=x&(-x))if(c[x])c[x]=0; return; } inline int ask(int x) { int sum=-INF; for(;x;x-=x&(-x))sum=max(sum,c[x]); return sum; } inline void CDQ(int l,int r) { int mid=(l+r)>>1; if(l<mid)CDQ(l,mid); if(r>mid+1)CDQ(mid+1,r); R int i=l,j=mid+1; for(j=mid+1;j<=r;j++) { if(k(j)==1) { for(;i<=mid&&x(i)<=x(j);i++) if(k(i)==0)add(y(i),x(i)+y(i)); int maxx=ask(y(j)); maxx==0?0:ans[id(j)]=min(ans[id(j)],y(j)+x(j)-maxx); } } for(j=l;j<i;j++)if(k(j)==0)remove1(y(j)); i=l,j=mid+1; for(R int k=l;k<=r;++k) { if(j>r||(i<=mid&&t[i]<t[j]))tmp[k]=t[i],++i; else tmp[k]=t[j],++j; } for(R int k=l;k<=r;++k)t[k]=tmp[k]; return; } int main() { //freopen("1.in","r",stdin); //freopen("1.out","w",stdout); n=read();m=read();num=n; for(R int i=1;i<=n;++i){x(i)=read()+1,y(i)=read()+1;maxx1=max(maxx1,y(i));maxx2=max(maxx2,x(i));} for(R int i=1;i<=m;++i) { int p=read(); x(++num)=read()+1;y(num)=read()+1; if(p==2)k(num)=1,id(num)=++cnt; maxx1=max(maxx1,y(num)); maxx2=max(maxx2,x(num)); } ++maxx1;++maxx2; for(R int i=1;i<=cnt;++i)ans[i]=INF; for(R int i=1;i<=num;++i)a[i]=t[i]; CDQ(1,num);//左下 for(R int i=1;i<=num;++i)t[i]=a[i],y(i)=maxx1-y(i); CDQ(1,num);//左上 for(R int i=1;i<=num;++i)t[i]=a[i],x(i)=maxx2-x(i); CDQ(1,num);//右下 for(R int i=1;i<=num;++i)t[i]=a[i],x(i)=maxx2-x(i),y(i)=maxx1-y(i); CDQ(1,num);//右上 for(R int i=1;i<=cnt;++i)put(ans[i]); return 0; }
当然 为了能够完美的A掉这道题 我决定去学一个更强的东西 KD-tree.
然后在抄了20min代码后觉定放弃 尧神说没多大卵用。。。
所以这篇博文到此啦，关键是不断优化的思想 少开O2