BZOJ 3235: [Ahoi2013]好方的蛇

标签(空格分隔): OI-BZOJ OI-DP OI-容斥原理


Time Limit: 10 Sec

Memory Limit: 64 MB


Description

有一天,可爱的蛇心花怒放,把自己变成了一个正方形!但是她改变的时候

被induce了导致改变出了些问题....

按照预设,她应该变成一个N*N的全黑正方形,但是这个正方形出现了一些白的格子...现在她的身体不幸出了些小反应,定义一个subsnake是一个至少有两格的全黑矩形。

现在蛇想让你帮忙求一下一共有多少对不相交的subsnake,答案模10007。

Input

第一行一个整数 N, 接下来N行,每行一个长度为N的字符串,如果是B,那么是黑的,如果是 W那么是白的。

Output

一行一个整数,表示答案

Sample Input

3

BBW

BBW

BWW

Sample Output

5

HINT

N<=1000


Solution####

dp计数题

首先可以用单调栈维护以某个点为右下角可能的矩形个数设为sum[i][j]

求出f[i][j]表示在[1-i][1-j]内的矩形个数,转移为

f[i][j]=f[i-1][j]+f[i][j-1]-f[i-1][j-1]+sum[i][j];

这样可以用f[i][n]表示下边界在i以内的矩形个数。

也可求出上边界为i的矩形个数

然后可以算出,可以被一条横线分开的矩形个数

同理可以算出,可以被一条竖线分开的矩形个数

发现这样2种情况会被统计多次:

1100

1100

0011

0011



0011

0011

1100

1100

对于第一种情况可以枚举右下正方形的左上端点然后用f统计

第二种类似。


Code####

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<bitset>
#include<vector>
using namespace std;
#define PA pair<int,int>
int read()
{
int s=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){s=(s<<1)+(s<<3)+ch-'0';ch=getchar();}
return s*f;
}
int n,mo=10007,ans;
bool p[1005][1005];
int f[1005][1005],g[1005][1005],u[1005],sum;
int s1[1005],s2[1005],s3[1005],tot;
char z[1005];
int main()
{
//freopen(".in","r",stdin);
//freopen(".out","w",stdout);
n=read();
for(int i=1;i<=n;i++)
{scanf("%s",z);
for(int j=1;j<=n;j++)
p[i][j]=(z[j-1]=='B');
}
for(int i=1;i<=n;i++)u[i]=0;
for(int i=1;i<=n;i++)
{for(int j=1;j<=n;j++)
u[j]=(p[i][j]?u[j]+1:0);
tot=sum=0;
for(int j=1;j<=n;j++)
{int k=0;
while(tot&&s1[tot]>u[j])k+=s2[tot],sum-=s3[tot--];
tot++;k++;
s1[tot]=u[j];s2[tot]=k;s3[tot]=u[j]*k;
sum+=s3[tot]-p[i][j];
f[i][j]=f[i-1][j]+f[i][j-1]-f[i-1][j-1]+sum;f[i][j]%=mo;
sum+=p[i][j];
}
} for(int i=1;i<=n;i++)u[i]=0;
for(int i=1;i<=n;i++)
{for(int j=1;j<=n;j++)
u[j]=(p[i][j]?u[j]+1:0);
tot=sum=0;
for(int j=n;j>=1;j--)
{int k=0;
while(tot&&s1[tot]>u[j])k+=s2[tot],sum-=s3[tot--];
tot++;k++;
s1[tot]=u[j];s2[tot]=k;s3[tot]=u[j]*k;
sum+=s3[tot]-p[i][j];
g[i][j]=g[i-1][j]+g[i][j+1]-g[i-1][j+1]+sum;g[i][j]%=mo;
sum+=p[i][j];
}
} for(int i=1;i<=n;i++)u[i]=0;
for(int i=n;i>=1;i--)
{for(int j=1;j<=n;j++)
u[j]=(p[i][j]?u[j]+1:0);
tot=sum=0;
for(int j=n;j>=1;j--)
{int k=0;
while(tot&&s1[tot]>u[j])k+=s2[tot],sum-=s3[tot--];
tot++;k++;
s1[tot]=u[j];s2[tot]=k;s3[tot]=u[j]*k;
sum+=s3[tot]-p[i][j];
ans+=sum*f[n][j-1]+sum*f[i-1][n]-sum*f[i-1][j-1];ans%=mo;
sum+=p[i][j];
}
} for(int i=1;i<=n;i++)u[i]=0;
for(int i=n;i>=1;i--)
{for(int j=1;j<=n;j++)
u[j]=(p[i][j]?u[j]+1:0);
tot=sum=0;
for(int j=1;j<=n;j++)
{int k=0;
while(tot&&s1[tot]>u[j])k+=s2[tot],sum-=s3[tot--];
tot++;k++;
s1[tot]=u[j];s2[tot]=k;s3[tot]=u[j]*k;
sum+=s3[tot]-p[i][j];
ans-=sum*g[i-1][j+1];ans%=mo;
sum+=p[i][j];
}
}
cout<<(ans+mo)%mo<<endl;
//fclose(stdin);
//fclose(stdout);
return 0;
}

BZOJ 3235: [Ahoi2013]好方的蛇的更多相关文章

  1. 【BZOJ 3235】 3235: [Ahoi2013]好方的蛇 (单调栈+容斥原理)

    3235: [Ahoi2013]好方的蛇 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 187  Solved: 95 Description 有一天, ...

  2. 3235: [Ahoi2013]好方的蛇

    3235: [Ahoi2013]好方的蛇 链接 分析: 可以求出以每个点为顶点的满足条件的矩形有多少个,单调栈求.设为sum. 然后对这个数组进行二维前缀和,可以求出每个矩阵内,以右下角.左下角为端点 ...

  3. BZOJ3235 [Ahoi2013]好方的蛇 【单调栈 + dp】

    题目链接 BZOJ3235 题解 求出每个点为顶点,分别求出左上,左下,右上,右下的矩形的个数\(g[i][j]\) 并预处理出\(f[i][j]\)表示点\((i,j)\)到四个角的矩形内合法矩形个 ...

  4. BZOJ 3233: [Ahoi2013]找硬币

    BZOJ 3233: [Ahoi2013]找硬币 标签(空格分隔): OI-BZOJ OI-DP Time Limit: 10 Sec Memory Limit: 64 MB Description ...

  5. BZOJ 3233: [Ahoi2013]找硬币( dp )

    dp(x)表示最大面值为x时需要的最少硬币数. 枚举x的质因数p,  dp(x) = min( dp(x/p) - (p-1) * sigma[a[i]/x] ). ----------------- ...

  6. BZOJ 3238: [Ahoi2013]差异 [后缀数组 单调栈]

    3238: [Ahoi2013]差异 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 2326  Solved: 1054[Submit][Status ...

  7. BZOJ 3237: [Ahoi2013]连通图

    3237: [Ahoi2013]连通图 Time Limit: 20 Sec  Memory Limit: 512 MBSubmit: 1161  Solved: 399[Submit][Status ...

  8. BZOJ 3236: [Ahoi2013]作业

    3236: [Ahoi2013]作业 Time Limit: 100 Sec  Memory Limit: 512 MBSubmit: 1393  Solved: 562[Submit][Status ...

  9. Bzoj 3236: [Ahoi2013]作业 莫队,分块

    3236: [Ahoi2013]作业 Time Limit: 100 Sec  Memory Limit: 512 MBSubmit: 1113  Solved: 428[Submit][Status ...

随机推荐

  1. hive 存储格式及压缩

    -- 设置参数 set hivevar:target_db_name=db_dw; use ${hivevar:target_db_name}; -- 创建textfile表 create table ...

  2. centos7-根据端口号查看所属应用

    [root@cent7-zuoys ~]# netstat -lnp | grep 8080 [root@cent7-zuoys ~]# ps 4735

  3. gulp不压缩打包layui

    从网上下载的layui都是压缩包,如何打包在一个文件且不压缩呢?如下方法: 1.https://gitee.com/sentsin/layui下载源码(本文的为2.4.5版本) 2.安装nodejs( ...

  4. CSS background 属性全家桶

    介绍我们都知道css的background属性是一个复合属性,可以简写成一行代码,也可以将每个属性分开来写. background 简写属性在一个声明中设置所有的背景属性.如:body{ backgr ...

  5. 1.1 Go安装与项目结构初始化

    软件安装安装包下载地址为:https://golang.org/dl/ 如果打不开可以: https://golang.google.cn/dl/ https://dl.google.com/go/g ...

  6. py---------面向对象进阶

    一.isinstance 和 issubclass isinstance(obj,cls)检查obj是否是类cls的对象,是则返回True class Foo(object): pass obj = ...

  7. Markdown emoji

    由于网页解析原因,部分符号无法正常显示 People emoji code emoji code :bowtie: :bowtie:

  8. Neutron命令测试3

    1.打开Ubuntu的/etc/network/interfaces文件 .默认内容如下 auto lo iface lo inet loopback 2.如果以DHCP方式配置网卡,则改为:auto ...

  9. Murano Weekly Meeting 2015.08.25

    Meeting time: 2015.August.25th 1:00~2:00 Chairperson:  Serg Melikyan, PTL from Mirantis Meeting summ ...

  10. php防止重复提交问题总结

    用户提交表单时可能因为网速的原因,或者网页被恶意刷新,致使同一条记录重复插入到数据库中,这是一个比较棘手的问题.我们可以从客户端和服务器端一起着手,设法避免同一表单的重复提交. 1.使用客户端脚本 提 ...