中南oj 2018年12月月赛Round2-Div1+Div2 题解（水题不写了）

想歪一波榜，结果B题读错题自闭了....被C题卡了，结果E题线段树给错过了.....（中南的人好喜欢出线段树啊）

B(2214): Sequence Magic

思路：我们分析一下：

 

假设上图的a-b插到下图的c后面（c>=a），我们要是知道下图的一个位置p，该怎么退出p在上图的位置呢？首先下图的p要是不在后两段区间内，倒推上图位置不变，（以下用的区间数值为上图的值）如果p在下图的中间段[ a , c ]，倒推回去，p就加上上图中间段的长，如果p在下图的第三段[ c+1,c+a-b+1 ]，倒推回去的话，p就减去上图的第三段的长。这种情况已经完全搞定了，还有一种c<a的情况靠你来自己推理咯。

#include<bits/stdc++.h>
using namespace std;
const int maxn=10005;
long long a[maxn],b[maxn],c[maxn];
int main()
{
	int n,m,k;
	cin>>n>>m>>k;
	for(int i=1;i<=m;i++)
	cin>>a[i]>>b[i]>>c[i];
	for(int j=1;j<=k;j++)
	{
		long long p;
		scanf("%lld",&p);
		for(int i=m;i;i--)
		{
			if(c[i]>=a[i])
			{
				if(p>=a[i]&&p<=c[i])p+=b[i]-a[i]+1;
				else if(p>c[i]&&p<=c[i]+b[i]-a[i]+1)p-=(c[i]-a[i]+1);
			}
			else
			{
				if(p>c[i]&&p<=c[i]+b[i]-a[i]+1)p+=a[i]-c[i]-1;
				else if(p>c[i]+b[i]-a[i]+1&&p<=b[i])p-=(b[i]-a[i]+1);
			}
		}
		if(p>n)p-=n;
		printf("%lld\n",p);
	}
}

C(2215): Maximum Submatrix

写超时了，暂时不会。//求这题的题解代码

E(2217): Data Structure

裸线段树，首先区间最大最小求和是基本操作，我主要讲讲怎么区间开根号以及求区间的方差。

其实区间开根号很简单，int 型最大2^32，一个数最多开5次就会变成1，因此每次开根号我都直接开到叶子节点就行，不用担心复杂度，每个叶子节点最多更新5次。

区间方差怎么求？设av为平均数，v为方差，s为区间的长度，对于一个节点，其左儿子方差为：，假设总区间的平均数已经变成了av-d，那么我怎么求这个节点的方差*这个节点的长度呢，对于左儿子的贡献变形一下：（下面的符号均是左儿子的符号），我再变形：s*v+2*d*sum-2*d*av*s+d*d*s，右儿子同理，最后把得到的值除以该节点的长度就是方差，这样就每个节点的方差都可以由两个儿子O(1)复杂度转化而来。

#include<bits/stdc++.h>
#define db double
using namespace std;
const int maxn=1e5+10;
struct node
{
	int s,mx,mn,sum;
	db av,v;
	node operator+(const node& t)const
	{
		node p;
		p.s=s+t.s;
		p.mx=max(mx,t.mx);
		p.mn=min(mn,t.mn);
		p.sum=sum+t.sum;
		p.av=(db(sum)+t.sum)/p.s;
		db d1=av-p.av;
		db d2=t.av-p.av;
		p.v=v*s+2*d1*sum-2*d1*av*s+d1*d1*s;
		p.v+=t.v*t.s+2*d2*t.sum-2*d2*t.av*t.s+d2*d2*t.s;
		p.v/=p.s;
		return p;
	}
}tree[maxn*4];
void build(int o,int l,int r)
{
	if(l==r)
	{
		scanf("%d",&tree[o].mx);
		tree[o].mn=tree[o].sum=tree[o].mx;
		tree[o].s=1;
		tree[o].av=(db)tree[o].mx;
		tree[o].v=0.0;
		return;
	}
	int m=(l+r)/2,ls=o*2,rs=o*2+1;
	build(ls,l,m);
	build(rs,m+1,r);
	tree[o]=tree[ls]+tree[rs]; 
}
node qu(int o,int l,int r,int ql,int qr)
{
	if(l>=ql&&r<=qr)return tree[o];
	int m=(l+r)/2,ls=o*2,rs=o*2+1;
	if(qr<=m)return qu(ls,l,m,ql,qr);
	else if(ql>m)return qu(rs,m+1,r,ql,qr);
	return qu(ls,l,m,ql,qr)+qu(rs,m+1,r,ql,qr);
}
void up(int o,int l,int r,int ql,int qr)
{
	if(tree[o].mx==1)return;
	int m=(l+r)/2,ls=o*2,rs=o*2+1;
	if(l>=ql&&r<=qr)
	{
		if(l==r)
		{
			tree[o].mn=tree[o].sum=tree[o].mx=(int)sqrt(tree[o].mx);
			tree[o].av=(db)tree[o].mx;
			return;
		}
		up(ls,l,m,ql,qr);
		up(rs,m+1,r,ql,qr);
	}
	else
	{
		if(ql<=m)up(ls,l,m,ql,qr);
		if(qr>m)up(rs,m+1,r,ql,qr);
	}
	tree[o]=tree[ls]+tree[rs];
}
int main()
{
	int n,m,op,l,r;
	cin>>n>>m;
	build(1,1,n);
	while(m--)
	{
		scanf("%d%d%d",&op,&l,&r);
		if(op!=5)
		{
			node p=qu(1,1,n,l,r);
			if(op==1)printf("%d\n",p.sum);
			else if(op==2)printf("%d\n",p.mx);
			else if(op==3)printf("%d\n",p.mn);
			else printf("%.2lf\n",p.v);
		}
		else
		up(1,1,n,l,r);
	}
}