1. 程式人生 > >憤怒的小鳥[狀壓DP]

憤怒的小鳥[狀壓DP]

傳送門

預處理所有可能的拋物線可以打哪些豬,存在g陣列中

f[i] 表示打了豬的狀態  

f[i|g[j]] = min(f[i|g[j]],f[i]+1)

#include<bits/stdc++.h>
#define N 20
#define M 1<<20
#define eps 0.00000001
using namespace std;
int f[M],g[N*N],cnt;
int T,n,m,vis[N];
double x[N],y[N];
bool check(double a,double b,double x,double y){
	if(fabs( a*x*x + b*x - y ) < eps) return true;
	return false;
}
int main(){
	scanf("%d",&T);
	while(T--){
		memset(f,127,sizeof(f)); f[0]=0;
		cnt=0; memset(g,0,sizeof(g));
		memset(vis,0,sizeof(vis));
		scanf("%d%d",&n,&m);
		for(int i=1;i<=n;i++){
			scanf("%lf%lf",&x[i],&y[i]);
		}
		
		for(int i=1;i<=n;i++){
			for(int j=i+1;j<=n;j++){
				double a = (y[i]*x[j] - y[j]*x[i])/((x[i]*x[j])*(x[i]-x[j]));
				double b = (y[i] - a*x[i]*x[i]) / x[i];
				if(a>-eps) continue;
				vis[i]=vis[j]=1; cnt++;
				for(int k=1;k<=n;k++){
					if(check(a,b,x[k],y[k])) g[cnt]|=1<<(k-1),vis[k]=1;
				}
			}
		}
		for(int i=1;i<=n;i++) if(!vis[i]) g[++cnt]=1<<(i-1);
		for(int i=0;i<=(1<<n)-1;i++)
			for(int j=1;j<=cnt;j++)	
				f[i|g[j]] = min(f[i|g[j]],f[i]+1);
		printf("%d\n",f[(1<<n)-1]);
	}
}