bzoj4262: Sum

Description

 

Input

第一行一个数 t，表示询问组数。

第一行一个数 t，表示询问组数。

接下来 t 行，每行四个数 l_1, r_1, l_2, r_2。

 

Output

一共 t 行，每行一个数 Sum。

Sample Input

4
1 3 5 7
2 4 6 8
1 1 9 9
9 9 1 1

Sample Output

9322587654
9025304064
1065645568
0

HINT

1<=t<=40000,1<=L1<R1<=10^5,1<=L2<=R2<=10^5

题解：

http://www.cnblogs.com/clrs97/p/4824806.html

code：

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long int64;
char ch;
bool ok;
void read(int &x){
    ok=0;
    for (ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;
    for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
    if (ok) x=-x;
}
const int maxn=100005;
const int mod=1E9;
const int inf=2147483647;
int n,q,l1,r1,l2,r2,stack[maxn],top;
int64 a[maxn],ans[maxn];
int cnt;
struct Quer{
    int x,l,r,id,op;
}quer[maxn];
bool cmp(const Quer &a,const Quer &b){return a.x<b.x;}
struct Data{
    int64 a,b,c,d;
    void init(){a=1,b=c=d=0;}    
};
Data operator+(const Data &x,const Data &y){return (Data){x.a*y.a,x.b*y.a+y.b,y.c*x.a+x.c,y.c*x.b+x.d+y.d};}
struct seg{
    #define ls k<<1
    #define rs (k<<1)+1
    int len[maxn<<2];
    int64 v[maxn<<2],s[maxn<<2];
    Data tag[maxn<<2];
    void addtag(int k,Data t){s[k]=t.c*v[k]+t.d*len[k]+s[k],v[k]=t.a*v[k]+t.b*len[k],tag[k]=tag[k]+t;}
    void pushdown(int k){addtag(ls,tag[k]),addtag(rs,tag[k]),tag[k].init();}
    void update(int k){v[k]=v[ls]+v[rs],s[k]=s[ls]+s[rs];}
    void build(int k,int l,int r){
        tag[k].init(),len[k]=r-l+1,v[k]=s[k]=0;
        if (l==r) return;
        int m=(l+r)>>1;
        build(ls,l,m),build(rs,m+1,r);
    }
    void cover(int k,int l,int r,int x,int y,Data t){
        if (l==x&&r==y){addtag(k,t);return;}
        int m=(l+r)>>1; pushdown(k);
        if (y<=m) cover(ls,l,m,x,y,t);
        else if (x<=m) cover(ls,l,m,x,m,t),cover(rs,m+1,r,m+1,y,t);
        else cover(rs,m+1,r,x,y,t);
        update(k);
    }
    void cover(int x,int y,int64 v){cover(1,1,n,x,y,(Data){0,v,0,0}),addtag(1,(Data){1,0,1,0});}
    int64 query(int k,int l,int r,int x,int y){
        if (l==x&&r==y) return s[k];
        int m=(l+r)>>1; pushdown(k);
        if (y<=m) return query(ls,l,m,x,y);
        else if (x<=m) return query(ls,l,m,x,m)+query(rs,m+1,r,m+1,y);
        else return query(rs,m+1,r,x,y);
    }
    int64 query(int x,int y){return query(1,1,n,x,y);}
}T;
int main(){
    read(q);
    for (int i=1;i<=q;i++){
        read(l1),read(r1),read(l2),read(r2),n=max(n,r1),n=max(n,r2);
        quer[++cnt]=(Quer){r2,l1,r1,i,1},quer[++cnt]=(Quer){l2-1,l1,r1,i,-1};
    }
    sort(quer+1,quer+cnt+1,cmp);
    for (int i=1,t1=1023,t2=1025;i<=n;i++,t1=1LL*t1*1023%mod,t2=1LL*t2*1025%mod) a[i]=t1^t2;
    T.build(1,1,n);
    a[0]=inf,stack[top=1]=0;
    for (int i=1,j=1;i<=n;i++){
        while (top&&a[stack[top]]<=a[i]) top--;
        T.cover(stack[top]+1,i,a[i]);
        while (quer[j].x<i) j++;
        for (;quer[j].x==i;j++) ans[quer[j].id]+=T.query(quer[j].l,quer[j].r)*quer[j].op;
        stack[++top]=i;    
    }
    T.build(1,1,n);
    a[0]=-inf,stack[top=1]=0;
    for (int i=1,j=1;i<=n;i++){
        while (top&&a[stack[top]]>=a[i]) top--;
        T.cover(stack[top]+1,i,a[i]);
        while (quer[j].x<i) j++;
        for (;quer[j].x==i;j++) ans[quer[j].id]-=T.query(quer[j].l,quer[j].r)*quer[j].op;
        stack[++top]=i;    
    }
    for (int i=1;i<=q;i++) printf("%lld
",ans[i]);
    return 0;
}