1. 程式人生 > >[NOIP模擬][好題]Problem

[NOIP模擬][好題]Problem

題目描述:
今天豆豆在做作業的時候遇到這麼一個問題:
給出 N 個正整數 a1..aN,再給出一個正整數 k ,現在可以進行如下操作:每次選擇一個大於 k 的正整數 ai,將 ai 減去 1 ,選擇 a[i-1] 或 a[i+1] 中的一個加上 1 。經過一定次數的操作後,最大能夠選出多長的一個連續子序列,使得這個子序列的每個數都不小於 k 。
總共有 M 次詢問,每次詢問給出的一個 k ,回答這個詢問。
輸入格式:
第一行兩個整數 N , m 代表陣列長度和詢問個數。
接下來一行,第 i 個正整數表示 ai 。
第三行 M 個正整數,第 i 個正整數表示第 i 次詢問的 k 。
輸出格式:


對於每個詢問,輸出一個整數表示問題的答案。
樣例輸入:
5 6
1 2 1 1 5
1 2 3 4 5 6
樣例輸出:
5 5 2 1 1 0
備註:
對於 20% 的資料,N≤10。
對於 40% 的資料,N≤1000。
對於 50% 的資料,N≤100000。
對於 100% 的資料,N≤1000000,M≤20,ai109k109
請選手注意自己的常數。
題目分析:
我們發現所謂的ai減一,ai1ai+1加一,其實就是數字轉移的過程。那麼對於一段序列只要它的總和大於等於個數*K(平均數大於K)就行。於是你維護字首和,n2列舉每一段的和來判斷,就可以得到40分。如果我們給原序列每一個都減去K,再維護字首和,那麼對於sum[i],在前面找到一個比它小的sum[j],所以j+1~i就是可以成立的(他們倆對減一下即大於0,且是j+1~i的和)。於是你線段樹維護一下,就可以拿50分了。
正解:對於上面的序列,我們可以把第一個數加入單調遞減棧,掃一遍依次把滿足要求的數加入棧。然後我們再從頭開始掃,對於一個sum[i],一個指標指向棧中的某一個數(初始時你可以指向0,或者棧的最後一個),如果sum[i]大於指標所指向的左邊數,你就向左邊移,直到不能移為止,此時指標指向的那個數小於你的sum[i],那麼這兩個數之間的區間是滿足要求的,計算一下,與當前算出最大的長度比較。如果小於指標所指向的那個數,指標右移,但是一旦出現右移就是不優秀的,因為你的sum[i]只比sum[i-1]右移了一個,而你的指標右移,在原序列也至少移了一個,所以不會更優秀,無須再右移。因為你最多隻會右移n次,而你的左移操作,其實是在抵消右移,也不會多於n次,所以整體是O(n)的。
附程式碼:

#include<iostream>
#include<cstring>
#include<string>
#include<cstdlib>
#include<cstdio>
#include<ctime>
#include<iomanip>
#include<cmath>
#include<cctype>
#include<algorithm>
using namespace std;

const int N=1e6+10;
int n,m,num,po,ans;
long
long sum[N],k,a[N]; struct node{ long long w; int po; }stack[N]; int readint() { char ch;int i=0,f=1; for(ch=getchar();(ch<'0'||ch>'9')&&ch!='-';ch=getchar()); if(ch=='-') {ch=getchar();f=-1;} for(;ch>='0'&&ch<='9';ch=getchar()) i=(i<<3)+(i<<1)+ch-'0'; return i*f; } int main() { //freopen("problem.in","r",stdin); //freopen("problem.out","w",stdout); n=readint();m=readint(); for(int i=1;i<=n;i++) a[i]=readint(); while(m--) { num=0; ans=0; k=readint(); for(int i=1;i<=n;i++) { sum[i]=sum[i-1]+a[i]-k; if(sum[i]<stack[num].w) { stack[++num].w=sum[i]; stack[num].po=i; } } po=0; for(int i=1;i<=n;i++) { if(sum[i]<stack[po].w) {po=min(po+1,num);continue;} while(sum[i]>=stack[po-1].w&&po>0) po--; ans=max(ans,i-stack[po].po); } printf("%d ",ans); } return 0; }