1.$B>0$

二分答案，然后掃描線，線段樹維護某個點作為左下角時的費用的最小值，支持區間加。

時間復雜度$O(n\log^2n)$。

?

2.$B=0$

枚舉左邊界，則最優右邊界可以通過雙指針求出。

用線段樹維護左右邊界之間最長的豎著的空的連續段的長度。

找到最大的連續段長度$\geq$左右邊界距離的位置，此時的距離就是答案。

時間復雜度$O(n\log n)$。

?

#include<cstdio> #include<algorithm> using namespace std; int n,m,B,P; inline void read(int&a){char c;while(!(((c=getchar())>='0')&&(c<='9')));a=c-'0';while(((c=getchar())>='0')&&(c<='9'))(a*=10)+=c-'0';} namespace BinarySearch{ const int N=30010,M=2100000; int i,j,l,r,mid,ans,ce,v[M],tag[M]; struct O{int a,b,c,d,p;}a[N]; struct E{int x,l,r,p;E(){}E(int _x,int _l,int _r,int _p){x=_x,l=_l,r=_r,p=_p;}}e[N*2]; inline bool cmp(const E&a,const E&b){return a.x<b.x;} void build(int x,int a,int b){v[x]=tag[x]=0;if(a==b)return;int mid=(a+b)>>1;build(x<<1,a,mid),build(x<<1|1,mid+1,b); } inline void add1(int x,int p){v[x]+=p;tag[x]+=p;} void add(int x,int a,int b,int c,int d,int p){if(c<=a&&b<=d){add1(x,p);return;}if(tag[x])add1(x<<1,tag[x]),add1(x<<1|1,tag[x]),tag[x]=0;int mid=(a+b)>>1;if(c<=mid)add(x<<1,a,mid,c,d,p);if(d>mid)add(x<<1|1,mid+1,b,c,d,p);v[x]=min(v[x<<1],v[x<<1|1]); } inline bool check(int mid){for(ce=0,i=1;i<=P;i++){int A=a[i].a-mid+1,B=a[i].b-mid+1,C=a[i].c,D=a[i].d;A=max(A,1),B=max(B,1),C=min(C,n-mid+1),D=min(D,m-mid+1);if(A<=C&&B<=D)e[++ce]=E(A,B,D,a[i].p),e[++ce]=E(C+1,B,D,-a[i].p);}sort(e+1,e+ce+1,cmp);build(1,1,m-mid+1);for(i=j=1;i<=n-mid+1;i++){while(j<=ce&&e[j].x<=i)add(1,1,m-mid+1,e[j].l,e[j].r,e[j].p),j++;if(v[1]<=B)return 1;}return 0; } void solve(){for(i=1;i<=P;i++)read(a[i].a),read(a[i].b),read(a[i].c),read(a[i].d),read(a[i].p);l=1,r=min(n,m);while(l<=r)if(check(mid=(l+r)>>1))l=(ans=mid)+1;else r=mid-1;printf("%d",ans); } } namespace Sweep{ const int N=1000010,M=2100000; int A,B,C,D,i,j,k,ans,gi[N],gd[N],v[N],w[N],nxt[N],ed; int tag[M],len[M],vl[M],vr[M],vm[M]; inline void add(int&x,int y,int z){v[++ed]=y;w[ed]=z;nxt[ed]=x;x=ed;} void build(int x,int a,int b){len[x]=vl[x]=vr[x]=vm[x]=b-a+1;if(a==b)return;int mid=(a+b)>>1;build(x<<1,a,mid),build(x<<1|1,mid+1,b); } inline void up(int x){if(tag[x]){vl[x]=vr[x]=vm[x]=0;return;}if(len[x]==1){vl[x]=vr[x]=vm[x]=1;return;}vl[x]=vl[x<<1];if(vl[x<<1]==len[x<<1])vl[x]+=vl[x<<1|1];vr[x]=vr[x<<1|1];if(vr[x<<1|1]==len[x<<1|1])vr[x]+=vr[x<<1];vm[x]=max(max(vm[x<<1],vm[x<<1|1]),vr[x<<1]+vl[x<<1|1]); } void add(int x,int a,int b,int c,int d,int p){if(c<=a&&b<=d){tag[x]+=p;up(x);return;}int mid=(a+b)>>1;if(c<=mid)add(x<<1,a,mid,c,d,p);if(d>mid)add(x<<1|1,mid+1,b,c,d,p);up(x); } void solve(){while(P--)read(A),read(B),read(C),read(D),read(k),add(gi[A],B,D),add(gd[C],B,D);build(1,1,m);for(i=1;i<=n;i++){for(k=gd[i-1];k;k=nxt[k])add(1,1,m,v[k],w[k],-1);while(vm[1]>=j-i+1&&j<=n)for(k=gi[++j];k;k=nxt[k])add(1,1,m,v[k],w[k],1);ans=max(ans,j-i);}printf("%d",ans); } } int main(){read(n),read(m),read(B),read(P);if(B)BinarySearch::solve();else Sweep::solve();return 0; }

　　

轉載于:https://www.cnblogs.com/clrs97/p/5648487.html

總結

以上是生活随笔為你收集整理的BZOJ1795 : [Ioi2008]Pyramid Base 金字塔地基的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。