BZOJ2752 洛谷2221 HAOI2012 高速公路

题目描述
Y901高速公路是一条重要的交通纽带，政府部门建设初期的投入以及使用期间的养护费用都不低，因此政府在这条高速公路上设立了许多收费站。
Y901高速公路是一条由N-1段路以及N个收费站组成的东西向的链，我们按照由西向东的顺序将收费站依次编号为1~N，从收费站i行驶到i+1(或从i+1行驶到i)需要收取Vi的费用。高速路刚建成时所有的路段都是免费的。
政府部门根据实际情况，会不定期地对连续路段的收费标准进行调整，根据政策涨价或降价。
无聊的小A同学总喜欢研究一些稀奇古怪的问题，他开车在这条高速路上行驶时想到了这样一个问题:对于给定的l,r(l

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<cstdlib>
using namespace std;
#define MAXN 100005
#define LL long long
struct Seg_Tree{ int l,r; LL s[4],v; }tre[MAXN<<2];
int n,m;
char s[5];
LL ans[5];
LL GCD(LL a,LL b){
    return b==0?a:GCD(b,a%b);
}
inline int read(int &x){
    x=0;int f=1; char c=getchar();
    while(c>'9'||c<'0'){ if(c=='-')f=-1; c=getchar(); }
    while(c>='0'&&c<='9'){ x=(x<<1)+(x<<3)+(c-'0');c=getchar(); }
    x*=f;
}
LL Get(int l,int r){
    return 1LL*r*(r+1)*(2*r+1)/6LL-1LL*(l-1)*l*(2*l-1)/6LL;
}
inline void UpDate(int u,LL v){
    int cnt=tre[u].r-tre[u].l+1;
    tre[u].s[1]+=1LL*cnt*v;
    tre[u].s[2]+=1LL*(tre[u].r+tre[u].l)*cnt/2LL*v;
    tre[u].s[3]+=1LL*Get(tre[u].l,tre[u].r)*v;
    tre[u].v+=v;
}
inline void Push_Up(int u){
    for(int i=1;i<=3;i++)
        tre[u].s[i]=tre[u<<1].s[i]+tre[u<<1|1].s[i];
}
inline void Push_Down(int u){
    UpDate(u<<1,tre[u].v);UpDate(u<<1|1,tre[u].v);
    tre[u].v=0;
}
void Add(int u,int l,int r,LL val){
    if(l<=tre[u].l&&tre[u].r<=r){ UpDate(u,val); return ; }
    if(tre[u].v!=0) Push_Down(u);
    int Mid=(tre[u].l+tre[u].r)>>1;
    if(l<=Mid) Add(u<<1,l,r,val);
    if(r>Mid) Add(u<<1|1,l,r,val);
    Push_Up(u);
}
LL Query(int u,int l,int r,int k){
    if(l<=tre[u].l&&tre[u].r<=r) return tre[u].s[k];
    LL ans=0;
    if(tre[u].v!=0) Push_Down(u);
    int Mid=(tre[u].l+tre[u].r)>>1;
    if(l<=Mid) ans+=Query(u<<1,l,r,k);
    if(r>Mid) ans+=Query(u<<1|1,l,r,k);
    return ans;
}
void Build(int u,int l,int r){
    tre[u].l=l;tre[u].r=r;
    tre[u].v=tre[u].s[0]=tre[u].s[1]=tre[u].s[2]=tre[u].s[3]=0;
    if(l==r) return;
    int Mid=(l+r)>>1;
    Build(u<<1,l,Mid);
    Build(u<<1|1,Mid+1,r);
}
int main(){
    scanf("%d",&n);scanf("%d",&m);
    Build(1,1,n-1);//第n个点没有用
    for(int i=1;i<=m;i++){
        scanf("%s",s);LL l,r,v;
        scanf("%lld%lld",&l,&r);
        if(s[0]=='C'){ scanf("%lld",&v);Add(1,(int)l,(int)r-1,v); }
        else{
            LL x=(r-l+1)*(r-l)/2LL;
            for(int i=1;i<=3;i++)
                ans[i]=Query(1,(int)l,(int)r-1,i);
            ans[0]=ans[1]*(r-l*r)+ans[2]*(l+r-1)-ans[3];
            LL g=GCD(ans[0],x);
            printf("%lld/%lld\n",ans[0]/g,x/g);
        }
    }
    return 0;
}

感谢来自cnblog CQzhangyu 的讲解

还有！
本机不过样例的AC：