轮廓线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 ...
随机推荐
- Win下常用命令大全
定时关机命令win+R输入 Shutdown -s //后面不加时间参数,默认的是倒计时30秒后关机Shutdown -s -t 0 //立刻关机shutdown -s -t 100 //100s后关 ...
- PHP连接数据库的方法
mysql可通过两种方式通过PHP和web相连,一种通过php的mysql相关函数,另一种通过php的ODBC相关函数. 相关函数如下: MYSQL函数 mysql_affected_rows: 得到 ...
- 多边形碰撞 -- SAT方法
检测凸多边形碰撞的一种简单的方法是SAT(Separating Axis Theorem),即分离轴定理. 原理:将多边形投影到一条向量上,看这两个多边形的投影是否重叠.如果不重叠,则认为这两个多边形 ...
- <meta>指定浏览器模式(browser mode)或文档模式(document mode)无效
这是前两天解决的一个故障,准确的说它不是一个SharePoint的问题,而是IE8浏览器或者说是HTML代码的问题,但我感觉还是挺有意思的,所以贴上来分享一下. 基础知识 简单的讲,就是IE浏览器中有 ...
- Java数据校验(Bean Validation / JSR303)
文档: http://beanvalidation.org/1.1/spec/ API : http://docs.jboss.org/hibernate/beanvalidation/spec/1. ...
- 掌握Thinkphp3.2.0----自动完成
自动完成,是在create数据的时候自动执行的一个或多个动作 自动完成也分为动态和静态两种 格式: //完成规则array( array(完成字段1,完成规则,[完成条件,附加规则]), array( ...
- adv钓鱼题
4 10 4 5 6 2 10 2 10 4 5 6 2 10 2 60 1 20 29 20 50 20 60 1 20 20 20 60 20 #include<iostream> u ...
- PGPool 配置错误定位 s_do_auth: expecting R got E
自从按照教程 http://www.pgpool.net/docs/latest/pgpool-zh_cn.html#hba配置好PGPool以后,每次启动 pgpool -c -n -D 都报 s_ ...
- 自定义组件-支持PNG图片的多态GraphicButton
按钮功能使用TButton也可以解决, 但是TButton是会获得焦点的, 很多时候我们要求按钮不获得焦点, 而Speedbutton又不支持PNG图片 所以按照TSpeedbutton的代码, 重新 ...
- django对数据查询结果进行排序的方法
在你的 Django 应用中,你或许希望根据某字段的值对检索结果排序,比如说,按字母顺序. 那么,使用 order_by() 这个方法就可以搞定了. ? 1 2 >>> Publis ...