bzoj4822: [Cqoi2017]老C的任务（扫描线+BIT/线段树）

　　裸题...

　　依旧是写了BIT和线段树两种（才不是写完线段树后才想起来可以写BIT呢

　　怎么卡常数都挺大...QAQ ccz和yy的写法好快哇%%%

BIT：

#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<map>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
struct poi{int fir,sec,trd,ty,pos,w;}a[maxn];
int n,m,x,y,z,cnt,N,x1,x2,y1,y2;
int b[maxn];
ll tree[maxn],ans[maxn];
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
inline bool cmp(poi a,poi b){return a.fir==b.fir?a.ty<b.ty:a.fir<b.fir;}
inline int lowbit(int x){return x&-x;}
inline void add(int x,int delta){for(;x<=N;x+=lowbit(x))tree[x]+=delta;}
inline ll query(int x){ll sum=;for(;x;x-=lowbit(x))sum+=tree[x];return sum;}
int main()
{
    read(n);read(m);cnt=n;
    for(int i=;i<=n;i++)
    {
        read(x);read(y);read(z);
        a[i].fir=x;a[i].sec=y;a[i].w=z;b[++N]=y;
    }
    for(int i=;i<=m;i++)
    {
        read(x1);read(y1);read(x2);read(y2);
        a[++cnt].fir=x1;a[cnt].sec=y1;a[cnt].trd=y2;a[cnt].ty=-;a[cnt].pos=i;
        a[++cnt]=a[cnt-];a[cnt].fir=x2;a[cnt].ty=;b[++N]=y1;b[++N]=y2;
    }
    sort(b+,b++N);
    for(int i=;i<=n;i++)a[i].sec=lower_bound(b+,b++N,a[i].sec)-b;
    for(int i=n+;i<=cnt;i+=)
    {
        a[i].sec=a[i+].sec=lower_bound(b+,b++N,a[i].sec)-b;
        a[i].trd=a[i+].trd=lower_bound(b+,b++N,a[i].trd)-b;
    }
    sort(a+,a++cnt,cmp);
    for(int i=;i<=cnt;i++)
    {
        if(a[i].ty)
        {
            ll sum=query(a[i].trd)-query(a[i].sec-);
            ans[a[i].pos]+=(a[i].ty==?sum:-sum);
        }
        else add(a[i].sec,a[i].w);
    }
    for(int i=;i<=m;i++)printf("%lld\n",ans[i]);
    return ;
}

线段树：

#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<map>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
struct poi{int fir,sec,trd,ty,pos,w;}a[maxn];
int n,m,x,y,z,cnt,N,x1,x2,y1,y2;
int b[maxn];
ll tree[maxn*],ans[maxn];
void read(int &k)
{
    int f=;k=;char c=getchar();
    while(c<''||c>'')c=='-'&&(f=-),c=getchar();
    while(c<=''&&c>='')k=k*+c-'',c=getchar();
    k*=f;
}
inline bool cmp(poi a,poi b){return a.fir==b.fir?a.ty<b.ty:a.fir<b.fir;}
inline void pushup(int x){tree[x]=tree[x<<]+tree[x<<|];}
inline void add(int x,int l,int r,int cx,int delta)
{
    if(l==r){tree[x]+=delta;return;}
    int mid=(l+r)>>;
    if(cx<=mid)add(x<<,l,mid,cx,delta);
    else add(x<<|,mid+,r,cx,delta);
    pushup(x);
}
inline ll query(int x,int l,int r,int cl,int cr)
{
    if(cl<=l&&r<=cr)return tree[x];
    int mid=(l+r)>>;ll ret=;
    if(cl<=mid)ret+=query(x<<,l,mid,cl,cr);
    if(cr>mid)ret+=query(x<<|,mid+,r,cl,cr);
    return ret;
}
int main()
{
    read(n);read(m);cnt=n;
    for(int i=;i<=n;i++)
    {
        read(x);read(y);read(z);
        a[i].fir=x;a[i].sec=y;a[i].w=z;b[++N]=y;
    }
    for(int i=;i<=m;i++)
    {
        read(x1);read(y1);read(x2);read(y2);
        a[++cnt].fir=x1;a[cnt].sec=y1;a[cnt].trd=y2;a[cnt].ty=-;a[cnt].pos=i;
        a[++cnt]=a[cnt-];a[cnt].fir=x2;a[cnt].ty=;b[++N]=y1;b[++N]=y2;
    }
    sort(b+,b++N);
    for(int i=;i<=n;i++)a[i].sec=lower_bound(b+,b++N,a[i].sec)-b;
    for(int i=n+;i<=cnt;i+=)
    {
        a[i].sec=a[i+].sec=lower_bound(b+,b++N,a[i].sec)-b;
        a[i].trd=a[i+].trd=lower_bound(b+,b++N,a[i].trd)-b;
    }
    sort(a+,a++cnt,cmp);
    for(int i=;i<=cnt;i++)
    {
        if(a[i].ty)
        {
            ll sum=query(,,N,a[i].sec,a[i].trd);
            ans[a[i].pos]+=(a[i].ty==?sum:-sum);
        }
        else add(,,N,a[i].sec,a[i].w);
    }
    for(int i=;i<=m;i++)printf("%lld\n",ans[i]);
    return ;
}

在yy大爷的指导下卡了一波常...艹上了#7，但是写得丑还是写得丑比yy大爷还是慢不少

#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<map>
#define ll long long
using namespace std;
const int maxn=,inf=1e9;
struct poi{int fir,sec,trd,ty,pos;}a[maxn];
int n,m,x,y,z,cnt,N,x1,x2,y1,y2;
int b[maxn],w[maxn];
ll tree[maxn],ans[maxn];
char buf[],*ptr=buf-;
int read()
{
    int x=,f=,c=*++ptr;
    while(c<)c=='-'&&(f=-),c=*++ptr;
    while(c>)x=x*+c-,c=*++ptr;
    return x*f;
}
inline bool cmp(poi a,poi b){return a.fir==b.fir?a.ty<b.ty:a.fir<b.fir;}
inline int lowbit(int x){return x&-x;}
inline void add(int x,int delta){for(;x<=N;x+=lowbit(x))tree[x]+=delta;}
inline ll query(int x){ll sum=;for(;x;x-=lowbit(x))sum+=tree[x];return sum;}
int main()
{
    fread(buf,,sizeof(buf),stdin)[buf]=;
    n=read();m=read();cnt=n;
    for(int i=;i<=n;i++)
    {
        x=read();y=read();w[i]=read();
        a[i].fir=x;a[i].sec=y;a[i].pos=i;b[++N]=y;
    }
    for(int i=;i<=m;i++)
    {
        x1=read();y1=read();x2=read();y2=read();
        a[++cnt].fir=x1;a[cnt].sec=y1;a[cnt].trd=y2;a[cnt].ty=-;a[cnt].pos=i;
        a[++cnt]=a[cnt-];a[cnt].fir=x2;a[cnt].ty=;
    }
    sort(b+,b++N);
    for(int i=;i<=n;i++)a[i].sec=lower_bound(b+,b++N,a[i].sec)-b;
    for(int i=n+;i<=cnt;i+=)
    {
        a[i].sec=a[i+].sec=lower_bound(b+,b++N,a[i].sec)-b;
        a[i].trd=a[i+].trd=upper_bound(b+,b++N,a[i].trd)-b-;
    }
    sort(a+,a++cnt,cmp);
    for(int i=;i<=cnt;i++)
    {
        if(a[i].ty)
        {
            ll sum=query(a[i].trd)-query(a[i].sec-);
            ans[a[i].pos]+=(a[i].ty==?sum:-sum);
        }
        else add(a[i].sec,w[a[i].pos]);
    }
    for(int i=;i<=m;i++)printf("%lld\n",ans[i]);
    return ;
}