程序博客网 > 阿里云空间多少钱一年

线段树&树状数组总结篇

来源：互联网 发布：阿里云空间多少钱一年 编辑：程序博客网 时间：2024/06/16 04:49

本文系转载，原创请访问：http://blog.csdn.net/zhou_yujia/article/details/51386549

1、 入门题：hdu1166敌兵布阵

单点更新区间求和

线段树：

[cpp] view plain copy print?
  1. using namespace std;  
  2. struct node  
  3. {  
  4.      int l,r,sum;  
  5. }num[maxn*4];  
  6. int m[maxn];  
  7. char s[10];  
  8. int t,n,a,b;  
  9. void build(int root,int l,int r)  
  10. {  
  11.      num[root].l=l;  
  12.      num[root].r=r;  
  13.      if(num[root].l==num[root].r)  
  14.      {  
  15.           num[root].sum=m[l];return;  
  16.      }  
  17.      int mid=(num[root].l+num[root].r)/2;  
  18.      build(root<<1,l,mid);  
  19.      build(root<<1|1,mid+1,r);  
  20.      num[root].sum=num[root<<1].sum+num[root<<1|1].sum;  
  21. }  
  22. void update(int root,int pos,int data)  
  23. {  
  24.      if(num[root].l==num[root].r)  
  25.      {  
  26.           num[root].sum=data;  
  27.           return;  
  28.      }  
  29.      int mid=(num[root].l+num[root].r)/2;  
  30.      if(pos<=mid) update(root<<1,pos,data);  
  31.      else update(root<<1|1,pos,data);  
  32.      num[root].sum=num[root<<1].sum+num[root<<1|1].sum;  
  33. }  
  34. int  query(int root,int L,int R)  
  35. {  
  36.      if(L<=num[root].l&&R>=num[root].r)  
  37.      {  
  38.           return num[root].sum;  
  39.      }  
  40.      int mid=(num[root].l+num[root].r)/2;  
  41.      int ans=0;  
  42.      if(L<=mid) ans+=query(root<<1,L,R);  
  43.      if(R>mid) ans+=query(root<<1|1,L,R);  
  44.      return ans;  
  45. }  
  46.   
  47. int main()  
  48. {  
  49.      freopen("data.in.txt","r",stdin);  
  50.      while(~scanf("%d",&t))  
  51.      {  
  52.           int cnt=1;  
  53.           while(t--)  
  54.           {  
  55.                scanf("%d",&n);  
  56.                for(int i=1;i<=n;i++)      scanf("%d",&m[i]);  
  57.                build(1,1,maxn);  
  58.                printf("Case %d: \n",cnt++);  
  59.                while(~scanf("%s",s))  
  60.                {  
  61.                     if(strcmp(s,"END")==0) break;  
  62.                     scanf("%d%d",&a,&b);  
  63.                     if(strcmp(s,"Add")==0)  
  64.                     {  
  65.                          m[a]+=b;  
  66.                          update(1,a,m[a]);  
  67.                     }  
  68.                     if(strcmp(s,"Sub")==0)  
  69.                     {  
  70.                          m[a]-=b;  
  71.                          update(1,a,m[a]);  
  72.                     }  
  73.                     if(strcmp(s,"Query")==0)  
  74.                     {  
  75.                          if(a>b) swap(a,b);  
  76.                          printf("%d\n",query(1,a,b));  
  77.                     }  
  78.                }  
  79.           }  
  80.      }  
  81.      return 0;  
  82. }  

树状数组

[cpp] view plain copy print?
  1. int lowbit(int i)  
  2. {  
  3.      return i&(-i);  
  4. }  
  5. void update(int i,int x)  
  6. {  
  7.      while(i<=n)  
  8.      {  
  9.           tree[i]+=x;  
  10.           i+=lowbit(i);  
  11.      }  
  12. }  
  13. int query(int n)  
  14. {  
  15.      int sum=0;  
  16.      while(n>0)  
  17.      {  
  18.           sum+=tree[n];  
  19.           n-=lowbit(n);  
  20.      }  
  21.      return sum;  
  22. }  
  23. int main()  
  24. {  
  25.     int t;scanf("%d",&t);  
  26.     for(int cas=1;cas<=T;cas++)  
  27.     {  
  28.          memset(tree,0,sizeof(tree));  
  29.          scanf("%d",&n);  
  30.          for(int i=1;i<=n;i++)  
  31.          {  
  32.               int x;scanf("%d",&x);  
  33.               update(i,x);  
  34.          }  
  35.          char str[10];  
  36.          while(~scanf("%s",str)&&strcmp(str,"END"))  
  37.          {  
  38.              inta,b;scanf("%d%d",&a,&b);  
  39.               if(str[0]=='Q')  
  40. printf("%d\n",Query(b)-Query(a-1));  
  41.             elseif(str[0]=='S')Update(a,-b);  
  42.               else Update(a,b);  
  43.          }  
  44.     }  
  45.     return 0;  
  46. }  

2、懒惰标记入门题：hdu1698 just a hook

区间更新求总和 ，不写懒惰标记也能过==看我的这篇

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include <cstdio>  
  3. using namespace std;  
  4. const int maxn=1e5+5;  
  5. #define l_t 2*root  
  6. #define r_t 2*root+1  
  7. struct node{  
  8.     int l,r,val,tag;  //tag:1 2 3 -1  
  9.     int mid(){  return (l+r)/2; }  
  10. }tree[maxn<<2];  
  11. void update(int a,int b,int root,int tg){  
  12.     if(a==tree[root].l&&b==tree[root].r){  
  13.          tree[root].tag=tg;  
  14.          tree[root].val=(tree[root].r-tree[root].l+1)*tg;  
  15.          return ;  
  16.     }  
  17.     if(tree[root].tag!=-1){  
  18.          tree[l_t].tag=tree[root].tag;  
  19.          tree[l_t].val=tree[l_t].tag*(tree[l_t].r-tree[l_t].l+1);  
  20.          tree[r_t].tag=tree[root].tag;  
  21.          tree[r_t].val=tree[r_t].tag*(tree[r_t].r-tree[r_t].l+1);  
  22.          tree[root].tag=-1;  
  23.     }  
  24.     int m=tree[root].mid();  
  25.     if(a>m) update(a,b,r_t,tg);  
  26.     else if(b<=m) update(a,b,l_t,tg);  
  27.     else {  
  28.         update(a,m,l_t,tg);  
  29.         update(m+1,b,r_t,tg);  
  30.     }  
  31.     tree[root].val=tree[l_t].val+tree[r_t].val;  
  32. }  
  33. void build(int l,int r,int root){  
  34.     tree[root].l=l;  
  35.     tree[root].r=r;  
  36.     tree[root].tag=-1;  
  37.     if(r==l){  
  38.         tree[root].val=1;  
  39.         return ;  
  40.     }  
  41.     int mid=(l+r)/2;  
  42.     build(l,mid,l_t);  
  43.     build(mid+1,r,r_t);  
  44.     tree[root].val=tree[l_t].val+tree[r_t].val;  
  45. }  
  46.   
  47. int main()  
  48. {  
  49.     //freopen("cin.txt","r",stdin);  
  50.     int ca;  
  51.     cin>>ca;  
  52.     for(int k=1;k<=ca;k++){  
  53.         int n,q;  
  54.         scanf("%d",&n);  
  55.         build(1,n,1);  
  56.         scanf("%d",&q);  
  57.         for(int i=0;i<q;i++){  
  58.             int x,y,z;  
  59.             scanf("%d%d%d",&x,&y,&z);  
  60.             update(x,y,1,z);  
  61.         }  
  62.         printf("Case %d: The total value of the hook is %d.\n",k,tree[1].val);  
  63.     }  
  64.     return 0;  
  65. }  

3、线段树求逆序数：hdu1394minimum inversion number

求依次把数组的第一个数放到最后，组成的新数组的逆序数的最小值。给定的数是0~n-1之间的。

整体的思路是先求出来把每个数插到最后相对比放在原始位置变动的多少，这个相对值+原始的逆序数值就是最终结果。关键是相对值如何求？由于每个数都是0~n-1之间，省了离散化这个步骤，我们知道n-1-a[i]是比这个数大的个数，也就是这个数后面有多少个比大的，（我们让每步的初始状态是读入的这个数是第一个）那么我们把这个刚刚读入的数放到最后，里外里就为逆序数贡献了n-1-2*a[i]，中间的不用管，反正没为逆序数总值做贡献。   可以看我这篇：感觉写的还没有我这个详细呢

hdu1394Minimum Inversion Number逆序数

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. using namespace std;  
  4. struct node  
  5. {  
  6.     int l,r,sum;  
  7. }tree[150000];  
  8. void build(int root,int l,int r)  
  9. {  
  10.     tree[root].l=l,tree[root].r=r;  
  11.     if(tree[root].l==tree[root].r)  
  12.     {  
  13.         tree[root].sum=0;  
  14.         return;  
  15.     }  
  16.     build(root<<1,l,(l+r)/2);  
  17.     build(root<<1|1,(l+r)/2+1,r);  
  18.     tree[root].sum=0;  
  19. }  
  20. void update(int root,int pos,int val)  
  21. {  
  22.     if(tree[root].l==tree[root].r)  
  23.     {  
  24.         tree[root].sum=val;  
  25.         return;  
  26.     }  
  27.     int mid=(tree[root].l+tree[root].r)/2;  
  28.     if(pos<=mid) update(root<<1,pos,val);  
  29.     else update(root<<1|1,pos,val);///这里写成了build 你脑子在想什么==  
  30.     tree[root].sum=min(tree[root<<1].sum,tree[root<<1|1].sum);  
  31. }  
  32. int main()  
  33. {  
  34.     //freopen("cin.txt","r",stdin);  
  35.     int n,a[12],ret;///  
  36.     while(~scanf("%d",&n))  
  37.     {  
  38.         build(1,1,n);  
  39.         ret=0;  
  40.         for(int i=0;i<n;i++) {  
  41.             scanf("%d",&a[i]);  
  42.             ret+=(n-1-2*a[i]);  
  43.             update(1,a[i],ret);  
  44.         }  
  45.         ret=0;  
  46.         for(int i=0;i<n-1;i++)  
  47.         {  
  48.             for(int j=i+1;j<n;j++)  
  49.             {  
  50.                 if(a[i]>a[j]) ret++;  
  51.             }  
  52.         }  
  53.         printf("%d\n",ret+tree[1].sum);  
  54.     }  
  55.     return 0;  
  56. }  

4、懒惰标记经典题：poj3667hotel

 题意：旅馆有人开房：给定区间长度，尽量从小号开始给连续的房间，询问最左侧的房间号、有人退房：给定左区间端点和区间长度

现在看来也没什么好说的，注意题中要求的顺序

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstring>  
  3. #include<cstdio>  
  4. using namespace std;  
  5. struct node  
  6. {  
  7.     int l,r,lsum,rsum,sum,loop;  
  8. }tree[400000];  
  9. int Find_Max(int a,int b ,int c)  
  10. {  
  11.     int maxn=a>b?a:b;  
  12.     maxn=maxn>c?maxn:c;  
  13.     return maxn;  
  14. }  
  15. void push_down(int t)  
  16. {  
  17.     if(tree[t].loop!=-1)  
  18.     {  
  19.         tree[t<<1].loop=tree[t<<1|1].loop=tree[t].loop;  
  20.         if(tree[t].loop)  
  21.         {  
  22.             tree[t<<1].lsum=tree[t<<1].rsum=tree[t<<1].sum=0;  
  23.             tree[t<<1|1].lsum=tree[t<<1|1].rsum=tree[t<<1|1].sum=0;  
  24.         }  
  25.         else //if(tree[t].loop==0)  
  26.         {  
  27.             tree[t<<1].lsum=tree[t<<1].rsum=tree[t<<1].sum=tree[t<<1].r-tree[t<<1].l+1;  
  28.             tree[t<<1|1].lsum=tree[t<<1|1].rsum=tree[t<<1|1].sum=tree[t<<1|1].r-tree[t<<1|1].l+1;  
  29.         }  
  30.         tree[t].loop=-1;  
  31.     }  
  32. }  
  33. void push_up(int l,int r,int t)  
  34. {  
  35.     tree[t].lsum=tree[t<<1].lsum;  
  36.     tree[t].rsum=tree[t<<1|1].rsum;  
  37.     int x=(l+r)/2;  
  38.     if(tree[t].lsum==x-l+1) tree[t].lsum+=tree[t<<1|1].lsum;  
  39.     if(tree[t].rsum==r-x) tree[t].rsum+=tree[t<<1].rsum;  
  40.     tree[t].sum=Find_Max(tree[t<<1].sum,tree[t<<1|1].sum,tree[t<<1].rsum+tree[t<<1|1].lsum);///  
  41. }  
  42. void build(int l,int r,int t)  
  43. {  
  44.     tree[t].l=l,tree[t].r=r;  
  45.     tree[t].sum=tree[t].lsum=tree[t].rsum=r-l+1;  
  46.     tree[t].loop=-1;  
  47.     if(l==r) return;  
  48.     int x=(l+r)/2;  
  49.     build(l,x,2*t);  
  50.     build(x+1,r,2*t+1);  
  51. }  
  52. void update_tree(int l,int r,int t,int cnt)  
  53. {  
  54.     if(tree[t].l==l&&tree[t].r==r)  
  55.     {  
  56.         tree[t].loop=cnt;  
  57.         if(cnt) tree[t].lsum=tree[t].rsum=tree[t].sum=0;  
  58.         else tree[t].lsum=tree[t].rsum=tree[t].sum=r-l+1;  
  59.         return;  
  60.     }  
  61.     push_down(t);  
  62.     int x=(tree[t].l+tree[t].r)/2;  
  63.     if(x>=r) update_tree(l,r,2*t,cnt);  
  64.     else if(x+1<=l) update_tree(l,r,2*t+1,cnt);  
  65.     else  
  66.     {  
  67.         update_tree(l,x,2*t,cnt);  
  68.         update_tree(x+1,r,t<<1|1,cnt);  
  69.     }  
  70.     push_up(tree[t].l,tree[t].r,t);  
  71. }  
  72. int Query(int l,int r,int t,int cnt)  
  73. {  
  74.     if(l==r) return 1;  
  75.     int x=(l+r)/2;  
  76.     push_down(t);  
  77.     if(tree[t<<1].sum>=cnt) return Query(l,x,t<<1,cnt);  
  78.     else if(tree[t<<1].rsum+tree[t<<1|1].lsum>=cnt)return x-tree[t<<1].rsum+1;  
  79.     else return Query(x+1,r,t<<1|1,cnt);  
  80. }  
  81. int main()  
  82. {  
  83.    // freopen("cin.txt","r",stdin);  
  84.     int i,j,n,m,x,y;  
  85.     while(~scanf("%d%d",&n,&m))  
  86.     {  
  87.         build(1,n,1);  
  88.         while(m--)  
  89.         {  
  90.             scanf("%d",&i);  
  91.             if(i==1)  
  92.             {  
  93.                 scanf("%d",&x);  
  94.                 if(tree[1].sum<x)  
  95.                 {  
  96.                     printf("0\n");  
  97.                     continue;  
  98.                 }  
  99.                 y=Query(1,n,1,x);  
  100.                 printf("%d\n",y);  
  101.                 update_tree(y,y+x-1,1,1);  
  102.             }  
  103.             else  
  104.             {  
  105.                 scanf("%d%d",&x,&y);  
  106.                 update_tree(x,x+y-1,1,0);  
  107.             }  
  108.         }  
  109.     }  
  110.     return 0;  
  111. }  

6、 hdu2795billbroad

展板上贴窄公告,优先选择最上最左的，输出行号。和刚刚那个题一样，也需要注意优先级的问题，更要注意长度是有限制的，不像上面那个是一维的

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. #define maxn 200000  
  5. using namespace std;  
  6. struct bill  
  7. {  
  8.      int l,r,ans;  
  9. }num[maxn*4];  
  10. int m;  
  11. int w,h,n;  
  12. void build(int L,int R,int root)  
  13. {  
  14.      num[root].l=L;  
  15.      num[root].r=R;  
  16.      num[root].ans=w;  
  17.      if(L==R) return;  
  18.      build(L,(L+R)/2,root<<1);  
  19.      build((L+R)/2+1,R,root<<1|1);  
  20.     //  
  21. }  
  22. int update(int root,int m)  
  23. {  
  24.      if(num[root].ans<m) return -1;  
  25.      if(num[root].l==num[root].r)  
  26.      {  
  27.           num[root].ans-=m;  
  28.           return num[root].l;  
  29.      }  
  30.      int ret;  
  31.      if(num[root<<1].ans>=m) ret=update(root<<1,m);  
  32.      else ret=update(root<<1|1,m);  
  33.      num[root].ans=max(num[root<<1].ans,num[root<<1|1].ans);  
  34.      return ret;  
  35. }  
  36. int main()  
  37. {  
  38.     //freopen("data.in.txt","r",stdin);  
  39.     while(~scanf("%d%d%d",&h,&w,&n))  
  40.     {  
  41.          build(1,min(n,h),1);  
  42.          for(int i=1;i<=n;i++)  
  43.          {  
  44.               scanf("%d",&m);  
  45.               printf("%d\n",update(1,m));  
  46.          }  
  47.     }  
  48.     return 0;  
  49. }  
===================================上面都是去年暑假做的水题，下面的是今年的练习，难度明显高了一个level======================================
7、树状数组、线段树求逆序数：poj2299Ultra-QuickSort 求交换相邻的数字使之有序的次数

简而言之就是交换次数==每个数前面比他大的总和 详细的代码见：

poj2299Ultra-QuickSort【树状数组求逆序数、离散化】、【归并排序模板】

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. #include<algorithm>  
  5. using namespace std;  
  6. #define maxn 500005  
  7. int n;  
  8. int c[maxn],b[maxn];  
  9. struct Node  
  10. {  
  11.     int index,v;  
  12. }node[maxn];  
  13. bool cmp(Node a,Node b)  
  14. {  
  15.     return a.v<b.v;  
  16. }  
  17. int lowbit(int x)  
  18. {  
  19.     return x&(-x);  
  20. }  
  21. void add(int i,int val)  
  22. {  
  23.     while(i<=n)  
  24.     {  
  25.         c[i]+=val;  
  26.         i+=lowbit(i);  
  27.     }  
  28. }  
  29. int sum(int i)  
  30. {  
  31.     int s=0;  
  32.     while(i>0)  
  33.     {  
  34.         s+=c[i];  
  35.         i-=lowbit(i);  
  36.     }  
  37.     return s;  
  38. }  
  39. int main()  
  40. {  
  41.     //freopen("cin.txt","r",stdin);  
  42.     while(~scanf("%d",&n)&&n)  
  43.     {  
  44.         for(int i=1;i<=n;i++)  
  45.         {  
  46.             scanf("%d",&node[i].v);  
  47.             node[i].index=i;  
  48.         }  
  49.         memset(b,0,sizeof(b));  
  50.         memset(c,0,sizeof(c));  
  51.         sort(node+1,node+1+n,cmp);  
  52.         b[node[1].index]=1;  
  53.         for(int i=2;i<=n;i++)  
  54.         {  
  55.             b[node[i].index]=i;  
  56.         }  
  57.         long long ans=0;  
  58.         for(int i=1;i<=n;i++)  
  59.         {  
  60.             add(b[i],1);  
  61.             ans+=(i-sum(b[i]));  
  62.         }  
  63.         printf("%I64d\n",ans);  
  64.     }  
  65.     return 0;  
  66. }  
线段树做法

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. #include<algorithm>  
  5. using namespace std;  
  6. #define maxn 500005  
  7. struct node  
  8. {  
  9.     int id,val;  
  10. }num[maxn];  
  11. bool cmp(node n1,node n2)  
  12. {  
  13.     return n1.val<n2.val;  
  14. }  
  15. struct Tree  
  16. {  
  17.     int l,r,tot;  
  18. }tree[maxn<<2];  
  19. int rank[maxn];  
  20. void build(int rt,int l,int r)  
  21. {  
  22.     tree[rt].l=l;tree[rt].r=r;  
  23.     tree[rt].tot=0;  
  24.     if(l==r)return;  
  25.     int mid=(l+r)/2;  
  26.     build(rt<<1,l,mid);  
  27.     build(rt<<1|1,mid+1,r);  
  28. }  
  29. void update(int rt,int x)  
  30. {  
  31.     tree[rt].tot++;  
  32.     if(tree[rt].l==tree[rt].r&&tree[rt].l==x)  
  33.         return;  
  34.     int mid=(tree[rt].l+tree[rt].r)/2;  
  35.     if(x<=mid)update(rt<<1,x);  
  36.     else update(rt<<1|1,x);  
  37. }  
  38. int query(int rt,int l,int r)  
  39. {  
  40.     int sum=0;  
  41.     if(tree[rt].l==l&&tree[rt].r==r)  
  42.         return tree[rt].tot;  
  43.     int mid=(tree[rt].l+tree[rt].r)/2;  
  44.     if(r<=mid) sum=query(rt<<1,l,r);  
  45.     else if(l>mid) sum=query(rt<<1|1,l,r);  
  46.     else  
  47.         sum=query(rt<<1,l,mid)+query(rt<<1|1,mid+1,r);  
  48.     return sum;  
  49. }  
  50. int main()  
  51. {  
  52.    // freopen("cin.txt","r",stdin);  
  53.     int n;  
  54.     while(~scanf("%d",&n)&&n)  
  55.     {  
  56.         for(int i=0;i<n;i++)  
  57.         {  
  58.             scanf("%d",&num[i].val);  
  59.             num[i].id=i;  
  60.         }  
  61.         sort(num,num+n,cmp);  
  62.         for(int i=0;i<n;i++) rank[num[i].id]=i+1;  
  63.         build(1,0,n+1);  
  64.         long long sum=0;  
  65.         for(int i=0;i<n;i++)  
  66.         {  
  67.             int x = rank[i];  
  68.             sum+=query(1,x+1,n+1);  
  69.             update(1,x);  
  70.         }  
  71.         printf("%I64d\n",sum);  
  72.     }  
  73.     return 0;  
  74. }  


、树状数组区间更新单点求值hdu4267点击我看博客

A Simple Problem with Integers

[cpp] view plain copy print?
  1. #include<cstdio>  
  2. #include<cstring>  
  3. using namespace std;  
  4. int n,q,c[12][12][50010],tmp[50010];  
  5. int lowbit(int i)  
  6. {  
  7.     return i&(-i);  
  8. }  
  9. void add(int t1,int t2,int i,int x)  
  10. {  
  11.     while(i>0)  
  12.     {  
  13.         c[t1][t2][i]+=x;  
  14.         i-=lowbit(i);  
  15.     }  
  16. }  
  17. int query(int t1,int t2,int x)  
  18. {  
  19.     int s=0;  
  20.     while(x<=n)  
  21.     {  
  22.         s+=c[t1][t2][x];  
  23.         x+=lowbit(x);  
  24.     }  
  25.     return s;  
  26. }  
  27. int main()  
  28. {  
  29.     //freopen("cin.txt","r",stdin);  
  30.     while(~scanf("%d",&n))  
  31.     {  
  32.         memset(c,0,sizeof(c));  
  33.         int a,b,k,e,m;  
  34.         for(int i=1;i<=n;i++)  
  35.         {  
  36.             scanf("%d",&tmp[i]);  
  37.         }  
  38.    //     for(int i=1;i<=n;i++) printf("%d  ",c[i]);  
  39.         scanf("%d",&q);  
  40.         while(q--)  
  41.         {  
  42.             scanf("%d",&m);  
  43.             if(m==1)  
  44.             {  
  45.                 scanf("%d%d%d%d",&a,&b,&k,&e);  
  46.                 int num=(b-a)/k;  
  47.                 int s=a%k;  
  48.                 add(k,s,a-1,-e);  
  49.                 add(k,s,b,e);  
  50.             }  
  51.             else  
  52.             {  
  53.                 scanf("%d",&a);  
  54.                 int sum=tmp[a];  
  55.                 for(int i=1;i<=10;i++)  
  56.                 {  
  57.                     sum+=query(i,a%i,a);  
  58.                 }  
  59.                 printf("%d\n",sum);  
  60.             }  
  61.         }  
  62.     }  
  63.     return 0;  
  64. }  

9.二维树状数组 单点更新 区间求值

poj1195Mobile phones【二维树状数组。单点更新/区间求和】

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. using namespace std;  
  5. int tree[1050][1050],n,t,x,y,a,l,b,r,cmd;  
  6. int lowbit(int i)  
  7. {  
  8.     return i&(-i);  
  9. }  
  10. void update(int x,int y,int a)  
  11. {  
  12.     for(int i=x;i<=n;)  
  13.     {  
  14.         for(int j=y;j<=n;)  
  15.         {  
  16.             tree[i][j]+=a;  
  17.             j+=lowbit(j);  
  18.         }  
  19.         i+=lowbit(i);  
  20.     }  
  21. }  
  22. int query(int x,int y)  
  23. {  
  24.     int s=0;  
  25.     for(int i=x;i>0;i-=lowbit(i))  
  26.     {  
  27.         for(int j=y;j>0;j-=lowbit(j))  
  28.         s+=tree[i][j];  
  29.     }  
  30.     return s;  
  31. }  
  32. int main()  
  33. {  
  34.    // freopen("cin.txt","r",stdin);  
  35.     scanf("%d%d",&cmd,&n);  
  36.     memset(tree,0,sizeof(tree));  
  37.     while(~scanf("%d",&cmd))  
  38.     {  
  39.         if(cmd==3) break;  
  40.         if(cmd==1)  
  41.         {  
  42.             scanf("%d%d%d",&x,&y,&a);  
  43.             x++;y++;  
  44.             update(x,y,a);  
  45.         }  
  46.         else  
  47.         {  
  48.             scanf("%d%d%d%d",&l,&b,&r,&t);  
  49.             int sum=0;  
  50.             l++;b++;r++;t++;  
  51.             sum=query(r,t)+query(l-1,b-1)-query(l-1,t)-query(r,b-1);  
  52.             printf("%d\n",sum);  
  53.         }  
  54.     }  
  55.     return 0;  
  56. }  

10、离散化区间更新线段树。贴布条

poj2528Mayor's posters【离散化线段树区间更新】

既然说了新贴的贴到原来的上面，那么新贴的染新色最后看有多少颜色就得了啊

hdu4027Can you answer these queries?【线段树区间更新区间求和】更新是区间内的值都开根号


[cpp] view plain copy print?
  1. #include <cstdio>  
  2. #include <algorithm>  
  3. #include<cstring>  
  4. using namespace std;  
  5. #define lson l , m , rt << 1  
  6. #define rson m + 1 , r , rt << 1 | 1  
  7. const int maxn = 11111;  
  8. bool hash[maxn];  
  9. int li[maxn] , ri[maxn];  
  10. int X[maxn*3];  
  11. int col[maxn<<4];  
  12. int cnt;  
  13. void push_down(int rt)  
  14. {  
  15.     if(col[rt]!=-1)  
  16.     {  
  17.         col[rt<<1]=col[rt<<1|1]=col[rt];  
  18.         col[rt]=-1;  
  19.     }  
  20. }  
  21. void update(int l,int r,int c,int L,int R,int rt)  
  22. {  
  23.     if(l<=L&&R<=r)  
  24.     {  
  25.         col[rt]=c;  
  26.         return;  
  27.     }  
  28.     push_down(rt);  
  29.     int mid=(L+R)/2;  
  30.     if(l<=mid)update(l,r,c,L,mid,rt<<1);  
  31.     if(r>mid) update(l,r,c,mid+1,R,rt<<1|1);///  
  32. }  
  33. void query(int l,int r,int rt)  
  34. {  
  35.     if(col[rt]!=-1)  
  36.     {  
  37.         if(!hash[col[rt]]) cnt++;  
  38.         hash[col[rt]]=1;  
  39.         return;  
  40.     }  
  41.     if(l==r) return;  
  42.     int m=(l+r)/2;  
  43.     query(l,m,rt<<1);  
  44.     query(m+1,r,rt<<1|1);  
  45. }  
  46. int Bin(int num,int R)  
  47. {  
  48.     int l=0,r=R-1,mid;  
  49.     while(l<=r)  
  50.     {  
  51.         mid=(l+r)/2;  
  52.         if(X[mid]==num) return mid;  
  53.         if(X[mid]<num) l=mid+1;  
  54.         else r=mid-1;  
  55.     }  
  56. }  
  57. int main()  
  58. {  
  59.    // freopen("cin.txt","r",stdin);  
  60.     int t,n;  
  61.     scanf("%d",&t);  
  62.     while(t--)  
  63.     {  
  64.         int nn=0,m=1;  
  65.         scanf("%d",&n);  
  66.         for(int i=0;i<n;i++) scanf("%d%d",&li[i],&ri[i]),X[nn++]=li[i],X[nn++]=ri[i];  
  67.         sort(X,X+nn);  
  68.         for(int i=1;i<nn;i++)  if(X[i]!=X[i-1]) X[m++]=X[i];  
  69.         for(int i=m-1;i>0;i--)if(X[i]!=X[i-1]+1)X[m++]=X[i-1]+1;  
  70.         sort(X,X+m);  
  71.         memset(col,-1,sizeof(col));  
  72.         for(int i=0;i<n;i++)  
  73.         {  
  74.             int l=Bin(li[i],m);  
  75.             int r=Bin(ri[i],m);  
  76.             update(l,r,i,0,m,1);  
  77.         }  
  78.         cnt=0;  
  79.         memset(hash , false , sizeof(hash));  
  80.         query(0,m,1);  
  81.         printf("%d\n",cnt);  
  82.     }  
  83.     return 0;  
  84. }  

11、线段树区间更新区间求和 区间更新是开根号，注意算到一就结束了啊

hdu4027Can you answer these queries?【线段树区间更新区间求和】

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. #include<cmath>  
  5. using namespace std;  
  6. struct node  
  7. {  
  8.     int l,r;  
  9.     long long sum;  
  10. }num[400005];  
  11. int n,m,t;  
  12. void build(int l,int r,int i)  
  13. {  
  14.     num[i].l=l;  
  15.     num[i].r=r;  
  16.     if(l==r)  
  17.     {  
  18.         scanf("%I64d",&num[i].sum);  
  19.         return;  
  20.     }  
  21.     int mid=(l+r)/2;  
  22.     build(l,mid,i<<1);  
  23.     build(mid+1,r,i<<1|1);  
  24.     num[i].sum=num[i<<1].sum+num[i<<1|1].sum;  
  25. }  
  26. void update(int l,int r,int rt)  
  27. {  
  28.     if(num[rt].sum==(num[rt].r-num[rt].l+1)&&l==num[rt].l&&r==num[rt].r)  
  29.      return;  
  30.     if(num[rt].l==num[rt].r)  
  31.     {  
  32.         num[rt].sum=(int)sqrt(num[rt].sum*1.0);  
  33.         return;  
  34.     }  
  35.     int mid=(num[rt].l+num[rt].r)/2;  
  36.     if(r<=mid) update(l,r,rt<<1);  
  37.     else if(l>mid) update(l,r,rt<<1|1);  
  38.     else  
  39.     {  
  40.         update(l,mid,rt<<1);  
  41.         update(mid+1,r,rt<<1|1);  
  42.     }  
  43.     num[rt].sum=num[rt<<1].sum+num[rt<<1|1].sum;  
  44. }  
  45. long long query(int l,int r,int rt)  
  46. {  
  47.     if(l==num[rt].l&&r==num[rt].r) return num[rt].sum;  
  48.     long long ans=0;  
  49.     int mid=(num[rt].l+num[rt].r)/2;  
  50.     if(r<=mid) ans+=query(l,r,rt<<1);  
  51.     else if(l>mid) ans+=query(l,r,rt<<1|1);  
  52.     else  
  53.     {  
  54.         ans+=query(l,mid,rt<<1);  
  55.         ans+=query(mid+1,r,rt<<1|1);  
  56.     }  
  57.     return ans;  
  58. }  
  59. int main()  
  60. {  
  61.    // freopen("cin.txt","r",stdin);  
  62.     int cas=1;  
  63.     while(~scanf("%d",&n))  
  64.     {  
  65.         printf("Case #%d:\n",cas++);  
  66.         build(1,n,1);  
  67.         scanf("%d",&m);  
  68.         while(m--)  
  69.         {  
  70.             int a,b;  
  71.             scanf("%d%d%d",&t,&a,&b);  
  72.             if(a>b) swap(a,b);  
  73.             if(t==0) update(a,b,1);  
  74.             else printf("%I64d\n",query(a,b,1));  
  75.         }  
  76.         puts("");  
  77.     }  
  78.     return 0;  
  79. }  

12、线段树求第k大、第k小、维护区间最值

hdu2852KiKi's K-Number【线段树第k小】

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. using namespace std;  
  5. struct node  
  6. {  
  7.     int l,r,tot;  
  8. }num[800000];  
  9. int tmp;  
  10. void build(int rt,int l,int r)  
  11. {  
  12.     num[rt].l=l;num[rt].r=r;  
  13.     if(l==r)  
  14.     {  
  15.         num[rt].tot=0;  
  16.         return;  
  17.     }  
  18.     int mid=(l+r)/2;  
  19.     build(rt<<1,l,mid);  
  20.     build(rt<<1|1,mid+1,r);  
  21.     num[rt].tot=0;  
  22. }  
  23. void update(int rt,int pos,bool flag)///1:add 0:del  
  24. {  
  25.     if(num[rt].l==num[rt].r)  
  26.     {  
  27.         if(flag) num[rt].tot++;  
  28.         else num[rt].tot--;  
  29.         return;///!!!  
  30.     }  
  31.     int mid=(num[rt].l+num[rt].r)/2;  
  32.     if(pos<=mid) update(rt<<1,pos,flag);  
  33.     else update(rt<<1|1,pos,flag);  
  34.     num[rt].tot=num[rt<<1].tot+num[rt<<1|1].tot;  
  35. }  
  36. int sum(int rt,int l,int r)  
  37. {  
  38.     if(num[rt].r<l||num[rt].l>r) return 0;  
  39.     if(num[rt].r<=r&&num[rt].l>=l) return num[rt].tot;  
  40.     int mid=(num[rt].l+num[rt].r)/2;  
  41.     return sum(rt<<1,l,r)+sum(rt<<1|1,l,r);  
  42. }  
  43. void query(int rt,int ans)  
  44. {  
  45.     if(ans>num[rt<<1].tot)  
  46.     {  
  47.         if(num[rt].l==num[rt].r)  
  48.         {  
  49.             tmp=num[rt].r;  
  50.             return;  
  51.         }  
  52.         query(rt<<1|1,ans-num[rt<<1].tot);  
  53.     }  
  54.     else query(rt<<1,ans);  
  55. }  
  56. int main()  
  57. {  
  58.   //  freopen("cin.txt","r",stdin);  
  59.     int q,a,b,m;  
  60.     while(~scanf("%d",&m))  
  61.     {  
  62.         build(1,1,100000);  
  63.         while(m--)  
  64.         {  
  65.             scanf("%d",&q);  
  66.             if(q==0)  
  67.             {  
  68.                 scanf("%d",&a);  
  69.                 update(1,a,1);  
  70.             }  
  71.             else if(q==1)  
  72.             {  
  73.                 scanf("%d",&a);  
  74.                 if(sum(1,1,a)==sum(1,1,a-1)) printf("No Elment!\n");  
  75.                 else update(1,a,0);  
  76.             }  
  77.             else  
  78.             {  
  79.                 scanf("%d%d",&a,&b);  
  80.                 {  
  81.                     int t=sum(1,1,a);  
  82.                     tmp=0;  
  83.                     query(1,t+b);  
  84.                     if(tmp==100000)printf("Not Find!\n");  
  85.                     else printf("%d\n",tmp);  
  86.                 }  
  87.             }  
  88.         }  
  89.     }  
  90.     return 0;  
  91. }  
hdu4006The kth great number【线段树第k大】


[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. using namespace std;  
  5. #define maxn 1000000  
  6. struct node  
  7. {  
  8.     int l,r,tot;  
  9. }num[maxn*8];  
  10. void build(int rt,int l,int r)  
  11. {  
  12.     num[rt].l=l;num[rt].r=r;  
  13.     if(l==r)  
  14.     {  
  15.         num[rt].tot=0;  
  16.         return;  
  17.     }  
  18.     int mid=(l+r)/2;  
  19.     build(rt<<1,l,mid);  
  20.     build(rt<<1|1,mid+1,r);  
  21.     num[rt].tot=0;  
  22. }  
  23. void update(int rt,int pos)  
  24. {  
  25.     if(num[rt].l==num[rt].r)  
  26.     {  
  27.         if(num[rt].l==pos) num[rt].tot++;  
  28.         return;  
  29.     }  
  30.     int mid=(num[rt].l+num[rt].r)/2;  
  31.     if(pos<=mid) update(rt<<1,pos);  
  32.     else update(rt<<1|1,pos);  
  33.     num[rt].tot=num[rt<<1].tot+num[rt<<1|1].tot;  
  34. }  
  35. int query(int rt,int k)  
  36. {  
  37.     if(k>num[rt<<1|1].tot)  
  38.     {  
  39.         if(num[rt].l==num[rt].r)  
  40.         {  
  41.             return num[rt].l;  
  42.         }  
  43.         query(rt<<1,k-num[rt<<1|1].tot);  
  44.     }  
  45.     else query(rt<<1|1,k);  
  46. }  
  47. int main()  
  48. {  
  49.    // freopen("cin.txt","r",stdin);  
  50.     int n,k,a;  
  51.     char st[3];  
  52.     while(~scanf("%d%d",&n,&k))  
  53.     {  
  54.         build(1,1,maxn);  
  55.         while(n--)  
  56.         {  
  57.             scanf("%s",st);  
  58.             if(st[0]=='I')  
  59.             {  
  60.                 scanf("%d",&a);  
  61.                 update(1,a);  
  62.             }  
  63.             else printf("%d\n",query(1,k));  
  64.         }  
  65.     }  
  66.     return 0;  
  67. }  

poj2823Sliding Window【线段树维护滚动区间最值】    也就是建好树之后区间查询就好

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. using namespace std;  
  5. struct Node  
  6. {  
  7.     int l,r,minn,maxn;  
  8. }num[4000000];  
  9. int cnt[1000000];  
  10. int min(int a,int b){if(a<b)return a;return b;}  
  11. int max(int a,int b){if(a>b)return a;return b;}  
  12. void build(int rt,int l,int r)  
  13. {  
  14.     num[rt].l=l;num[rt].r=r;  
  15.     if(l==r)  
  16.     {  
  17.         num[rt].maxn=cnt[l];  
  18.         num[rt].minn=cnt[l];  
  19.         return;  
  20.     }  
  21.     int mid=(l+r)/2;  
  22.     build(rt<<1,l,mid);  
  23.     build(rt<<1|1,mid+1,r);  
  24.     num[rt].maxn=max(num[rt<<1].maxn,num[rt<<1|1].maxn);  
  25.     num[rt].minn=min(num[rt<<1].minn,num[rt<<1|1].minn);  
  26. }  
  27. int query1(int rt,int l,int r)  
  28. {  
  29.     if(num[rt].l==l&&num[rt].r==r)  
  30.         return num[rt].minn;  
  31.     int mid=(num[rt].l+num[rt].r)/2;  
  32.     if(r<=mid) return query1(rt<<1,l,r);  
  33.     else if(l>mid) return query1(rt<<1|1,l,r);  
  34.     else  
  35.     return min(query1(rt<<1,l,mid),query1(rt<<1|1,mid+1,r));  
  36. }  
  37. int query2(int rt,int l,int r)  
  38. {  
  39.     if(num[rt].l==l&&num[rt].r==r)  
  40.         return num[rt].maxn;  
  41.     int mid=(num[rt].l+num[rt].r)/2;  
  42.     if(r<=mid) return query2(rt<<1,l,r);  
  43.     else if(l>mid) return query2(rt<<1|1,l,r);  
  44.     else  
  45.     return max(query2(rt<<1,l,mid),query2(rt<<1|1,mid+1,r));  
  46. }  
  47. int main()  
  48. {  
  49.    // freopen("cin.txt","r",stdin);  
  50.     int n,k;  
  51.     while(~scanf("%d%d",&n,&k))  
  52.     {  
  53.         for(int i=1;i<=n;i++)scanf("%d",&cnt[i]);  
  54.         build(1,1,n);  
  55.         for(int i=1;i<=n-k;i++)  
  56.             printf("%d ",query1(1,i,i+k-1));  
  57.         printf("%d\n",query1(1,n-k+1,n));  
  58.         for(int i=1;i<=n-k;i++)  
  59.             printf("%d ",query2(1,i,i+k-1));  
  60.         printf("%d\n",query2(1,n-k+1,n));  
  61.     }  
  62.     return 0;  
  63. }  


13、单点更新，区间求最值 

poj2886Who Gets the Most Candies?【线段树单点更新】

除了涉及反素数，其他只是把删掉的小孩当前位清空而已

[cpp] view plain copy print?
  1. #include <iostream>  
  2. #include<cstdio>  
  3. #include<cstring>  
  4. using namespace std;  
  5. #define maxn 500005  
  6. int a[37]={1,2,4,6,12,24,36,48,60,120,180,240,360,720,840,1260,1680,2520,5040,7560,10080,15120,20160,25200,27720,45360,50400,  
  7.            55440,83160,110880,166320,221760,277200,332640,498960,500001};  
  8. int b[37]={1,2,3,4,6,8,9,10,12,16,18,20,24,30,32,36,40,48,60,64,72,80,84,90,96,100,108,120,128,144,160,168,180,192,200,1314521};  
  9. char name[maxn][15];  
  10. int pos[maxn];  
  11. struct node  
  12. {  
  13.     int l,r,tot;  
  14. }num[maxn*4];  
  15. void build(int rt,int l,int r)  
  16. {  
  17.     int mid=(l+r)/2;  
  18.     num[rt].l=l;num[rt].r=r;  
  19.     if(l==r)  
  20.     {  
  21.         num[rt].tot=1;  
  22.         return;  
  23.     }  
  24.    // printf("rt=%d tot=%d\n",rt,num[rt].tot);  
  25.     build(rt<<1,l,mid);  
  26.     build(rt<<1|1,mid+1,r);  
  27.     num[rt].tot=num[rt<<1].tot+num[rt<<1|1].tot;  
  28. }  
  29. int update(int rt,int k)  
  30. {  
  31.     num[rt].tot--;  
  32.     if(num[rt].l==num[rt].r)  
  33.         return num[rt].l;  
  34.     if(k<=num[rt<<1].tot) update(rt<<1,k);  
  35.     else update(rt<<1|1,k-num[rt<<1].tot);  
  36. }  
  37. int main()  
  38. {  
  39.     //freopen("cin.txt","r",stdin);  
  40.     int n,k;  
  41.     while(~scanf("%d%d",&n,&k))  
  42.     {  
  43.         for(int i=1;i<=n;i++)scanf("%s%d",name[i],&pos[i]);  
  44.        // for(int i=1;i<=n;i++)printf("%s% d\n",name[i],pos[i]);  
  45.         build(1,1,n);  
  46.         int cnt=0;  
  47.         while(a[cnt]<=n) cnt++;  
  48.         cnt--;  
  49.         int ps=0;  
  50.         pos[ps]=0;  
  51.       //  printf("tot=%d\n",num[1].tot);  
  52.         for(int i=0;i<a[cnt];i++)  
  53.         {  
  54.             if(pos[ps]>0)  
  55.                 k = ((k + pos[ps] - 2) % num[1].tot + num[1].tot) % num[1].tot + 1;  
  56.             else k = ((k + pos[ps] - 1) % num[1].tot + num[1].tot) % num[1].tot + 1;  
  57.            // cout<<k<<endl;  
  58.             ps=update(1,k);//cout<<ps<<endl;  
  59.         }  
  60.         printf("%s %d\n",name[ps],b[cnt]);  
  61.     }  
  62.     return 0;  
  63. }  

14、

POJ 2155 Matrix 【二维线段树模板题】

[cpp] view plain copy print?
  1. #include <stdio.h>  
  2. #include <string.h>  
  3. #include <iostream>  
  4.   
  5. using namespace std;  
  6. const int MAXN = 1010;  
  7. struct Nodey  
  8. {  
  9.     int l,r;  
  10.     int val;  
  11. };  
  12. int n;  
  13. int locx[MAXN],locy[MAXN];  
  14. struct Nodex  
  15. {  
  16.     int l,r;  
  17.     Nodey sty[MAXN*3];  
  18.     void build(int i,int _l,int _r)  
  19.     {  
  20.         sty[i].l = _l;  
  21.         sty[i].r = _r;  
  22.         sty[i].val = 0;  
  23.         if(_l == _r)  
  24.         {  
  25.             locy[_l] = i;  
  26.             return;  
  27.         }  
  28.         int mid = (_l + _r)>>1;  
  29.         build(i<<1,_l,mid);  
  30.         build((i<<1)|1,mid+1,_r);  
  31.     }  
  32.     void add(int i,int _l,int _r,int val)  
  33.     {  
  34.         if(sty[i].l == _l && sty[i].r == _r)  
  35.         {  
  36.             sty[i].val += val;  
  37.             return;  
  38.         }  
  39.         int mid = (sty[i].l + sty[i].r)>>1;  
  40.         if(_r <= mid)add(i<<1,_l,_r,val);  
  41.         else if(_l > mid)add((i<<1)|1,_l,_r,val);  
  42.         else  
  43.         {  
  44.             add(i<<1,_l,mid,val);  
  45.             add((i<<1)|1,mid+1,_r,val);  
  46.         }  
  47.     }  
  48. }stx[MAXN*3];  
  49. void build(int i,int l,int r)  
  50. {  
  51.     stx[i].l = l;  
  52.     stx[i].r = r;  
  53.     stx[i].build(1,1,n);  
  54.     if(l == r)  
  55.     {  
  56.         locx[l] = i;  
  57.         return;  
  58.     }  
  59.     int mid = (l+r)>>1;  
  60.     build(i<<1,l,mid);  
  61.     build((i<<1)|1,mid+1,r);  
  62. }  
  63. void add(int i,int x1,int x2,int y1,int y2,int val)  
  64. {  
  65.     if(stx[i].l == x1 && stx[i].r == x2)  
  66.     {  
  67.         stx[i].add(1,y1,y2,val);  
  68.         return;  
  69.     }  
  70.     int mid = (stx[i].l + stx[i].r)/2;  
  71.     if(x2 <= mid)add(i<<1,x1,x2,y1,y2,val);  
  72.     else if(x1 > mid)add((i<<1)|1,x1,x2,y1,y2,val);  
  73.     else  
  74.     {  
  75.         add(i<<1,x1,mid,y1,y2,val);  
  76.         add((i<<1)|1,mid+1,x2,y1,y2,val);  
  77.     }  
  78. }  
  79. int sum(int x,int y)  
  80. {  
  81.     int ret = 0;  
  82.     for(int i = locx[x];i;i >>= 1)  
  83.         for(int j = locy[y];j;j >>= 1)  
  84.             ret += stx[i].sty[j].val;  
  85.     return ret;  
  86. }  
  87.   
  88. int main()  
  89. {  
  90.     //freopen("in.txt","r",stdin);  
  91.     //freopen("out.txt","w",stdout);  
  92.     int T;  
  93.     scanf("%d",&T);  
  94.     while(T--)  
  95.     {  
  96.         int q;  
  97.         scanf("%d%d",&n,&q);  
  98.         build(1,1,n);  
  99.         char op[10];  
  100.         int x1,x2,y1,y2;  
  101.         while(q--)  
  102.         {  
  103.             scanf("%s",op);  
  104.             if(op[0] == 'C')  
  105.             {  
  106.                 scanf("%d%d%d%d",&x1,&y1,&x2,&y2);  
  107.                 add(1,x1,x2,y1,y2,1);  
  108.             }  
  109.             else  
  110.             {  
  111.                 scanf("%d%d",&x1,&y1);  
  112.                 if(sum(x1,y1)%2 == 0)printf("0\n");  
  113.                 else printf("1\n");  
  114.             }  
  115.         }  
  116.         if(T)printf("\n");  
  117.     }  
  118.     return 0;  
  119. }  
15、

hdu4325 Flowers【树状数组区间更新单点求值 离散化】

[cpp] view plain copy print?
  1. #include <cstdio>  
  2. #include <algorithm>  
  3. #include<cstring>  
  4. using namespace std;  
  5. #define maxn 200008  
  6. int n,tree[maxn],X[maxn],li[maxn],ri[maxn],nn,m,mm,ma;  
  7. int lowbit(int i)  
  8. {  
  9.     return i&(-i);  
  10. }  
  11. void update(int i,int x)  
  12. {  
  13.     while(i<=ma)  
  14.     {  
  15.         tree[i]=tree[i]+x;  
  16.         i=i+lowbit(i);  
  17.     }  
  18. }  
  19. long long query(int n)  
  20. {  
  21.     long long sum=0;  
  22.     while(n>0)  
  23.     {  
  24.         sum+=tree[n];  
  25.         n=n-lowbit(n);  
  26.     }  
  27.     return sum;  
  28. }  
  29. int Bin(int num,int R)  
  30. {  
  31.     int l=0,r=R-1,mid;  
  32.     while(l<=r)  
  33.     {  
  34.         mid=(l+r)/2;  
  35.         if(X[mid]==num) return mid;  
  36.         if(X[mid]<num) l=mid+1;  
  37.         else r=mid-1;  
  38.     }  
  39. }  
  40. int q[maxn];  
  41. int main()  
  42. {  
  43.   //  freopen("cin.txt","r",stdin);  
  44.     int t,cas=1;  
  45.     scanf("%d",&t);  
  46.     while(t--)  
  47.     {  
  48.         nn=0,m=1,mm;  
  49.         memset(X,0,sizeof(X));  
  50.         scanf("%d%d",&n,&mm);  
  51.         for(int i=0;i<n;i++) scanf("%d%d",&li[i],&ri[i]),X[nn++]=li[i],X[nn++]=ri[i];  
  52.         for(int i=0;i<mm;i++)scanf("%d",&q[i]),X[nn++]=q[i];  
  53.         sort(X,X+nn);  
  54.        // for(int i=0;i<nn;i++)printf("x=%d  ",X[i]);puts("");  
  55.         for(int i=1;i<nn;i++)  if(X[i]!=X[i-1]) X[m++]=X[i];  
  56.        // for(int i=m-1;i>0;i--)if(X[i]-X[i-1]>=1)X[m++]=X[i-1]+1;  
  57.         sort(X,X+m);  
  58.         ma=m+1;  
  59.       //  for(int i=0;i<m;i++)printf("x=%d  ",X[i]);puts("");  
  60.       //  printf("m=%d\n",m);  
  61.         memset(tree,0,sizeof(tree));  
  62.         for(int i=0;i<n;i++)  
  63.         {  
  64.             int l=Bin(li[i],m);  
  65.             int r=Bin(ri[i],m);  
  66.          //   printf("l=%d r=%d\n",l+1,r+1);  
  67.             update(l+1,1);  
  68.             update(r+2,-1);  
  69.           //  update(li[i],1);update(ri[i]+1,-1);  
  70.         }  
  71.         printf("Case #%d:\n",cas++);  
  72.         for(int i=0;i<mm;i++)  
  73.         {  
  74.             int qu=Bin(q[i],m);  
  75.            // printf("qu=%d\n",qu+1);  
  76.             printf("%I64d\n",query(qu+1));  
  77.         }  
  78.   
  79.     }  
  80.     return 0;  
  81. }  


16 扫描线

hdu 1828 Picture【扫描线求周长模板题】

POJ1151 (HDU 1542) Atlantis【矩形面积并，线段树+离散化+扫描线模板】 


阅读全文
0 0
阿里云空间多少钱一年 阿里云空间多少钱一年
原创粉丝点击
热门IT博客
c语言 bcdc语言 bcd
淘宝网店铺购买淘宝网店铺购买
恐怖游戏知乎
linux下搭建lamp环境
java web高级编程豆瓣
数据交换平台产品
java共享内存
photoshop mac版 2015
傲剑降龙升级数据
qq邮箱有mac版本
ui和程序员漫画
apm飞控用的什么单片机
无损卡尔曼滤波知乎
济南儿童编程教育
mac如何打开隐藏文件夹
sql 随机选取
淘宝hot
我的世界怎么开枪械js
德意志斯坦 知乎
机房装修 知乎
热门问题 老师的惩罚 人脸识别 我在镇武司摸鱼那些年 重生之率土为王 我在大康的咸鱼生活 盘龙之生命进化 天生仙种 凡人之先天五行 春回大明朝 姑娘不必设防，我是瞎子 虎牌保险柜初始密码 密码牌 密码牌怎么看 密码牌图解 密码箱怎么改密码 密码箱改密码 密码箱怎么恢复出厂设置 密码箱快速打开图解 密码箱怎么设置密码 密码箱密码 密码箱打不开 密码箱初始密码 旅行箱密码锁 旅行箱密码锁设置图解 密码箱初始密码怎么改 旅行箱密码 密码箱密码正确打不开 怎样开密码箱的锁 密码箱开锁公司 密码箱卡住了打不开怎么办 铝合金密码箱 男士密码箱 小密码箱 密码箱密码忘了 密码箱的密码怎么改 密码箱的图片 密码箱如何设置密码 密码箱怎么锁 密码箱图片 密码箱密码怎么设置 密码箱改密码怎么改 密码箱排行榜 开密码箱的方法 密码箱怎么开 密码箱锁了 密码箱怎么用 密码箱多少钱 小型密码箱 密码箱什么牌子的好 怎么开密码箱的锁 怎么开密码箱

程序博客网，程序员的互联网技术博客家园。csdn论坛精品 msdn技术资料都在这里