Kiana最近沉迷于一款神奇的游戏无法自拔。
简单来说，这款游戏是在一个平面上进行的。
有一架弹弓位于(0,0)处，每次Kiana可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如的曲线，其中a,b是Kiana指定的参数，且必须满足a<0。
当小鸟落回地面(即x轴）时，它就会瞬间消失。
在游戏的某个关卡里，平面的第一象限中有n只绿色的小猪，其中第i只小猪所在的坐标为(xi,yi)。
如果某只小鸟的飞行轨迹经过了(xi,yi),那么第i只小猪就会被消灭掉，同时小鸟将会沿着原先的轨迹继续飞行；
如果一只小鸟的飞行轨迹没有经过(xi,yi)，那么这只小鸟飞行的全过程就不会对第i只小猪产生任何影响。
例如，若两只小猪分别位于(1,3)和(3,3)，Kiana可以选择发射一只飞行轨迹为 的小鸟，这样两只小猪就会被这只小鸟一起消灭。
而这个游戏的目的，就是通过发射小鸟消灭所有的小猪。
这款神奇游戏的每个关卡对Kiana来说都很难，所以Kiana还输入了一些神秘的指令，使得自己能更轻松地完成这个游戏。这些指令将在【输入格式】中详述。
假设这款游戏一共有T个关卡，现在Kiana想知道，对于每一个关卡，至少需要发射多少只小鸟才能消灭所有的小猪。由于她不会算，所以希望由你告诉她。
输入
第一行包含一个正整数T，表示游戏的关卡总数。
下面依次输入这T个关卡的信息。每个关卡第一行包含两个非负整数n,m，分别表示该关卡中的小猪数量和Kiana输入的神秘指令类型。接下来的n行中，第i行包含两个正实数(xi,yi)，表示第i只小猪坐标为(xi,yi)。数据保证同一个关卡中不存在两只坐标完全相同的小猪。
如果m=0，表示Kiana输入了一个没有任何作用的指令。
首先m没用
然后这题数据那么小，显然状压
状压无非就是处理打掉和没打掉
那么我们可以预处理每两个小鸟的轨迹，ab可以数学推导出来，然后计算这个轨迹经过了哪一些点，计算出来之后把经过点的状态压出来，我认为有必要贴出来2的幂次方
2^1=2
2^2=4
2^3=8
2^4=16
2^5=32
2^6=64
2^7=128
2^8=256
2^9=512
2^10=1024
2^11=2048
2^12=4096
2^13=8192
2^14=16384
2^15=32768
2^16=65536
2^17=131072
2^18=262144
2^19=524288
2^20=1048576
2^21=2097152
2^22=4194304
2^23=8388608
2^24=16777216
2^25=33554432
2^26=67108864
2^27=134217728
2^28=268435456
2^29=536870912
2^30=1073741824
2^31=2147483648
2^32=4294967296
然后就暴力枚举全集就完了
重点！！！精度问题，注意在比较相等的时候精度一定要开的足够小！！！！我就是这么被卡了！！！！请开1e-6
同时还有一个问题，如果有些点只能原点到那里没有其他任何点经过（就是这个点与其他点组合成的二次函数系数一定是正数），最后直接一个个单独处理就行了，非常简单
记住特判n=1
代码：
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 using namespace std; 5 double x[20],y[20],a,b; 6 int now,s[262150]; 7 double abs(double x){ 8 if(x<0)return -x; 9 return x; 10 } 11 int main() { 12 int T; 13 cin>>T; 14 while(T--) { 15 memset(s,0x3f3f3f3f,sizeof(s)); 16 int n,m; 17 cin>>n>>m; 18 int tot=(1<<n)-1; 19 for(int i=1; i<=n; i++) { 20 cin>>x[i]>>y[i]; 21 } 22 if(n==1){ 23 cout<<1<<endl; 24 continue; 25 } 26 s[0]=0; 27 for(int i=1; i<n; i++) { 28 for(int j=i+1; j<=n; j++) {//进行两两枚举 29 if(x[i]==x[j])continue; 30 now=0; 31 b=(y[j]*x[i]*x[i]-x[j]*x[j]*y[i])*1.0/(x[i]*x[i]*x[j]-x[i]*x[j]*x[j]); 32 a=(y[i]-b*x[i])*1.0/(x[i]*x[i]*1.0); 33 if(a>=0.0)continue; 34 for(int k=1; k<=n; k++) { 35 if(abs((a*x[k]*x[k]+b*x[k])*1.0-y[k]*1.0)<=0.000001) { 36 now+=1<<(k-1); 37 } 38 } 39 for(int k=0; k<=tot; k++) { 40 s[k|now]=min(s[k|now],s[k]+1); 41 } 42 } 43 } 44 for(int i=1;i<=n;i++){ 45 int now=1<<(i-1); 46 for(int j=0;j<=tot;j++){ 47 s[j|now]=min(s[j|now],s[j]+1); 48 } 49 } 50 cout<<s[tot]<<endl; 51 } 52 }