CF 958E2. Guard Duty (medium)

这道题是昨天linkfqy dalao上课讲的一道题

当时他讲的时候就想到了一种玄学的搞法，然后不敢相信自己切掉了

没想到后来CHJ dalao也想到了这种算法，然后发现是对的

后来10min就切掉了

看来以后要活跃一点了

这里主要讲两种算法（想WQS二分这种高深的算法就不讲了）

首先讲题意：一个人它只能在某个时间点去跟看k个他的客人对话，或是结束他和当前这个客人的对话。但是它不能在两个连续的时间段进行讲话，问它最少需要多少时间来和k个客人对话。

1.堆

我们先按时间排序，处理出每一个时间段的长度

然后就建模转化成：n-1个点中取k个并且每两个都不相邻的最少代价和。

这就和一道题很像了，看Luogu P1484&sol

用堆水一水即可

CODE

#include<cstdio>
#include<algorithm>
#include<queue>
using namespace std;
const int N=500005;
int k,n,x[N],a[N],pre[N],nxt[N],ans;
bool vis[N];
struct data
{
    int x,num;
    bool operator <(const data &s) const
    {
        return s.x<x;
    }
};
priority_queue<data> small;
inline char tc(void)
{
    static char fl[100000],*A=fl,*B=fl;
    return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
    x=0; char ch=tc();
    while (ch<'0'||ch>'9') ch=tc();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
int main()
{
    //freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
    register int i; read(k); read(n);
    for (i=1;i<=n;++i)
    read(x[i]); 
    sort(x+1,x+n+1);
    for (i=2;i<=n;++i)
    {
        a[i-1]=x[i]-x[i-1];
        pre[i-1]=i-2; nxt[i-1]=i;
        small.push((data){a[i-1],i-1});
    }
    a[0]=a[n]=1e9; pre[n]=n-1; nxt[0]=1;
    while (k--)
    {
        while (vis[small.top().num]) small.pop();
        data p=small.top(); small.pop(); ans+=p.x;
        p.x=a[p.num]=a[pre[p.num]]+a[nxt[p.num]]-a[p.num];
        vis[pre[p.num]]=vis[nxt[p.num]]=1;
        pre[p.num]=pre[pre[p.num]]; nxt[pre[p.num]]=p.num; 
        nxt[p.num]=nxt[nxt[p.num]]; pre[nxt[p.num]]=p.num;
        small.push(p);
    }
    printf("%d",ans);
    return 0;
}

2.DP

首先这种题目就是看一眼就想到DP的吧

我们令f[i][j][0/1]表示前i个点中取j个并且第i个点取or不取（1表示取，2表示不取）的最少代价

则很容易写出转移：

f[i][j][1]=f[i-1][j-1][0]+a[i]
f[i][j][0]=min(f[i-1][j][0],f[i-1][j][1])

首先我们发现i的那一维可以滚动存储优化掉

然后空间就够了，但是时间为O(nk)

~~然后坚信CF的机子是无敌的满怀信仰地交上去说不定就过了~~

然后这里我们就要讲一个结论了：

无论怎么取，都只可能取到3*k小的数

Why？很简单，如果所有的数左右都取，算上它自己，也就3个，扩展开就是3*k

所以我们就把时间优化到了O(k^2)

CODE

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=500005,K=5005;
int x[N],a[N],n,k,f[2][K][2],pre,num,tot;
inline char tc(void)
{
    static char fl[100000],*A=fl,*B=fl;
    return A==B&&(B=(A=fl)+fread(fl,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{
    x=0; char ch=tc();
    while (ch<'0'||ch>'9') ch=tc();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=tc();
}
inline int min(int a,int b)
{
    return a<b?a:b;
}
int main()
{
    //freopen("CODE.in","r",stdin); freopen("CODE.out","w",stdout);
    register int i,j;
    read(k); read(n);
    for (i=1;i<=n;++i)
    read(x[i]); sort(x+1,x+n+1);
    for (i=2;i<=n;++i)
    x[i-1]=a[i-1]=x[i]-x[i-1];
    sort(x+1,x+n); num=3*k<n?x[3*k]:x[n-1]; 
    memset(f[0],63,sizeof(f[0])); f[0][0][0]=f[0][0][1]=0;
    int now=0,lst=1,pre=-1;
    for (i=1;i<n;++i)
    {
        if (a[i]>num) continue;
        now^=1; lst^=1; memset(f[now],63,sizeof(f[now])); 
        ++tot;
        if (i-1!=pre) 
        {
            for (j=1;j<=min(tot,k);++j)
            {
                f[now][0][0]=f[now][0][1]=0;
                f[now][j][1]=min(f[lst][j-1][0],f[lst][j-1][1])+a[i];
                f[now][j][0]=min(f[lst][j][0],f[lst][j][1]);
            } 
        } else
        {
            for (j=1;j<=min(tot,k);++j)
            {
                f[now][0][0]=f[now][0][1]=0;
                f[now][j][1]=f[lst][j-1][0]+a[i];
                f[now][j][0]=min(f[lst][j][0],f[lst][j][1]);
            }
        } pre=i;
    }
    printf("%d",min(f[tot&1][k][0],f[tot&1][k][1]));
    return 0;
}

CF 958E2. Guard Duty (medium)

1.堆

2.DP

猜你喜欢