BZOJ3636: 教义问答手册

Description

“汉中沃野如关中，四五百里烟蒙蒙。黄云连天夏麦熟，水稻漠漠吹秋风。”——摘自 黄裳《汉中行》
“泉岭精神不朽，汉中诸球永生。”——摘自《泉岭精神创立者语录》
“把神犇烤一烤，味道会更好。”——摘自《xhr语录》
“秀恩爱有利于身心健康！”——摘自《泉岭精神集大成者语录》
“楼上说的对！”——摘自《泉岭精神信徒语录合集》
“不会做积分，怎么找妹子！”——摘自《xhr语录》
“切实保护耕地以放置更多的哨戒炮。”——摘自《泉岭精神信徒语录合集》
“就算两个包子一起吃掉，也不能阻止我修筑梯田。”——摘自《泉岭精神创立者语录》
“我来自泉岭，他来自汉中，我们半道而逢。”——摘自《泉岭精神集大成者语录》
【问题描述】
    作为泉岭精神的缔造者、信奉者、捍卫者、传承者，Pear决定印制一些教义问答手册，以满足泉岭精神日益增多的信徒。Pear收集了一些有关的诗选、语录，其中部分内容摘录在了【题目背景】里。这些语录是按出现的时间排好序的——Pear很喜欢这样的作风，于是决定在按时间排好序的基础上，选择部分语录，制作成若干本教义问答手册。
    一共有N条语录。Pear决定从中选出某一段时间内的所有语录，在此基础上印制大小为L的若干本教义问答手册。Pear对印制的手册有如下要求：
    1.每本手册必须包含这个区间内连续的恰好L条语录。
    2.不同手册包含的语录不能相同。
    3.每条语录有一个“主题相关程度”，这个数可正可负。Pear希望所有手册的语录的“主题相关程度”之和尽可能大。
    例如，对于区间[3,15]和L=3，一种选择方法是：[4,6]+[9,11]+[12,14]。这三个区间长度都恰好为L，且互不重叠。
    Pear并没有决定选哪段时间的语录，因此他有Q次询问。每次询问，给出两个数[l,r]表示候选语录的范围是第l条到第r条。你能回答出每个询问的最大“主题相关程度”之和么？

Input

第一行两个正整数N，L，含义如上所述。注意对于所有询问，L都是一样的。
    第二行N个整数，绝对值<=10000。第i个数表示第i条语录的“主题相关程度”。
    接下来Q行，每行两个正整数l和r，表示询问区间。

Output

输出Q行，每行表示这组询问的答案。注意，这个答案可以是0，如果区间负数过于多的话。

Sample Input

15 3
3 1 5 -2 3 -2 -2 2 2 2 0 3 2 -1 0
9
8 10
10 10
9 11
2 14
5 14
5 13
12 13
7 13
2 10

Sample Output

6
0
4
17
11
11
0
11
12

HINT

【数据范围】

对于10%的数据，N=1000,Q=1000,L<=50

对于另外20%的数据，N=100000,Q=100000,L<=5

对于另外20%的数据，N=100000,Q=100000,L<=10

对于100%的数据，N=100000,Q=100000,L<=50

Source

2014年国家集训队十五人互测

考虑线段树分治，对于每个询问向覆盖它的节点打上标记。

然后一起处理时枚举跨中点的一块的左端点，然后左右DP，对覆盖区间的询问更新答案。

数组好像要开到4Q*logN。

时间复杂度为O(NLlogN)。

#include<cstdio>
#include<cctype>
#include<queue>
#include<cmath>
#include<cstring>
#include<algorithm>
#define rep(i,s,t) for(int i=s;i<=t;i++)
#define dwn(i,s,t) for(int i=s;i>=t;i--)
#define ren for(int i=first[x];i;i=next[i])
using namespace std;
const int BufferSize=<<;
char buffer[BufferSize],*head,*tail;
inline char Getchar() {
    if(head==tail) {
        int l=fread(buffer,,BufferSize,stdin);
        tail=(head=buffer)+l;
    }
    return *head++;
}
inline int read() {
    int x=,f=;char c=Getchar();
    for(;!isdigit(c);c=Getchar()) if(c=='-') f=-;
    for(;isdigit(c);c=Getchar()) x=x*+c-'';
    return x*f;
}
const int maxn=;
const int maxnode=;
int n,L,m,A[maxn],ans[maxn];
int first[maxn<<],next[maxnode],num[maxnode],ql[maxnode],qr[maxnode],ToT;
void AddQuery(int pos,int l,int r,int id) {
    num[++ToT]=id;ql[ToT]=l;qr[ToT]=r;next[ToT]=first[pos];first[pos]=ToT;
}
int f[maxn],g[maxn],S[maxn];
void solve(int x,int l,int r) {
    if(r-l+<L) return;
    int mid=l+r>>,lc=x<<,rc=lc|;
    ren {
        if(qr[i]<=mid) AddQuery(lc,ql[i],qr[i],num[i]);
        else if(ql[i]>mid) AddQuery(rc,ql[i],qr[i],num[i]);
        else if(ql[i]!=l||qr[i]!=r) AddQuery(lc,ql[i],mid,num[i]),AddQuery(rc,mid+,qr[i],num[i]);
    }
    rep(p,max(l,mid-L+),mid) {
        f[p-]=g[p]=;
        rep(i,p,r) {
            if(i-p+<L) f[i]=;
            else f[i]=max(f[i-],f[i-L]+(S[i]-S[i-L]));
        }
        dwn(i,p-,l) {
            if(p-i<L) g[i]=;
            else g[i]=max(g[i+],g[i+L]+(S[L+i-]-S[i-]));
        }
        ren if(qr[i]>=p&&ql[i]<=p) ans[num[i]]=max(ans[num[i]],f[qr[i]]+(ql[i]<p?g[ql[i]]:));
    }
    if(l<r) solve(lc,l,mid),solve(rc,mid+,r);
}
int main() {
    n=read();L=read();
    rep(i,,n) S[i]=S[i-]+(A[i]=read());
    m=read();
    rep(i,,m) {
        int l=read(),r=read();
        if(l<=r) AddQuery(,l,r,i);
    }
    solve(,,n);
    rep(i,,m) printf("%d\n",ans[i]);
    return ;
}