[Balkan 2007] Mokia

[题目链接]

https://www.lydsy.com/JudgeOnline/problem.php?id=1176

[算法]

CDQ分治 + 树状数组即可

时间复杂度 : O(Nlog^2N)

[代码]

#include<bits/stdc++.h>
using namespace std;
const int N = ;
const int M = ;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;

struct Query
{
        int pos , x , y , value , type , id;
} q[N * ] , t1[N * ] , t2[N * ];

int s , w , m , k;
int c[M] , ans[N];

template <typename T> inline void chkmax(T &x,T y) { x = max(x,y); }
template <typename T> inline void chkmin(T &x,T y) { x = min(x,y); }
template <typename T> inline void read(T &x)
{
    T f = ; x = ;
    char c = getchar();
    for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
    for (; isdigit(c); c = getchar()) x = (x << ) + (x << ) + c - '';
    x *= f;
}
inline bool cmp(Query a , Query b)
{
        if (a.x != b.x) return a.x < b.x;
        else if (a.y != b.y) return a.y < b.y;
        else return a.type < b.type;
}
inline int lowbit(int x)
{
        return x & (-x);
}
inline void modify(int x , int val)
{
        for (int i = x; i <= w; i += lowbit(i))
                c[i] += val;
}
inline int query(int x)
{
        int ret = ;
        for (int i = x; i; i -= lowbit(i))
                ret += c[i];
        return ret;
}
inline void cdq(int l , int r)
{
        int mid = (l + r) >> ;
        if (l == r) return;
        for (int i = l; i <= r; i++)
        {
                if (q[i].type ==  && q[i].pos <= mid) modify(q[i].y , q[i].value);
                else if (q[i].type ==  && q[i].pos > mid) ans[q[i].id] += q[i].value * query(q[i].y);
        }
        for (int i = l; i <= r; i++)
        {
                if (q[i].type ==  && q[i].pos <= mid)
                        modify(q[i].y , -q[i].value);
        }
        int l1 =  , l2 = ;
        for (int i = l; i <= r; i++)
                if (q[i].pos <= mid) t1[++l1] = q[i];
                else t2[++l2] = q[i];
        for (int i = ; i <= l1; i++) q[l + i - ] = t1[i];
        for (int i = ; i <= l2; i++) q[l + l1 + i - ] = t2[i];
        cdq(l , mid);
        cdq(mid +  , r);
}

int main()
{

        read(s); read(w);
        while (true)
        {
                int type;
                read(type);
                if (type == ) break;
                if (type == )
                {
                        int x , y , a;
                        read(x); read(y); read(a);
                        q[++m] = (Query){m , x , y , a ,  , k};
                } else
                {
                        int X1 , Y1 , X2 , Y2;
                        read(X1); read(Y1); read(X2); read(Y2);
                        ans[++k] = s * (X2 - X1 + ) * (Y2 - Y1 + );
                        q[++m] = (Query){m , X2 , Y2 ,  ,  , k};
                        q[++m] = (Query){m , X1 -  , Y2 , - ,  , k};
                        q[++m] = (Query){m , X2 , Y1 -  , - ,  , k};
                        q[++m] = (Query){m , X1 -  , Y1 -  ,  ,  , k};
                }
        }
        sort(q +  , q + m +  , cmp);
        cdq( , m);
        for (int i = ; i <= k; i++) printf("%d\n" , ans[i]);

        return ;

}