【TMD模拟赛】上低音号 链表

这道题一看有两个出发现点，一枚举点去找边界，想了一会就Pass了...，二是枚举相框，我们最起码枚举两个边界，然后发现平行边界更好处理，然而仍然只有30分，这个时候就来到了链表的神奇应用，我们枚举上界u，下界d在u的基础之上从下往上枚举，我们每次枚举上界的开始就把上界以下的点建成链表（它的形状大概是在从左到右的基础上对于同一列的从上倒下，就是蛇形），然后让下届去逼近并结算答案，十分巧妙。

关于链表：大概有单向链表，双向链表，以及循环链表，他们作为数据结构的应用十分狭窄，只有在特定情境下才有大用处，所以在oi层面只有双向链表有一定用处，他的删除在记录位置的前提下是O(1)的，查询O(n)，插入的话主要是找到合适的位置（前驱），剩下的就是O(1)了。

#include <cstdio>
#include <cstring>
#include <algorithm>
typedef long long LL;
const int N=;
int pre[N],next[N],mem[N][N],last[N];
int line,column,n,k;
LL ans;
bool ex[N];
struct Point{int x,y;}p[N];
inline bool comp(Point a,Point b){return a.y<b.y||(a.y==b.y&&a.x<b.x);}
int main(){
  scanf("%d%d%d%d",&line,&column,&n,&k);
  for(int i=;i<=n;++i)
    scanf("%d%d",&p[i].x,&p[i].y);
  std::sort(p+,p+n+,comp);
  for(int i=;i<=n;++i)
    mem[p[i].x][++mem[p[i].x][]]=i;
  for(int u=;u<=n;++u){
    int pr=;memset(ex,,sizeof(ex));
    for(int i=u;i<=n;++i)
      for(int j=;j<=mem[i][];++j)
        ex[mem[i][j]]=true,last[mem[i][j]]=line;
    for(int i=;i<=n;++i)
      if(ex[i])pre[i]=pr,next[pr]=i,pr=i;
    pre[n+]=pr,next[pr]=n+;
    for(int i=;i<=k;++i)
      pre[n+i]=n+i-,next[n+i-]=n+i;
    for(int d=n;d>=u;--d){
      for(int i=;i<=mem[d][];++i){
        int nk=mem[d][i];
        for(int j=;j<=k;++j)nk=next[nk];
        if(next[mem[d][i]]<=n&&p[next[mem[d][i]]].x!=d)
          ans+=(LL)(last[next[mem[d][i]]]-d+)*
          (p[next[mem[d][i]]].y-p[pre[next[mem[d][i]]]].y)*
          (nk>n?:column-p[nk].y+),
          last[next[mem[d][i]]]=d-;
        nk=pre[nk],
        ans+=(LL)(last[mem[d][i]]-d+)*
        (p[mem[d][i]].y-p[pre[mem[d][i]]].y)*
        (nk>n?:column-p[nk].y+);
        for(int j=,now=mem[d][i];j<k&&now;++j){
          now=pre[now],nk=pre[nk];
          if(now&&p[now].x!=d&&last[now]>=d)
            ans+=(LL)(last[now]-d+)*
            (p[now].y-p[pre[now]].y)*
            (nk>n?:column-p[nk].y+),
            last[now]=d-;
        }
      }
      for(int i=;i<=mem[d][];++i)
        next[pre[mem[d][i]]]=next[mem[d][i]],
        pre[next[mem[d][i]]]=pre[mem[d][i]];
    }
  }
  printf("%lld",ans);
}