记录编号 452397 评测结果 AAAAAAAAAA
题目名称 [HZOI 2016] 数列操作d 最终得分 100
用户昵称 Gravatarswttc 是否通过 通过
代码语言 C++ 运行时间 6.613 s
提交时间 2017-09-19 15:31:27 内存使用 137.64 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<cmath>
#define mid (l+r>>1)
#define ls (o<<1)
#define rs ((o<<1)+1)
#define av 0.5 

using namespace std;

const int MOD=1e9+7;
long long n,m,dep[3000010],ha[3000010],ans,lazc[3000010];
long double lazy[3000010],sum[3000010];
void buildt(int l,int r,int o)
{
	if(l==r){
		return;
	}
	buildt(l,mid,ls);
	buildt(mid+1,r,rs);
	dep[o]=dep[ls]+1;
	return;
}

void pushdown(int l,int r,int o)
{
	if(dep[o]==0||lazy[o]==0) return;
	lazc[ls]+=lazc[o];
	sum[ls]+=(long double)(mid-l+1)*(lazy[o]-(long double)lazc[o]*pow(2,dep[o]-1)*av);
	lazy[ls]+=lazy[o]-(long double)lazc[o]*pow(2,dep[o]-1)*av;
	lazc[rs]+=lazc[o];
	sum[rs]+=(long double)(r-mid)*(lazy[o]+(long double)lazc[o]*pow(2,dep[o]-1)*av);
	lazy[rs]+=lazy[o]+(long double)lazc[o]*pow(2,dep[o]-1)*av;
	lazy[o]=0.00;
	lazc[o]=0;
}

void change(int l,int r,int o,int ll,int rr,int x)
{
	if(l>=ll&&r<=rr){
		if(l==r){
			sum[o]+=(long double)(l-ll)*(long double)x;
		}
		else{
			sum[o]+=(long double)(l-ll+r-ll)/2.0*(long double)(r-l+1)*(long double)x;
			lazy[o]+=(long double)(l-ll+r-ll)/2.0*(long double)x;
			lazc[o]+=x;
		}
		return;
	}
	pushdown(l,r,o);
	if(mid>=ll) change(l,mid,ls,ll,rr,x);
	if(mid<rr) change(mid+1,r,rs,ll,rr,x);
	sum[o]=sum[ls]+sum[rs];
	return;
}

void query(int l,int r,int o,int ll,int rr)
{	 
	if(l>=ll&&r<=rr){
		ans+=(long long)sum[o];
		return;
	}
	pushdown(l,r,o);
	if(mid>=ll) query(l,mid,ls,ll,rr);
	if(mid<rr) query(mid+1,r,rs,ll,rr);
	return;
}

int main()
{
	freopen("segment.in","r",stdin);
	freopen("segment.out","w",stdout);
	scanf("%lld%lld",&n,&m);
	long long temp=1;
	for(int i=1;i<=n;i++){
		if(i>temp){
			temp*=2;
		}
		ha[i]=temp;
	}
	buildt(1,ha[n],1);
	for(int i=1;i<=m;i++){
		char c;
		cin>>c;
		if(c=='1'){
			int a,b,c;
			scanf("%d%d%d",&a,&b,&c);
			change(1,ha[n],1,a,b,c);
		}
		else{
			int a,b;
			scanf("%d%d",&a,&b);
			ans=0;
			query(1,ha[n],1,a,b);
			printf("%lld\n",ans%MOD);
		}
	}//system("pause");
	return 0;
}