Educational Codeforces Round 98 (Rated for Div. 2) E. Two Editorials 细节题

题目连接
真细节题… 选了一个很复杂的做法。

枚举第一个作者$A$选的框，枚举选手$i$。你会发现另一个作者$B$选的框对于选手$i$来说，会有一些右端点连续的区间比$A$更优。那这样就可以知道作者$B$选哪些框更优并且也能用二次差分统计出优多少。
搞完所有选手后，找出作者$B$选哪个框最优即可。

二次差分的模板：

void add(int l,int r,int op,int d){    //区间[l,r] 加上op为首项，d为公差的等差数列                  
  tw[l]+=d;
  tw[r+1]-=d;
  cat[l]+=op-d;
  cat[r+1]-=op-d;
  cat[r+1]-=d*(r+1-l);
}

分类讨论的细节还希望读者可以自己多想想，理清楚。。。。。（因为我也是因为细节问题 debug了很久…）

细节见代码：

//#pragma GCC optimize(2)
//#pragma GCC optimize(3)
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N = 2e3 + 10;
#define fi first
#define se second
#define pb push_back
#define wzh(x) cerr<<#x<<'='<<x<<endl;
int n,m,k;
struct uzi{
  int l,r;
}p[N];
int a[N][N],L[N][N],R[N][N];
int cat[N],tw[N];
int F[N],G[N];
void add(int l,int r,int op,int d){                     
  tw[l]+=d;
  tw[r+1]-=d;
  cat[l]+=op-d;
  cat[r+1]-=op-d;
  cat[r+1]-=d*(r+1-l);
}
void solve(int x,int y,int z){//多的
  int mx=min(k,p[x].r-p[x].l+1);
  if(L[x][y+1]<L[x][mx]){//左边等差
    int l=L[x][mx-1];
    add(L[x][y+1],l,1,1);
  }else if(a[F[x]][x]>y&&a[F[x]][x]!=mx){
  	add(F[x],L[x][mx]-1,a[F[x]][x]-y,1);
  }
  add(L[x][mx],R[x][mx],mx-y,0);//中间一段恒等 
  if(G[x]>R[x][mx]&&y+1<mx&&R[x][y+1]>R[x][mx]){//第三段等差
    int dx=R[x][mx-1];
    int dy=R[x][y+1];
  	add(dx,dy,mx-y-1,-1);  
  }else if(a[G[x]][x]>y && a[G[x]][x]!=mx){
  	add(R[x][mx]+1,G[x],mx-y-1,-1);
  }
}
void solve_(int x){
	int l=F[x],r=G[x];
    int mx=min(k,p[x].r-p[x].l+1);
	int dx=L[x][mx],dy=R[x][mx];
	add(dx,dy,mx,0);
	if(l<dx)add(l,dx-1,a[l][x],1);
	if(r>dy)add(dy+1,r,mx-1,-1);
}
int main() {
  ios::sync_with_stdio(false);
  cin>>n>>m>>k;
  for(int i=1;i<=m;i++)cin>>p[i].l>>p[i].r;
  for(int i=k;i<=n;i++){
    int x=i-k+1,y=i;
    for(int j=1;j<=m;j++){
      int l=p[j].l,r=p[j].r;
      if(r<x||l>y)continue;
      l=max(l,x);
      r=min(r,y);
      a[i][j]=r-l+1;
    }
  }
  for(int i=1;i<=m;i++){
    for(int j=k;j<=n;j++){
      if(a[j][i]){
        if(!F[i])F[i]=j;
        G[i]=j;
        if(!L[i][a[j][i]])L[i][a[j][i]]=j;
        R[i][a[j][i]]=j;
      }
  
    }
  }
  int ok=0;
  for(int i=k;i<=n;i++){
    int res=0;
    for(int j=1;j<=m;j++)res+=a[i][j];
    for(int j=1;j<=m;j++){
      int now=a[i][j];
      if(now==k)continue;
      if(now==p[j].r-p[j].l+1)continue;
      if(now)solve(j,now,i);
      else solve_(j);
    }
    int pp=0;
    for(int j=1;j<=n;j++){
      tw[j]+=tw[j-1];
      cat[j]+=cat[j-1]+tw[j];
      pp=max(pp,cat[j]);
    }
    for(int j=0;j<=n+1;j++)cat[j]=tw[j]=0;
    res+=pp;
    ok=max(ok,res);
  }
  cout<<ok<<'\n';
  return 0;
}