【LGR-052】洛谷9月月赛II（加赛）

题解：

没打。。。

ab题满世界都过了应该没什么意思

c题是个比较有意思的思维题（先看了题解才会的。。。）

我们考虑这么一件事情

没钥匙的人出门后

门一定是开着的

他进来的时候，门一定是开着的

其他时候，一定门能关着比较好

但是这个东西并不是很好维护，因为可能重复算

我们考虑重新排个序

我们把a[]在b[]前面的那个排在他前面

这样我们只需要记录一下上一个取了没有就可以了

感觉说出来不是很好理解这个东西。。。

d题还是个有点难度的数据结构题

一件比较显然的事情是我们对a的操作等价于是要

找当前x修改后是不是前缀最大值，以及他能延伸到哪（这个权值线段树显然可以维护）

然后就变成了区间覆盖成一个值了

然后 刚开始就想错了。。。

刚开始想的是去考虑这一段的值是怎么通过原先值变化的

其实这样反而不好做

直接去暴力统计现在的值

那么我们要知道的其实就是 当前区间>k的数的个数以及<k的数的乘积

这玩意我好像只会树套树或者cdq分治+排序+线段树。。。(写起来挺烦的啊）

于是复杂度就是nlog^2n了

不知道正解是怎么玩的。。。

代码：

c:

#include <bits/stdc++.h>
using namespace std;
#define rint register int
#define IL inline
#define rep(i,h,t) for (rint i=h;i<=t;i++)
#define dep(i,t,h) for (rint i=t;i>=h;i--)
const int N=;
const int INF=1e9;
struct re{
  int a,b,c,pos;
}a[N],b[N],c1[N],c2[N],p[N];
int pre[N],pos[N],x1[N],x2[N],f[N][N][];
bool tt[N];
bool cmp(re x,re y)
{
  return x.b<y.b;
}
int main()
{
  freopen("1.in","r",stdin);
  freopen("1.out","w",stdout);
  ios::sync_with_stdio(false);
  int n,m;
  cin>>n>>m;
  rep(i,,n)
  {
    cin>>a[i].a>>a[i].b; a[i].pos=i;
    b[i*-].a=i*-; b[i*-].b=a[i].a;
    b[i*].a=i*; b[i*].b=a[i].b;
  }
  sort(a+,a+n+,cmp);
  sort(b+,b+*n+,cmp);
  rep(i,,n) pos[a[i].pos]=i;
  rep(i,,n*)
  {
    if (b[i].a%==&&b[i-].a%==)
      pre[b[i].a/]=(b[i-].a+)/;
    if (b[i].a%==) x2[b[i].a/]=b[i].b-b[i-].b;
    if (b[i].a%==) x1[(b[i].a+)/]=b[i+].b-b[i].b;
  }
  int cnt=;
  rep(i,,n)
    if (!tt[i])
    {
      int x=i;
      p[++cnt].a=a[x].a; p[cnt].b=a[x].b; p[cnt].c=; tt[i]=;
      if (pos[pre[a[x].pos]]!=x)
        c1[cnt].a=x1[a[x].pos],c1[cnt].b=x2[a[x].pos];
      else
        c1[cnt].a=x1[a[x].pos],c1[cnt].b=;
      while (pos[pre[a[x].pos]]&&pos[pre[a[x].pos]]!=x)
      {
        x=pos[pre[a[x].pos]];
        tt[x]=;
        p[++cnt].a=a[x].a; p[cnt].b=a[x].b; p[cnt].c=;
        c1[cnt].a=x1[a[x].pos]; c1[cnt].b=x2[a[x].pos];
      }
    }
  rep(i,,n)
    rep(j,,m)
    {
      f[i][j][]=f[i][j][]=INF;
      if (p[i].c)
      {
        if (j) f[i][j][]=min(f[i][j][],min(f[i-][j-][],f[i-][j-][]));
        f[i][j][]=min(f[i][j][],f[i-][j][]+c1[i].b);
        f[i][j][]=min(f[i][j][],f[i-][j][]+c1[i].b+c1[i].a);
      } else
      {
        if (j) f[i][j][]=min(f[i][j][],min(f[i-][j-][],f[i-][j-][]));
        f[i][j][]=min(f[i][j][],f[i-][j][]+c1[i].b+c1[i].a);
        f[i][j][]=min(f[i][j][],f[i-][j][]+c1[i].b+c1[i].a);
      }
    }
  cout<<min(f[n][m][],f[n][m][]);
  return ;
}