二进制的力量

状态压缩DP

愤怒的小鸟

第一次接触状态压缩DP是在NOIP2016的愤怒的小鸟,当时菜得连题目都没看懂,不过现在回过头来看还是挺简单的,那么我们再来看看这道题吧。

题意&数据范围看这

考虑预处理出两个点构成的抛物线,因为经过原点,所以对于二次函数

ax2+bx+c
因此已知两个点 (x1,y1),(x2,w2)可得出

a=(y2/x2+y1/x1)/(x2-x1)

b=(y2-ax22)/x2

所以思路也就自然而然的来了,枚举两个点,求出它们所构成的抛物线,再枚举其它的点,看这条抛物线经过另外的哪些点。

然后问题又来了,该如何存储这条抛物线上的信息?

思路1:储存这条抛物线上的个数

这应该是最好想到的一条思路了,但也很快可以否决掉,会发现只存个数的话会无法判断是否打完。

思路2:对每条抛物线开一个数组

这样做的正确性无法否定,可是你的空间呢?

思路3:二进制信息存储

有的时候暴力的思想也是很重要的,这往往预示着正解。

既然我们只要表示某个点的存在或否,我们何不用2进制的一位来表示呢?

这便是状态压缩的核心思想了,对于简单的状态存储,重新开数组太过浪费,可以用一个二进制数来代替数组。

假设我们现在有5个点,那么一条抛物线的初始状态是什么也没有,用一个五位的二进制数表示即 00000

现在假设我们发现第3个点在这条抛物线上,即要让它变成00100 要怎么办呢。

不难想到我们只要把 1 左移 2 位,再or上去就行了。

因此每当我们发现一个编号为 i 的点,假设表示抛物线的数为f,那么我们可以写出如下运算式:

f=f|2i-1

这样就可以预处理出所有的抛物线了,

然后有了抛物线,状态转移方程也不难写出来。

设dpi表示当前状态最小需要的最小猪,则:
dpi|f=min { dpi+1}
其中f为抛物线。

#include<cstdio>
#include<cstring>
#include<cmath>
#define maxm 1000+10
#define maxn 1000+10
#define min(a,b) (a<b?a:b)
using namespace std;
int f[(1<<20)],get[20][20];
double x[20],y[20];
template<class T>
inline bool read(T&n,char ch=getchar(),int sign=1){
if(ch==EOF)return 0;for(n=0;(ch<'0'||ch>'9')&&ch!='-';ch=getchar());
if(ch=='-')sign=-1,ch=getchar();for(;ch>='0'&&ch<='9';ch=getchar())
(n*=10)+=(ch-'0');n*=sign;return 1;
}
inline bool equ(double a,double b){
return abs(a-b)<10e-8;
}
int main(){
#ifdef Files
freopen("f.in","r",stdin);
freopen("f.out","w",stdout);
#endif
#ifdef die
for(int i=1;i<=19;i++) printf("%d ",un[i]);
#endif
int T;read(T);
for(int i=1;i<=T;i++){
int n,m,tot=0;read(n),read(m);
memset(get,0,sizeof(get));
for(int i=1;i<=n;i++) scanf("%lf%lf",&x[i],&y[i]);
for(int i=1;i<=n;i++)
for(int j=i+1;j<=n;j++){
if(equ(x[i],x[j]))continue;
double a=(y[j]/x[j]-y[i]/x[i])/(x[j]-x[i]);
double b=(y[j]-a*x[j]*x[j])/x[j];
if(a>-1e-8) continue;
for(int k=1;k<=n;k++)
if(equ(a*x[k]*x[k]+b*x[k],y[k])) get[i][j]|=1<<(k-1);
}
memset(f,127,sizeof(f));f[0]=0;
for(int i=0;i<=(1<<n)-1;i++)
for(int j=1;j<=n;j++)
if(!((1<<j-1)&i)){
for(int k=j+1;k<=n;k++)
f[i|get[j][k]]=min(f[i]+1,f[i|get[j][k]]);
f[i|1<<(j-1)]=min(f[i]+1,f[i|1<<(j-1)]);
}
printf("%d\n",f[(1<<n)-1]);
}
}

【dp】状压dp的更多相关文章

  1. 【BZOJ】1076 [SCOI2008]奖励关 期望DP+状压DP

    [题意]n种宝物,k关游戏,每关游戏给出一种宝物,可捡可不捡.每种宝物有一个价值(有负数).每个宝物有前提宝物列表,必须在前面的关卡取得列表宝物才能捡起这个宝物,求期望收益.k<=100,n&l ...

  2. CCF 201312-4 有趣的数 (数位DP, 状压DP, 组合数学+暴力枚举, 推公式, 矩阵快速幂)

    问题描述 我们把一个数称为有趣的,当且仅当: 1. 它的数字只包含0, 1, 2, 3,且这四个数字都出现过至少一次. 2. 所有的0都出现在所有的1之前,而所有的2都出现在所有的3之前. 3. 最高 ...

  3. hdu 4352 "XHXJ's LIS"(数位DP+状压DP+LIS)

    传送门 参考博文: [1]:http://www.voidcn.com/article/p-ehojgauy-ot.html 题解: 将数字num字符串化: 求[L,R]区间最长上升子序列长度为 K ...

  4. [转]状态压缩dp(状压dp)

    状态压缩动态规划(简称状压dp)是另一类非常典型的动态规划,通常使用在NP问题的小规模求解中,虽然是指数级别的复杂度,但速度比搜索快,其思想非常值得借鉴. 为了更好的理解状压dp,首先介绍位运算相关的 ...

  5. 状态压缩dp 状压dp 详解

    说到状压dp,一般和二进制少不了关系(还常和博弈论结合起来考,这个坑我挖了还没填qwq),二进制是个好东西啊,所以二进制的各种运算是前置知识,不了解的话走下面链接进百度百科 https://baike ...

  6. 洛谷 P3343 - [ZJOI2015]地震后的幻想乡(朴素状压 DP/状压 DP+微积分)

    题面传送门 鸽子 tzc 竟然来补题解了,奇迹奇迹( 神仙题 %%%%%%%%%%%% 解法 1: 首先一件很明显的事情是这个最小值可以通过类似 Kruskal 求最小生成树的方法求得.我们将所有边按 ...

  7. 51nod 1673 树有几多愁(链表维护树形DP+状压DP)

    题意 lyk有一棵树,它想给这棵树重标号. 重标号后,这棵树的所有叶子节点的值为它到根的路径上的编号最小的点的编号. 这棵树的烦恼值为所有叶子节点的值的乘积. lyk想让这棵树的烦恼值最大,你只需输出 ...

  8. BZOJ3836 [Poi2014]Tourism 【树形dp +状压dp】

    题目链接 BZOJ3836 题解 显然这是个\(NP\)完全问题,此题的解决全仗任意两点间不存在节点数超过10的简单路径的性质 这意味着什么呢? \(dfs\)树深度不超过\(10\) \(10\)很 ...

  9. SCUT - 254 - 欧洲爆破 - 概率dp - 状压dp

    https://scut.online/p/254 思路很清晰,写起来很恶心. #include<bits/stdc++.h> using namespace std; #define l ...

  10. bzoj1076: [SCOI2008]奖励关(期望dp+状压dp)

    1076: [SCOI2008]奖励关 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 2989  Solved: 1557[Submit][Statu ...

随机推荐

  1. Govern EventBus - 历经多年生产环境验证的事件驱动架构框架

    Govern EventBus Govern EventBus 是一个历经四年生产环境验证的事件驱动架构框架, 通过事件总线机制来治理微服务间的远程过程调用. 使用本地事务来支持微服务内强一致性,事件 ...

  2. liunx磁盘挂载操作

    目标: 欲挂载目录:/home欲挂载磁盘:/dev/sdd 1.查看磁盘使用情况 [root@localhost ~]# df -h 文件系统 容量 已用 可用 已用% 挂载点 devtmpfs 47 ...

  3. 【转载】CentOS 7 系统区域(语言)和键盘设置

    CentOS 7 系统区域(语言)和键盘设置   即使是在window中,平常说的语言设置这一项也是归类为系统区域,CentOS可以通过修改/etc/locale.conf配置文件或使用localec ...

  4. Linux上使用iSCSI概述

    iSCSI简介 1. scsi和iscsi SCSI技术是存储设备最基本的标准协议,通常需要设备互相靠近并用SCSI总线连接,因此受到物理环境的限制 iSCSI(Internet Small Comp ...

  5. centos7安装google-chrome

    完整的安装步骤:https://www.tecmint.com/install-google-chrome-on-redhat-centos-fedora-linux/ 1.简单安装测试版:sudo ...

  6. python3 使用random函数批量产生注册邮箱

    '''你是一个高级测试工程师,现在要做性能测试,需要你写一个函数,批量生成一些注册使用的账号. 1.产生的账号是以@163.com结尾,长度由用户输,产生多少条也由用户输入,2.用户名不能重复,用户名 ...

  7. 安装beanstalkd队列问题——No package beanstalkd available

    CentOS7.4安装beanstalkd 时无可用源 No package beanstalkd availableError:Nothing to do 可从以下获取:wget /etc/yum. ...

  8. tomcat---starup.bat点击窗口自动关闭

  9. Spring中的注解@

    @ResponseBody作用:@ResponseBody注解用于将Controller的方法返回的对象,通过springmvc提供的HttpMessageConverter接口转换为指定格式的数据如 ...

  10. [LeetCode] 1744. 你能在你最喜欢的那天吃到你最喜欢的糖果吗?

    都儿童节了,为什么要折磨一个几百个月大的孩子? 把题意读懂挺难的.不过读懂后基本也就知道怎么做了.恶心的是int类型可能会越界,要用long类型(很难想到).这题不好 [1744. 你能在你最喜欢的那 ...