憤怒的小鳥[狀壓DP]
阿新 • • 發佈:2018-12-19
預處理所有可能的拋物線可以打哪些豬,存在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]); } }