BZOJ_2734_[HNOI2012]集合选数_构造+状压DP
BZOJ_2734_[HNOI2012]集合选数_构造+状压DP
题意:《集合论与图论》这门课程有一道作业题,要求同学们求出{1, 2, 3, 4, 5}的所有满足以 下条件的子集:若 x 在该子集中,则 2x 和 3x 不能在该子集中。同学们不喜欢这种具有枚举性 质的题目,于是把它变成了以下问题:对于任意一个正整数 n≤100000,如何求出{1, 2,..., n} 的满足上述约束条件的子集的个数(只需输出对 1,000,000,001 取模的结果),现在这个问题就 交给你了。
分析:
我们构造出一个矩阵
$\begin{matrix}
1&2^03^1&2^03^2\\
2^13^0&2^13^1&2^13^2\\
2^23^0&2^23^1&2^23^2\\
\end{matrix}
$
发现矩阵的相邻两个格子的数不能同时取
状压DP一下
要把所有不在矩阵中的数当作1重新构造,比如5,7等等
每个矩阵的结果乘起来就是答案
代码:
#include <stdio.h>
#include <string.h>
#include <algorithm>
using namespace std;
#define LL long long
LL p=1000000001,A,f[18][1<<12];
int vis[100050],s[18],mat[18][18];
LL ans=1;
void build(int x){
int n=1,m=1,now=x;
while(now*3<=A)m++,now*=3;
now=x;
while(now*2<=A)n++,now<<=1;
int mask=(1<<m)-1;
memset(s,0,sizeof(s));
memset(f,0,sizeof(f));
memset(mat,0,sizeof(mat));
mat[1][1]=x;vis[x]=1;
for(int i=2;i<=m;i++){
mat[1][i]=mat[1][i-1]*3;
vis[mat[1][i]]=1;
}
s[1]=mask;
for(int i=2;i<=n;i++){
mat[i][1]=mat[i-1][1]*2;
vis[mat[i][1]]=1;
for(int j=2;j<=m;j++){
mat[i][j]=mat[i-1][j]*2;
if(mat[i][j]>A){
s[i]=mask^((1<<m-j+1)-1);
break;
}
vis[mat[i][j]]=1;
}
if(!s[i])s[i]=mask;
}
f[0][0]=1;
s[0]=mask;
for(int i=0;i<n;i++){
for(int j=0;j<=mask;j++){
if((j|s[i])!=s[i])continue;
if(j&(j<<1))continue;
for(int k=0;k<=mask;k++){
if((k|s[i+1])!=s[i+1])continue;
if(k&(k<<1))continue;
if(j&k)continue;
f[i+1][k]+=f[i][j];
f[i+1][k]%=p;
}
}
}
LL re=0;
for(int i=0;i<=mask;i++)re+=f[n][i],re%=p;
ans=re*ans%p; /*for(int i=1;i<=n;i++){
for(int j=1;j<=m;j++){
printf("%d ",mat[i][j]);
}
puts("");
}*/ /*for(int i=1;i<=n;i++){
printf("%d\n",s[i]);
}*/
}
int main(){
scanf("%lld",&A);
for(int i=1;i<=A;i++){
if(!vis[i])build(i);
}
printf("%lld",ans);
}
BZOJ_2734_[HNOI2012]集合选数_构造+状压DP的更多相关文章
- [HNOI2012]集合选数(构造,状态压缩,DP)
神仙题. 莫名其妙的就试一试把所有数放进一个类似矩阵的东西里面. 首先把 \(1\) 放到左上角,然后在每个数的右边放它的 \(3\) 倍(大于 \(n\) 就不用放了),下面放它的 \(2\) 倍( ...
- [HNOI2012]集合选数 --- 状压DP
[HNOI2012]集合选数 题目描述 <集合论与图论>这门课程有一道作业题,要求同学们求出\({1,2,3,4,5}\)的所有满足以 下条件的子集:若 x 在该子集中,则 2x 和 3x ...
- bzoj 2734: [HNOI2012]集合选数 状压DP
2734: [HNOI2012]集合选数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 560 Solved: 321[Submit][Status ...
- $HNOI2012\ $ 集合选数 状压$dp$
\(Des\) 求对于正整数\(n\leq 1e5\),{\(1,2,3,...,n\)}的满足约束条件:"若\(x\)在该子集中,则\(2x\)和\(3x\)不在该子集中."的子 ...
- 洛谷$P3226\ [HNOI2012]$集合选数 状压$dp$
正解:$dp$ 解题报告: 传送门$QwQ$ 考虑列一个横坐标为比值为2的等比数列,纵坐标为比值为3的等比数列的表格.发现每个数要选就等价于它的上下左右不能选. 于是就是个状压$dp$板子了$QwQ$ ...
- 2734: [HNOI2012]集合选数
2734: [HNOI2012]集合选数 链接 分析: 转化一下题意. 1 3 9 27... 2 6 18 54... 4 12 36 108... 8 24 72 216... ... 写成这样的 ...
- 2734: [HNOI2012]集合选数 - BZOJ
Description <集合论与图论>这门课程有一道作业题,要求同学们求出{1, 2, 3, 4, 5}的所有满足以 下条件的子集:若 x 在该子集中,则 2x 和 3x 不能在该子集中 ...
- [HNOI2012]集合选数(状压DP+构造)
题目要求若出现x,则不能出现2x,3x 所以我们考虑构造一个矩阵 \(1\ 2\ 4 \ 8--\) \(3\ 6\ 12\ 24--\) \(9\ 18\ 36--\) \(--\) 不难发现,对于 ...
- BZOJ 2734: [HNOI2012]集合选数 [DP 状压 转化]
传送门 题意:对于任意一个正整数 n≤100000,如何求出{1, 2,..., n} 的满足若 x 在该子集中,则 2x 和 3x 不能在该子集中的子集的个数(只需输出对 1,000,000,001 ...
随机推荐
- javascript随机一个1-9的数字
window.onload=function(){ var oTxt=document.getElementById('txt'); for(i=1;i<=200;i ...
- 《JUnit实战(第2版)》读书笔记
第1章 JUnit起步 主要了解JUnit,如何安装.运行JUnit 要点 JUnit4不需要像JUnit3那样extends TestCase类 Junit4基本都是用注解(该书都翻译为注释,但我喜 ...
- 修改was数据源
本机的RAD运行的工程可以通过修改jpa中的persistence中的jni修改数据源: 对于通过was控制台部署的ear需要在was控制台:资源--jdbc 修改数据源
- 基础概念:Oracle数据库、实例、用户、表空间、表之间的关系
基础概念:Oracle数据库.实例.用户.表空间.表之间的关系 数据库: Oracle数据库是数据的物理存储.这就包括(数据文件ORA或者DBF.控制文件.联机日志.参数文件).其实Oracle数据库 ...
- 手动编译Flume
1.源码下载: 我用的是1.6版,因为加了kafka-sink,下载地址 http://www.apache.org/dyn/closer.cgi/flume/1.6.0/apache-flume-1 ...
- SpringBoot集成jsp
一.springBoot集成jsp: 1.修改pom文件 <!--集成jsp所需jar包--> <!--jsp页面使用jstl标签--> <dependency> ...
- golang基础- ElasticSearch搜索引擎、kibana可视化工具、向ES输出数据
转载自:https://blog.csdn.net/u013210620/article/details/78647366 安装ElasticSearch ElasticSearch是一个基于Luce ...
- 从前端开发看HTTP协议的应用
一.Chrome Developer Network Tab Cheome Developer作为现在前端开发者最常用的开发调试工具,其具有前端可以涉及到的各方面的强大功能,为我们的开发和定位问题提供 ...
- [Luogu 4135] 作诗
Description 神犇SJY虐完HEOI之后给傻×LYD出了一题: SHY是T国的公主,平时的一大爱好是作诗. 由于时间紧迫,SHY作完诗之后还要虐OI,于是SHY找来一篇长度为N的文章,阅读M ...
- 初探Apache Beam
文章作者:luxianghao 文章来源:http://www.cnblogs.com/luxianghao/p/9010748.html 转载请注明,谢谢合作. 免责声明:文章内容仅代表个人观点, ...