2017 暑假艾教集训 day9(整体二分 + cdq分治 cdq真是我女神！！！)

最新推荐文章于 2019-08-01 08:12:51 发布

德鸭

最新推荐文章于 2019-08-01 08:12:51 发布

阅读量268

点赞数

本文链接：https://blog.csdn.net/qq_36908995/article/details/77018659

版权

POJ 2104 求第k小

做法：整体二分加树状数组。

离线把询问和插入做整体二分，每次把s，t扫一般如果是插入就在线段树上相应的位置插入（val，1）。

对于每个询问查询在l，r 中用多少个数

如果加上累加的数小于k那么答案应该在右区域把该询问扔到右区域并把查询的结果累加（左区间本来有查询结果个前k小，直接进入右边会少算）。

如果大于 k 的话那么答案应该在左端扔进左端

在一系列操作做完后还要在树状数组中相应的拆点本轮插入的点！！！

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
const int maxn=200006;
const int inf=0x3f3f3f3f;
int n,m;
int c[maxn];

int lowbit(int x)
{
    return x&-x;
}
void updata(int pos,int val)
{
    for(int i=pos;i<=n;i+=lowbit(i)) c[i]+=val;
}
int query(int p)
{
    int sum=0;
    for(int i=p;i>0;i-=lowbit(i)) sum+=c[i];
    return sum;
}


int ans[maxn];
int tmp[maxn<<1];
struct node
{
    int l,r,k;
    int cur,val,index;
    int kind;
}q[maxn<<1],q1[maxn],q2[maxn];

void init()
{
    memset(tmp,0,sizeof(tmp));
    memset(c,0,sizeof(c));
    for(int i=1;i<=n+m;++i)
    {
        q[i].cur = q1[i].cur =q2[i].cur=0;
        q[i].val = q1[i].val =q2[i].val=0;
    }
}

void work(int s,int t,int l,int r)
{
    if(s>t) return;
    if(l==r)
    {
        for(int i=s;i<=t;++i)
        {
            if(q[i].kind==2)
            {
                ans[q[i].index] = l;
            }
        }
        return;
    }
    int mid=(l+r)>>1;
    int n1=0 , n2=0 , flag1=0, flag2=0 ;
    for(int i=s;i<=t;++i)
    {
        if(q[i].kind==1)
        {
            if(q[i].val<=mid)  { updata(q[i].index,1); q1[n1++] = q[i]; }
            else q2[n2++] = q[i];
        }

        else if(q[i].kind==2)
        {
            tmp[i] = query(q[i].r) - query(q[i].l-1);
            if(tmp[i]+q[i].cur >= q[i].k)
            {
                q1[n1++] = q[i]; flag1=1;
            }
            else
            {
               q[i].cur+= tmp[i] ;  q2[n2++] = q[i]   ; flag2=1;
            }
        }
    }

    for(int i=s;i<=t;++i)
    {
        if(q[i].kind==1 && q[i].val<=mid) updata(q[i].index,-1);
    }

    for(int i=0;i<n1;++i) q[i+s] = q1[i];
    for(int i=0;i<n2;++i) q[i+s+n1]=q2[i];

    if(flag1)   work(s,s+n1-1,l,mid);
    if(flag2)   work(s+n1,t,mid+1,r);

}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        init();
        for(int i=1;i<=n;++i)
        {
            scanf("%d",&q[i].val);
            q[i].kind=1;
            q[i].index=i;
        }
        for(int i=n+1;i<=n+m;++i)
        {
            scanf("%d%d%d",&q[i].l , &q[i].r ,&q[i].k);
            q[i].kind=2;
            q[i].index=i-n;
        }
        work(1,n+m,-inf,inf);
        for(int i=1;i<=m;++i)
        {
            printf("%d\n",ans[i]);
        }
    }
    return 0;
}

Bzoj 4237（Bzoj好多好题，但要权限号无语了）

做法：祖传分治+单调栈+二分。注意二分，二分写炸了看了别人的题解。

首先把稻田从中间分为两边，分别是上边和下边。用cdq分治可以算出只有上边或下边的情况。我们只需要算出一个点在上边，一个点在下边的情况。

首先不难发现我们要保证 x 的递增性。同时对于上边的稻草人应该有一个从下到上的数据结构进行维护，而下边的稻草人需要从上往下进行维护。用个单调栈即可。

对于每一个上面的点算出下面相对于的值累加即可。

很容易发现他应该满足 x值在前一个点和当前点之间。这些满足条件的下边点数就是答案。

应该x是单调递增的，所以只需要二分就可以了（二分要小心！！！）；

#include <bits/stdc++.h>
using namespace std;
const int maxn=200005;
typedef long long ll;
int n;
struct point
{
   int x,y;
}a[maxn],b[maxn];
int q[maxn],p[maxn];
bool cmp(const point a , const point b)
{
    return a.y<b.y;
}
ll ans=0;

int bs(int x,int l,int r)
{
    while(l+1<r)
    {
        int mid=(l+r)>>1;
        if(a[q[mid]].x < x) l=mid;
        else r=mid;
    }
    return l;
}
void aqx_bs(int l,int r)
{
    if(l==r) return;
    int mid=(l + r)>>1;
    aqx_bs(l,mid);
    aqx_bs(mid+1,r);
    int top1=0 ,top2=0;
    for(int i=mid+1,j=l;i<=r;++i)
    {
        while(top1 && a[i].y < a[ p[top1]].y ) top1--;
        p[++top1]=i;
        for(; j<=mid && a[j].x < a[i].x ; j++)
        {
           while(top2 && a[j].y>a[q[top2]].y ) top2--;
           q[++top2]=j;
        }
        ans += top2 - bs(a[p[top1-1]].x,0,top2+1);
    }
    int lx = l , ly=mid+1;
    for(int i=l;i<=r;++i) b[i] = ((lx<=mid &&a[lx].x<a[ly].x) || (ly>r)) ? a[lx++]:a[ly++];
    for(int i=l;i<=r;++i) a[i]=b[i];
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;++i)
    {
       scanf("%d%d",&a[i].x,&a[i].y);
    }
    sort(a+1,a+1+n,cmp);
    ans=0;
    aqx_bs(1,n);
    printf("%lld\n",ans);
    return 0;
}

HDU 5324（这道题刚开始拿线段树杠了四个小时，然后一个小时的树状数组就过了，哎）谢谢各位大佬的博客哇

做法：三维偏序问题，自然想到cdq的分治对吧。但这道题有些不一样，它有两个条件 1.要输出路径 2.当答案相同时输出字典序最小的（输出答案的时候正序扫到第一个满足条件的点即可）

那么问题关键就是怎么维护它的路径，其实这类题我们都是处理过的，我们需要维护每个答案的pre ，然后从后往前走。

为了满足这一条件，我们把题目的所有条件都倒过来做一遍分治，然后输出时找它的pre即可；

在做分治的时候不要改变原来的信息，应该新开一个数组来做。

1.每次在做的时候因为 id 是从大到小的所以先做右区间，然后处理右区间和左区间的关系，最后处理左区间。和刚才一样，每一遍做完后还要把树状数组还原！！！

对于每一个dp 只需要找1-L中最大的dp即可。

#include <bits/stdc++.h>
using namespace std;
const int maxn=1e5+100;
typedef long long ll;
struct node
{
    int l,r,id;
    bool flag;
    node(int a=0,int b=0,int c=0,bool d=0) :l(a), r(b) ,id(c),flag(d){}

}xhd[maxn<<1],p[maxn<<1];
int n,book[maxn<<1],tree[maxn<<1],dp[maxn<<1];
bool cmp(const node a, const node b)
{
    if(a.r != b.r) return a.r > b.r ;
    if(a.l != b.l) return a.l < b.l;
    return a.id > b.id;
}
int lowbit(int x)
{
    return x& (-x);
}
void updata(int pos,int val)
{
    for(int i=pos;i<maxn;i+=lowbit(i))
    tree[i] = max(tree[i],val);
}
int query(int pos)
{
    int ans=0;
    for(int i=pos;i>0;i-=lowbit(i))
        ans=max(ans,tree[i]);
    return ans;
}
void clr(int pos)
{
    for(int i=pos;i<maxn;i+=lowbit(i))
        tree[i]=0;
}
void cdq(int begin,int end)
{
    if(begin>=end) return;
    int mid=(begin + end)>>1;
    cdq(mid+1,end);
    int tot=0;
    for(int i=begin;i<=mid;++i)
    {
        p[++tot] = node(xhd[i].l,xhd[i].r,xhd[i].id,1);
    }
    for(int i=mid+1;i<=end;++i)
    {
        p[++tot] = node(xhd[i].l,xhd[i].r,xhd[i].id,0);
    }
    sort(p+1,p+tot+1,cmp);
    for(int i=1;i<=tot;++i)
    {
        if(p[i].flag==0)
        {
            updata(p[i].l,dp[p[i].id]);
        }
        else
        {
            dp[p[i].id] = max(dp[p[i].id],query(p[i].l)+1);
        }
    }
    for(int i=1;i<=tot;++i)
    {
        if(p[i].flag==0) clr(p[i].l);
    }
    cdq(begin,mid);
}

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        memset(tree,0,sizeof(tree));
        int cnt=0;
        for(int i=1;i<=n;++i)    scanf("%d",&xhd[i].l);
        for(int i=1;i<=n;++i)    scanf("%d",&xhd[i].r);
        for(int i=1;i<=n;++i)
        {
            dp[i]=1;
            xhd[i].id=i;
        }
        for(int i=1;i<=n;++i)
        {
            book[cnt++] = xhd[i].l;
            book[cnt++] = xhd[i].r;
        }
        sort(book,book+cnt);
        cnt = unique(book,book+cnt) -book;
        for(int i=1;i<=n;++i)
        {
            xhd[i].l = lower_bound(book , book+cnt , xhd[i].l) -book +1;
            xhd[i].r = lower_bound(book , book+cnt , xhd[i].r) -book +1;
        }
        cdq(1,n);

        int temp= -1;
        for(int i=1;i<=n;++i) temp=max(dp[i],temp);
        printf("%d\n",temp);

        int pre=-1;
        for(int i=1;i<=n;++i)
        {
            if(dp[i]==temp && (pre==-1 ||(xhd[i].l<=xhd[pre].l && xhd[i].r >= xhd[pre].r)))
            {
                temp--;
                printf("%d%c",i,(temp==0? '\n':' '));
                pre=i;
            }
        }
    }
    return 0;
}

德鸭

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
2017 暑假艾教集训 day9(整体二分 + cdq分治 cdq真是我女神！！！)

POJ 2104 求第k小做法：整体二分加树状数组。离线把询问和插入做整体二分，每次把s，t扫一般如果是插入就在线段树上相应的位置插入（val，1）。对于每个询问查询在l，r 中用多少个数如果加上累加的数小于k那么答案应该在右区域把该询问扔到右区域并把查询的结果累加（左区间本来有查询结果个前k小，直接进入右边会少算）。如果大于 k 的话那么答案应该在左端扔进左端
复制链接

扫一扫