「CCO 2020」购物计划

「CCO 2020」购物计划

核心思想：堆+调整临近

\(x_i=y_i=1\)

这个限制相当于每一组内的权值排名可以确定，设组内为\(A_{i,j}(j\ge 1)\)

那么我们一个方案的选择可以用\(M\)个指针\(P_i\)表示，和为\(\displaystyle \sum A_{i,P_i}\)

考虑用调整的方式解决这个问题，大体思路上，我们可以记录当前移动指针\(P_p\)

每次可以选择移动\(P_p\)，或者某一个\(P_{i},i>p\)

如果直接进行，每次移动的指针数量是\(O(m)\)级别的，显然不可行

考虑优化一下进行，每次只能选取\(i=p+1\)

此时，编号小的会先被移动

为了保证答案单调性，我们需要将\(A_{i,2}-A_{i,1}\)较小的组先移动

同时，并不是每一个组都会被移动，因此转移还要支持一个特殊回撤操作，来撤回当前组的指针

也就是说，若\(P_p=2\)，可以选择把\(P_p\)回撤为\(1\)，然后将\(P_{p+1}\)改为\(2\)

由此，每个点状态可以选择：

1.移动自己

2.移动下一个

3.若\(P_p=2\)，回撤自己，同时移动下一个

这样的调整法，可以保证每一个状态恰好有一个前驱，且转移过程中值不断变大

由此可以\(O(k)\)状态数进行调整，用堆维护，复杂度为\(O(k\log k)\)

\[\ \]

组内调整

一个组内会选择若干个数\(A_{b_i},i\in [1,c]\)

初始最小值，显然满足\(b_i=i\)

类似的，我们记录当前指针\(p\)，前驱指针\(l\)，后继指针\(r\)

显然\(p\)要往后移，且不能达到\(r\)，因此决策只有两种

1.移动前驱\(l\)，并将当前指针变为前驱

2.移动自己\(p\)

这样的调整是固定个数的，因此，一开始就把\(c\in[x_i,y_i]\)的所有情况插入即可

\[\ \]

最后，将两部分一同进行，每次组间调整时，通过组内调整查询答案

总的组内和组间调整次数均为\(O(k)\)，状态数分别不超过\(2k,3k\)

复杂度为\(O(k\log k)\)

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define pb push_back
#define rep(i,a,b) for(int i=a,i##end=b;i<=i##end;++i)
char IO;
int rd(){
    int s=0;
    while(!isdigit(IO=getchar()));
    do s=(s<<1)+(s<<3)+(IO^'0');
    while(isdigit(IO=getchar()));
    return s;
}
enum{N=200010,LINF=1ll<<60};

int n,m,k,I[N]; ll L;
struct Group{
    vector <int> A;
    int l,r,c;
    struct Node{
        int l,p,r; ll s;
        bool operator < (const Node &__) const { return s>__.s; }
    };
    priority_queue <Node> que;
    vector <ll> V;
    void Init(){
        l=rd(),r=rd(),sort(A.begin(),A.end()),c=A.size();
        if(c<l) {
            rep(i,1,k) puts("-1");
            exit(0);
        }
        r=min(r,c);
        rep(i,0,l-1) L+=A[i];
        ll s=0;
        if(!l) V.pb(0);
        rep(i,0,c-1) {
            if(i>=r) break;
            s+=A[i];
            if(i>=l-1) que.push((Node){i-1,i,c-1,s});
        }
    }
    void Next(){
        if(que.empty()) return V.pb(LINF);
        Node t=que.top(); que.pop();
        V.pb(t.s);
        if(t.p<t.r) que.push((Node){t.l,t.p+1,t.r,t.s-A[t.p]+A[t.p+1]}); // Move current point
        if(~t.l && t.l<t.p-1) que.push((Node){t.l-1,t.l+1,t.p-1,t.s-A[t.l]+A[t.l+1]}); // Move previous point 
    }
    // get kth sum
    ll operator [] (const int &k){
        while((int)V.size()<k) Next();
        return V[k-1];
    }
} S[N];

struct Node{
    int x,y; ll s;
    bool operator < (const Node &__) const { return s>__.s; }
};
priority_queue <Node> que;

int main(){
    n=rd(),m=rd(),k=rd();
    rep(i,1,n) { int x=rd(); S[x].A.pb(rd()); }
    rep(i,1,m) S[i].Init();
    printf("%lld\n",L),k--;
    rep(i,1,m) I[i]=i;
    sort(I+1,I+m+1,[&](int x,int y){ return S[x][2]-S[x][1]<S[y][2]-S[y][1]; });
    que.push((Node){1,2,L-S[I[1]][1]+S[I[1]][2]});
    while(k) {
        Node t=que.top(); que.pop();
        if(t.s>=LINF) break;
        k--,printf("%lld\n",t.s);
        int i=I[t.x],j=I[t.x+1];
        que.push((Node){t.x,t.y+1,t.s-S[i][t.y]+S[i][t.y+1]});// Move current point
        if(j) que.push((Node){t.x+1,2,t.s-S[j][1]+S[j][2]}); // Move next point
        if(t.y==2 && j) que.push((Node){t.x+1,2,t.s-S[i][2]+S[i][1]-S[j][1]+S[j][2]}); 
        // Back current point ,and move next point
    }
    while(k--) puts("-1");
}