[USACO18JAN]Sprinklers
作者:互联网
一个矩形要符合什么条件
右上角的右上有点,左下角的左下有点
所以每列的选择高度为一个区间,小于后缀最大值大于前缀最小值(不管是作为右上角还是左下角)
然后对于一个:
求完全在这个区域里的矩形个数
从上往下考虑,处理左下角在这一行,右上角在上面的方案数
发现每行的点对于之后的贡献分别是0+1+2+3。。。
全局变量tot维护决策点个数
删除的时候注意删除对应决策点即可。
// luogu-judger-enable-o2 #include<bits/stdc++.h> #define reg register int #define il inline #define fi first #define se second #define mk(a,b) make_pair(a,b) #define numb (ch^'0') using namespace std; typedef long long ll; template<class T>il void rd(T &x){ char ch;x=0;bool fl=false; while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true); for(x=numb;isdigit(ch=getchar());x=x*10+numb); (fl==true)&&(x=-x); } template<class T>il void output(T x){if(x/10)output(x/10);putchar(x%10+'0');} template<class T>il void ot(T x){if(x<0) putchar('-'),x=-x;output(x);putchar(' ');} template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');} namespace Miracle{ const int N=1e5+5; const int mod=1e9+7; int n; int u[N],d[N]; ll a[N],sum[N]; int h[N]; vector<int>dele[N]; int main(){ rd(n); int x,y; for(reg i=1;i<=n;++i){ rd(x);rd(y); ++x;++y;h[x]=y; } int mx=0; for(reg i=n;i>=1;--i){ mx=max(mx,h[i]); u[i]=mx; } int mi=n+233; for(reg i=1;i<=n;++i){ mi=min(mi,h[i]); d[i]=mi;dele[d[i]-1].push_back(i); } for(reg i=1;i<=n;++i){ ++a[d[i]];--a[u[i]+1]; } for(reg i=1;i<=n;++i){ a[i]+=a[i-1]; } for(reg i=n;i>=1;--i){ sum[i]=sum[i+1]+a[i]; } ll ans=0,tot=0; ll len=0; for(reg i=n;i>=1;--i){ for(reg j=0;j<(int)dele[i].size();++j){ int now=dele[i][j]; len+=u[now]-d[now]+1; tot-=sum[i+1]-len; } ans=(ans+tot)%mod; tot=tot+(a[i]-1)*a[i]/2; } cout<<ans; return 0; } } signed main(){ Miracle::main(); return 0; } /* Author: *Miracle* Date: 2019/4/13 19:58:12 */View Code
标签:USACO18JAN,ch,10,int,void,reg,Sprinklers,define 来源: https://www.cnblogs.com/Miracevin/p/10718670.html