poj 2823 Sliding Window (单调队列)
Sliding Window
看这个问题：An array of size n ≤ 106 is given to you. There is a sliding window of size k which is moving from the very left of the array to the very right. You can only see the k numbers in the window. Each time the sliding window moves rightwards by one position.Your task is to determine the maximum and minimum values in the sliding window at each position.
也就是有一个数列a，要求你求数列b和c，b[i]是a[i]…a[i+w-1]中的最小值，c[i]是最大值。如果a是1,3,-1,-3,5,3,6,7，则b为-1,-3,-3,-3,3,3，c为3,3,5,5,6,7。
这个问题相当于一个数据流（数列a）在不断地到来，而数据是不断过期的，相当于我们只能保存有限的数据（sliding window中的数据，此题中就是窗口的宽度w），对于到来的查询（此题中查询是每时刻都有的），我们要返回当前滑动窗口中的最大值最小值。注意，元素是不断过期的。
解决这个问题可以使用一种叫做单调队列的数据结构，它维护这样一种队列：
a)从队头到队尾，元素在我们所关注的指标下是递减的（严格递减，而不是非递增），比如查询如果每次问的是窗口内的最小值，那么队列中元素从左至右就应该递增，如果每次问的是窗口内的最大值，则应该递减，依此类推。这是为了保证每次查询只需要取队头元素。
b)从队头到队尾，元素对应的时刻（此题中是该元素在数列a中的下标）是递增的，但不要求连续，这是为了保证最左面的元素总是最先过期，且每当有新元素来临的时候一定是插入队尾。
满足以上两点的队列就是单调队列，首先，只有第一个元素的序列一定是单调队列。
那么怎么维护这个单调队列呢？无非是处理插入和查询两个操作。
对于插入，由于性质b，因此来的新元素插入到队列的最后就能维持b)继续成立。但是为了维护a)的成立，即元素在我们关注的指标下递减，从队尾插入新元素的时候可能要删除队尾的一些元素，具体说来就是，找到第一个大于（在所关注指标下）新元素的元素，删除其后所有元素，并将新元素插于其后。因为所有被删除的元素都比新元素要小，而且比新元素要旧，因此在以后的任何查询中都不可能成为答案，所以可以放心删除。
对于查询，由于性质b，因此所有该时刻过期的元素一定都集中在队头，因此利用查询的时机删除队头所有过期的元素，在不含过期元素后，队头得元素就是查询的答案（性质a），将其返回即可。
由于每个元素都进队出队一次，因此摊销复杂度为O(n)。
POJ2823就是上面描述的那道题。
http://www.cnblogs.com/Jason-Damon/archive/2012/04/19/2457889.html
讲得很好.我认为重点的地方加颜色了.
如果要找最大值,那么比新加入的元素小且旧的元素是不可能成为答案的,可以放心删去!
理解了这句话就理解了单调队列...
上面那个博客将得是挺好,只是代码写得略丑.
自己按照思想写了个还看得过去的
1
2 /*************************************************************************
3 > File Name: code/2015summer/单调队列/B.cpp
4 > Author: 111qqz
5 > Email: rkz2013@126.com
6 > Created Time: 2015年07月29日 星期三 20时18分16秒
7 ************************************************************************/
8
9 #include<iostream>
10 #include<iomanip>
11 #include<cstdio>
12 #include<algorithm>
13 #include<cmath>
14 #include<cstring>
15 #include<string>
16 #include<map>
17 #include<set>
18 #include<queue>
19 #include<vector>
20 #include<stack>
21 #define y0 abc111qqz
22 #define y1 hust111qqz
23 #define yn hez111qqz
24 #define j1 cute111qqz
25 #define tm crazy111qqz
26 #define lr dying111qqz
27 using namespace std;
28 #define REP(i, n) for (int i=0;i<int(n);++i)
29 typedef long long LL;
30 typedef unsigned long long ULL;
31 const int inf = 0x7fffffff;
32 const int N=1E6+7;
33 int a[N],mxans[N],q[N],indx[N],mians[N];
34 int n,k;
35 int head,tail;
36
37 int main()
38 {
39 cin>>n>>k;
40 for ( int i = 0 ; i < n ; i++ )
41 {
42 scanf("%d",&a[i]);
43 }
44 head = 1 ; tail = 0;
45 for ( int i = 0; i < n ; i++)
46 {
47 if (indx[head]+k-1<i)
48 {
49 head++;
50 }
51 while (head<=tail&&q[tail]<=a[i]) tail--;
52 tail++;
53 q[tail]=a[i];
54 indx[tail]=i;
55 if (i>=k-1)
56 {
57 mxans[i-k+1]=q[head];
58 }
59 }
60 head = 1 ; tail = 0;
61 memset(q,0,sizeof(q));
62 for ( int i = 0 ; i < n ; i++)
63 {
64 if (indx[head]+k-1<i)
65 {
66 head++;
67 }
68 while (head<=tail&&q[tail]>=a[i]) tail--;
69 tail++;
70 q[tail]=a[i];
71 indx[tail]=i;
72 if (i>=k-1)
73 {
74 mians[i-k+1]=q[head];
75 }
76 }
77 for ( int i = 0 ; i < n-k+1 ; i++ )
78 {
79 if (i==0)
80 cout<<mians[i];
81 else cout<<" "<<mians[i];
82 }
83 cout<<endl;
84 for ( int i = 0 ; i < n-k+1 ; i++ )
85 {
86 if (i==0)
87 cout<<mxans[i];
88 else cout<<" "<<mxans[i];
89 }
90
91 return 0;
92 }
93
94
更新一个双端队列实现的（为毛交g++会Tle，交c++才能过2333）：
1
2 /* ***********************************************
3 Author :111qqz
4 Created Time :2016年08月05日 星期五 00时09分33秒
5 File Name :code/poj/2823.cpp
6 ************************************************ */
7
8 #include <cstdio>
9 #include <cstring>
10 #include <iostream>
11 #include <algorithm>
12 #include <vector>
13 #include <queue>
14 #include <set>
15 #include <map>
16 #include <string>
17 #include <cmath>
18 #include <cstdlib>
19 #include <ctime>
20 #include <deque>
21 #define fst first
22 #define sec second
23 #define lson l,m,rt<<1
24 #define rson m+1,r,rt<<1|1
25 #define ms(a,x) memset(a,x,sizeof(a))
26 typedef long long LL;
27 #define pi pair < int ,int >
28 #define MP make_pair
29
30 using namespace std;
31 const double eps = 1E-8;
32 const int dx4[4]={1,0,0,-1};
33 const int dy4[4]={0,-1,1,0};
34 const int inf = 0x3f3f3f3f;
35 const int N=1E6+7;
36 int n,k;
37 int a[N];
38 int mn[N],mx[N];
39 deque<int>dq;
40 int main()
41 {
42 #ifndef ONLINE_JUDGE
43 freopen("code/in.txt","r",stdin);
44 #endif
45 scanf("%d%d",&n,&k);
46 for ( int i = 1 ; i <= n ; i++) scanf("%d",&a[i]);
47 for ( int i = 1 ; i <= n ; i++)
48 {
49 while (!dq.empty()&&dq.front()<i-k+1)dq.pop_front();
50 while (!dq.empty()&&a[dq.back()]>=a[i]) dq.pop_back();
51 dq.push_back(i);
52 if (i>=k)
53 mn[i-k+1] = a[dq.front()];
54 }
55
56 dq.clear();
57
58 for ( int i = 1 ; i <= n ; i++)
59 {
60 while (!dq.empty()&&dq.front()<i-k+1) dq.pop_front();
61 while (!dq.empty()&&a[dq.back()]<=a[i]) dq.pop_back();
62 dq.push_back(i);
63 if (i>=k)
64 mx[i-k+1] = a[dq.front()];
65 }
66
67 for ( int i = 1 ; i <n-k+1 ; i++)
68 printf("%d ",mn[i]);
69 printf("%d\n",mn[n-k+1]);
70 for ( int i = 1 ; i < n-k+1 ; i++)
71 printf("%d ",mx[i]);
72 printf("%d\n",mx[n-k+1]);
73
74 #ifndef ONLINE_JUDGE
75 fclose(stdin);
76 #endif
77 return 0;
78 }
79
80