P4242【NOIP2016 DAY2】愤怒的小鸟

输入格式



输出格式

对每个关卡依次输出一行答案。

输出的每一行包含一个正整数，表示相应的关卡中，消灭所有小猪最少需要的小鸟数量。

题解

状压dp

先两两枚举鸟坐标求出抛物线表达式并枚举可击杀数量，将其视作一个商品存入，然后背包dp。

注意精度 用rep来限制精度查

代码

#include<stdio.h>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
#define maxn 20
#define maxn1 270001
#define maxn2 500
#define inf 1e9
const double rep=1e-7;
struct node{
double x,y;
};
node pos[maxn];
int v[maxn2],f[maxn1];
int t,n,m;
int full;
int num,mc,cnt;
void work()
{
int i,j,m;
double a,b;
cnt=0;
for(i=1;i<=n;i++)
{
for(j=i+1;j<=n;j++)
{

double x1,x2,y1,y2;
x1=pos[i].x;y1=pos[i].y;
x2=pos[j].x;y2=pos[j].y;
a=(y1*x2-y2*x1)/(x2*x1)/(x1-x2);
if(a>=0) continue;
b=y1/x1-a*x1;
int val=0;
cnt++;
for(m=1;m<=n;m++){
if(fabs(a*pos[m].x*pos[m].x+b*pos[m].x-pos[m].y)<=rep){
val|=(1<<(m-1));
}

}

v[cnt]=val;

}
}
for(i=1;i<=n;i++){
v[++cnt]=1<<(i-1);
}
}
void dp()
{
int i,j,k;
int ful=((1<<n)-1);
for(i=1;i<=ful;i++) f[i]=inf;
f[0]=0;
for(i=1;i<=cnt;i++)
{
for(j=0;j<=ful;j++)
{
if(f[j]!=inf){
f[j|v[i]]=min(f[j|v[i]],f[j]+1);
}
}
}
}
int main()
{
int i,j;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++) scanf("%lf%lf",&pos[i].x,&pos[i].y);
work();
dp();
cout<<f[(1<<n)-1]<<endl;
}
}