noip愤怒的小鸟&&vijos2008

vijos题目传送门

题目有点长

就不搬过来了

写之前一直觉得很难

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

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

了解了一波精度问题

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

所以会产生精度问题

一般取个eps=1e-10

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

刚开始写了a<=eps

很明显是错误的

我竟然没发现

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

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

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
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;
}