bzoj4604: The kth maximum number

Description

神犇Aleph在陪蒟蒻Bob玩一个游戏。神犇Aleph随手在地面上画了一个巨大无比的二维平面，然后在其上做一些微
小的贡献或教蒟蒻Bob做人。由于神犇Aleph是队爷，所以有时他会施用”队爷光环”这一魔法，在二维平面上的某坐
标整点处添加一个权值为v的贡献；又由于神犇Aleph是神犇，所以有时他会施用”嘲讽”这一技能，询问蒟蒻Bob在
矩形区域(x1, y1), (x2, y2)(x1≤x2且y1≤y2，包括边界)中，神犇Aleph所做的第k大的贡献是多少。由于神犇Al
eph是Au爷，所以他不会在同一个坐标整点处做两次或两次以上的贡献。现在神犇Aleph希望蒟蒻Bob回答他的每次
询问。然而蒟蒻Bob傻傻不会做，于是来求助您，宇宙第一神犇，请您来回答神犇Aleph的每次询问。
Input

输入的第一行为两个正整数N, Q，表示横纵坐标的范围和神犇Aleph的操作次数（包括贡献次数和询问次数）。
接下来Q行，每行代表神犇Aleph的一个操作，操作格式如下：
首先第一个数字type，表示操作种类。type=1表示贡献，type=2表示询问。
若type=1，接下来会有三个正整数x, y, v，表示在坐标整点(x, y)添加一个贡献v。(1≤x, y≤N, 1≤v≤10^9)
若type=2，接下来会有五个正整数x1, y1, x2, y2, k，表示询问矩形区域(x1, y1), (x2, y2)中第k大的贡献。
(1≤x1≤x2≤N，1≤y1≤y2≤N，1≤k≤Q)
初始时平面上不存在贡献。
本题共有7组测试数据。对于所有的数据，N≤500,000。
Q的范围见下表：
测试点1-2 Q=1,000
测试点3-7 Q=50,000
Output

对于每个询问(type=2的操作)，回答第k大的贡献。若不存在第k大的贡献，请输出”NAIVE!ORZzyz.”（输出不含双引号）。
Sample Input

10 7

1 1 1 1

1 2 2 3

1 4 1 2

1 3 4 4

2 1 1 4 1 3

2 2 2 3 5 4

2 2 1 4 4 2
Sample Output

NAIVE!ORZzyz.

NAIVE!ORZzyz.

3

题解

反正可以离线
我们写一个整体二分套CDQ分治就好了
不只是第一次写，还是效率变低的缘故，居然写了一个晚上QAQ
一开始还误打误撞地过了样例TAT，然而不仅bug多，找的还是k小，这样例也是没谁了
不知道为什么，跑得特别慢

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<cstring>
using namespace std;
const int N=500005*4;
int n,q;
struct qq
{
    int x1,y1,x2,y2;
    int k;
    bool tf;//true:增加操作   false:询问操作 
    int id;//哪一个询问
    int f;//权   CDQ的时候用 
}s[N];
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;
}
int maxx=0;
int mymax (int x,int y){return x>y?x:y;}
int mymin (int x,int y){return x,y?x:y;}
int ans[N];
qq k[N];int num;
int tt[N];
int f[N];//树状数组 
bool cmp (qq a,qq b){return a.x1<b.x1;}
int lb (int x){return x&(-x);}
void change (int x,int y)
{
    while (x<=n)
    {
        f[x]+=y;
        x+=lb(x);
    }
}
int get (int x)
{
    int lalal=0;
    while (x>=1)
    {
        lalal=lalal+f[x];
        x-=lb(x);
    }
    return lalal;
}
void cdq (int l,int r)
{
    if (l>=r) return ;
    int mid=(l+r)>>1;
    cdq(l,mid);cdq(mid+1,r);
    sort(k+l,k+1+mid,cmp);sort(k+1+mid,k+1+r,cmp);
    int L=l;
    for (int u=mid+1;u<=r;u++)//更新答案
    {
        while (L<=mid)
        {
            if (k[L].x1>k[u].x1) break;
            if (k[L].tf==true) change(k[L].y1,1);
            L++;
        }
        if (k[u].tf==true) continue;
        tt[k[u].id]=tt[k[u].id]+k[u].f*get(k[u].y1);
    }
    for (int u=l;u<L;u++)
        if (k[u].tf)    
            change(k[u].y1,-1);
}
qq aa[N],bb[N];//分治 
void solve (int l,int r,int L,int R)//二分的范围   区间 
{
    if (l==r)
    {
        for (int u=L;u<=R;u++)  
        if (s[u].tf==false)
        {
            ans[s[u].id]=l;
        }
        return ;
    }
    int mid=(l+r)>>1;   
    for (int u=L;u<=R;u++) tt[u]=0;
    num=0;
    for (int u=L;u<=R;u++)
    {
        if (s[u].tf==true&&s[u].k>mid)//修改操作
        { 
            k[++num]=s[u];
        }
        else if (s[u].tf==false)
        {
            k[++num].f=1;  k[num].x1=s[u].x2;   k[num].y1=s[u].y2;   k[num].id=u;    k[num].tf=false;
            k[++num].f=-1; k[num].x1=s[u].x1-1; k[num].y1=s[u].y2;   k[num].id=u;     k[num].tf=false;
            k[++num].f=-1; k[num].x1=s[u].x2;   k[num].y1=s[u].y1-1; k[num].id=u;    k[num].tf=false;
            k[++num].f=1;  k[num].x1=s[u].x1-1; k[num].y1=s[u].y1-1; k[num].id=u;    k[num].tf=false;
        }
    }
    cdq(1,num);
    int a_num=0,b_num=0;
    for (int u=L;u<=R;u++)
    {
        if (s[u].tf==true)
        {
            if (s[u].k>mid)
                aa[++a_num]=s[u];
            else
                bb[++b_num]=s[u];
        }
        else
        {
            if (tt[u]>=s[u].k)
            {
                aa[++a_num]=s[u];
            }
            else
            {
                bb[++b_num]=s[u];
                bb[b_num].k-=tt[u];
            }
        }
    }
    bool tf=false,tf1=false;
    for (int u=1;u<=a_num;u++)
    {
        s[L+u-1]=aa[u];
        if (aa[u].tf==false) tf=true;
    }
    for (int u=1;u<=b_num;u++)
    {
        s[L+a_num+u-1]=bb[u];
        if (bb[u].tf==false) tf1=true;
    }
    if (tf==true)  solve(mid+1,r,L,L+a_num-1);
    if (tf1==true) solve(l,mid,L+a_num,R);
}
bool ooo[N];
int main()
{
    memset(f,0,sizeof(f));
    n=read();q=read();
    for (int u=1;u<=q;u++)
    {
        int op=read();
        if (op==1)//增加操作
        {
            s[u].x1=read();
            s[u].y1=read();
            s[u].k=read();
            ooo[u]=true;
            s[u].tf=true;
            maxx=mymax(maxx,s[u].k);
        }
        else
        {
            ooo[u]=false;
            s[u].x1=read(); 
            s[u].y1=read();
            s[u].x2=read();
            s[u].y2=read();
            s[u].k=read();
            s[u].tf=false;
            s[u].id=u;
        }
    }
    memset(ans,0,sizeof(ans));
    solve(0,maxx,1,q);
    for (int u=1;u<=q;u++)
        if (ooo[u]==false)
        {
            if (ans[u]==0) printf("NAIVE!ORZzyz.\n");
            else printf("%d\n",ans[u]);
        }
    return 0;
}