愤怒的小鸟

Description:
给你\(n<=18\)个小猪,发射的小鸟轨迹为抛物线,求最小用多少个小鸟可以将小猪全部干掉
看到n很小,我想到了搜索,于是我用\(dfs\)枚举出,每个抛物线打掉的小猪集合然后判断他的合法性,结果TLE成了50分,mmp,瞄了一眼题解,看到他是枚举小猪,来确定抛物线,感觉妙了很多,于是我写了如下的
code:
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const double eps=0.0000001;
int t,n,m,ans;
double x[20],y[20],a[20],b[20];
bool vis[20];
inline double fabs(double x){
if(x<0)return -x;
else return x;
}
inline void dfs(int pos,int num,int cnt){
//当前决策第几个位置,构造了多少个抛物线,剩下几个独立的
if(pos>n){
ans=min(ans,num+cnt);
return ;
}
bool flag=0;
for(int i=1;i<=num;++i){//枚举是否可以被前面的抛物线覆盖
double xx=a[i]*x[pos]*x[pos]+b[i]*x[pos];
double yy=y[pos];
if(fabs(fabs(xx)-fabs(yy))<=eps&&(xx*yy>0)){
flag=1;
vis[pos]=1;//不独立
dfs(pos+1,num,cnt-1);
vis[pos]=0;
break;
}
}
if(flag)return ;
for(int i=1;i<pos;++i){
if(vis[i])continue;//不独立
double a1=x[i];
double b1=y[i];
double a2=x[pos];
double b2=y[pos];
double aa=(b1*a2-b2*a1)/(a1*a1*a2-a1*a2*a2);
double bb=(a1*a1*b2-a2*a2*b1)/(a1*a1*a2-a2*a2*a1);
if(aa>=0)continue;//不可行
a[num+1]=aa;
b[num+1]=bb;//可行
vis[i]=vis[pos]=1;//不独立了
dfs(pos+1,num+1,cnt-2);
a[num+1]=0;
b[num+1]=0;
vis[i]=vis[pos]=0;//回溯
}
dfs(pos+1,num,cnt);//自己独立
}
int main(){
scanf("%d",&t);
for(int i=1;i<=t;++i){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i)scanf("%lf%lf",&x[i],&y[i]);
ans=0x3f3f3f3f;
dfs(1,0,n);
cout<<ans<<endl;
}
}
80分还是TLE
于是乎,又瞄了一眼题解,他加了一个最优性剪枝\(num+cnt>=ans,return ;\)
于是我也加了一个,这下可惨了,一下WA成了40
注意看我的代码,\(dfs(pos+1,num,cnt-1)\text{&&}dfs(pos+1,num+1,cnt-2)\text{&&}dfs(pos+1,num,cnt)\)
\(num+cnt\)的总和是变小了,如果使用最优性剪枝,有可能将最优值剪掉
仔细比较题解和我的代码,他并不是把所有的小猪刚开始都变成独立的,这样对于一个新小猪,他的code

\(dfs(pos+1,num,cnt)\text{&&} dfs(pos+1,num+1,cnt-1)\text{&&}dfs(pos+1,num,cnt+1)\)

\(num+cnt\)的总和单调不降,可以使用最优性剪枝
我是上来就把所有小猪看成独立,而他是将小猪后放进去
code:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
using namespace std;
const double eps=1e-8;
int t,n,m,ans;
double x[20],y[20],a[20],b[20],xx[20],yy[20];
inline void dfs(int pos,int num,int cnt){
if(num+cnt>=ans)return ;//最优性剪枝
if(pos>n){
ans=num+cnt;
return ;
}//边界
bool flag=0;
for(int i=1;i<=num;++i){//枚举是否可以被前面的抛物线覆盖
double xx=a[i]*x[pos]*x[pos]+b[i]*x[pos];
double yy=y[pos];
if(fabs(xx-yy)<eps){
dfs(pos+1,num,cnt);
flag=1;
break;
}
}
if(!flag){
for(int i=1;i<=cnt;++i){
double a1=xx[i];
double b1=yy[i];
double a2=x[pos];
double b2=y[pos];
if(fabs(a1-a2)<=eps)continue;
double aa=(b1*a2-b2*a1)/(a1*a1*a2-a1*a2*a2);
double bb=(a1*a1*b2-a2*a2*b1)/(a1*a1*a2-a2*a2*a1);
if(aa>=0)continue;//不可行
a[num+1]=aa;
b[num+1]=bb;
double va=xx[i];
double vb=yy[i];
for(int j=i;j<cnt;++j){
xx[j]=xx[j+1];
yy[j]=yy[j+1];
}
dfs(pos+1,num+1,cnt-1);
for(int j=cnt;j>i;j--)
{
xx[j]=xx[j-1];
yy[j]=yy[j-1];
}
xx[i]=va;
yy[i]=vb;
}
xx[cnt+1]=x[pos];
yy[cnt+1]=y[pos];
dfs(pos+1,num,cnt+1);//自己独立
}
}
int main(){
scanf("%d",&t);
for(int i=1;i<=t;++i){
scanf("%d%d",&n,&m);
for(int i=1;i<=n;++i)scanf("%lf%lf",&x[i],&y[i]);
ans=0x3f3f3f3f;
dfs(1,0,0);
cout<<ans<<endl;
}
}

搜索:状态要定好,剪枝要想好

最新文章

  1. 3.Java基础之Date对象
  2. 每日Scrum(7)
  3. 给日志添加&ldquo;复制&rdquo;效果
  4. mvc5 知识点01
  5. text输入框中按下enter键时阻止刷新页面
  6. FreeMarker 乱码解决方案 生成静态html文件
  7. DataDictionaryTool 一款生成数据库字典工具支持mysql和oracle
  8. 【Linux】环境变量设置
  9. Linux中利用crontab创建计划任务
  10. C# winform调用东软动态库的问题
  11. jQuery 学习笔记一
  12. 【Solution】API测试工具,访问方式
  13. DOM访问和处理HTML文档的标准方法
  14. 一套高可用、易伸缩、高并发的IM群聊架构方案设计实践
  15. C# Excel转换成Json工具(含源码)
  16. nginx-1.12.1编译参数详情
  17. 【cocos2dx中Node类getParent和getChildByTag()】学习体会
  18. git编译
  19. 使用 shinydashboard
  20. 【targeting学习笔记】Display Advertising Targeting

热门文章

  1. drupal7 怎样将一个date字段加入上日期插件效果
  2. hdu5391Zball in Tina Town
  3. Java - Thinking in Java 第2章 一切都是对象
  4. 自己封装js组件 - 初级
  5. cpc,tank
  6. poj_3667线段树区间合并
  7. thinkphp5项目--个人博客(八)
  8. BZOJ 1711 网络流
  9. CentOS7下安装二进制MYSQL8
  10. CUDA笔记(十)