51nod 1672

原题链接

找出和最大的区间交。

所有区间按左端点排序，建立线段树，枚举交区间的左端点，并不断加入右端点，每次加入右端点时，对应区间的值皆加一，找出值不小于K的最右端点，区间和可以通过预处理前缀和得出。

为保证左端点被K个区间覆盖的存在，需要预先加入K个区间。

#include<map>
#include<set>
#include<list>
#include<cmath>
#include<ctime>
#include<queue>
#include<stack>
#include<cstdio>
#include<string>
#include<bitset>
#include<vector>
#include<utility>
#include<numeric>
#include<cstring>
#include<iterator>
#include<iostream>
#include<algorithm>
#include<functional>
using namespace std;
#define fi first
#define se second
#define MP make_pair
#define ll long long
#define PR pair<ll,ll>
#define LS (root<<1)
#define RS ((root<<1)|1)
#define LSON LS,l,mid
#define RSON RS,(mid+1),r
#define MID mid=((l+r)/2)
#define mm(a,b) memset(a,b,sizeof(a))
#define rep(a,b,c) for(int a=b;a<c;a++)
#define INF 0x3f3f3f3f3f3f3f3f
#define eps 1e-9
#define Pi acos(-1.0)
#define MAXN 100005
#define MAXE 3000050
#define Inf 9999999
#define debug puts("**************");
template<class T> T f_max(T a, T b){ return a > b ? a : b; }
template<class T> T f_min(T a, T b){ return a < b ? a : b; }
template<class T> T f_abs(T a){ return a > 0 ? a : -a; }
template<class T> T gcd(T a, T b){ return b ? gcd(b, a%b) : a; }
template<class T> T lcm(T a, T b){ return a / gcd(a, b)*b; }
const int mod=1000007;
struct seg{
    int l, r;
    bool operator <(const seg &a)const{
        return l<a.l;
    }
}S[MAXN];
int a[MAXN];
ll sum[MAXN];
int cover[MAXN<<2];
void build(int root,int l,int r){
    cover[root]=0;
    if(l==r)
        return;
    int mid=(l+r)>>1;
    build(root<<1,l,mid);
    build(root<<1|1,mid+1,r);
}
void push_up(int root){
    cover[root]=cover[root<<1]+cover[root<<1|1];
}
void update(int root,int l,int r,int x){
    if(l==r){
        cover[root]++;
        return;
    }
    int mid=(l+r)>>1;
    if(x>mid)
        update(root<<1|1,mid+1,r,x);
    else
        update(root<<1,l,mid,x);
    push_up(root);
}
int query(int root,int l,int r,int k){ //find the Kth rightmost end
    if(l==r)
        return l;
    int mid=(l+r)>>1;
    if(cover[root<<1|1]>=k)
        return query(root<<1|1,mid+1,r,k);
    else
        return query(root<<1,l,mid,k-cover[root<<1|1]);
}
int main(){
    int n, k, m;

    while(~scanf("%d%d%d",&n, &k, &m)){
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            sum[i]=sum[i-1]+a[i];
        }
        for(int i=1;i<=m;i++){
            scanf("%d%d",&S[i].l,&S[i].r);
        }
        sort(S+1,S+1+m);
//        for(int i=1;i<=m;i++)
//            printf("l=%d r=%d\n",S[i].l,S[i].r);
        build(1,1,n);
        for(int i=1;i<k;i++)
            update(1,1,n,S[i].r);
        ll ans=0;
        for(int i=k;i<=m;i++){
            update(1,1,n,S[i].r);
            int pos=query(1,1,n,k);
//            printf("pos=%d\n",pos);
            if(pos>=S[i].l)
                ans=max(ans,sum[pos]-sum[S[i].l-1]);
        }
        printf("%lld\n",ans);
    }

    return 0;
}