【CQOI2012】局部极小值

Description

  有一个\(n\)行\(m\)列的整数矩阵,其中\(1\)到\(nm\)之间的每个整数恰好出现一次。如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点)都小,我们说这个格子是局部极小值。

  给出所有局部极小值的位置,你的任务是判断有多少个可能的矩阵。

答案对\(123456789\)取模。

\(1\leq n\leq 4,1\leq m\leq 7\)。

样例输入:

1 3

.X.

样例输出:

2

我们可以发现一张图上的最小值点最多\(8\)个。于是我们找到所有的最小值点以及与这些点的所有相邻点。

我们考虑从\(1\)到\(nm\)将每个数字填入一个格子中。显然一个格子\((x,y)\)只有在\((x,y)\)周围的最小值点都已经填充了数字的时候它才能填充。

我们设\(size_{S}\)表示最小值填充的情况为\(S\)的时候有多少非最小值点可以填充。

我们考虑从\(nm\)到\(1\)倒着枚举每个数\(i\),然后在枚举集合\(S\)表示\(i+1\)到\(nm\)中已经填充了\(S\)集合中的最小值。然后用组合数搞一搞就好了。

但是我们会发现一个问题,我们算出的答案中的最小值点可能不止题目中给定的点,这种情况是要舍去的。所以我们就\(dfs\)枚举所有可能的最小值点的情况,然后容斥就行了。

代码:

#include<bits/stdc++.h>
#define ll long long
#define N 10
#define M 10 using namespace std;
inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;} const ll mod=12345678;
int n,m;
char mp[N][N];
#define pr pair<int,int>
#define mp(a,b) make_pair(a,b) vector<pr>mn;
int q[15];
ll f[1<<8],g[1<<8];
int num;
ll fac[N*M];
ll C[N*M][N*M];
int sta[N][M];
int dx[]={1,1,0,-1,-1,-1,0,1},dy[]={0,-1,-1,-1,0,1,1,1};
int size[1<<8];
int bin[1<<8]; bool vis[N][M]; ll DP() {
memset(sta,0,sizeof(sta));
memset(f,0,sizeof(f));
memset(size,0,sizeof(size));
int cnt=0;
for(int i=0;i<mn.size();i++) {
int x=mn[i].first,y=mn[i].second;
for(int j=0;j<8;j++) {
int a=x+dx[j],b=y+dy[j];
sta[a][b]|=1<<cnt;
}
cnt++;
}
if(!cnt) return fac[n*m];
for(int T=0;T<1<<cnt;T++) {
for(int i=1;i<=n;i++) {
for(int j=1;j<=m;j++) {
if(sta[i][j]&&(sta[i][j]&T)==sta[i][j]) {
size[T]++;
}
}
}
}
ll ans=0;
f[0]=1;
for(int i=n*m;i>=1;i--) {
memset(g,0,sizeof(g));
for(int T=0;T<1<<cnt;T++) {
if(!f[T]) continue ;
int res=((1<<cnt)-1)^T;
for(int j=0;j<cnt;j++) {
if(T>>j&1) continue ;
int now=size[res]-size[res^(1<<j)];
int emp=n*m-i-bin[T]-(size[(1<<cnt)-1]-size[res]);
if(emp>=now) {
(g[T|(1<<j)]+=C[emp][now]*fac[now]%mod*f[T])%=mod;
}
}
}
for(int T=0;T<1<<cnt;T++) (f[T]+=g[T])%=mod;
}
ans=f[(1<<cnt)-1];
return ans*fac[n*m-size[(1<<cnt)-1]-bin[(1<<cnt)-1]]%mod;
} ll dfs(int x,int y,ll flag) {
if(x>n) {
ll tem=DP();
return DP()*flag%mod;
}
if(y>m) return dfs(x+1,1,flag);
else {
ll ans=0;
if(!vis[x][y]) {
int tag=1;
for(int i=0;i<8;i++) {
int a=x+dx[i],b=y+dy[i];
if(vis[a][b]) tag=0;
}
if(tag) {
vis[x][y]=1;
mn.push_back(mp(x,y));
(ans+=dfs(x,y+1,flag*(mod-1)%mod))%=mod;
vis[x][y]=0;
mn.pop_back();
}
}
(ans+=dfs(x,y+1,flag))%=mod;
return ans;
}
} void out(int S) {
for(int i=0;i<mn.size();i++) cout<<(S>>i&1);
cout<<"\n";
} int main() {
n=Get(),m=Get();
fac[0]=1;
for(int i=1;i<=n*m;i++) fac[i]=fac[i-1]*i%mod;
for(int i=0;i<=n*m;i++)
for(int j=0;j<=i;j++)
C[i][j]=(!j||i==j)?1:(C[i-1][j-1]+C[i-1][j])%mod;
for(int S=0;S<1<<8;S++)
for(int i=0;i<8;i++) bin[S]+=(S>>i&1); for(int i=1;i<=n;i++) {
scanf("%s",mp[i]+1);
}
int flag=0;
for(int i=1;i<=n;i++) {
for(int j=1;j<=m;j++) {
if(mp[i][j]=='X') {
mn.push_back(mp(i,j));
vis[i][j]=1;
if(mp[i-1][j]=='X'||mp[i][j-1]=='X'||mp[i-1][j-1]=='X'||mp[i-1][j+1]=='X') {
flag=1;
}
}
}
}
if(flag) return cout<<0,0;
cout<<dfs(1,1,1);
return 0;
}

【CQOI2012】局部极小值的更多相关文章

  1. bzoj 2669 [cqoi2012]局部极小值 DP+容斥

    2669: [cqoi2012]局部极小值 Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 838  Solved: 444[Submit][Status ...

  2. bzoj2669[cqoi2012]局部极小值 容斥+状压dp

    2669: [cqoi2012]局部极小值 Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 774  Solved: 411[Submit][Status ...

  3. [BZOJ2669] [cqoi2012]局部极小值

    [BZOJ2669] [cqoi2012]局部极小值 Description 有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次.如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点) ...

  4. 【BZOJ 2669】 2669: [cqoi2012]局部极小值 (状压DP+容斥原理)

    2669: [cqoi2012]局部极小值 Time Limit: 3 Sec  Memory Limit: 128 MBSubmit: 667  Solved: 350 Description 有一 ...

  5. P3160 [CQOI2012]局部极小值

    题目 P3160 [CQOI2012]局部极小值 一眼就是状压,接下来就不知道了\(qwq\) 做法 我们能手玩出局部小值最多差不多是\(8,9\)个的样子,\(dp_{i,j}\)为填满\(1~i\ ...

  6. P3160 [CQOI2012]局部极小值 题解(状压DP+容斥)

    题目链接 P3160 [CQOI2012]局部极小值 双倍经验,双倍快乐 解题思路 存下来每个坑(极小值点)的位置,以这个序号进行状态压缩. 显然,\(4*7\)的数据范围让极小值点在8个以内(以下示 ...

  7. BZOJ2669 [cqoi2012]局部极小值 状压DP 容斥原理

    欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目传送门 - BZOJ2669 题意概括 有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次.如果一个格子比所 ...

  8. 【bzoj2669】[cqoi2012]局部极小值 容斥原理+状压dp

    题目描述 有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次.如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点)都小,我们说这个格子是局部极小值. 给出所有局部极小值的位置,你的任 ...

  9. BZOJ 2669 CQOI2012 局部极小值 状压dp+容斥原理

    题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=2669 题意概述:实际上原题意很简洁了我就不写了吧.... 二话不说先观察一下性质,首先棋盘 ...

  10. [CQOI2012]局部极小值

    题目链接 注意到\(4\times 7\)的矩阵的局部极小值最多只有8个,可以状压. 设\(f[i][sta]\)表示从小到大填数,当前填到\(i\),极小值的填充状态为\(sta\)的方案数. 考虑 ...

随机推荐

  1. Startup在不同环境中的处理

    ASP.NET Core引进了在多种环境中对控制应用程序行为的进一步支持,例如开发环境(Development Environment).预发布环境(Staging Environment),和生产环 ...

  2. <a>标签的特殊和文本的样式

    a是特殊的,要改变a里面的颜色,必须直接给a设置,给a的父级设置不行 属性继承:明明是父级上的的设置样式,结果后代标签也跟着发生变化,这就叫做属性继承. Html 标记语言, 不是编程语言.说白了就是 ...

  3. PATH环境变量

    PATH是环境变量,要大写 那几个目录是你放置linux命令的目录,输入命令后系统会去PATH中寻找是否存在该命令 查看当前环境变量: echo $PATH 也可以用set命令看一下 设置: expo ...

  4. Spring Cloud Feign 使用方法与性能优化

    1. feign自定义Configuration和root 容器有效隔离. 用@Configuration注解 不能在主@ComponentScan (or @SpringBootApplicatio ...

  5. cloudera manager 安装配置

    前面cloudera manager 环境准备和安装我参考的是: https://blog.csdn.net/m0_38017084/article/details/82218559 这篇博客,写的非 ...

  6. NPOI 读取Excel文件

    private void buttonExcel_Click(object sender, EventArgs e) { FileStream fs = null; List<ISheet> ...

  7. react-router-dom v^4路由、带参路由的配置

    首先安装路由 npm install --save react-router-dom 新建一个router.js文件 然后我们的router.js代码如下↓ import React from 're ...

  8. 通过git上传本地代码到github仓库

    最近呢,武汉天气燥热,在公司没啥事,就自己写了一下小demo. 作为一个菜鸟,只在github上扒过别人的代码,还没自己上传过,就试了一下,遇到了一些坑,记录一下. 前提是电脑上安装了git,没有安装 ...

  9. Thrift的C++服务端(线程池和非阻塞)模式

    非阻塞模式 #include "RpcServiceHandler.h" #include <thrift/concurrency/ThreadManager.h> # ...

  10. JavaScript中的原型链和继承

    理解原型链 在 JavaScript 的世界中,函数是一等公民. 上面这句话在很多地方都看到过.用我自己的话来理解就是:函数既当爹又当妈."当爹"是因为我们用函数去处理各种&quo ...