二进制的力量

状态压缩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. JS String总结

    String常用总结 1.字符 length属性:表示字符串包含多少16位码元 charAt():方法返回给定索引位置的字符 charCodeAt() :可以查看指定码元的字符编码 String.fr ...

  2. C# 多线程技术

    这节讲一下多线程(Thread)技术. 在讲线程之前,先区分一下程序,进程,线程三者的区别,大体上说,一个程序可以分为多个进程,一个进程至少由一个线程去执行,它们是层层包含的关系.我们写的程序,就是一 ...

  3. Ping命令浅析

    Ping ​ Ping基于ICMP协议. ​ Ping可以分为 内网Ping 和 外网Ping ​ 下面以内网Ping为例,使用的软件是eNSP和WireShark ​ ​ Step1.创建拓扑,PC ...

  4. Azure Storage 利用 azCopy 复制迁移数据

    一,引言 前两天遇到了Azure Blob Storage 需要迁移到另外的一个 Azure Blob Storage 中.手动下载.上传已经无法满足了,得另寻一种方式了 AzCopy.Azure 为 ...

  5. Java集合详解(二):ArrayList原理解析

    概述 本文是基于jdk8_271版本进行分析的. ArrayList是Java集合中出场率最多的一个类.底层是基于数组实现,根据元素的增加而动态扩容,可以理解为它是加强版的数组.ArrayList允许 ...

  6. Sentinel导航

    简介 最近都在弄微服务的东西,现在来记录下收获.我从一知半解到现在能从0搭建使用最大的感触有两点 1.微服务各大组件的版本很多,网上很多博客内容不一定适合你的版本,很多时候苦苦琢磨都是无用功 2.网上 ...

  7. zimbra启用SMTP认证

    zmprov modifyServer {{ you domain }} zimbraMtaTlsAuthOnly FALSE zmcontrol restart 查看对应配置 zmprov getS ...

  8. [刷题] 198 House Robber

    要求 你是一个小偷,每个房子中有价值不同的宝物,但若偷连续的两栋房子,就会触发报警系统,求最多可偷价值多少的宝物 示例 [3,4,1,2],返回6[3,(4),1,(2)] [4,3,1,2],返回6 ...

  9. 【山外笔记-SVN命令】svn命令详解

    本文打印版文件下载地址 [山外笔记-SVN命令]svn命令详解-打印版.pdf 一.命令简介 svn命令用于Subversion命令行客户端,执行svn相关的操作. 二.命令语法 1.svn语法: ( ...

  10. nginx 的三种虚拟主机配置方法

    nginx三种虚拟主机配置的方法. 基于端口 在生产环境中一般使用端口或者域名. [root@web01 /etc/nginx/conf.d]# cat web01.conf server { lis ...