【noip模拟】局部最小值
TimeLimit: 1000ms MemoryLimit: 256MB
Description
有一个n行m列的整数矩阵,其中1到n×m之间的每个整数恰好出现一次。如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点)都小,我们说这个格子是局部极小值。
给出所有局部极小值的位置,你的任务是判断有多少个可能的矩阵。
Input
输入第一行包含两个整数n和m(1<=n<=4, 1<=m<=7),即行数和列数。以下n行每行m个字符,其中“X”表示局部极小值,“.”表示非局部极小值。
Output
输出仅一行,为可能的矩阵总数除以998244353的余数。
Sample Input
【样例输入1】
2 4
.X..
...X
【样例输入2】
4 2
X.
..
..
X.
【样例输入3】
1 2
XX
Sample Output
【样例输出1】
2100
【样例输出2】
2520
【样例输出3】
0
HINT
对于10%的数据,$n≤3,m≤3$
对于20%的数据,$n≤3,m≤4$
对于45%的数据,$n≤4,m≤4$
对于60%的数据,$n≤4,m≤5$
对于80%的数据,$n≤4,m≤6$
对于100%的数据,$n≤4,m≤7$
[吐槽]
这题的话。。在模拟题里面出现了。。两遍啊哈哈。。
第一次做的时候啃得挺。。辛苦。。不过不得不说这个状压还是很好玩的
[题解]
首先看看这个局部最小值有什么可以利用的地方。。
画一下图就会发现好像最多只有8个
这个数字很小所以就可以。。。
考虑状压
考虑用$f_{i,j}$表示填完了$1$到$i$这$i$个数,局部最小值的填充情况为$j$的方案数
用$cnt_i$表示局部最小值的填充情况为$i$时,还有多少个位置是可以填数的
那么就可以得到:
$f_{i,j} = f_{i-1,j}*(cnt_j - (i -1)) + \sum\limits_{k\in j} f_{i-j,k}$
为啥是$cnt_j - (i - 1) $呢
因为本来有$cnt_j$个位置可以填,但是现在已经有$i - 1$个位置已经确定下来了,并不能填所以要减掉
考虑cnt怎么算
首先我们枚举一下状态$ i (0 <= i <= (1<<tot) -1 )$ (tot为X的个数)
因为现在已经填充的局部最小值有哪些是确定了的,那就意味着状态中非填充的局部最小值是不能填的
然后因为我们接下来填的数是剩下的数中最小的,所以这个数肯定不能填在局部最小值的相邻位置
那么就很显然是枚举一下$i$中没有填充的局部最小值,然后把这些地方减掉就好啦
然后就是$k$怎么枚举
我们可以每次取lowbit,统计完之后就把最后一位去掉(实现起来就是$k$ & $(k - 1) $)就ok了
最后还有一个小小的问题
题目要求只能有这些标为X的地方是局部最小值,其他的地方不能是
而我们的这种奇妙填数方法可能会导致有别的格子不小心变成了局部最小值了,这是不合法的
所以就应该要把这些方案减掉
具体实现其实很简单粗暴,直接枚举有哪些位置可能变成局部最小值
用同样的方法求出方案数,然后减掉就好(其实就是一个dfs巨无敌粗暴)
然后就很愉快滴搞完啦ovo
#include<iostream>
#include<cstdio>
#include<cstring>
#define MOD 998244353
#define ll long long
using namespace std;
const int dx[]={,-,,,,-,-,,};
const int dy[]={,,,,-,-,,-,};
struct xxx
{
int x,y;
}X[];
char a[][];
ll b[][],cnt[<<],f[][<<];
int n,m,tot;
ll ans;
int dfs(int x,int y,int op);
ll calc();
bool ok(int x,int y);
bool check(); int main()
{
scanf("%d%d\n",&n,&m);
for (int i=;i<=n;++i)
{
for (int j=;j<=m;++j)
scanf("%c",&a[i][j]);
scanf("\n");
}
if (!check()) {printf("0\n");return ;}
ans=;
dfs(,,);
printf("%lld\n",ans);
} int dfs(int x,int y,int op)
{
int mark=op?-:;
if (x==n+)
{ans=(ans+MOD+mark*calc())%MOD;return ;}
if (a[x][y]=='.')
{
a[x][y]='X';
if (check())
{
if (y==m) dfs(x+,,op^);
else dfs(x,y+,op^);
}
a[x][y]='.';
}
if (y==m) dfs(x+,,op);
else dfs(x,y+,op);
} ll calc()
{
tot=;
for (int i=;i<=n;++i)
for (int j=;j<=m;++j)
if (a[i][j]=='X')
X[++tot].x=i,X[tot].y=j;
int x,y;
memset(cnt,,sizeof(cnt));
for (int i=;i<<<tot;++i)
{
for (int j=;j<=n;++j)
for (int k=;k<=m;++k)
b[j][k]=;
for (int j=;j<tot;++j)
{
if (<<j&i) continue;
x=X[j+].x,y=X[j+].y;
b[x][y]=;
for (int k=;k<=;++k)
if (ok(x+dx[k],y+dy[k]))
b[x+dx[k]][y+dy[k]]=;
}
for (int j=;j<=n;++j)
for (int k=;k<=m;++k)
cnt[i]=(cnt[i]+b[j][k])%MOD;
}
//for (int i=0;i<1<<tot;++i) printf("%d ",cnt[i]);
//printf("\n");
for (int i=;i<=n*m;++i)
for (int j=;j<<<tot;++j)
f[i][j]=;
f[][]=;
for (int i=;i<=n*m;++i)
{
for (int j=;j<<<tot;++j)
{
f[i][j]=((ll)f[i-][j]*((cnt[j]-(i-))%MOD))%MOD;
for (int k=j;k;k&=(k-))
{
x=k&(-k);
f[i][j]=((ll)f[i][j]+f[i-][j-x])%MOD;
}
}
}
return f[n*m][(<<tot)-];
} bool ok(int x,int y)
{
if (x<||y<||x>n||y>m) return false;
return true;
} bool check()
{
int x,y;
for (int i=;i<=n;++i)
for (int j=;j<=m;++j)
{
if (a[i][j]!='X') continue;
for (int k=;k<=;++k)
{
x=i+dx[k];
y=j+dy[k];
if (ok(x,y)&&a[x][y]=='X') return false;
}
}
return true;
}
挫挫滴代码
【noip模拟】局部最小值的更多相关文章
- 关于过拟合、局部最小值、以及Poor Generalization的思考
Poor Generalization 这可能是实际中遇到的最多问题. 比如FC网络为什么效果比CNN差那么多啊,是不是陷入局部最小值啊?是不是过拟合啊?是不是欠拟合啊? 在操场跑步的时候,又从SVM ...
- 2018.9.22 NOIP模拟赛
*注意:这套题目应版权方要求,不得公示题面. 从这里开始 Problem A 妹子 Problem B 旅程 Problem C 老大 因为业务水平下滑太严重,去和高一考NOIP模拟,sad... P ...
- NOIP模拟赛-2018.11.7
NOIP模拟赛 如果用命令行编译程序可以发现没加头文件之类的错误. 如果用命令行编译程序可以发现没加头文件之类的错误. 如果用命令行编译程序可以发现没加头文件之类的错误. 编译之前另存一份,听说如果敲 ...
- NOI.AC NOIP模拟赛 第四场 补记
NOI.AC NOIP模拟赛 第四场 补记 子图 题目大意: 一张\(n(n\le5\times10^5)\)个点,\(m(m\le5\times10^5)\)条边的无向图.删去第\(i\)条边需要\ ...
- NOIP模拟题汇总(加厚版)
\(NOIP\)模拟题汇总(加厚版) T1 string 描述 有一个仅由 '0' 和 '1' 组成的字符串 \(A\),可以对其执行下列两个操作: 删除 \(A\)中的第一个字符: 若 \(A\)中 ...
- 9.9 NOIP模拟题
9.9 NOIP模拟题 T1 两个圆的面积求并 /* 计算圆的面积并 多个圆要用辛普森积分解决 这里只有两个,模拟计算就好 两圆相交时,面积并等于中间两个扇形面积减去两个三角形面积 余弦定理求角度,算 ...
- 20190908 NOIP 模拟40
考试过程: 刚看完题,发现T1是个类lis 问题,但要求$O(nlogn)$,应该是个数据结构优化dp,T2应该是个数据结构,T3是个字符串?没有匹配,不会是后缀数组吧,这是NOIP模拟啊,可能是个d ...
- NOIP模拟17.9.22
NOIP模拟17.9.22 前进![问题描述]数轴的原点上有一只青蛙.青蛙要跳到数轴上≥
- noip模拟30[毛毛毛探探探]
\(noip模拟30\;solutions\) 所以说,这次被初中的大神给爆了????? 其实真的不甘心,这次考场上的遗憾太多,浪费的时间过多,心情非常不好 用这篇题解来结束这场让人伤心的考试吧 \( ...
随机推荐
- 2017年StackOverflow上最好的20个Python问题
1.Python的 .. (点号 点号) 是什么语法? 答案地址:https://stackoverflow.com/questions/43487811/what-is-python-dot-dot ...
- CentOS Crontab(定时任务)
安装crontab: yum install crontabs 说明: service crond start //启动服务 service crond stop //关闭服务 service cro ...
- MySQL之表的数据类型
一 介绍 存储引擎决定了表的类型,而表内存放的数据也要有不同的类型,每种数据类型都有自己的宽度,但宽度是可选的 详细参考: http://www.runoob.com/mysql/mysql-data ...
- CodeForces - 796D Police Stations bfs
思路:删除尽量多的边使得所有点都能在限制距离之内到达一个警局,删除边会形成多棵子树,最多只能k棵.其实就是以每个警局为根结点,把整棵树划分为以警局为根结点的k棵树,说明要删除的边的数量就是k-1条,即 ...
- 【Learning】 莫比乌斯反演
莫比乌斯反演 对于两个定义域为非负整数的函数\(F(n)\)和\(f(n)\) 若满足:\(F(n)=\sum\limits_{d|n}f(d)\),则反演得到\(f(n)=\sum\limi ...
- 笔记︱范数正则化L0、L1、L2-岭回归&Lasso回归(稀疏与特征工程)
机器学习中的范数规则化之(一)L0.L1与L2范数 博客的学习笔记,对一些要点进行摘录.规则化也有其他名称,比如统计学术中比较多的叫做增加惩罚项:还有现在比较多的正则化. -------------- ...
- R语言︱决策树族——随机森林算法
每每以为攀得众山小,可.每每又切实来到起点,大牛们,缓缓脚步来俺笔记葩分享一下吧,please~ --------------------------- 笔者寄语:有一篇<有监督学习选择深度学习 ...
- 【html5】html5本地简单存储
html5本地简单存储 HTML5 提供了四种在客户端存储数据的新方法,即 localStorage .sessionStorage.globalStorage.Web Sql Database. 前 ...
- Java利用递归算法统计1-6的数组排列组合数
Java利用递归算法统计1-6的数组排列组合数 1.设计源码 /** * @Title:ArrayCombination.java * @Package:com.you.data * @Descrip ...
- AM335x(TQ335x)学习笔记——GPIO按键驱动移植
还是按照S5PV210的学习顺序来,我们首先解决按键问题.TQ335x有六个用户按键,分别是上.下.左.右.Enter和ESC.开始我想到的是跟学习S5PV210时一样,编写输入子系统驱动解决按键问题 ...