天天看點

[bzoj 4241]曆史研究

IOI國曆史研究的第一人——JOI教授,最近獲得了一份被認為是古代IOI國的住民寫下的日記。JOI教授為了通過這份日記來研究古代IOI國的生活,開始着手調查日記中記載的事件。

日記中記錄了連續N天發生的時間,大約每天發生一件。

事件有種類之分。第i天(1<=i<=N)發生的事件的種類用一個整數Xi表示,Xi越大,事件的規模就越大。

JOI教授決定用如下的方法分析這些日記:

1. 選擇日記中連續的一些天作為分析的時間段

2. 事件種類t的重要度為t*(這段時間内重要度為t的事件數)

3. 計算出所有事件種類的重要度,輸出其中的最大值 現在你被要求制作一個幫助教授分析的程式,每次給出分析的區間,你需要輸出重要度的最大值。

這道題用莫隊添加很容易,但是删除很難,是以我們就用一個叫復原莫隊的東西。

我們定義l為目前詢問的左端點所在的塊的下一個塊的左端點,r為目前詢問的左端點所在的塊的右端點(一開始),由于排序之後,如果詢問的左端點所在的塊相同,那麼它們的詢問的右端點就會遞增。

那麼我們r就可以單調遞增,記錄答案為s2,l就向左改變到詢問的左端點,記錄最終答案為s1,然後l就向右移動回原來的位置,消除l向左的影響,s1=s2。那這道題就解決了。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdlib>
using namespace std;
inline long long Max(long long a,long long b){return a>b?a:b;}
struct node
{
    int l,r,id;
    long long s;
    node(){l=s=;}
}q[100010],b[];
int n,m,cnt; 
int a[],belong[],ex[];
long long s1,s2,sum[],num[];
void fk()
{
    cnt=sqrt(n);
    for(int i=;i<=n;i++)belong[i]=(i-)/cnt+;
}
bool cmp1(node a,node b)
{
    if(belong[a.l]==belong[b.l])
    {
        if(a.r>b.r)return false;
        if(a.r<b.r)return true;
        return ;
    }
    if(belong[a.l]>belong[b.l])return false;
    if(belong[a.l]<belong[b.l])return true;
}
bool cmp2(node a,node b)
{
    if(a.id>b.id)return false;
    if(a.id<b.id)return true;
    return ;
}
bool cmp3(node a,node b)
{
    if(a.s>b.s)return false;
    if(a.s<b.s)return true;
    return ;
}
void add(int x)
{
    sum[x]+=ex[x];
    s1=Max(s1,sum[x]);
}
int main()
{
    int ss=;
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)scanf("%d",&b[i].s),b[i].id=i;
    sort(b+,b+n+,cmp3);
    for(int i=;i<=n;i++)
    {
        if(b[i].s!=b[i-].s)ss++;
        a[b[i].id]=ss;
        ex[ss]=b[i].s;
    }
    for(int i=;i<=m;i++)
    {
        scanf("%d%d",&q[i].l,&q[i].r);q[i].id=i;
    }
    fk();
    sort(q+,q+m+,cmp1);
    int l,r;
    for(int i=;i<=m;i++)
    {
        if(belong[q[i].l]!=belong[q[i-1].l])
        {
            memset(sum,,sizeof(sum));
            s1=s2=;
            l=belong[q[i].l]*cnt+;r=l-;
        }
        if(belong[q[i].l]==belong[q[i].r])
        {
            long long ans=;      
            for(int j=q[i].l;j<=q[i].r;j++)
            {
                num[a[j]]+=ex[a[j]];ans=Max(ans,num[a[j]]);
            }
            q[i].s=ans;
            for(int j=q[i].l;j<=q[i].r;j++)num[a[j]]-=ex[a[j]];
        }
        else
        {
            for(int j=r+;j<=q[i].r;j++)add(a[j]);
            s2=s1;
            for(int j=l-;j>=q[i].l;j--)add(a[j]);
            q[i].s=s1;
            for(int j=q[i].l;j<=l-;j++)sum[a[j]]-=ex[a[j]];
            s1=s2;
            r=q[i].r;
        }
    }
    sort(q+,q+m+,cmp2);
    for(int i=;i<=m;i++)printf("%lld\n",q[i].s);
    return ;
}
           

繼續閱讀