zoukankan      html  css  js  c++  java
  • 板子2

      线段树(单点更新)

    /*  gyt
           Live up to every day            */
    #include<cstdio>
    #include<cmath>
    #include<iostream>
    #include<algorithm>
    #include<vector>
    #include<stack>
    #include<cstring>
    #include<queue>
    #include<set>
    #include<string>
    #include<map>
    #include <time.h>
    #define PI acos(-1)
    using namespace std;
    typedef long long ll;
    typedef double db;
    const int maxn = 500002;
    const ll maxm = 1e7;
    const ll base = 2333;
    const int INF = 1<<30;
    const db eps = 1e-8;
    const ll mod = 1e9+13;
    struct nodee{
        int left, right, mid, num;
    }node[maxn];
    int ca=1;
    
    void build(int left, int right, int n) {
        node[n].left=left;
        node[n].right=right;
        node[n].mid=(left+right)/2;
        node[n].num=0;
        if (left+1==right)  return;
        build(left, (left+right)/2, 2*n);
        build((left+right)/2, right, 2*n+1);
    }
    void updata(int pos, int value, int n) {
        node[n].num+=value;
        if (node[n].left+1==node[n].right)
            return;
        if (pos<node[n].mid)
            updata(pos, value, 2*n);
        else
            updata(pos, value, 2*n+1);
    }
    int query(int left, int right, int n) {
        if (node[n].left==left && node[n].right==right) {
            return node[n].num;
        }
        if (left<node[n].mid) {
            if(right <= node[n].mid ){
                return query(left,right,2*n);
            }
            else{
                return query(left,node[n].mid,2*n) + query(node[n].mid,right,2*n + 1);
            }
        }
        else
            return query(left,right,2*n+1);
    }
    void solve() {
        int n;  scanf("%d", &n);
        build(1, n+1, 1);
        for (int i=1; i<=n; i++) {
            int x;  scanf("%d", &x);
            updata(i, x, 1);
        }
        char c[10];
        printf("Case %d:
    ", ca++);
        while(1) {
            scanf("%s", c);
            if (strcmp(c, "End")==0)  break;
            int x, y;  scanf("%d%d", &x, &y);
            if (strcmp(c, "Query")==0) {
                int ans=query(x, y+1, 1);
                printf("%d
    ", ans);
            }
             else if(strcmp(c,"Add") == 0){
                updata(x, y, 1);
            }
            else if(strcmp(c,"Sub") == 0){
                updata(x, -y, 1);
            }
        }
    }
    int main() {
        int t = 1;
       // freopen("in.txt","r",stdin);
      //  freopen("gcd.out","w",stdout);
        scanf("%d", &t);
        while(t--)
            solve();
        return 0;
    }

       dfs序+线段树

        将dfs序作为建树的标准

    /*  gyt
           Live up to every day            */
    #pragma comment(linker, "/STACK:1024000000,1024000000")
    #include<cstdio>
    #include<cmath>
    #include<iostream>
    #include<algorithm>
    #include<vector>
    #include<stack>
    #include<cstring>
    #include<queue>
    #include<set>
    #include<string>
    #include<map>
    #include <time.h>
    #define PI acos(-1)
    using namespace std;
    typedef long long ll;
    typedef double db;
    #define LL(x) (x<<1)
    #define RR(x) (x<<1|1)
    const int maxn = 100000+5;
    const ll maxm = 1e7;
    const ll base = 2333;
    const ll inf = 1e17;
    const db eps = 1e-8;
    const ll mod = 1e9+13;
    int a[maxn], n, m, cnt, head[maxn];
    int in[maxn], out[maxn], xx, idx[maxn];
    ll dis[maxn];  int ca=1;
    struct Edge{
        int v, next;
    }e[maxn*2];
    struct TREE{
        int l, r;  ll val, maxx, lazy;
        void fun(ll tmp) {
            lazy+=tmp;
            val+=(r-l+1)*tmp;
            maxx+=tmp;
        }
    }tree[maxn*4];
    
    void init() {
        xx=cnt=0;
        memset(head, -1, sizeof(head));
    }
    void add(int u, int v) {
        e[cnt].v=v, e[cnt].next=head[u];
        head[u]=cnt++;
    }
    void dfs(int u, int fa) {
        in[u]=++xx;  idx[xx]=u;
        for (int i=head[u]; ~i; i=e[i].next) {
            int v=e[i].v;
            if (v==fa)  continue;
            dis[v]=dis[u]+a[v];
            dfs(v, u);
        }
        out[u]=xx;
    }
    void pushup(int id) {
        tree[id].val=tree[id*2].val+tree[id*2+1].val;
        tree[id].maxx=max(tree[id*2].maxx, tree[id*2+1].maxx);
    }
    void pushdown(int id) {
        if (tree[id].lazy) {
            tree[id*2].fun(tree[id].lazy);
            tree[id*2+1].fun(tree[id].lazy);
            tree[id].lazy=0;
        }
    }
    void build(int id, int l, int r) {
        tree[id].l=l, tree[id].r=r, tree[id].lazy=0;
        if (l==r) {
            tree[id].val=tree[id].maxx=dis[idx[l]];
        }
        else {
            int mid=(l+r)/2;
            build(id*2, l, mid);
            build(id*2+1, mid+1, r);
            pushup(id);
        }
    }
    void updata(int id, int st, int ed, int val) {
        int l=tree[id].l, r=tree[id].r;
        if (st<=l&&r<=ed)  tree[id].fun(val);
        else {
            pushdown(id);
            int mid=(l+r)/2;
            if (st<=mid)  updata(id*2, st, ed, val);
            if (ed>mid)  updata(id*2+1, st, ed, val);
            pushup(id);
        }
    }
    ll quary(int id, int st, int ed) {
        int l=tree[id].l, r=tree[id].r;
        if (st<=l&&r<=ed)  return tree[id].maxx;
        else {
            pushdown(id);
            int mid=(l+r)/2;
            ll sum1=-inf, sum2=-inf;
            if (st<=mid)  sum1=quary(id*2, st, ed);
            if (ed>mid)  sum2=quary(id*2+1, st, ed);
            pushup(id);
            return max(sum1, sum2);
        }
    }
    void solve() {
        scanf("%d%d", &n, &m);
        init();
        for (int i=0; i<n-1; i++) {
            int u, v;  scanf("%d%d", &u, &v);
            add(u, v);  add(v, u);
        }
        for (int i=0; i<n; i++)  scanf("%d", a+i);
        dis[0]=a[0];
        dfs(0, -1);
        build(1, 1, n);
        printf("Case #%d:
    ", ca++);
        for (int i=0; i<m; i++) {
            int op, x;  scanf("%d%d", &op, &x);
            if (!op) {
                int y;  scanf("%d", &y);
                updata(1, in[x], out[x], y-a[x]);
                a[x]=y;
            }
            else {
                printf("%I64d
    ", quary(1, in[x], out[x]));
            }
        }
    }
    int main() {
        int t = 1;
        //freopen("in.txt","r",stdin);
      //  freopen("gcd.out","w",stdout);
        scanf("%d", &t);
        while(t--)
            solve();
        return 0;
    }
  • 相关阅读:
    NHibernate介绍
    dwr配置文件dwr.xml详解
    架构设计师与SOA
    SOA是什么
    JDK常用命令
    在WPF的WebBrowser控件中抑制脚本错误
    通过编程计算一个游戏的胜率
    在C#中模拟大数乘法
    解决HttpWebRequest首次连接特别慢的问题
    布隆过滤器(Bloom Filter)
  • 原文地址:https://www.cnblogs.com/gggyt/p/7249559.html
Copyright © 2011-2022 走看看