【简解】C2CRNI - Crni
【题目大意】
给定一个N行N列的矩阵,每个格子要么为白色要么为黑色。黑矩形为所涵单元格数大于等于2且所涵单元格均为黑色的矩表。要解决的问题是在给定的矩形中找出两个没有共公部分的黑矩形,输出所有方案数,由于数较大,输出它对10007的余数。
传送门。
【分析】
听说是道很老的套路题,考试考了这道。。。
有个比较明显的东西,就是我们要求的两个黑矩形没有公共部分,那么必然会有一条竖线(或横线)把这两个黑矩形分开,我们就可以枚举这条竖线(横线)来统计答案。
为了避免计算重复,我们需要一边固定,另一边统计总和。并且,还应减去既能被一条竖线分开又能被一条横线分开的两个黑矩形(因为它被算了两次嘛)。
大体思路就是这样,然后就是实现。
比较容易想到的就是计算以每个格子为左下角(右下角、左上角、右上角)的黑矩形有多少个。
考试的时候太弱了,只会写\(O(n^4)\)的暴力。。。(它还挂了)
简单说一点暴力吧,就是\(O(n^2)\)枚举每个格子,然后\(O(n^2)\)枚举矩形的长度,可以用二维前缀和\(O(1)\)来判断枚举的矩形是否合法,然后就\(O(n^4)\)预处理,再\(O(n^2)\)统计贡献。
可以看到,预处理是瓶颈,我们应考虑如何优化预处理的复杂度。
对每一行拆开考虑,我们会发现比较像直方图里的最大矩形,然后就做啊。
我们以计算以每个格子为右下角的黑矩形有多少个为例。
设\(f[i][j]\)就是以\((i,j)\)为右下角的黑矩形有多少个,我们会发现,它可以直接继承它前面出现的第一个小于它高度的位置的答案,那中间那些比它还高的矩形呢,没错,只能扩展它自己的高度。
即\(f[i][j]=f[i][k]+1+h[j]*(j-k)-1\)
其中\(k\)表示从\(j\)向前走第一次出现的\(h[k]<h[j]\)的位置,那么直接可以继承\(f[i][k]\),然后加一,然后对于那块大的矩形里的所有点均可作为当前矩形的左上角,除了自己本身,因为\(1\times 1\)的黑矩形不能算入答案,所以会减一。
考虑到\(k\)的计算,想到单调栈优化,可以做到\(O(1)\)转移。
所以我们做到了\(O(n^2)\)预处理,然后就行了。
细节有点多,还有要注意时刻取模。
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
const int P = 10007;
const int N = 1000 + 5;
inline int read(){
int f = 1, x = 0; char ch;
do { ch = getchar(); if (ch == '-') f = -1; } while (ch < '0' || ch > '9');
do {x = (x << 3) + (x << 1) + ch - '0'; ch = getchar(); } while (ch >= '0' && ch <= '9');
return f * x;
}
inline void hand_in() {
freopen("crni.in", "r", stdin);
freopen("crni.out", "w", stdout);
}
int n, mp[N][N], res;
int f[N][N], g[N][N], h[N];
char ch[N][N]; int sta[N], top;
int main(){
// hand_in();
n = read();
for (int i = 1;i <= n; ++i) scanf("%s", ch[i] + 1);
for (int i = 1;i <= n; ++i) {
for (int j = 1;j <= n; ++j) {
mp[i][j] = (ch[i][j] == 'C');
}
}
// 左上
memset(h, 0, sizeof h);
for (int i = 1;i <= n; ++i) {
top = 0;
sta[top] = 0;
for (int j = 1;j <= n; ++j) {
if (mp[i][j] == 1) h[j] ++;
else h[j] = 0;
if (h[j] == 0) {
top = 0;
sta[top] = j;
continue;
}
int size;
while (top && h[sta[top]] >= h[j]) top --;
size = h[j] * (j - sta[top]) - 1;
f[i][j] = f[i][sta[top]] + size + (top > 0);
f[i][j] %= P;
sta[++top] = j;
}
}
//右下
memset(h, 0, sizeof h);
for (int i = n;i >= 1; --i) {
top = 0;
sta[top] = n + 1;
for (int j = n;j >= 1; --j) {
if (mp[i][j] == 1) h[j] ++;
else h[j] = 0;
if (h[j] == 0) {
top = 0;
sta[top] = j;
continue;
}
int size;
while (top && h[sta[top]] >= h[j]) top --;
size = h[j] * (sta[top] - j) - 1;
g[i][j] = g[i][sta[top]] + size + (top > 0);
g[i][j] %= P;
sta[++top] = j;
}
}
int suma = 0, sumb = 0;
for (int i = 1;i < n; ++i) {
for (int j = 1;j <= n; ++j) {
suma += f[i][j], sumb += g[i + 1][j];
suma %= P, sumb %= P;
}
res += suma * sumb;
res %= P;
sumb = 0;
}
suma = 0, sumb = 0;
for (int j = 1;j < n; ++j) {
for (int i = 1;i <= n; ++i) {
suma += f[i][j], sumb += g[i][j + 1];
suma %= P, sumb %= P;
}
res += suma * sumb;
res %= P;
sumb = 0;
}
for (int i = 1;i <= n; ++i) {
for (int j = 1;j <= n; ++j) {
f[i][j] += f[i - 1][j] + f[i][j - 1] - f[i - 1][j - 1] + P;
f[i][j] %= P;
}
}
for (int i = 1;i < n; ++i) {
for (int j = 1;j < n; ++j) {
res -= (f[i][j] * g[i + 1][j + 1] + P) % P;
while (res < 0) res += P;
}
}
//右上
memset(f, 0, sizeof f);
memset(h, 0, sizeof h);
for (int i = 1;i <= n; ++i) {
top = 0;
sta[top] = n + 1;
for (int j = n;j >= 1; --j) {
if (mp[i][j] == 1) h[j] ++;
else h[j] = 0;
if (h[j] == 0) {
top = 0;
sta[top] = j;
continue;
}
int size;
while (top && h[sta[top]] >= h[j]) top --;
size = h[j] * (sta[top] - j) - 1;
f[i][j] = f[i][sta[top]] + size + (top > 0);
f[i][j] %= P;
sta[++top] = j;
}
}
//左下
memset(g, 0, sizeof g);
memset(h, 0, sizeof h);
for (int i = n;i >= 1; --i) {
top = 0;
sta[top] = 0;
for (int j = 1;j <= n; ++j) {
if (mp[i][j] == 1) h[j] ++;
else h[j] = 0;
if (h[j] == 0) {
top = 0;
sta[top] = j;
continue;
}
int size;
while (top && h[sta[top]] >= h[j]) top --;
size = h[j] * (j - sta[top]) - 1;
g[i][j] = g[i][sta[top]] + size + (top > 0);
g[i][j] %= P;
sta[++top] = j;
}
}
for (int i = 1;i <= n; ++i) {
for (int j = n;j >= 1; --j) {
f[i][j] += f[i - 1][j] + f[i][j + 1] - f[i - 1][j + 1] + P;
f[i][j] %= P;
}
}
for (int i = 1;i < n; ++i) {
for (int j = 2;j <= n; ++j) {
res -= (f[i][j] * g[i + 1][j - 1] + P) % P;
while (res < 0) res += P;
}
}
printf("%d", res);
return 0;
}
【简解】C2CRNI - Crni的更多相关文章
- python ConfigParser、shutil、subprocess、ElementTree模块简解
ConfigParser 模块 一.ConfigParser简介ConfigParser 是用来读取配置文件的包.配置文件的格式如下:中括号“[ ]”内包含的为section.section 下面为类 ...
- 【题解】C2Crni - Crni [COCI2010] [SP7884]
[题解]C2Crni - Crni [COCI2010] [SP7884] 传送门:\(\text{C2Crni - Crni}\) \(\text{[COCI2010]}\) \(\text{[SP ...
- AC题目简解-数据结构
A - Japan POJ 3067 要两条路有交叉,(x1,y1)(x2,y2)那么需要满足:(x1-x2)*(y1-y2)<0判断出这是求逆序的问题 树状数组求逆序,先通过自定义的比较器实 ...
- UE4 RHI与Render模块简解
UE4中的RHI指的是Render hardware interface,作用像Ogre里的RenderSystem,针对Dx11,Dx12,Opengl等等平台抽象出相同的接口,我们能方便能使用相同 ...
- zabbix基本监控各指标简解
监控项目及使用模板 监控http和https: Template App HTTP Service Template App HTTPS Service 监控cpu,内存,网络等: Templ ...
- (转)FFMPEG解码H264拼帧简解
http://blog.csdn.net/ikevin/article/details/7649095 H264的I帧通常 0x00 0x00 0x00 0x01 0x67 开始,到下一个帧头开始之前 ...
- Spring ApplicationContext 简解
ApplicationContext是对BeanFactory的扩展,实现BeanFactory的所有功能,并添加了事件传播,国际化,资源文件处理等. configure locations:(C ...
- HTTP协议简解
1.什么是http协议 http协议: 浏览器客户端 与 服务器端 之间数据传输的规范 2.查看http协议的工具 1)使用火狐的firebug插件(右键->查看元素->网络) 2)使用 ...
- linux netstat 命令简解
Netstat 简介: Netstat是在内核中访问网络及相关信息的程序,它能提供TCP连接,TCP和UDP监听,进程内存管理的相关报告.常见参数-a (all)显示所有选项,默认不显示LISTEN相 ...
随机推荐
- Shell 中eval的用法
test.sh:pipe="|"eval ls $pipe wc -l 输出bogon:Desktop macname$ ./test.sh 45 test.sh:eval ech ...
- yum本地安装rrdtool报错,提示版本冲突的解决办法
[root@mysql-mon40 mm]# yum localinstall rrdtool-1.4.7-1.el6.rfx.x86_64.rpm -y 错误: Multilib version p ...
- Linux里使用rz和sz命令
lrzsz是一个unix通信套件提供的X,Y,和ZModem文件传输协议,官网:http://freecode.com/projects/lrzsz/ windows 需要向centos服务器上传文件 ...
- GCC编译UDF和gdb调试UDF
Fluent版本:19.0 前面我们介绍过使用VS来编译和调试UDF,其实我们也可以用GCC来编译UDF.gdb调试UDF.本次介绍的方法更具有通用性,也适用于Linux下Fluent的UDF编译和调 ...
- html内获取当前文件路径,页面获取当前路径
function getRealPath(){ var curWwwPath = window.document.location.href; var pathName = window.docume ...
- 组合模式( Composite Pattern)
参考文档:http://blog.csdn.net/ai92/article/details/298336 定义: 组合多个对象形成树形结构以表示“整体-部分”的结构层次. 设计动机: 这幅图片我们都 ...
- Siam R-CNN: Visual Tracking by Re-Detection
Siam R-CNN: Visual Tracking by Re-Detection 2019-12-02 22:21:48 Paper:https://128.84.21.199/abs/1911 ...
- 【WPF】修改ListBox的Item的样式
<ListBox> <!-- 数据 --> <ListBoxItem>AAAA</ListBoxItem> <ListBoxItem>BB& ...
- Anaconda + PyCharm + Pytorch
Anaconda 1. 下载Anaconda https://www.anaconda.com/download/ 2. 安装 3. 添加环境变量 Path - C:\Users\Godzilla ...
- phpspreadsheet 中文文档(三) 计算引擎
2019年10月11日13:59:52 使用PhpSpreadsheet计算引擎 执行公式计算 由于PhpSpreadsheet表示内存中的电子表格,因此它还提供公式计算功能.单元格可以是值类型(包含 ...