BNU 12680 Jupiter Atacks! (線段樹)

來源:互聯網
上載者:User

題意:給定一個區間,,兩種操作,一種是將x的值賦為y,另一種為詢問區間(x,y)的sum值,此處sum值的運算為:

例如: 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 <iostream>#include <algorithm>#include <cmath>#include<functional>#include <cstdio>#include <cstdlib>#include <cstring>#include <string>#include <vector>#include <set>#include <queue>#include <stack>#include <climits>//形如INT_MAX一類的#define MAX 100050#define INF 0x7FFFFFFF#define L(x) x<<1#define R(x) x<<1|1using 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<n; i++) {            getchar();            scanf("%c%d%d", &op,&x,&y);            if(op == 'E') {                update(x,1,y - a[x]);                a[x] = y;            } else {                ans = 0;                //test();                query(x,y,1);                printf("%lld\n",ans);            }        }        puts("-");    }    return 0;}


相關文章

聯繫我們

該頁面正文內容均來源於網絡整理,並不代表阿里雲官方的觀點,該頁面所提到的產品和服務也與阿里云無關,如果該頁面內容對您造成了困擾,歡迎寫郵件給我們,收到郵件我們將在5個工作日內處理。

如果您發現本社區中有涉嫌抄襲的內容,歡迎發送郵件至: info-contact@alibabacloud.com 進行舉報並提供相關證據,工作人員會在 5 個工作天內聯絡您,一經查實,本站將立刻刪除涉嫌侵權內容。

A Free Trial That Lets You Build Big!

Start building with 50+ products and up to 12 months usage for Elastic Compute Service

  • Sales Support

    1 on 1 presale consultation

  • After-Sales Support

    24/7 Technical Support 6 Free Tickets per Quarter Faster Response

  • Alibaba Cloud offers highly flexible support services tailored to meet your exact needs.