POJ 2104 K-th Number

时空隧道

题意：

给出一个序列，每次查询区间第k小

分析：

整体二分入门题？

代码：

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 //by NeighThorn

 #define inf 0x3f3f3f3f

 using namespace std;

 const int maxn=+,maxm=+;

 int n,m,a[maxn],ans[maxm],tr[maxn];

 struct M{

     int x,y,k,id,flag;

     M(int a=,int b=,int c=,int d=,int e=){

         x=a,y=b,k=c,id=d,flag=e;

     }

 }q[maxm+maxn],q1[maxm+maxn],q2[maxm+maxn];

 inline void add(int x,int y){

     for(;x<=n;x+=x&(-x))

         tr[x]+=y;

 }

 inline int query(int x){

     int res=;

     for(;x;x-=x&(-x))

         res+=tr[x];

     return res;

 }

 inline void solve(int L,int R,int l,int r){

     if(L>R)

         return;

     if(l==r){

         for(int i=L;i<=R;i++)

             if(q[i].flag==)

                 ans[q[i].id]=l;

         return;

     }

     int mid=(l+r)>>,l1=,l2=;

     for(int i=L;i<=R;i++){

         if(q[i].flag==){

             if(q[i].x<=mid)

                 add(q[i].id,),q1[l1++]=q[i];

             else

                 q2[l2++]=q[i];

         }

         else{

             int lala=query(q[i].y)-query(q[i].x-);

             if(lala>=q[i].k)

                 q1[l1++]=q[i];

             else

                 q[i].k-=lala,q2[l2++]=q[i];

         }

     }

     for(int i=;i<l1;i++)

         if(q1[i].flag==)

             add(q1[i].id,-);

     memcpy(q+L,q1,sizeof(q[])*l1);

     memcpy(q+L+l1,q2,sizeof(q[])*l2);

     solve(L,L+l1-,l,mid);solve(L+l1,R,mid+,r);

 }

 signed main(void){

     while(scanf("%d%d",&n,&m)!=EOF){

         memset(tr,,sizeof(tr));

         for(int i=;i<=n;i++){

             scanf("%d",&a[i]);

             q[i]=(M){a[i],,,i,};

         }

         for(int i=;i<=m;i++)

             scanf("%d%d%d",&q[i+n].x,&q[i+n].y,&q[i+n].k),q[i+n].id=i,q[i+n].flag=;

         solve(,n+m,-inf,inf);

         for(int i=;i<=m;i++)

             printf("%d\n",ans[i]);

     }

     return ;

 }

BZOJ 1901: Zju2112 Dynamic Rankings

时空隧道

分析：

和上一题一样，就是把修改操作拆成删除和插入操作...

代码：

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 //by NeighThorn

 #define inf 0x3f3f3f3f

 using namespace std;

 const int maxn=+;

 int n,m,cnt,tot,tr[maxn],pre[maxn],ans[maxn];

 char ch[];

 struct M{

     int x,y,k,id,flag;

 }q[maxn*],q1[maxn*],q2[maxn*];

 inline int read(void){

     char ch=getchar();int f=,x=;

     while(!(ch>=''&&ch<='')){

         if(ch=='-')

             f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

         x=x*+ch-'',ch=getchar();

     return x;

 }

 inline void insert(int x,int y){

     for(;x<=n;x+=x&(-x))

         tr[x]+=y;

 }

 inline int query(int x){

     int res=;

     for(;x;x-=x&(-x))

         res+=tr[x];

     return res;

 }

 inline void solve(int L,int R,int l,int r){

     if(L>R)

         return;

     if(l==r){

         for(int i=L;i<=R;i++)

             if(q[i].flag==)

                 ans[q[i].id]=l;

         return;

     }

     int mid=(l+r)>>,l1=,l2=;

     for(int i=L;i<=R;i++){

         if(q[i].flag==){

             if(q[i].x<=mid)

                 insert(q[i].id,q[i].y),q1[l1++]=q[i];

             else

                 q2[l2++]=q[i];

         }

         else{

             int lala=query(q[i].y)-query(q[i].x-);

             if(lala>=q[i].k)

                 q1[l1++]=q[i];

             else

                 q[i].k-=lala,q2[l2++]=q[i];

         }

     }

     for(int i=;i<l1;i++)

         if(q1[i].flag==)

             insert(q1[i].id,-q1[i].y);

     memcpy(q+L,q1,sizeof(q[])*l1);

     memcpy(q+L+l1,q2,sizeof(q[])*l2);

     solve(L,L+l1-,l,mid);solve(L+l1,R,mid+,r);

 }

 signed main(void){

     n=read(),m=read();tot=n;

     for(int i=;i<=n;i++)

         pre[i]=read(),q[i].x=pre[i],q[i].y=,q[i].k=,q[i].id=i,q[i].flag=;

     for(int i=,x,y;i<=m;i++){

         scanf("%s",ch);tot++;

         if(ch[]=='Q')

             q[tot].x=read(),q[tot].y=read(),q[tot].k=read(),q[tot].id=++cnt,q[tot].flag=;

         else

             x=read(),y=read(),q[tot].x=pre[x],q[tot].y=-,q[tot].k=,q[tot].id=x,q[tot].flag=,

             pre[x]=y,q[++tot].x=y,q[tot].y=,q[tot].k=,q[tot].id=x,q[tot].flag=;

     }

     solve(,tot,,inf);

     for(int i=;i<=cnt;i++)

         printf("%d\n",ans[i]);

     return ;

 }//Cap ou pas cap. Cap.

BZOJ 3110: [Zjoi2013]K大数查询

时空隧道

分析：

还是一样的...只是把树状数组改成了线段树...

然而...虽然题水...但是我WA了半页...快读写错了...没有读入负数

代码：

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 //by NeighThorn

 #define int long long

 using namespace std;

 const int maxn=+;

 int n,m,cnt,ans[maxn];

 struct M{

     int x,y,id,flag;

     long long k;

 }q[maxn],q1[maxn],q2[maxn];

 struct Tree{

     int l,r;

     long long sum,lazy;

 }tree[maxn<<];

 inline long long longread(void){

     char ch=getchar();int f=;long long x=;

     while(!(ch>=''&&ch<='')){

         if(ch=='-')

             f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

         x=x*+ch-'',ch=getchar();

     return f*x;

 }

 inline int intread(void){

     char ch=getchar();int f=,x=;

     while(!(ch>=''&&ch<='')){

         if(ch=='-')

             f=-;

         ch=getchar();

     }

     while(ch>=''&&ch<='')

         x=x*+ch-'',ch=getchar();

     return f*x;

 }

 inline void build(int l,int r,int tr){

     tree[tr].l=l,tree[tr].r=r,tree[tr].sum=,tree[tr].lazy=;

     if(l==r)

         return;

     int mid=(l+r)>>;

     build(l,mid,tr<<),build(mid+,r,tr<<|);

 }

 inline void change(int l,int r,int val,int tr){//cout<<l<<" "<<r<<"   "<<tree[tr].l<<" "<<tree[tr].r<<" "<<tr<<"   ";

     if(tree[tr].l==l&&tree[tr].r==r){

         tree[tr].sum+=val*(r-l+);

         tree[tr].lazy+=val;return;

     }

     int mid=(tree[tr].l+tree[tr].r)>>;//cout<<mid<<endl;

     if(tree[tr].lazy)

         change(tree[tr].l,mid,tree[tr].lazy,tr<<),change(mid+,tree[tr].r,tree[tr].lazy,tr<<|),tree[tr].lazy=;

     if(r<=mid)

         change(l,r,val,tr<<);

     else if(l>mid)

         change(l,r,val,tr<<|);

     else

         change(l,mid,val,tr<<),change(mid+,r,val,tr<<|);

     tree[tr].sum=tree[tr<<].sum+tree[tr<<|].sum;

 }

 inline long long query(int l,int r,int tr){//cout<<l<<" "<<r<<" "<<tree[tr].l<<" "<<tree[tr].r<<" "<<tr<<" ";

     if(tree[tr].l==l&&tree[tr].r==r)

         return tree[tr].sum;

     int mid=(tree[tr].l+tree[tr].r)>>;//cout<<mid<<endl;

     if(tree[tr].lazy)

         change(tree[tr].l,mid,tree[tr].lazy,tr<<),change(mid+,tree[tr].r,tree[tr].lazy,tr<<|),tree[tr].lazy=;

     if(r<=mid)

         return query(l,r,tr<<);

     else if(l>mid)

         return query(l,r,tr<<|);

     else

         return query(l,mid,tr<<)+query(mid+,r,tr<<|);

 }

 inline void solve(int L,int R,int l,int r){//cout<<L<<" "<<R<<" "<<l<<" "<<r<<endl;

     if(L>R)

         return;

     if(l==r){

         for(int i=L;i<=R;i++)

             if(q[i].flag==)

                 ans[q[i].id]=l;

         return;

     }

     int mid=(l+r)>>,l1=,l2=;

     for(int i=L;i<=R;i++){//cout<<i<<endl;

         if(q[i].flag==){//cout<<"***"<<endl;

             if(q[i].k<=mid)

                 change(q[i].x,q[i].y,,),q1[l1++]=q[i];

             else

                 q2[l2++]=q[i];

         }

         else{//cout<<"&&&"<<endl;

             int lala=query(q[i].x,q[i].y,);//cout<<q[i].x<<" "<<q[i].y<<" "<<lala<<" "<<tree[1].sum<<endl;

             if(lala>=q[i].k)

                 q1[l1++]=q[i];

             else

                 q[i].k-=lala,q2[l2++]=q[i];

         }

     }//cout<<"lala"<<endl;

     for(int i=;i<l1;i++)

         if(q1[i].flag==)

             change(q1[i].x,q1[i].y,-,);

 //  memcpy(q+L,q1,sizeof(q[0])*l1);

 //  memcpy(q+L+l1,q2,sizeof(q[0])*l2);

     for(int i=;i<l1;i++)

         q[L+i]=q1[i];

     for(int i=;i<l2;i++)

         q[L+l1+i]=q2[i];

     solve(L,L+l1-,l,mid),solve(L+l1,R,mid+,r);

 }

 signed main(void){

 //  freopen("sequence10.in","r",stdin);

 //  freopen("out.out","w",stdout);

     n=intread();m=intread();cnt=;

 //  scanf("%lld%lld",&n,&m);

     for(int i=,l,s,x,y;i<=m;i++){

         l=intread();

 //      scanf("%lld",&l);

         if(l==)

             q[i].x=intread(),q[i].y=intread(),q[i].k=longread(),q[i].k=n-q[i].k+,q[i].flag=;

         else

             q[i].x=intread(),q[i].y=intread(),q[i].k=longread(),q[i].id=++cnt,q[i].flag=;

     }build(,n,);solve(,m,,n);

     for(int i=;i<=cnt;i++)

         printf("%lld\n",n-ans[i]+);

     return ;

 }//Cap ou pas cap. Cap.

BZOJ 2527: [Poi2011]Meteors

时空隧道

还是板子题...感觉我要废了TAT...

代码：

 #include<algorithm>

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 #include<vector>

 //by NeighThorn

 using namespace std;

 const int maxn=+;

 int n,m,k,cnt,ans[maxn],need[maxn];

 vector<int> v[maxn];

 unsigned long long tr[maxn];

 struct M{

     int x,y,k,id,flag;

 }q[maxn<<],q1[maxn<<],q2[maxn<<];

 inline void insert(int x,int y){

     for(;x<=m;x+=x&(-x))

         tr[x]+=y;

 }

 inline unsigned long long query(int x){

     unsigned long long res=;

     for(;x;x-=x&(-x))

         res+=tr[x];

     return res;

 }

 inline void solve(int L,int R,int l,int r){

     if(L>R)

         return;

     if(l==r){

         for(int i=L;i<=R;i++)

             if(q[i].flag==)

                 ans[q[i].id]=l;

         return;

     }

     int mid=(l+r)>>,l1=,l2=;

     for(int i=L;i<=R;i++){

         if(q[i].flag==){

             if(q[i].id<=mid)

                 insert(q[i].x,q[i].k),insert(q[i].y+,-q[i].k),q1[l1++]=q[i];

             else

                 q2[l2++]=q[i];

         }

         else{

             unsigned long long lala=;

             for(int j=;j<v[q[i].id].size();j++)

                 lala+=query(v[q[i].id][j]);

             if(lala>=q[i].k)

                 q1[l1++]=q[i];

             else

                 q[i].k-=lala,q2[l2++]=q[i];

         }

     }

     for(int i=;i<l1;i++)

         if(q1[i].flag==)

             insert(q1[i].x,-q1[i].k),insert(q1[i].y+,q1[i].k);

     for(int i=;i<l1;i++)

         q[L+i]=q1[i];

     for(int i=;i<l2;i++)

         q[L+l1+i]=q2[i];

     solve(L,L+l1-,l,mid),solve(L+l1,R,mid+,r);

 }

 signed main(void){

     scanf("%d%d",&n,&m);

     for(int i=,x;i<=m;i++)

         scanf("%d",&x),v[x].push_back(i);

     for(int i=;i<=n;i++)

         scanf("%d",&need[i]);

     scanf("%d",&k);

     for(int i=,s,x,y;i<=k;i++){

         scanf("%d%d%d",&x,&y,&s);

         if(x<=y)

             q[++cnt].x=x,q[cnt].y=y,q[cnt].k=s,q[cnt].flag=,q[cnt].id=i;

         else

             q[++cnt].x=x,q[cnt].y=m,q[cnt].k=s,q[cnt].flag=,q[cnt].id=i,

             q[++cnt].x=,q[cnt].y=y,q[cnt].k=s,q[cnt].flag=,q[cnt].id=i;

     }

     for(int i=;i<=n;i++)

         q[++cnt].k=need[i],q[cnt].id=i,q[cnt].flag=;

     solve(,cnt,,k+);

     for(int i=;i<=n;i++){

         if(ans[i]==k+)

             puts("NIE");

         else

             printf("%d\n",ans[i]);

     }

     return ;

 }//Cap ou pas cap. Cap.

By NeighThorn

整体二分QAQ的更多相关文章

  1. BZOJ2527 [Poi2011]Meteors 整体二分 树状数组

    原文链接http://www.cnblogs.com/zhouzhendong/p/8686460.html 题目传送门 - BZOJ2527 题意 有$n$个国家. 太空里有$m$个太空站排成一个圆 ...

  2. BZOJ3110 [Zjoi2013]K大数查询 树套树 线段树 整体二分 树状数组

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ3110 题意概括 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位 ...

  3. 一篇自己都看不懂的CDQ分治&整体二分学习笔记

    作为一个永不咕咕咕的博主,我来更笔记辣qaq CDQ分治 CDQ分治的思想还是比较简单的.它的基本流程是: \(1.\)将所有修改操作和查询操作按照时间顺序并在一起,形成一段序列.显然,会影响查询操作 ...

  4. 洛谷P3242 接水果 [HNOI2015] 整体二分

    正解:整体二分+树状数组 解题报告: 传送门! 题目还是大概解释下?虽然其实是看得懂的来着,,, 大概就是说给一棵树.给定一些询问,每个询问都是说在两个点之间的路径上的子路径的第k大是什么 然后看到这 ...

  5. 洛谷P3250 网络 [HNOI2016] 整体二分

    正解:整体二分+树状数组 解题报告: 传送门! 亲这里的建议是用整体二分呢 dbq最近看sd淘宝说话体看多了有点脑抽,,, 首先考虑如果是单组询问怎么做昂QAQ 考虑二分答案 对于所有比mid小的操作 ...

  6. BZOJ.2738.矩阵乘法(整体二分 二维树状数组)

    题目链接 BZOJ 洛谷 整体二分.把求序列第K小的树状数组改成二维树状数组就行了. 初始答案区间有点大,离散化一下. 因为这题是一开始给点,之后询问,so可以先处理该区间值在l~mid的修改,再处理 ...

  7. [学习笔记] CDQ分治&整体二分

    突然诈尸.png 这两个东西好像都是离线骗分大法... 不过其实这两个东西并不是一样的... 虽然代码长得比较像 CDQ分治 基本思想 其实CDQ分治的基本思想挺简单的... 大概思路就是长这样的: ...

  8. BZOJ3110 K大数查询 【线段树 + 整体二分 或 树套树(非正解)】

    Description 有N个位置,M个操作.操作有两种,每次操作如果是1 a b c的形式表示在第a个位置到第b个位置,每个位置加入一个数c 如果是2 a b c形式,表示询问从第a个位置到第b个位 ...

  9. P5163 WD与地图（整体二分+权值线段树）

    传送门 细节要人命.jpg 这题思路太新奇了--首先不难发现可以倒着做变成加边,但是它还需要我们资瓷加边的同时维护强连通分量.显然加边之后暴力跑是不行的 然后有一个想法,对于一条边\((u,v)\), ...

随机推荐

  1. java head space/ java.lang.OutOfMemoryError: Java heap space内存溢出

    上一篇JMX/JConsole调试本地还可以在centos6.5 服务器上进行监控有个问题端口只开放22那么设置的9998端口 你怎么都连不上怎么监控?(如果大神知道还望指点,个人见解) 线上项目出现 ...

  2. 23种设计模式--责任链模式-Chain of Responsibility Pattern

    一.责任链模式的介绍 责任链模式用简单点的话来说,将责任一步一步传下去,这就是责任,想到这个我们可以相当击鼓传花,这个是为了方便记忆,另外就是我们在项目中经常用到的审批流程等这一类的场景时我们就可以考 ...

  3. 利用bootstrap的carousel.js实现轮播图动画

    前期准备: 1.jquery.js. 2.bootstrap的carousel.js. 3.bootstrap.css. 如果大家不知道在哪下载,可以联系小颖,小颖把这些js和css可以发送给你. 一 ...

  4. Power BI官方视频(3) Power BI Desktop 8月份更新功能概述

    Power BI Desktop 8月24日发布了更新版本.现将更新内容翻译整理如下,可以根据后面提供的链接下载最新版本使用. 1.主要功能更新 1.1 数据钻取支持在线版 以前的desktop中进行 ...

  5. 算法与数据结构(八) AOV网的关键路径

    上篇博客我们介绍了AOV网的拓扑序列,请参考<数据结构(七) AOV网的拓扑排序(Swift面向对象版)>.拓扑序列中包括项目的每个结点,沿着拓扑序列将项目进行下去是肯定可以将项目完成的, ...

  6. 简约之美Jodd-http--深入源码理解http协议

    Jodd 是一个开源的 Java 工具集, 包含一些实用的工具类和小型框架.简单,却很强大! jodd-http是一个轻巧的HTTP客户端.现在我们以一个简单的示例从源码层看看是如何实现的? Http ...

  7. 获取打开的Word文档

    using Word = Microsoft.Office.Interop.Word; int _getApplicationErrorCount=0; bool _isMsOffice = true ...

  8. beans.xml

    <?xml version="1.0" encoding="UTF-8"?> <beans xmlns="http://www.sp ...

  9. （翻译）FIFO In Hardware

    翻译一些自己觉得有价值的材料,工作中碰到英语大多数是读,基本没有写或者翻的,翻得不好不到位的敬请指摘. 同时也附原文以供参考. http://electronics.stackexchange.com ...

  10. SQL中字符串拼接

    1. 概述 在SQL语句中经常需要进行字符串拼接,以sqlserver,oracle,mysql三种数据库为例,因为这三种数据库具有代表性. sqlserver: select '123'+'456' ...