【BZOJ-2732】集合选数 状压DP (思路题)
2734: [HNOI2012]集合选数
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 1070 Solved: 623
[Submit][Status][Discuss]
Description
《集合论与图论》这门课程有一道作业题,要求同学们求出{1, 2, 3, 4, 5}的所有满足以 下条件的子集:若 x 在该子集中,则 2x 和 3x 不能在该子集中。同学们不喜欢这种具有枚举性 质的题目,于是把它变成了以下问题:对于任意一个正整数 n≤100000,如何求出{1, 2,..., n} 的满足上述约束条件的子集的个数(只需输出对 1,000,000,001 取模的结果),现在这个问题就 交给你了。
Input
只有一行,其中有一个正整数 n,30%的数据满足 n≤20。
Output
仅包含一个正整数,表示{1, 2,..., n}有多少个满足上述约束条件 的子集。
Sample Input
Sample Output
【样例解释】
有8 个集合满足要求,分别是空集,{1},{1,4},{2},{2,3},{3},{3,4},{4}。
HINT
Source
Solution
一开始真没想到..
就是写一个矩阵
$$\begin{bmatrix}x & 3x& 9x& ...&\\ 2x & 6x& 18x& ...&\\ 4x& 12x& 36x& ...&\\ ...& ...& ...& ...& \end{bmatrix}$$
然后我们发现实际上就是需要求矩阵中,不允许出现相邻的两个数。
因为这个矩阵的列数是$log_{3}$,行数是$log_{2}$的,显然可以状压..
但是有些数并不会出现在这个矩阵中,所以构造完一个后,找下一个没出现过的x重新构造,乘法原理结合起来即可。
Code
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define LL long long
inline int read()
{
int x=,f=; char ch=getchar();
while (ch<'' || ch>'') {if (ch=='-') f=-; ch=getchar();}
while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}
return x*f;
}
#define MAXN 100010
#define P 1000000001
int a[][],N,M,end[];
LL f[][],bin[],ANS=;
bool visit[MAXN];
inline void Pre(int x)
{
memset(a,,sizeof(a));
memset(end,,sizeof(end)); end[]=;
a[][]=x;
for (int i=; ; i++)
if ((a[i-][]<<)<=N) a[i][]=a[i-][]<<; else {M=i-; break;}
for (int i=; i<=M; i++)
for (int j=; ; j++)
if (a[i][j-]*<=N) a[i][j]=a[i][j-]*; else {end[i]=j-; break;}
// puts("======================================");
// for (int i=1; i<=M; i++,puts(""))
// for (int j=1; a[i][j]; j++)
// printf("%d ",a[i][j]);
// puts("======================================");
for (int i=; i<=M; i++)
for (int j=; j<=end[i]; j++)
visit[a[i][j]]=;
// for (int i=1; i<=M; i++) printf("%d ",end[i]); puts("");
}
inline LL DP(int x)
{
Pre(x); memset(f,,sizeof(f)); f[][]=;
for (int i=; i<=M; i++)
for (int j=; j<bin[end[i]]; j++)
if (f[i][j])
for (int k=; k<bin[end[i+]]; k++)
(f[i+][k]=(!(j&k) && !(k&(k>>)))? (f[i][j]+f[i+][k]) : f[i+][k])%=P;
return f[M+][];
}
int main()
{
N=read();
bin[]=; for (int i=; i<=; i++) bin[i]=bin[i-]<<;
for (int i=; i<=N; i++) if (!visit[i]) (ANS*=DP(i))%=P;
printf("%lld\n",ANS);
return ;
}
断断续续写了2节课....一开始边界搞得有问题特别不科学.....这是弱智,
【BZOJ-2732】集合选数 状压DP (思路题)的更多相关文章
- BZOJ 2734 [HNOI2012]集合选数 (状压DP、时间复杂度分析)
题目链接 https://www.lydsy.com/JudgeOnline/problem.php?id=2734 题解 嗯早就想写的题,昨天因为某些不可告人的原因(大雾)把这题写了,今天再来写题解 ...
- bzoj 2734: [HNOI2012]集合选数 状压DP
2734: [HNOI2012]集合选数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 560 Solved: 321[Submit][Status ...
- [HNOI2012]集合选数 --- 状压DP
[HNOI2012]集合选数 题目描述 <集合论与图论>这门课程有一道作业题,要求同学们求出\({1,2,3,4,5}\)的所有满足以 下条件的子集:若 x 在该子集中,则 2x 和 3x ...
- 【BZOJ-2734】集合选数 状压DP (思路题)
2734: [HNOI2012]集合选数 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1070 Solved: 623[Submit][Statu ...
- 洛谷$P3226\ [HNOI2012]$集合选数 状压$dp$
正解:$dp$ 解题报告: 传送门$QwQ$ 考虑列一个横坐标为比值为2的等比数列,纵坐标为比值为3的等比数列的表格.发现每个数要选就等价于它的上下左右不能选. 于是就是个状压$dp$板子了$QwQ$ ...
- $HNOI2012\ $ 集合选数 状压$dp$
\(Des\) 求对于正整数\(n\leq 1e5\),{\(1,2,3,...,n\)}的满足约束条件:"若\(x\)在该子集中,则\(2x\)和\(3x\)不在该子集中."的子 ...
- bzoj 2734 [HNOI2012]集合选数 状压DP+预处理
这道题很神啊…… 神爆了…… 思路大家应该看别的博客已经知道了,但大部分用的插头DP.我加了预处理,没用插头DP,一行一行来,速度还挺快. #include <cstdio> #inclu ...
- BZOJ.4145.[AMPPZ2014]The Prices(状压DP)
BZOJ 比较裸的状压DP. 刚开始写麻烦惹... \(f[i][s]\)表示考虑了前\(i\)家商店,所买物品状态为\(s\)的最小花费. 可以写求一遍一定去\(i\)商店的\(f[i]\)(\(f ...
- bzoj 2734 集合选数
Written with StackEdit. Description <集合论与图论>这门课程有一道作业题,要求同学们求出\(\{1, 2, 3, 4, 5\}\)的所有满足以 下条件的 ...
随机推荐
- 《java JDK7 学习笔记》之对象封装
1.构造函数实现对象初始化流程的封装.方法封装了操作对象的流程.java中还可以使用private封装对象私有数据成员.封装的目的主要就是隐藏对象细节,将对象当做黑箱子进行操作. 2.在java命名规 ...
- LVM 类型的 Storage Pool - 每天5分钟玩转 OpenStack(8)
LVM 类型的 Storage Pool 不仅一个文件可以分配给客户机作为虚拟磁盘,宿主机上 VG 中的 LV 也可以作为虚拟磁盘分配给虚拟机使用. 不过,LV 由于没有磁盘的 MBR 引导记录,不能 ...
- 国内经典BI系统架构分析
谈起商业智能BI,也许大家并不陌生,但你是否了解国内的各类BI系统架构? 自国内商业智能发展以来,就系统结构方面已经历了多次优化性的变革.目前国内商业智能BI系统的经典架构的模式包括数据层.业务层和应 ...
- Web报表工具FineReport的JS API开发(二)
上次介绍FineReport的JS API中的第一类开发--FR,这次就来介绍一下FS和contentWindow类的开发. 1 FS FS是数据决策系统中的js接口,比如说FS.tabPane.ad ...
- 理解 OpenStack 高可用(HA)(3):Neutron 分布式虚拟路由(Neutron Distributed Virtual Routing)
本系列会分析OpenStack 的高可用性(HA)概念和解决方案: (1)OpenStack 高可用方案概述 (2)Neutron L3 Agent HA - VRRP (虚拟路由冗余协议) (3)N ...
- 用CSS开启硬件加速来提高网站性能
国外一篇文章,有点意思,转载过来,准备尝试下~ 中文地址:http://www.cnblogs.com/rubylouvre/p/3471490.html 原文地址:http://blog.teamt ...
- JSONResult 封装
import java.util.List; import com.fasterxml.jackson.databind.JsonNode; import com.fasterxml.jackson. ...
- String类的equals是如何进行字符串比较的
先看一段代码: public static void main(String[] args) { String a = new String("abc"); String b = ...
- Hibernate 分组查询 子查询 原生SQL
分组查询: 使用group by关键字对数据分组,使用having关键字对分组数据设定约束条件,从而完成对数据分组和统计 1.1 聚合函数:常被用来实现数据统计功能 ① count() 统计记录条数 ...
- 第1章 重构,第一个案例(3):运用多态取代switch
3. 运用多态取代与价格相关的条件逻辑 3.1 switch和“常客积分”代码的再次搬迁 (1)switch:最好不要在另一个对象的属性上运用switch语句 switch(getMovie().ge ...