预处理所有可能的抛物线可以打哪些猪,存在g数组中
f[i] 表示打了猪的状态
#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]);
}
}