• bzoj3110[Zjoi2013]K大数查询


    刚刚学习了整体二分,练习一下。

    太长时间没写线段树了,结果在这道题里,线段树没有写update函数,怎么调都调不出来,无奈只好翻别人程序,改了一下就过了。

    还有一个非常麻烦的事情是边界问题,程序里有一些需要打小于,有一些打小于等于,解决这个问题一种方法是对程序逻辑良好的整体把握,另一种方法是大力对拍,把可能情况测一遍就行了。

    样例太水,不要相信样例。

    #include<cstdio>
    #include<algorithm>
    using namespace std;
    #define LL  long long
    #define FILE "dealing"
    #define up(i,j,n) for(int i=j;i<=n;i++)
    LL read(){
    	LL x=0,f=1,ch=getchar();
    	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    	while(ch>='0'&&ch<='9')x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
    	return x*f;
    }
    const LL maxn=800000,mod=1000000007,inf=10000000000000LL;
    bool cmin(LL& a,LL b){return a>b?a=b,true:false;}
    bool cmax(LL& a,LL b){return a<b?a=b,true:false;}
    LL n,m;
    struct node{LL x,y,op,ans,k,id;}a[maxn],q1[maxn],q2[maxn];
    bool cmp(node a,node b){return a.id<b.id;}
    LL key,left,right;
    LL delet[maxn],s[maxn],siz[maxn],flag[maxn];
    void add(LL x,LL del){delet[x]+=del,s[x]+=siz[x]*del;}
    void qing(LL x){delet[x]=s[x]=0;flag[x]=1;}
    void updata(LL x){s[x]=s[x<<1]+s[x<<1|1];}
    void pushdown(LL x){
    	if(flag[x])qing(x<<1),qing(x<<1|1),flag[x]=0;
    	if(delet[x])add(x<<1,delet[x]),add(x<<1|1,delet[x]),delet[x]=0;
    }
    void change(LL l,LL r,LL o){
    	if(l>right||r<left)return;
    	if(l>=left&&r<=right){add(o,key);return;}
    	LL mid=(l+r)>>1;
    	pushdown(o);
    	change(l,mid,o<<1);
    	change(mid+1,r,o<<1|1);
    	updata(o);
    }
    void Change(LL L,LL R,LL d){
    	left=L,right=R,key=d;
    	change(1,n,1);
    }
    LL query(LL l,LL r,LL o){
    	if(r<left||l>right)return 0;
    	if(l>=left&&r<=right)return s[o];
    	LL mid=(l+r)>>1;
    	pushdown(o);
    	return query(l,mid,o<<1)+query(mid+1,r,o<<1|1);
    }
    LL Query(LL L,LL R){
    	left=L,right=R;
    	return query(1,n,1);
    }
    void build(LL l,LL r,LL o){
    	if(l==r){siz[o]=1;return;}
    	LL mid=(l+r)>>1;
    	build(l,mid,o<<1);
    	build(mid+1,r,o<<1|1);
    	siz[o]=siz[o<<1]+siz[o<<1|1];
    }
    LL sum[maxn];
    void divide(LL L,LL R,LL l,LL r){
    	if(l>r)return;
    	if(L==R){
    		up(i,l,r)if(a[i].op==2)a[i].ans=L;
    		return;
    	}
    	qing(1);
    	LL mid=(L+R)>>1;
    	up(i,l,r){
    		if(a[i].op==1&&a[i].k>mid)Change(a[i].x,a[i].y,1);
    		else if(a[i].op==2)sum[i]=Query(a[i].x,a[i].y);
    	}
    	LL top1=0,top2=0;
    	up(i,l,r){
    		if(a[i].op==2){
    			if(sum[i]<=a[i].k){
    				a[i].k-=sum[i];
    				q1[++top1]=a[i];
    			}else q2[++top2]=a[i];
    		}
    		else {
    			if(a[i].k>mid)q2[++top2]=a[i];
    			else q1[++top1]=a[i];
    		}
    	}
    	up(i,l,l+top1-1)a[i]=q1[i-l+1];
    	up(i,l+top1,r)a[i]=q2[i-top1-l+1];
    	divide(L,mid,l,l+top1-1);
    	divide(mid+1,R,l+top1,r);
    }
    int main(){
    	freopen(FILE".in","r",stdin);
    	freopen(FILE".out","w",stdout);
    	n=read(),m=read();LL Min=inf,Max=-inf;
    	up(i,1,m){
    		a[i].op=read(),a[i].x=read(),a[i].y=read(),a[i].k=read();a[i].id=i;
    		if(a[i].op==2)a[i].k-=1;
    		if(a[i].op==1)cmin(Min,a[i].k),cmax(Max,a[i].k);
    	}
    	build(1,n,1);
    	divide(Min,Max,1,m);
    	sort(a+1,a+m+1,cmp);
    	up(i,1,m)if(a[i].op==2)printf("%lld
    ",a[i].ans);
    	return 0;
    }
    

      

  • 相关阅读:
    无法在WEB服务器上启动调试
    Zedgraph悬停时显示内容闪烁的解决
    用ZedGraph控件作图圆
    34.node.js之Url & QueryString模块
    33.Node.js 文件系统fs
    32.Node.js中的常用工具类util
    31.Node.js 常用工具 util
    30.Node.js 全局对象
    28.Node.js 函数和匿名函数
    27.Node.js模块系统
  • 原文地址:https://www.cnblogs.com/chadinblog/p/6439190.html
Copyright © 2020-2023  润新知