洛谷1950 长方形 （单调栈）

懒得放题目链接了
qwq

（果然我是菜的真实，单调栈都不会，gg）

首先我们看到这个题。
应该会想到就是直接枚举行，然后计算当前行的答案。

那现在，对于每一行来说，如果我们能够维护出\(h[j]\)表示第\(j\)列的最近的一个不合法的位置。
那么实际上就是求一堆矩形的并的一个图形中。
有多少个矩形。

首先考虑暴力，我们可以直接枚举每一列，然后枚举他前面的列进行计算，这个复杂度是\(O(n^3)\)的。
那我们应该怎么去优化他呢？

这时候就需要单调栈了！

我们用单调栈维护一个单调上升的序列。

通过枚举每个列，算以\((i,j)\)这个点为右下角的\(ans\)。

然后我们通过维护一个\(now\)来表示已经处理过的值，每次更新
然后\(ans+=now\)

具体应该怎么做呢？
对于每次加入，我们尝试弹出栈里面的元素，然后删除他和他之前那个元素的位置之间的贡献（也就是他们两个之间的距离乘上栈顶元素的高度！）
（其实就相当于中间一段矩阵都是不合法的，没法选）

最后加入的时候，记得加入当前列与栈顶之间的贡献（就是当前列和栈顶元素之间的距离乘上当前列的\(i-h[j]\)）
表示中间一段矩形都是可以算的。‘

这里运用了一个很巧妙的思路就是每次只处理相邻两个元素的贡献

具体实现还是看代码吧

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
using namespace std;
inline int read()
{
   int x=0,f=1;char ch=getchar();
   while (!isdigit(ch)){if (ch=='-') f=-1;ch=getchar();}
   while (isdigit(ch)){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}
   return x*f;
}
const int maxn = 1e6+1e2;
struct Node
{
  int pos,val;  
};
Node a[maxn];
int n,m;
char s[3010][3010];
int h[maxn];
long long ans,now;
int main()
{
  n=read(),m=read();
  for (int i=1;i<=n;i++) scanf("%s",s[i]+1);
  for (int i=1;i<=n;i++)
  {
     //if (s[i][j]=='*') h[j]=i;
     int top=0;
     now=0;
     for (int j=1;j<=m;j++)
     {
        if (s[i][j]=='*') h[j]=i; //维护每一列的最近的一个不合法点在哪里 
        while (top>=1 && a[top].val<=h[j]) //之所以是>=，是因为h表示点的位置，越小反而越高  
        {
            now-=(a[top].pos-a[top-1].pos)*(i-h[a[top].pos]); //每次减去这一部分的答案 
            top--;
        }
        now+=(j-a[top].pos)*(i-h[j]);///计算当前元素和栈顶之间对ans的贡献 
        a[++top].val=h[j];
        a[top].pos=j;
        ans+=now;
     }
  }
  cout<<ans;
  return 0;
}

转载于:https://www.cnblogs.com/yimmortal/p/10153485.html