【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. JQuery官方学习资料(译):避免与其他库的冲突

    避免与其他库的冲突     JQuery库和它所有的插件都是包含在jQuery命名空间下的,作为一般规则,全局对象被存储在jQuery的命名空间内,所以你不会得到jQuery和任何其他库之间的冲突(例 ...

  2. 从零开始学安全(二十六)●利用Nmap目标的本版进行探测

    通过对对方电脑的服务探测 对本版较低的服务 或者无补丁的服务 可以直入侵 版本探测 version  后边就是版本

  3. 禅道导入bugfree 3.0的数据

    禅道项目导入bugfree功能只支持到2.0, 官方不提供3.0的导入,只好自己写了一个.因为bugfree 3.0换人开发了,表结构和禅道差别很大,所以,这个工具不是完全转换,一些History表内 ...

  4. 买or不买?如何测试博彩公司赔率是否合理?

    世界杯期间,烧烤店.酒吧都热闹起来了,柔柔我的朋友圈也热闹起来了,有酱紫的: 还有酱紫的: 然后还有酱紫的: 酱紫的: 当然天台也是一如既然的热闹: 似乎人人都在输钱,那真正的赢家在哪里呢?博彩业的真 ...

  5. SpringBoot数据库集成-Mybatis

    一.java web开发环境搭建 网上有很多教程,参考教程:http://www.cnblogs.com/Leo_wl/p/4752875.html 二.Spring boot搭建 1.Intelli ...

  6. 【Java】Properties文件的解析

    public abstract class ReadProperties { public ReadProperties() {} /** * 回调函数,由调用者处理 * @param key * @ ...

  7. MySQL 查看执行计划

    MySQL 使用 explain + sql 语句查看 执行计划,该执行计划不一定完全正确但是可以参考. EXPLAIN SELECT * FROM user WHERE nid = 3; selec ...

  8. 洛谷P4578 [FJOI2018]所罗门王的宝藏(dfs)

    题意 题目链接 Sol 对于每个询问\(x, y, c\) 从在\((x, y)\)之间连一条边权为\(c\)的双向边,然后就是解\(K\)个二元方程. 随便带个数进去找找环就行了 #include& ...

  9. 详解LSTM

    https://blog.csdn.net/class_brick/article/details/79311148 今天的内容有: LSTM 思路 LSTM 的前向计算 LSTM 的反向传播 关于调 ...

  10. 腾讯.NET&PHP面试题

    在整个面试过程中,作为面试者的你,角色就是小怪兽,面试官的角色则是奥特曼,更不幸的是,作为小怪兽的你是孤身一人,而奥特曼却往往有好几个助攻,你总是被虐得不要不要的~ 作为复读一年才考上专科的我,遗憾的 ...