轮廓线DP POJ3254 && BZOJ 1087
补了一发轮廓线DP,发现完全没有必要从右往左设置状态,自然一点:
5 6 7 8 9
1 2 3 4
如此设置轮廓线标号,转移的时候直接把当前j位改成0或者1就行了。注意多记录些信息对简化代码是很有帮助的,尤其对于我这种代码经常错的一塌糊涂的人来说。。
呆马:
POJ 3254 Corn Fields
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#define inf 1000000007
#define oo 100000000
#define maxn 15
#define maxm 5200 using namespace std; int n,m,mb;
int a[maxn][maxn];
int f[maxn][maxn][];
int bit[][maxn]; bool legal(int x, int y,int mask)
{
if (y==) if (bit[mask][y]&bit[mask][y+]) return ;
for (int i=;i<y-;i++) if (bit[mask][i]&bit[mask][i+]) return ;
for (int i=y;i<m;i++) if (bit[mask][i]&bit[mask][i+]) return ;
return ;
} int main()
{
scanf("%d%d",&n,&m);
memset(a,,sizeof(a));
for (int i=;i<=n;i++)
for (int j=;j<=m;j++)
scanf("%d",&a[i][j]);
mb=(<<m)-;
for (int mask=;mask<=mb;mask++)
{
int tmp=mask;
for (int i=;i<=m;i++)
{
bit[mask][i]=tmp%;
tmp/=;
}
}
memset(f,,sizeof(f));
f[][m][]=;
for (int i=;i<=n;i++)
for (int j=;j<=m;j++)
for (int mask=;mask<=mb;mask++)
if (legal(i,j,mask))
{
int tmp=<<(j-);
if (j==)
{
f[i][j][mask&(~tmp)]+=f[i-][m][mask];
f[i][j][mask&(~tmp)]%=oo;
if (!bit[mask][j] && a[i][j])
{
f[i][j][mask|tmp]+=f[i-][m][mask];
f[i][j][mask|tmp]%=oo;
}
}
else
{
f[i][j][mask&(~tmp)]+=f[i][j-][mask];
f[i][j][mask&(~tmp)]%=oo;
if (!bit[mask][j] && a[i][j])
{
f[i][j][mask|tmp]+=f[i][j-][mask];
f[i][j][mask|tmp]%=oo;
}
}
}
int ans=;
for (int mask=;mask<=mb;mask++)
if (legal(n+,,mask))
{
ans+=f[n][m][mask];
ans%=oo;
}
printf("%d\n",ans);
return ;
}
Corn Fields
BZOJ 1087 Kings
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <vector>
#define inf 1000000007
#define oo 100000000
#define maxn 11
#define maxm 1050 using namespace std; int n,m,mb;
long long f[maxn][maxn][][maxm];
int bit[maxm][maxn]; bool legal(int x, int y,int mask)
{
if (y==)
{
if (bit[mask][n]&bit[mask][n+]) return ;
if (bit[mask][n-]&bit[mask][n+]) return ;
}
else
{
if (bit[mask][y-]&bit[mask][n+]) return ;
if (bit[mask][y-]&bit[mask][n+]) return ;
}
for (int i=;i<n;i++) if (bit[mask][i]&bit[mask][i+]) return ;
return ;
} int main()
{
scanf("%d%d",&n,&m);
mb=(<<(n+))-;
memset(bit,,sizeof(bit));
for (int mask=;mask<=mb;mask++)
{
int tmp=mask;
for (int i=;i<=n+;i++)
{
bit[mask][i]=tmp%;
tmp/=;
}
}
memset(f,,sizeof(f));
f[][n][][]=;
for (int i=;i<=n;i++)
for (int j=;j<=n;j++)
for (int k=;k<=m;k++)
for (int mask=;mask<=mb;mask++)
if (legal(i,j,mask))
{
int tmp=<<(j-);
int ost=mask&(~tmp),nst=mask|tmp;
if (((mask>>(j-))&)!=((mask>>n)&)) ost^=(<<n),nst^=(<<n);
if (j==)
{
f[i][j][k][ost]+=f[i-][n][k][mask];
if (k && !bit[mask][j] && !bit[mask][j+] && !bit[mask][j-])
f[i][j][k][nst]+=f[i-][n][k-][mask];
}
else
{
f[i][j][k][ost]+=f[i][j-][k][mask];
if (k && !bit[mask][j] && !bit[mask][j+] && !bit[mask][j-] && !bit[mask][n+])
f[i][j][k][nst]+=f[i][j-][k-][mask];
}
}
long long ans=;
for (int mask=;mask<=mb;mask++)
if (legal(n+,,mask))
{
ans+=f[n][n][m][mask];
}
printf("%lld\n",ans);
return ;
}
Kings
轮廓线DP POJ3254 && BZOJ 1087的更多相关文章
- BZOJ.4572.[SCOI2016]围棋(轮廓线DP)
BZOJ 洛谷 \(Description\) 给定\(n,m,c\).\(Q\)次询问,每次询问给定\(2*c\)的模板串,求它在多少个\(n*m\)的棋盘中出现过.棋盘的每个格子有三种状态. \( ...
- HDU4804 Campus Design 轮廓线dp
跟上面那篇轮廓线dp是一样的,但是多了两个条件,一个是在原图上可能有些点是不能放的(即障碍),所以转移的时候要多一个判断color[i][j]是不是等于1什么的,另外一个是我们可以有多的1*1的骨牌, ...
- POJ2411 Mondriaan's Dream 轮廓线dp
第一道轮廓线dp,因为不会轮廓线dp我们在南京区域赛的时候没有拿到银,可见知识点的欠缺是我薄弱的环节. 题目就是要你用1*2的多米诺骨排填充一个大小n*m(n,m<=11)的棋盘,问填满它有多少 ...
- UVA - 11270 轮廓线DP
其实这题还能用状压DP解决,可是时间达到2000ms只能过掉POJ2411.状压DP解法详见状压DP解POJ2411 贴上POJ2411AC代码 : 2000ms 时间复杂度h*w*(2^w)*(2^ ...
- [UOJ422][集训队作业2018]小Z的礼物——轮廓线DP+min-max容斥
题目链接: [集训队作业2018]小Z的礼物 题目要求的就是最后一个喜欢的物品的期望得到时间. 根据$min-max$容斥可以知道$E(max(S))=\sum\limits_{T\subseteq ...
- 【UOJ#422】【集训队作业2018】小Z的礼物(min-max容斥,轮廓线dp)
[UOJ#422][集训队作业2018]小Z的礼物(min-max容斥,轮廓线dp) 题面 UOJ 题解 毒瘤xzy,怎么能搬这种题当做WC模拟题QwQ 一开始开错题了,根本就不会做. 后来发现是每次 ...
- 4572: [Scoi2016]围棋 轮廓线DP KMP
国际惯例的题面:这种题目显然DP了,看到M这么小显然要状压.然后就是具体怎么DP的问题.首先我们可以暴力状压上一行状态,然后逐行转移.复杂度n*3^m+3^(m*2),显然过不去. 考虑状态的特殊性, ...
- Solid Dominoes Tilings (轮廓线dp打表 + 容器)
第一步先打一个表,就是利用轮廓线DP去打一个没有管有没有分界线组合数量的表 #include<bits/stdc++.h> using namespace std; ; <<; ...
- HDU 4802 && HDU 4803 贪心,高精 && HDU 4804 轮廓线dp && HDU 4805 计算几何 && HDU 4811 (13南京区域赛现场赛 题目重演A,B,C,D,J)
A.GPA(HDU4802): 给你一些字符串对应的权重,求加权平均,如果是N,P不计入统计 GPA Time Limit: 2000/1000 MS (Java/Others) Memory ...
随机推荐
- EMF学习,为了实现可扩展可自定义的模型验证 - 各种实现方法学习
自: http://blog.csdn.net/javaman_chen/article/details/6057033 http://www.ibm.com/developerworks/cn/op ...
- JMeter--集合点设置(转)
集合点:简单来理解一下,虽然我们的"性能测试"理解为"多用户并发测试",但真正的并发是不存在的,为了更真实的实现并发这感念,我们可以在需要压力的地方设置集合点, ...
- 利用jemalloc优化mysql
一.下载安装jemalloc #wget http://www.canonware.com/download/jemalloc/jemalloc-3.6.0.tar.bz2 #tar jxvf jem ...
- 给Source Insight做个外挂系列之一--发现Source Insight
一提到外挂程序,大家肯定都不陌生,QQ就有很多个版本的去广告外挂,很多游戏也有用于扩展功能或者作弊的工具,其中很多也是以外挂的形式提供的.外挂和插件的区别在于插件通常依赖于程序的支持,如果程序不支持插 ...
- Python之路----------基础 一(数据类型、变量、基本语法、流程控制)
一. 数据类型与变量 1.数据类型 整数 #Python在程序中的表示方法和数学上的写法一模一样,-1,0,1都是整数. 浮点数 1 #浮点数就是小数. 字符串 1 #在Python中字符串是以 ...
- MySQL 联合索引详解
MySQL 联合索引详解 联合索引又叫复合索引.对于复合索引:Mysql从左到右的使用索引中的字段,一个查询可以只使用索引中的一部份,但只能是最左侧部分.例如索引是key index (a,b,c ...
- yii框架的安装
在安装YII2.0首先我们需要下载一个composer,何为composer,简单解释为PHP的组件工具, Composer是PHP中用来管理依赖(dependency)关系的工具,你可以在自己的项目 ...
- 解决Office互操作错误"检索COML类工厂中 CLSID为 {xxx}的组件时失败,原因是出现以下错误: 80070005"
Excel为例(其他如Word也适用)文件数据导入时报出以下错误: 检索COML类工厂中 CLSID为 {00024500-0000-0000-C000-000000000046}的组件时失败,原因是 ...
- DNS CNAME的一些细节
1, 概述 DNS中的CNAME可以减轻运维压力,使得已有的DNS配置具有一定的灵活性和可扩展性.本文对CNAME中的一些细节做阐述, 使DNS服务器的运维人员和开发人员能合理地使用CNAME. 2, ...
- Spark&Hadoop:scala编写spark任务jar包,运行无法识别main函数,怎么办?
昨晚和同事一起看一个scala写的程序,程序都写完了,且在idea上debug运行是ok的.但我们不能调试的方式部署在客户机器上,于是打包吧.打包时,我们是采用把外部引入的五个包(spark-asse ...