【BZOJ-2669】局部极小值 状压DP + 容斥原理
2669: [cqoi2012]局部极小值
Time Limit: 3 Sec Memory Limit: 128 MB
Submit: 561 Solved: 293
[Submit][Status][Discuss]
Description
Input
Output
Sample Input
X.
..
.X
Sample Output
HINT
Source
Solution
这道题有点劲!自己没想出来,于是看的论文 传送门 下面引用论文里的题解
对于一个合法的数填写方案,其中的数放置的顺序对其是没有影响的,于是我们可以从1开始填数,并且一个一个地填进格子。如果采取这样的做法,那么所有的“X”必然要在其周边所有的格子填数之前就填好一个数,而"X"有多少呢?很显然最多只有8个而已。这时我们就可以想到这样的一个状态压缩方式:opt[i][j](j是一个二进制表达)表示的是i及其以后的数还没有填进格子,被填写了数的“X”集合状态为j的情况下的方案数。
如上图4*7的矩阵中,红色的"X"表示已经填写数的"X",红色的格子表示已经填写数的非"X"格子,那么可以表述成这样的状态opt[8][num](8表示已经填写了7个数,下一个填写8,num是011010的表示,含义是第2、3、5个"X"已经填写了数了)
如果我们转移的话就会有两种情况:
第一种情况就是把i填进一个"X"中,这个显然只要枚举一下放哪一个"X",然后把这个"X"加入j表示的集合里就可以了。
如上图,下一步我们填写"X"是可以随意的,因为只要存在解,任意的"X"都是互不影响的。当前的状态为f[8][num1](num1为011010的表示),可以推导到f[9][num2](num2为111010、011110、011011的表示)。
第二种情况就是把i填进一个非"X"中,这样的选择就有很多了。对于全图我们一共有n*m个格子,若没有填进去数的"X"格子以及其周边的格子共有tot个,显然这tot个格子都是不能填i的(因为填进的是一个非"X",并且一个没有填进去数的"X"格子其周边因为都要比它小,所以这两者都不可以填i),又因为已经填写了1到i-1所有的数,所以剩下能填的选择数就是n*m-tot-(i-1)。
如上图,所有的蓝色区域都是无法填写i的,而下一步能填写的格子就只有白色的格子,即4*7-17-7=4个格子。
由于这样的处理方式,尤其是第二种转移可能会导致非"X"点变为最小值,所以还需要使用容斥原理来解决。
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define LL long long
#define P 12345678
int N,M,ANS,bin[],f[][],cnt[];
char mp[][];
int dx[]={-,-,-,,,,,,,},dy[]={-,,,-,,,-,,,};
inline bool OK(int x,int y) {return x>= && x<=N && y>= && y<=M;}
#define Pa pair<int,int>
Pa stack[]; int top;
bool visit[][];
inline void PreWork()
{
top=;
for (int i=; i<=N; i++)
for (int j=; j<=M; j++)
if (mp[i][j]=='X') stack[top++]=make_pair(i,j);
for (int i=; i<bin[top]; i++)
{
cnt[i]=; memset(visit,,sizeof(visit));
for (int j=; j<top; j++) if (~i&bin[j]) visit[stack[j].first][stack[j].second]=;
for (int j=; j<=N; j++)
for (int k=; k<=M; k++)
if (!visit[j][k])
{
bool flag=;
for (int d=,x,y; d<= && flag; d++)
x=j+dx[d],y=k+dy[d],flag=!visit[x][y];
cnt[i]+=flag;
}
}
}
inline int DP()
{
PreWork();
memset(f,,sizeof(f));
f[][]=;
for (int i=; i<=N*M; i++)
for (int j=; j<bin[top]; j++)
{
for (int k=; k<top; k++)
if (j&bin[k]) (f[i][j]+=f[i-][j^bin[k]])%=P;
(f[i][j]+=(LL)f[i-][j]*(cnt[j]-(i-))%P)%=P;
}
return f[N*M][bin[top]-];
}
inline void DFS(int dep,int x,int y)
{
if (y==M+) {DFS(dep,x+,); return;}
if (x==N+) {(ANS+=(LL)DP()*(dep&? -:)%P)%=P; return;}
DFS(dep,x,y+);
bool flag=;
for (int i=; i<= && flag; i++)
if (mp[x+dx[i]][y+dy[i]]=='X') flag=;
if (flag) mp[x][y]='X',DFS(dep+,x,y+),mp[x][y]='.';
}
int main()
{
bin[]=; for (int i=; i<=; i++) bin[i]=bin[i-]<<;
scanf("%d%d",&N,&M);
for (int i=; i<=N; i++) scanf("%s",mp[i]+);
DFS(,,);
printf("%d\n",(ANS+P)%P);
return ;
}
菜鸡.jpg
【BZOJ-2669】局部极小值 状压DP + 容斥原理的更多相关文章
- BZOJ2669 [cqoi2012]局部极小值 状压DP 容斥原理
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ2669 题意概括 有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次.如果一个格子比所 ...
- BZOJ 2669 CQOI2012 局部极小值 状压dp+容斥原理
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2669 题意概述:实际上原题意很简洁了我就不写了吧.... 二话不说先观察一下性质,首先棋盘 ...
- 【uoj#37/bzoj3812】[清华集训2014]主旋律 状压dp+容斥原理
题目描述 求一张有向图的强连通生成子图的数目对 $10^9+7$ 取模的结果. 题解 状压dp+容斥原理 设 $f[i]$ 表示点集 $i$ 强连通生成子图的数目,容易想到使用总方案数 $2^{sum ...
- 【bzoj2560】串珠子 状压dp+容斥原理
题目描述 有 $n$ 个点,点 $i$ 和点 $j$ 之间可以连 $0\sim c_{i,j}$ 条无向边.求连成一张无向连通图的方案数模 $10^9+7$ .两个方案不同,当且仅当:存在点对 $(i ...
- BZOJ 2064: 分裂( 状压dp )
n1+n2次一定可以满足..然后假如之前土地集合S1的子集subs1和之后土地集合S2的子集subs2相等的话...那么就少了2个+操作...所以最后答案就是n1+n2-少掉的最多操作数, 由状压dp ...
- 4455: [Zjoi2016]小星星|状压DP|容斥原理
OrzSDOIR1ak的晨神 能够考虑状压DP枚举子集,求出仅仅保证连通性不保证一一相应的状态下的方案数,然后容斥一下就是终于的答案 #include<algorithm> #includ ...
- 【BZOJ 2669】 2669: [cqoi2012]局部极小值 (状压DP+容斥原理)
2669: [cqoi2012]局部极小值 Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 667 Solved: 350 Description 有一 ...
- bzoj2669 [cqoi2012]局部极小值 状压DP+容斥
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=2669 题解 可以发现一个 \(4\times 7\) 的矩阵中,有局部最小值的点最多有 \(2 ...
- HDU5838 Mountain(状压DP + 容斥原理)
题目 Source http://acm.hdu.edu.cn/showproblem.php?pid=5838 Description Zhu found a map which is a N∗M ...
随机推荐
- linux版基金看板
程序员的吊丝们,还在害怕上班时偷偷看基金被老板发现吗?今天你们的福利来了,专属程序员吊丝一族的礼物,linux版基金看板. 优点: 1.自定义设置关注基金 2.linux系统,让别人可以以为你一直都在 ...
- crontab不执行perl脚本分析
在新装的Linux服务器上部署了一个作业监控磁盘空间并提前告警,在shell脚本里面调用了一个perl脚本发送告警邮件.结果出现了一个很奇怪的现象:如果手工执行该脚本/home/oracle/scri ...
- ORA-02292: integrity constraint (xxxx) violated - child record found
在更新表的主键字段或DELETE数据时,如果遇到ORA-02292: integrity constraint (xxxx) violated - child record found 这个是因为主外 ...
- C语言核心之数组和指针详解
指针 相信大家对下面的代码不陌生: int i=2; int *p; p=&i;这是最简单的指针应用,也是最基本的用法.再来熟悉一下什么是指针:首先指针是一个变量,它保存的并不是平常的数据,而 ...
- (转)Div和Table的区别
原文:http://www.cnblogs.com/lovebear/archive/2012/04/18/2456081.html Div与table的区别 1:速度和加载方式方面的区别 div 和 ...
- 解惑spring嵌套事物
工作中一直对spring中的事物管理都是最简单的配置 但是spring中的事物传播性配置 还有很多种,有时候经常疑惑service调用service的问题,今天的论坛上看到一篇写的非常详细的文字.记录 ...
- 【转帖】分享一个迅为4412开发板OTG烧录批处理文件
平台:iTOP-4412开发板 Bat 功能: 1.可以分条的执行烧录,不需要每次烧录都去复制命令 2.可以批量烧录 开发板系统烧录批处理文件,请将此文件放置在fastboot程序同目录下,下载地址: ...
- Nodejs以后台服务启动
1: 从网上查找 LINUX中我们可以使用这种简单的方式让node.js在后台运行: nohup node your_app.js & 经多次实验一直没有成功 2:使用 forever ...
- hdu2082 找单词 (母函数)
找单词 题意: 中文题,考虑是不是要写个英文题意..(可惜英语水平不够 囧rz) (题于文末) 知识点: 母函数(生成函数): 生成函数有普通型生成函数和指数型生成函数 ...
- 基于移动端Reactive Native轮播组件的应用与开发详解
总结下这段时间学习reactive native的一些东西,我们来认识一下,被炒得这么火的rn,究竟是个什么东西,以及如何去搭建自己的demo. reactive native是什么 由facebo ...