题目描述
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: