题解 P2831 【愤怒的小鸟】
我的天,这题是真的卡精度......
主要是精度很不好处理,经本蒟蒻测验,精度在\(10^{-6}\)会比较好优雅
【分析】
对于这种某个变量特别小\((\leq 31)\)的题目,本蒟蒻第一反应就是状压
对于某个抛物线,一定要打到起码一个小猪(不然不如不要这一条抛物线)
有人觉得最少会打掉两只小猪的,可以仔细想一下,万一\(a \geq 0\)呢......
好的,我们继续
那么,我们可以这么考虑,枚举每一只小猪的坐标
首先,有一条抛物线是只过它的
其次,再枚举其他小猪,算出抛物线方程(见下方),若\(a<0\)就记录这条抛物线,否则可以直接跳出
记录完这条抛物线后,我们枚举其他的小猪,查看是否在线上
这里要注意精度问题,如果计算出来的\(y\)与题目所给的\(y\)偏差值不超过\(10^{-6}\),就直接视为同一个值(毕竟计算是有一定的精度问题)
那么对于接下来,我已经预处理过了所有的小鸟的轨迹(抛物线),只需要在状压方程中直接算就可以了
在打\(0\)只小猪的时候,需要用\(0\)只小鸟,于是有:
\(dp[0]=0\)
假设当前状态为\(i\),抛物线为第\(j\)条,抛物线打掉的小猪状态为\(para[j]\),那么有:
\(dp[i|para[j]]=min(dp[i|para[j]],dp[i]+1)\)
接下来我们说一下怎么求\(a\)和\(b\):
假设枚举到两个小猪,坐标分别为\((x_1,y_1)\)和\((x_2,y_2)\),那么就对应地会有:
\(\begin{cases} y_1=ax_1^2+bx_1\\ \\y_2=ax_2^2+bx_2\end{cases}\)
这里有一个很明显的矩阵关系:
\(\left[\begin{matrix}\ y_1\ \\ \\y_2\end{matrix}\right]=\left[\begin{matrix}\ x_1^2&x_1\ \\ \\x_2^2&x_2\end{matrix}\right]\times\left[\begin{matrix}\ a\ \\ \\b\end{matrix}\right]\)
于是有:
\(\left[\begin{matrix}\ a\ \\ \\b\end{matrix}\right]=\left[\begin{matrix}\ x_1^2&x_1\ \\ \\x_2^2&x_2\end{matrix}\right]^{-1}\times\left[\begin{matrix}\ y_1\ \\ \\y_2\end{matrix}\right]\)
又因为:
\(\left|\begin{matrix}\ x_1^2&x_1\\ \\x_2^2&x_2\end{matrix}\right|=x_1^2x_2-x_1x_2^2=(x_1-x_2)x_1x_2\)
所以有:
\(\left[\begin{matrix}\ x_1^2&x_1\ \\ \\ x_2^2&x_2\end{matrix}\right]^{-1}=\frac{1}{(x_1-x_2)x_1x_2}\left[\begin{matrix}\ x_2&-x_1\ \\ \\-x_2^2&x_1^2\end{matrix}\right]\)
(逆矩阵的求法)
所以有:
\(\left[\begin{matrix}\ a\ \\ \\b\end{matrix}\right]={1\over(x_1-x_2)x_1x_2}\left[\begin{matrix}\ x_2&-x_1\ \\ \\-x_2^2&x_1^2\end{matrix}\right]\times\left[\begin{matrix}\ y_1\ \\ \\y_2\end{matrix}\right]\)
即:
\(\begin{cases}a={1\over(x_1-x_2)x_1x_2}\times(x_2y_1-x_1y_2)\\ \\b={1\over(x_1-x_2)x_1x_2}\times(x_1^2y_2-x_2^2y_1)\end{cases}\)
还有一件事
我们对于抛物线,在后面枚举状态的时候是不需要知道除了能打的小猪以外的所有信息
所以直接维护这个信息就可以了,其他都没有必要维护
这个众位大犇可以直接一下本蒟蒻的代码
【代码】
那本蒟蒻就放代码了:
#include<cstdio>
#include<cstring>
using namespace std;
#define min(a,b) ((a<b)?a:b)
inline void built(double &a,double &b,double x1,double y1,double x2,double y2){
a=(x2*y1-x1*y2)/(x1*x2*(x1-x2));
b=(x1*x1*y2-x2*x2*y1)/(x1*x2*(x1-x2));
}//计算a,b
inline bool inc(double a,double b,double x,double y){
double abs=a*x*x+b*x-y;
if(abs<0) abs=-abs;
return abs<=0.000001;
}//判定某个小猪是否在抛物线上
inline int read(){
int ans=0;char c=getchar();bool neg=0;
while((c<'0')|(c>'9')) { neg^=!(c^'-'); c=getchar(); }
while((c>='0')&(c<='9')) { ans=(ans<<3)+(ans<<1)+c-'0'; c=getchar(); }
return neg?-ans:ans;
}//无聊的读入优化
int n,para[200],dp[1<<18],countpara;
inline void pre(){
memset(dp,0x3f,sizeof(dp));
countpara=0;
double x[18]={0},y[18]={0};
n=read();read();
for(int i=0;i<n;i++) scanf("%lf %lf",&x[i],&y[i]);
//dp定义为无限大,抛物线的条数清空,读入
for(int i=0;i<n;i++){
para[countpara++]=(1<<i);
//只打一只小猪的
for(int j=i+1,vis=0;j<n;j++)//定义vis表示打到过的小猪,避免重复枚举
if((vis>>j)&1) continue;
else{
double a,b;
built(a,b,x[i],y[i],x[j],y[j]);
if(a>=0) continue;
para[countpara]=(1<<i);
for(int k=j;k<n;k++)//枚举小猪,查看是否在线上
if(inc(a,b,x[k],y[k])){
vis|=(1<<k);
para[countpara]|=(1<<k);
}
countpara++;
}
}
}
inline int ans(){//状压
dp[0]=0;
for(int i=0;i<(1<<n);i++)
for(int j=0;j<countpara;j++)
dp[i|para[j]]=min(dp[i|para[j]],dp[i]+1);
return dp[(1<<n)-1];
}
int main(){
int t=read();
while(t--){
pre();//先皮一下
printf("%d\n",ans());
}
return 0;
}
最后安利一下 本蒟蒻的博客
题解 P2831 【愤怒的小鸟】的更多相关文章
- 【题解】P2831 愤怒的小鸟 - 状压dp
P2831愤怒的小鸟 题目描述 \(Kiana\) 最近沉迷于一款神奇的游戏无法自拔. 简单来说,这款游戏是在一个平面上进行的. 有一架弹弓位于 \((0,0)\) 处,每次 \(Kiana\) 可以 ...
- 洛谷P2831 愤怒的小鸟
洛谷P2831 愤怒的小鸟 原题链接 题解 首先简单数学公式送上. \(ax_1^2+bx_1=y_1\) \(ax_2^2+bx_2=y_2\) \(ax_1^2x_2+bx_1x_2=y_1x_2 ...
- P2831 愤怒的小鸟
P2831 愤怒的小鸟 从 \((0, 0)\) 发射一只鸟, 轨迹满足抛物线, 问最少几只鸟可以打完 \(n <= 18\) 只猪 错误日志: 处理抛物线数组没有初始化 Solution 数据 ...
- P2831 愤怒的小鸟(状压dp)
P2831 愤怒的小鸟 我们先预处理出每个猪两两之间(设为$u,v$)和原点三点确定的抛物线(当两只猪横坐标相等时显然无解) 处理出$u,v$确定的抛物线一共可以经过多少点,记为$lines[u][v ...
- 洛谷 P2831 愤怒的小鸟
P2831 愤怒的小鸟 题目描述 Kiana 最近沉迷于一款神奇的游戏无法自拔. 简单来说,这款游戏是在一个平面上进行的. 有一架弹弓位于 (0,0)(0,0) 处,每次 Kiana 可以用它向第一象 ...
- P2831 愤怒的小鸟——状压
P2831 愤怒的小鸟 抛物线过原点,只要再找两个就能确定抛物线: 处理出两两之间的抛物线能过哪些点,状态压缩: 但是直接枚举每一条抛物线常数太大会T,所以我们需要预处理一个low_bit表示当前状态 ...
- Luogu P2831 愤怒的小鸟(状压+记忆化搜索)
P2831 愤怒的小鸟 题意 题目描述 Kiana最近沉迷于一款神奇的游戏无法自拔. 简单来说,这款游戏是在一个平面上进行的. 有一架弹弓位于\((0,0)\)处,每次Kiana可以用它向第一象限发射 ...
- [Luogu P2831] 愤怒的小鸟 (状压DP)
题面: 传送门:https://www.luogu.org/problemnew/show/P2831 Solution 首先,我们可以先康一康题目的数据范围:n<=18,应该是状压或者是搜索. ...
- 【luogu P2831 愤怒的小鸟】 题解
题目链接:https://www.luogu.org/problemnew/show/P2831 写点做题总结:dp,搜索,重在设计状态,状态设的好,转移起来也方便. 对于一条抛物线,三点确定.(0, ...
随机推荐
- 123-PHP类构造函数
<?php class ren{ //定义人类 private $name; //定义成员属性 public function __construct($name){ //定义构造函数 $thi ...
- git 在企业里的基本操作
拖下来码云上的代码: git add . 若把单个文件加入到暂存区,则用git add 某文件 若把所有文件加入到暂存区,则使用git add . git commit -m"提交" ...
- flink初识及安装flink standalone集群
flink architecture 1.可以看出,flink可以运行在本地,也可以类似spark一样on yarn或者standalone模式(与spark standalone也很相似),此外fl ...
- 吴裕雄--天生自然C++语言学习笔记:C++ 数据结构
C/C++ 数组允许定义可存储相同类型数据项的变量,但是结构是 C++ 中另一种用户自定义的可用的数据类型,它允许存储不同类型的数据项. 结构用于表示一条记录,假设想要跟踪图书馆中书本的动态,可能需要 ...
- slave_net_timeout, MASTER_HEARTBEAT_PERIOD, MASTER_CONNECT_RETRY,以及 MASTER_RETRY_COUNT设置和查看
在主从复制中,有几个参数是非常重要的,包括slave_net_timeout, MASTER_HEARTBEAT_PERIOD, MASTER_CONNECT_RETRY,以及 MASTER_RETR ...
- Ctrl + c 强制退出
linux 使用cat 命令进入一个很大的文件,一直在不停地显示内容,让它停止的方法,可以用Ctrl + c 强制退出
- [NCTF2019]Fake XML cookbook
0x00 知识点 XXE攻击 附上链接: https://xz.aliyun.com/t/6887 XXE(XML External Entity Injection)全称为XML外部实体注入 XML ...
- HDU - 4430 Yukari's Birthday(二分+枚举)
题意:已知有n个蜡烛,过生日在蛋糕上摆蜡烛,将蜡烛围成同心圆,每圈个数为ki,蛋糕中心最多可摆一个蜡烛,求圈数r和看,条件为r*k尽可能小的情况下,r尽可能小. 分析:n最大为1012,k最少为2,假 ...
- BZOJ:2186: [Sdoi2008]沙拉公主的困惑
问题:可能逆元不存在吗? 题解: Gcd(a,b)==Gcd(b,a-b); 从数据范围可以看出应该求M!的欧拉函数: 然后通过Gcd转化过去 一开始没想到 #include<iostream& ...
- Pickle incompatibility of numpy arrays between Python 2 and 3
Traceback (most recent call last): File "mnist.py", line 7, in <module> train_set, v ...