其他分享
首页 > 其他分享> > 6358. 【NOIP2019模拟2019.9.15】小ω的仙人掌

6358. 【NOIP2019模拟2019.9.15】小ω的仙人掌

作者:互联网

题目

题目大意

给你一串二元组\((a_i,b_i)\)的数列。
求最小的区间\([l,r]\)长度,满足\([l,r]\)中的每个二元组选或不选,使得\(\sum a_i=w\)且\(\sum b_i\leq k\)


思考历程

想了好久,想来想去都是一个背包……
最终决定打暴力……


正解

先说说GMH大爷的神奇解法。
首先是二分答案\(ans\),转化成判定问题。然后在数列中每\(ans\)个点设置一个观测点。
以每个观测点为中心,向左和向右背包,然后合并。

然而正解并不需要一个\(\log\)
考虑双指针,就是记一个当前的最佳答案\(ans\),后面的区间长度都要小于\(ans\)。脑补一下这个过程,其实这就是一个队列,只需要支持左边出右边入的队列。
但是背包问题不满足可减性。于是就有个非常骚的解法:
把这个队列用两个栈来代替,栈顶分别为队头和队尾。
加入的时候,就在第二个栈的栈顶加入;弹出的时候,就直接弹出第一个栈的栈顶。
如果第一个栈为空,那就将第二个栈里的东西倒过来放到第一个栈中,然后暴力重构。
每个元素只会暴力重构一次,所以不会时间超限。


代码

using namespace std;
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cassert>
#define N 10010
#define maxW 5010
inline int input(){
    char ch=getchar();
    while (ch<'0' || '9'<ch)
        ch=getchar();
    int x=0;
    do{
        x=x*10+ch-'0';
        ch=getchar();
    }
    while ('0'<=ch && ch<='9');
    return x;
}
int n,W,K,a[N],b[N];
int f[N][maxW];
int st1[N],top1,st2[N],top2;
inline void update(int &a,int b){a>b?a=b:0;}
inline bool ok(int j){
    int jj=st1[top1];
    for (int k=0;k<=W;++k)
        if (f[jj][k]+f[j][W-k]<=K)
            return 1;
    return 0;
}
int main(){
    freopen("cactus.in","r",stdin);
    freopen("cactus.out","w",stdout);
    n=input(),W=input(),K=input();
    for (int i=1;i<=n;++i)
        a[i]=input(),b[i]=input();
    int ans=n+1;
    f[0][0]=0;
    for (int i=1;i<=W;++i)
        f[0][i]=K+1;
    for (int i=1,j=1;i<=n;++i){
        if (ok(st2[top2]))
            ans=j-i;
        for (;j<=n && j-i+1<ans;++j){
            st2[++top2]=j;
            int lst=st2[top2-1];
            memcpy(f[j],f[lst],sizeof(int)*(W+1));
            for (int k=0;k+a[j]<=W;++k)
                update(f[j][k+a[j]],f[lst][k]+b[j]);
            if (ok(j))
                ans=j-i+1;
        }
        if (!top1){
            for (int j=top2;j>=1;--j)
                st1[++top1]=st2[j];
            top2=0;
            for (int j=1;j<top1;++j){
                int now=st1[j],lst=st1[j-1];
                memcpy(f[now],f[lst],sizeof(int)*(W+1));
                for (int k=0;k+a[now]<=W;++k)
                    update(f[now][k+a[now]],f[lst][k]+b[now]);
            }
        }
        top1--;
    }
    if (ans==n+1)
        printf("-1\n");
    else
        printf("%d\n",ans);
    return 0;
}

总结

还有这么骚的栈操作……
这告诉我们有时候维护队列的东西可以用两个栈来搞。

标签:背包,15,2019.9,int,栈顶,NOIP2019,队列,ans,include
来源: https://www.cnblogs.com/jz-597/p/11536140.html