poj 2823 Sliding Window（樹狀數組）

小記：求區間最值是可以用樹狀數組解的，但是這題我交c++最開始9600+ms過了，然後G++過不了，然後就優化，最後C++優化到了5600+ms過了，但是G++還是過不了

這個鬱悶啊，一看discuss都是這樣， 向來緣淺，奈何情深。 我不想放棄，遂看了一下別人的優化， 一看嚇一跳， 優化IO，791ms過了。。。。

然後我改了我的代碼的輸入，提交c++要7000+ms，提交G++2000+ms過了。

思路：參考這篇：點擊

代碼：

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <stdlib.h>

#define mst(a,b) memset(a,b,sizeof(a))
#define REP(a,b,c) for( a = b; a < c; ++a)
#define eps 10e-8

#define MAX_ 1000010


int a[MAX_];

int idx[MAX_], idy[MAX_];
int mmx[MAX_], mmi[MAX_];
int n, m;

inline int lowbit(int x){return x&(-x);}

inline int max(int a, int b){
    return a<b?b:a;
}
inline int min(int a, int b){
    return a<b?a:b;
}

void init(){
    int i, j;
    for( i = 1; i <= n; ++i){
        idx[i] = a[i];
        idy[i] = a[i];
        for( j = 1; j < lowbit(i); j<<=1){
            idx[i] = max(idx[i],idx[i-j]);
            idy[i] = min(idy[i],idy[i-j]);
        }
    }
}
inline int query1(int l,int r){
    int ans = a[r];

    while(1){
        ans = min(a[r], ans);
        if(l == r)break;
        for(r--; r-l >= lowbit(r); r -= lowbit(r)){
            ans = min(ans,idy[r]);
        }
    }
    return ans;
}


inline int query(int l,int r){
    int ans = a[r];

    while(1){
        ans = max(a[r],ans);
        if(l == r)break;
        for(r--; r-l >= lowbit(r); r -= lowbit(r)){
            ans = max(ans,idx[r]);
        }
    }
    return ans;
}
inline void put(int x){
    if(x< 0){
        putchar('-');
        x = -x;
    }
    if(x == 0){
        putchar('0');
        return;
    }
    char s[20];
    int bas = 0;
    for(;x;x/=10)s[bas++] = x%10+'0';
    for(;bas--;)putchar(s[bas]);
    return;
}

int main(){
	int T,  s, t, tmp, i, cnt;
	while(~scanf("%d%d", &n, &m)){

	    REP(i, 1, n+1){
            scanf("%d", &a[i]);
	    }
	    init();

        cnt = query1(1, m);
        put(cnt);
        REP(i, 2, n+1){
           if(i + m - 1<= n){
                    if(cnt == a[i-1]){
                        cnt = query1(i, i+m-1);
                    }else{
                        cnt = min(cnt, a[i+m-1]);
                    }
                    putchar(' ');
                    put(cnt);
            }
            else break;
        }
        cnt = query(1, m);
        putchar('\n');
        put(cnt);
        REP(i, 2, n+1){
            if(i + m - 1<= n){
                    if(cnt == a[i-1]){
                        cnt = query(i, i+m-1);
                    }else{
                        cnt = max(cnt, a[i+m-1]);
                    }
                    putchar(' ');
                    put(cnt);
                    //printf(" %d", cnt);
            }
            else break;
        }

        putchar('\n');


	}
	return 0;
}