多颗线段树+间隔点组成的区间求和 hdu4288 coder

2019-04-13 21:59发布

站内文章 / 模拟电子
12591 0 1697
题意： hdu4288
有3种操作：1.往集合里加入元素 2.删除集合里的元素 3.对集合里下标模5等于3的元素求和
加入和删除元素就是线段树里的单点更新，但第三种操作就不那么显然了。由于需要求和的元素都是模5等于3的等间隔的点，当我们对一个节点进行更新的时候，它的左子结点的满足条件的下标在这个节点肯定也满足条件，但右子节点模5等于3的下标并不是这个节点的模5等于3的下标，因为当我们计算右子节点的时候下标从1开始计算，而在当前节点的时候右子节点的下标是从左子结点的最后一个下标开始算的，所以我们用一个数组cnt来记录每个节点所在区间有多少点，每个节点有一个sum【5】数组记录它的模5的每个值的下标代表的值的和。这样当前节点rt的sum【rt】【i】=sum【左子结点】【i】+sum【右子节点】【i-左子结点数的个数】。
代码： #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; #define PB push_back #define MP make_pair #define REP(i,x,n) for(int i=x;i<(n);++i) #define FOR(i,l,h) for(int i=(l);i<=(h);++i) #define FORD(i,h,l) for(int i=(h);i>=(l);--i) #define SZ(X) ((int)(X).size()) #define ALL(X) (X).begin(), (X).end() #define RI(X) scanf("%d", &(X)) #define RII(X, Y) scanf("%d%d", &(X), &(Y)) #define RIII(X, Y, Z) scanf("%d%d%d", &(X), &(Y), &(Z)) #define DRI(X) int (X); scanf("%d", &X) #define DRII(X, Y) int X, Y; scanf("%d%d", &X, &Y) #define DRIII(X, Y, Z) int X, Y, Z; scanf("%d%d%d", &X, &Y, &Z) #define OI(X) printf("%d",X); #define RS(X) scanf("%s", (X)) #define MS0(X) memset((X), 0, sizeof((X))) #define MS1(X) memset((X), -1, sizeof((X))) #define LEN(X) strlen(X) #define F first #define S second #define Swap(a, b) (a ^= b, b ^= a, a ^= b) #define Dpoint strcut node{int x,y} #define cmpd int cmp(const int &a,const int &b){return a>b;} /*#ifdef HOME freopen("in.txt","r",stdin); #endif*/ const int MOD = 1e9+7; typedef vector VI; typedef vector VS; typedef vector VD; typedef long long LL; typedef pair PII; //#define HOME int Scan() { int res = 0, ch, flag = 0; if((ch = getchar()) == '-') //判断正负 flag = 1; else if(ch >= '0' && ch <= '9') //得到完整的数 res = ch - '0'; while((ch = getchar()) >= '0' && ch <= '9' ) res = res * 10 + ch - '0'; return flag ? -res : res; } /*----------------PLEASE-----DO-----NOT-----HACK-----ME--------------------*/ char op[100000+10][20]; int X[100000+10]; long long int sum[400000+10][5]; int cnt[400000+10]; int a[100000+10]; void pushup(int rt) { cnt[rt]=cnt[rt<<1]+cnt[rt<<1|1]; REP(i,0,5) { sum[rt][i]=sum[rt<<1][i]; } int offset=cnt[rt<<1]; REP(i,0,5) { sum[rt][(i+offset)%5]+=sum[rt<<1|1][i]; } } void update(int l,int r,int rt,int num,int o) {if(l==r) { cnt[rt]=o; sum[rt][1]=o*(X[l]); return; } int m=(l+r)>>1; if(num<=m) update(l,m,rt<<1,num,o); else update(m+1,r,rt<<1|1,num,o); pushup(rt); } int main() {int n; while(RI(n)!=EOF) { int nn=0; REP(i,0,n) { scanf("%s",op[i]); if(op[i][0]!='s') { int x; RI(x); a[i]=x; X[nn++]=x; } } sort(X,X+nn); int k=unique(X,X+nn)-X; MS0(sum); MS0(cnt); REP(i,0,n) { if(op[i][0]!='s') {int pos=lower_bound(X,X+k,a[i])-X; if(op[i][0]=='a') update(0,k-1,1,pos,1); else update(0,k-1,1,pos,0);} else printf("%I64d ",sum[1][3]); } } return 0; }

热门文章