5798 - Jupiter Atacks!

http://livearchive.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=3809

单点更新+区间查询+维护区间公式和

这题算是昨天开始刷线段树做过的比较难的了吧....首先要把维护量想出来, 再把合并区间的方式想出来, 然后写的时候 还有很多要注意的地方......


注: 思考的时候一定要明确一点, 让每个区间都满足一个性质, 合并结果的时候, 面对的都是满足这个性质的子区间. 合并出来的也是满足这个性质的.


三个地方被卡了(T_T.........):

1/ 在query的时候, 左右区间合并结果也要按公式算.而不是直接相加.

2/ 在query的时候, 如果并不是左右区间都占, 则 不能按 公式算.....写的时候要分开, 或者 判断.

3/ 在query的时候, 右区间长度并不是 (r-mid) , 而是 (R-mid)........R-mid 才是查询 的 真实的长度...


代码:

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#include<string>
#include<vector>
#include<map>
#include<algorithm>
using namespace std;
inline int Rint() { int x; scanf("%d", &x); return x; }
inline int max(int x, int y) { return (x>y)? x: y; }
inline int min(int x, int y) { return (x<y)? x: y; }
#define FOR(i, a, b) for(int i=(a); i<=(b); i++)
#define FORD(i,a,b) for(int i=(a);i>=(b);i--)
#define REP(x) for(int i=0; i<(x); i++)
typedef long long int64;
#define INF (1<<30)
const double eps = 1e-8;
#define bug(s) cout<<#s<<"="<<s<<" "

//	单点更新, 区间查询,	维护 ....一个恶心的式子的求和

#define MAXN 100002
int a[MAXN<<2];
int fac[MAXN];
int B, P, L, N;
void pushup(int l, int r, int e)
{
	int mid =  (l+r)>>1;
	int cnt = (r-mid);
	//int fac = 1;
	//FOR(k, 1, cnt) fac*=B;
	//a[e] = (int)(((int64)a[e>>1]*fac+a[e<<1|1])%P);
	a[e] = (int)(((int64)a[e<<1]*fac[cnt]+a[e<<1|1])%P);
	//bug(a[e<<1]);bug(fac[cnt]);bug(a[e<<1|1])<<endl;
}
//void build()
void update(int p, int v, int l, int r, int e)
{
	if(l==r)
		a[e] = v;
	else
	{
		int mid = (l+r)>>1;
		if(p<=mid)
			update(p, v, l, mid, e<<1);
		else
			update(p, v, mid+1,  r,  e<<1|1);
		pushup(l, r, e);
	}
}
int query(int L, int R, int l, int r, int e)
{
	if(L<=l && r<=R)
		return a[e];
	else
	{
		int mid = (l+r)>>1;
		int flag1=0, flag2=0;
		int ret = 0;
		if(L<=mid)
		{
			ret=query(L, R, l, mid, e<<1);			//这里也是合并子区间!!!! 也要按 那个公式, 而不是简单 相加
			flag1=1;
		}
		//ret=(int)(((int64)ret*fac[r-mid])%P);												//如果是单个区间就不用合并!!!!!		很tricky啊 
		int ret2 = 0;
		if(mid+1<=R)
		{
			ret2=query(L, R, mid+1, r, e<<1|1);
			flag2 = 1;
		}

		if(flag1 && flag2)
			//ret = (int)(((int64)ret*fac[r-mid]+ret2)%P);
			ret = (int)(((int64)ret*fac[R-mid]+ret2)%P);		//是R-mid不是r-mid.... R是真正查询的区间.....恩
		else
			ret+=ret2;

		//	也 可以这样, 直接 判断是左右都占还是单单占一边.
		//int64  res, res1,res2, tmp;
		//if(R<=mid)
		//	res=query(L,R,l,mid,e<<1);
		//else if(mid<L)
		//	res=query(L,R,mid+1,r,e<<1|1);
		//else
		//{
		//	res1=query(L,mid,l,mid,e<<1);
		//	res2=query(mid+1,R,mid+1,r,e<<1|1);
		//	int Len=R-mid;
		//	res=res1*fac[Len]+res2;
		//	if(res>=P)res%=P;
		//}
		//return (int)res;
		return ret;
	}
}
bool read()
{
	scanf("%d%d%d%d", &B, &P, &L, &N);
	return B+P+L+N != 0;
}
int main()
{
	while(read())
	{
		//bug(L)<<endl;
		memset(a, 0, sizeof(a));
		fac[0] = 1;
		FOR(i, 1, L)
		{
			fac[i]=(int)(((int64)fac[i-1]*B)%P);
			//bug(fac[i])<<endl;
		}
		REP(N)
		{
			char  ch[2];
			scanf("%s",ch); char op=ch[0];
			int x = Rint();
			int y =Rint();
			if(op == 'E')
				update(x, y, 1, L, 1);
			else
				printf("%d\n", query(x, y, 1, L, 1));
		}
		puts("-");
	}
}


发布了179 篇原创文章 · 获赞 5 · 访问量 20万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 大白 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览