憤怒的小鳥(angry bird )
阿新 • • 發佈:2020-10-15
又是一款當年風靡一時的遊戲所改編出的題,但是跟遊戲有所不同的是,當小鳥打到豬之後,小鳥並不會改變路線,依舊照的原來的拋物線運動軌跡運動,所以我們這裡不難想到要用到狀態壓縮,而且求資料範圍給的很小也不難猜出需要$dp + $搜尋。
#include<bits/stdc++.h> using namespace std; #define scy(x) freopen(x".in","r",stdin); freopen(x".out","w",stdout); int t,n,m,l[23][23],low[1<<20],dp[1<<20]; const double eps=1e-8; double x[23],y[23]; inline int read(){ int x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9'){ if(ch=='-') f=-1; ch=getchar(); } while(ch>='0'&&ch<='9'){ x=(x<<1)+(x<<3)+(ch^48); ch=getchar(); } return x*f; } void equation(double &x,double &y,double a1,double b1,double c1,double a2,double b2,double c2){//二次函式係數a,b y=(a1*c2-a2*c1)/(a1*b2-a2*b1); x=(c1-b1*y)/a1; } int main(){ // scy("in"); for(int i=0;i<(1<<18);i++){ int j=1; for(;j<=18 && i&(1<<(j-1));j++); low[i]=j; } t=read(); while(t--){ memset(l,0,sizeof(l)); memset(dp,0x3f,sizeof(dp)); dp[0]=0; n=read(),m=read(); for(int i=1;i<=n;i++) scanf("%lf%lf",x+i,y+i); for(int i=1;i<=n;i++){ for(int j=1;j<=n;j++){ if(fabs(x[i]-x[j])<eps) continue; double a,b; equation(a,b,x[i]*x[i],x[i],y[i],x[j]*x[j],x[j],y[j]); if(a>-eps) continue; for(int k=1;k<=n;k++){ if(fabs(a*x[k]*x[k]+b*x[k]-y[k])<eps) l[i][j]|=(1<<(k-1)); } } } for(int i=0;i<(1<<n);i++){//狀壓開始 int j=low[i]; dp[i|1<<(j-1)]=min(dp[i|(1<<(j-1))],dp[i]+1); for(int k=1;k<=n;k++) { dp[i|l[j][k]]=min(dp[i|l[j][k]],dp[i]+1); } } printf("%d\n",dp[(1<<n)-1]); } }