博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
noip愤怒的小鸟&&vijos2008
阅读量:6900 次
发布时间:2019-06-27

本文共 1300 字,大约阅读时间需要 4 分钟。

题目有点长

就不搬过来了

写之前一直觉得很难

写了之后才发现也不是特别难

其实就是状压dp,预处理出取某两个猪可以消除哪些其他的猪,然后对于每个状态,找到第一个没有消除的猪(贪心,因为迟早要杀死它),枚举它和哪个配对,复杂度O(2^n * n*n)。

了解了一波精度问题

不是所有的小数都是能用二进制表示的

所以会产生精度问题

一般取个eps=1e-10

然后注意判a<0时要用a<-eps

刚开始写了a<=eps

很明显是错误的

我竟然没发现

然后注意有一种情况是单独加一个点的

因为可能没有其他的点可以跟它凑一个抛物线

#include
#include
#include
#include
struct node{ double x,y;}e[20];int g[20][20];int f[300000];const double eps=1e-10;int main(){ int t,n,m; double a,b; scanf("%d",&t); while(t--) { memset(g,0,sizeof(g));memset(f,127,sizeof(f)); scanf("%d %d",&n,&m); for(int i=1;i<=n;i++) scanf("%lf %lf",&e[i].x,&e[i].y); for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(i!=j) { a=(e[i].y*e[j].x-e[j].y*e[i].x)/(e[i].x*e[j].x*(e[i].x-e[j].x)); b=(e[j].x*e[j].x*e[i].y-e[i].x*e[i].x*e[j].y)/(e[i].x*e[j].x*(e[j].x-e[i].x)); if(a<-eps) { for(int k=1;k<=n;k++) if(fabs(a*e[k].x*e[k].x+b*e[k].x-e[k].y)<=eps) g[i][j]|=1<<(k-1); } } f[0]=0;int top=(1<<(n))-1; for(int i=0;i<=top;i++) { int j=1; while(i&(1<<(j-1))) j++; if(j>n) continue; f[i|(1<<(j-1))]=std::min(f[i|(1<<(j-1))],f[i]+1); for(int k=j+1;k<=n;k++) f[i|g[j][k]]=std::min(f[i|g[j][k]],f[i]+1); } printf("%d\n",f[top]); } return 0;}

转载于:https://www.cnblogs.com/Brian551/p/7353018.html

你可能感兴趣的文章
python oop面向对象笔记
查看>>
python numpy模块使用笔记(更新)
查看>>
vue-cli构建项目 npm run build后应该怎么运行在本地查看效果
查看>>
unix平台下I/O聚集和分离的一种方案
查看>>
1081. Binary Lexicographic Sequence(找规律)
查看>>
Postman笔记 - 入门好文
查看>>
通过游戏来学习编程
查看>>
周记(飞船一号
查看>>
openssl初步使用
查看>>
Vue项目碰到"‘webpack-dev-server’不是内部或外部命令,也不是可运行的程序或批处理文件"报错...
查看>>
模拟任务资源管理器的小程序
查看>>
通过一个例子,总结下检测数组属性的N种方法
查看>>
Samba结合AD实现域帐号认证的文件服务器
查看>>
laravel的Eloquent中的get()和Query/Builder中的get()
查看>>
bzoj 2286(洛谷 2495) [Sdoi2011]消耗战——虚树
查看>>
51nod 1301 集合异或和——异或dp
查看>>
weblogic.jdbc.wrapper.Blob_oracle_sql_BLOB cannot be cast to oracle.sql.BLOB 解决方法
查看>>
表格边框设置
查看>>
问题 K: A/B Problem
查看>>
Django实战(7):改造ProductList界面
查看>>