【NOIP2016提高A组模拟8.15】Garden
题目
分析
其实原题就是【cqoi2012】【bzoj2669】局部极小值。
有一个n行m列的整数矩阵,其中1到nm之间的每个整数恰好出现一次。如果一个格子比所有相邻格子(相邻是指有公共边或公共顶点)都小,我们说这个格子是局部极小值。
给出所有局部极小值的位置,你的任务是判断有多少个可能的矩阵。
发现,X的位置最多有8个,那我们考虑状压dp。
我们从小到大把数填进去,用\(f_{i,j}\)表示,把第i个数填进去后,每个X是否被填了数,用二进制数j表示。
预处理出\(rest_j\)表示填充状态为j时共有多少位置是可以填充的(包括已填充的局部极小值位置)
转移:
\]
但是有些不是为X的位置有可能也是局部极小值,那么我们用容斥,每次把一下有可能出现局部极小值的地方改为X,当额外增加的X的个数为奇数,ans就减去dp得所得的答案,否则ans加上dp得所得的答案。
其中dp方面这篇论文(第5页到第8页)讲的非常清楚
#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
const int mo=12345678;
const int N=10;
int z[9][2]=
{
{-1,-1},
{-1,0},
{-1,1},
{0,-1},
{0,1},
{1,-1},
{1,0},
{1,1},
{0,0}
};
int a[N][N],T,n,m,ans,f[N*N][2000],mi[10],sign[N][2],tot,rest[2000];
char c[N][N];
bool bz[N][N];
int val()
{
tot=0;
int state=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
if(c[i][j]=='X')
{
state+=mi[tot];
sign[++tot][0]=i;
sign[tot][1]=j;
}
}
for(int i=0;i<=state;i++)
{
memset(bz,true,sizeof(bz));
rest[i]=0;
for(int j=1;j<=tot;j++)
if((mi[j-1]&i)==0)
{
for(int k=0;k<=8;k++)
{
bz[sign[j][0]+z[k][0]][sign[j][1]+z[k][1]]=false;
}
}
for(int j=1;j<=n;j++)
for(int k=1;k<=m;k++)
{
if(bz[j][k])
rest[i]++;
}
}
f[0][0]=1;
for(int i=1;i<=n*m;i++)
for(int j=0;j<=state;j++)
{
f[i][j]=0;
(f[i][j]+=f[i-1][j]*(rest[j]-i+1)%mo)%=mo;
for(int k=1;k<=tot;k++)
{
if(mi[k-1]&j)
{
(f[i][j]+=f[i-1][j-mi[k-1]])%=mo;
}
}
}
return f[n*m][mi[tot]-1];
}
int dg(int x,int y,int z1)
{
if(x>n)
{
(ans+=(val()*(z1%2?1:-1))%mo)%mo;
return 0;
}
int xx=x,yy=y+1;
if(yy>m)
{
yy=1;
xx++;
}
dg(xx,yy,z1);
bool q=true;
for(int i=0;i<=8;i++)
{
if(c[x+z[i][0]][y+z[i][1]]=='X')
{
q=false;
break;
}
}
if(q)
{
c[x][y]='X';
dg(xx,yy,z1+1);
c[x][y]='.';
}
}
int main()
{
mi[0]=1;
for(int i=1;i<=9;i++)
mi[i]=mi[i-1]*2;
scanf("%d",&T);
while(T--)
{
tot=0;
ans=0;
memset(c,0,sizeof(c));
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
c[i][j]=getchar();
while(c[i][j]!='X' && c[i][j]!='.')
c[i][j]=getchar();
if(c[i][j]=='.')
tot++;
}
}
if(tot==n*m)
{
printf("0\n");
continue;
}
for(int i=1;i<=n && ans!=-1;i++)
{
for(int j=1;j<=m;j++)
{
bool q=true;
if(c[i][j]=='X')
for(int k=0;k<=7;k++)
{
if(c[i+z[k][0]][j+z[k][1]]=='X')
{
q=false;
ans=-1;
break;
}
}
if(!q)
{
ans=-1;
break;
}
}
}
if(ans==-1)
{
printf("0\n");
}
if(ans!=-1)
{
dg(1,1,1);
printf("%d\n",(ans%mo+mo)%mo);
}
}
}
【NOIP2016提高A组模拟8.15】Garden的更多相关文章
- NOIP2016提高A组模拟10.15总结
第一题,就是将原有的式子一步步简化,不过有点麻烦,搞了很久. 第二题,枚举上下边界,维护一个单调队列,二分. 比赛上没有想到,只打了个暴力,坑了80分. 第三题,贪心,最后的十多分钟才想到,没有打出来 ...
- 【NOIP2016提高A组模拟10.15】打膈膜
题目 分析 贪心, 先将怪物按生命值从小到大排序(显然按这个顺序打是最优的) 枚举可以发对少次群体攻击, 首先将所有的群体攻击发出去, 然后一个一个怪物打,当当前怪物生命值大于2,如果还有魔法值就放重 ...
- 【NOIP2016提高A组模拟10.15】最大化
题目 分析 枚举两个纵坐标i.j,接着表示枚举区域的上下边界, 设对于每个横坐标区域的前缀和和为\(s_l\),枚举k, 显然当\(s_k>s_l\)时,以(i,k)为左上角,(j,k)为右下角 ...
- 【NOIP2016提高A组模拟10.15】算循环
题目 分析 一步步删掉循环, 首先,原式是\[\sum_{i=1}^n\sum_{j=1}^m\sum_{k=i}^n\sum_{l=j}^m\sum_{p=i}^k\sum_{q=j}^l1\] 删 ...
- 【NOIP2016提高A组模拟9.15】Map
题目 分析 发现,当原图是一棵树的时候,那么新建一条边后,就会变成环套树, 而环内的所有点对都是安全点对,如果环中有k个点,答案就是\(k(k-1)\) 联想到,当把原图做一遍tarjan缩点,每个环 ...
- 【NOIP2016提高A组模拟9.15】Osu
题目 分析 考虑二分答案, 二分小数显然是不可取的,那么我们将所有可能的答案求出来,记录在一个数组上,排个序(C++调用函数很容易超时,手打快排,时间复杂度约为\(O(>8*10^7)\),但相 ...
- 【NOIP2016提高A组模拟9.15】Math
题目 分析 因为\((-1)^2=1\), 所以我们只用看\(\sum_{j=1}^md(i·j)\)的值模2的值就可以了. 易证,一个数x,只有当x是完全平方数时,d(x)才为奇数,否则为偶数. 那 ...
- 【NOIP2016提高A组模拟8.15】Throw
题目 分析 首先对于一个状态(a,b,c),假定a<=b<=c: 现在考虑一下这个状态,的转移方案: \[1,中间向两边跳(a,b,c)-->(a*2-b,a,c).(a,b,c)- ...
- 【NOIP2016提高A组模拟8.15】Password
题目 分析 首先我们知道,原A序列其实表示一个矩阵,而这个矩阵的对角线上的数字就是答案B序列. 接着\(a.b>=gcd(a,b)\),所以序列A中的最大的数就是ans[1],第二大的数就是an ...
随机推荐
- Day02:对象和类(下)
构造方法 什么是构造方法 在类实例化对象时运行的一种特殊的方法 为什么需要构造方法 构造方法适合为对象的属性赋初值 编写构造方法 public 类名(){ //方法体 } 构造方法不写返回值 构造方法 ...
- C#学习笔记一(概念,对象与类型,继承)
一.基础 1.CLR为公共语言运行库,类似于JVM 2..NET Framwork是一个独立发布的程序包,其包含了CLR,类库及相关的语言编辑器等工具,类似于JDK,除了C#,还有其他几种语言在CLR ...
- xpath 的用法
爬虫利器xpath xpath术语 节点(node)
- angular [NgClass] [NgStyle],NgIf,[ngSwitch][ngSwitchCase]
[NgClass] CSS 类会根据表达式求值结果进行更新,更新逻辑取决于结果的类型: string - 会把列在字符串中的 CSS 类(空格分隔)添加进来, Array - 会把数组中的各个元素作 ...
- input type=file accept中文件格式限制
原文链接:https://blog.csdn.net/usuallyuser/article/details/83060341 accept="application/msexcel,app ...
- 自带的simple认证
参考: hdfs权限: 官网http://hadoop.apache.org/docs/r1.0.4/cn/hdfs_permissions_guide.html hdfs权限: http://dyi ...
- 【转帖】国产x86处理器KX-6000发布
国产最先进x86处理器KX-6000发布:8核3.0GHz 力压酷睿i5处理器 https://www.cnbeta.com/articles/tech/858981.htm 全网所有的网页都写错了 ...
- 【中间件】Kafka 学习 01
KafKa 博客教程-1 博客教程-2 kafka简介 kafka起源 Kafka是由LinkedIn开发并开源的分布式消息系统,2012年捐赠给Apache基金会,采用Scala语言,运行在JVM中 ...
- Java最新学习线路(基础,源码,项目,实战)
如需获取以下学习资源请关注公众号:Java编程指南 我们为自学者编程的或初学java的小伙伴们准备了一整套完整的学习资源和文章,还有我自己在自学路上的一些总结和学习线路,希望能帮到小伙伴们,如果有什么 ...
- Java第二周总结报告
第二周的学习,开始正式实践进行Java的学习. 本周做了什么? 了解的Java的一些基本知识,如Java变量,数据类型和运算符等.Java变量对不同的数据类型最好采用不同的命名规则,合理的命名有利于提 ...