rmq模板

[cpp]  view plain  copy

1. #include <stdio.h>  
2. #include <string.h>  
3. #include <math.h>  
4. #define M 50010  
5. int mx[M][17],mi[M][17],d[M];  
6. int n,t;  
7. int min(int a,int b)  
8. {  
9.     return a<b?a:b;  
10. }  
11. int max(int a,int b)  
12. {  
13.     return a>b?a:b;  
14. }  
15. void init()  
16. {  
17.     int i,j,m;  
18.     for(i=1;i<=n;i++)  
19.     mi[i][0]=mx[i][0]=d[i];  
20.     m=(int)floor(log(1.0*n)/log(2.0));  
21.     for(i=1;i<=m;i++)  
22.     {  
23.         for(j=n;j>=1;j--)  
24.         {  
25.             mx[j][i]=mx[j][i-1];  
26.             if(j+(1<<(i-1))<=n)  
27.             mx[j][i]=max(mx[j][i],mx[j+(1<<(i-1))][i-1]);  
28.             mi[j][i]=mi[j][i-1];  
29.             if(j+(1<<(i-1))<=n)  
30.             mi[j][i]=min(mi[j][i],mi[j+(1<<(i-1))][i-1]);  
31.         }  
32.     }  
33. }  
34. int rmq_max(int l,int r)  
35. {  
36.     int m=(int)floor(log(1.0*(r-l+1))/log(2.0));  
37.     return max(mx[l][m],mx[r-(1<<m)+1][m]);  
38. }  
39. int rmq_min(int l,int r)  
40. {  
41.     int m=(int)floor(log(1.0*(r-l+1))/log(2.0));  
42.     return min(mi[l][m],mi[r-(1<<m)+1][m]);  
43. }  
44. int main()  
45. {  
46.     int i,l,r;  
47.     while(scanf("%d%d",&n,&t)!=EOF)  
48.     {  
49.         for(i=1;i<=n;i++)  
50.         scanf("%d",&d[i]);  
51.         init();  
52.         while(t--)  
53.         {  
54.             scanf("%d%d",&l,&r);  
55.             int ans=rmq_max(l,r)-rmq_min(l,r);  
56.             printf("%d\n",ans);  
57.         }  
58.     }  
59.     return 0;  
60. }  

源博客地址：http://blog.csdn.net/yuan_lo/article/details/8744343