BNU 12680 Jupiter Atacks! (线段树) - c++编程基础

题意：给定一个区间，，两种操作，一种是将x的值赋为y，另一种为询问区间（x，y）的sum值，此处sum值的运算为：

H(fi,..., fj) = Bkfj-k ( mod P)

例如： H(f2,..., f5) = B^0 * f5 + B^1 * f4 + B^2 * f3 + B^3 * f2 ( mod P)

这题隐藏的好深~~~ 用线段树的话，用va域存储树内所有的H（tr[x].l，tr[x].r），更新时根据更新点和区间右端点乘以个（B ^ (tr[x].r - pos)）,求和时根据求和区间右端点和线段树区间的右端点乘以个（B ^( r - tr[x].r））。

需要注意的是，取余过程中可能会产生负数，需要归正。

#include   
#include   
#include   
#include  
#include   
#include   
#include   
#include   
#include   
#include   
#include   
#include   
#include //形如INT_MAX一类的  
#define MAX 100050  
#define INF 0x7FFFFFFF  
#define L(x) x<<1  
#define R(x) x<<1|1  
using namespace std;  
  
int b,p,l,n;  
long long B[MAX],ans;  
int a[MAX];  
struct node {  
    int l,r,mid;  
    long long va;  
}tr[MAX * 4];  
  
void init() {  
    memset(a,0,sizeof(a));  
    B[0] = 1;  
    for(int i=1; i<=l; i++) {  
        B[i] = (B[i - 1] * b) % p;  
    }  
}  
  
void build(int l,int r,int x) {  
    tr[x].l = l;  
    tr[x].r = r;  
    tr[x].mid = (l + r) >> 1;  
    tr[x].va = 0;  
    if(l == r) return ;  
    build(l,tr[x].mid,L(x)) ;  
    build(tr[x].mid + 1, r, R(x));  
}  
  
void update(int l,int x,int va) {  
    tr[x].va = (tr[x].va + B[tr[x].r - l] * va % p ) % p;  
    if(tr[x].va < 0) tr[x].va += p;  
    if(tr[x].l == tr[x].r) return ;  
    if(l >




 tr[x].mid) update(l,R(x),va);  
    else update(l,L(x),va);  
}  
  
void query(int l,int r,int x) {  
    if(l <= tr[x].l && r >= tr[x].r) {  
        ans = (ans + B[r - tr[x].r] * tr[x].va % p) % p;  
        if(ans < 0) ans += p;  
        return ;  
    }  
    if(r > tr[x].mid) query(l,r,R(x));  
    if(l <= tr[x].mid) query(l,r,L(x));  
}  
  
void test() {  
    for(int i=0; i<=3 * l; i++) {  
        printf("%d %d %lld\n",tr[i].l,tr[i].r,tr[i].va);  
    }  
}  
  
int main() {  
    while(scanf("%d%d%d%d",&b,&p,&l,&n)) {  
        if(b == 0 && p == 0 && l == 0 && n == 0) break;  
        char op;  
        int x,y;  
        init();  
        build(1,l,1);  
        for(int i=0; i