BZOJ 2752 HAOI2012高速公路

最新推荐文章于 2019-05-20 22:25:18 发布

Rayment_cc

最新推荐文章于 2019-05-20 22:25:18 发布

阅读量199

点赞数

分类专栏：线段树 BZOJ 文章标签：线段树

本文链接：https://blog.csdn.net/As_A_Kid/article/details/80370557

版权

BZOJ 同时被 2 个专栏收录

89 篇文章 0 订阅

订阅专栏

线段树

19 篇文章 0 订阅

订阅专栏

Problem

BZOJ

Solution

前置技能： $\sum_{i=l}^r i=(l+r)(r-l+1)/2$ , $\sum_{i=1}^n i^2=n*(n+1)*(2*n+1)/6$

一道期望的题目，直接推公式。
我们知道，对于一个长度为n的路段，第i段路计入贡献的方案一共有(i-l+1)(r-i)种，即左右两边的端点随意搭配。
先只考虑分子，暴力拆公式

\sum i = l r - 1 (i - l + 1) (r - i) v i

$\sum_{i=l}^{r-1} (i-l+1)(r-i)v_i$

= \sum i = l r - 1 (- i 2 + (l + r - 1) i - l r + r) v i

$=\sum_{i=l}^{r-1} (-i^2+(l+r-1)i-lr+r)v_i$

= - \sum i = l r - 1 i 2 v i + (l + r - 1) \sum i = l r - 1 i v i + (r - l r) \sum i = l r - 1 v i

$=-\sum_{i=l}^{r-1} i^2v_i+(l+r-1)\sum_{i=l}^{r-1}iv_i+(r-lr)\sum_{i=l}^{r-1}v_i$
所以我们的线段树就需要维护这么三个玩意，

i2vi,ivi,vi i 2 v i , i v i , v i $i^2v_i,iv_i,v_i$ ，然后区间更新的时候再手推一下公式就好了。

Code

#include <cstdio>
#define pushup(x) s1[x]=s1[x<<1]+s1[x<<1|1],s2[x]=s2[x<<1]+s2[x<<1|1],s3[x]=s3[x<<1]+s3[x<<1|1]
using namespace std;
typedef long long ll;
const int maxn=100010;
int n,m;
ll a=-1,b,c,aa,ab;
char op[5];
ll add[maxn<<2],s1[maxn<<2],s2[maxn<<2],s3[maxn<<2];
template <typename Tp> inline void read(Tp &x)
{
    x=0;int f=0;char ch=getchar();
    while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();
    if(ch=='-') f=1,ch=getchar();
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    if(f) x=-x;
}
inline ll gcd(ll a,ll b){return b?gcd(b,a%b):a;}
inline ll sum1(int l,int r){return (ll)(l+r)*(r-l+1)/2;}
inline ll sum2(int n){return (ll)n*(n+1)*(n<<1|1)/6;}
void pushdown(int l,int m,int r,int rt)
{
    s1[rt<<1]+=(ll)add[rt]*(m-l+1);s2[rt<<1]+=sum1(l,m)*add[rt];
    s3[rt<<1]+=(sum2(m)-sum2(l-1))*add[rt];add[rt<<1]+=add[rt];
    s1[rt<<1|1]+=(ll)add[rt]*(r-m);s2[rt<<1|1]+=sum1(m+1,r)*add[rt];
    s3[rt<<1|1]+=(sum2(r)-sum2(m))*add[rt];add[rt<<1|1]+=add[rt];
    add[rt]=0;
}
void update(int l,int r,int L,int R,int val,int rt)
{
    if(L<=l&&r<=R)
    {
        s1[rt]+=(ll)val*(r-l+1);s2[rt]+=sum1(l,r)*val;
        s3[rt]+=(sum2(r)-sum2(l-1))*val;add[rt]+=val;
        return ;
    }
    int m=(l+r)>>1;
    if(l!=r&&add[rt]) pushdown(l,m,r,rt);
    if(L<=m) update(l,m,L,R,val,rt<<1);
    if(m<R) update(m+1,r,L,R,val,rt<<1|1);
    pushup(rt);
}
void query(int l,int r,int L,int R,int rt)
{
    if(L<=l&&r<=R){aa+=a*s3[rt]+b*s2[rt]+c*s1[rt];return ;}
    int m=(l+r)>>1;
    if(l!=r&&add[rt]) pushdown(l,m,r,rt);
    if(L<=m) query(l,m,L,R,rt<<1);
    if(m<R) query(m+1,r,L,R,rt<<1|1);
}
int main()
{
    #ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
    #endif
    int l,r,val;ll d;
    read(n);read(m);
    for(int i=1;i<=m;i++)
    {
        scanf("%s",op);
        if(op[0]=='C')
        {
            read(l);read(r);read(val);
            update(1,n,l,r-1,val,1);
        }
        else if(op[0]=='Q')
        {
            read(l);read(r);
            b=(l+r-1);c=r-(ll)l*r;aa=0;
            query(1,n,l,r-1,1);ab=(ll)(r-l+1)*(r-l)/2;
            d=gcd(aa,ab);aa/=d;ab/=d;
            printf("%lld/%lld\n",aa,ab);
        }
    }
    return 0;
}