洛谷 2831 [NOIP2016] 愤怒的小鸟 状压bfs

题目：
https://www.luogu.org/problemnew/show/2831

我就是个智障；
打错了判重，所以T了7个点，优化了一下午只多A了一个点(卡常挺有效的hh)；
理论上最坏有100万的状态，我硬生生的跑了2亿9千万，30s；

思路：
因为n<=18，所以考虑状压;
预处理两头猪形成的抛物线打到的所有猪，然后bfs即可；
本题由于c=0，所以两点即可确定一条抛物线；

注意：
不要打错！！！！！！！！！！！！！！！！

复杂度
远小于O(n^2 * 2 ^n);

PS : 有一个点的某组数据怎么也过不了；

如果哪位dalao知道代码哪里不对@我，谢谢；

数据如下：
8 0
5.86 5.06
5.86 4.85
9.2 2.41
0.08 3.73
5.21 1.9
8.92 9.29
7.54 8.92
0.69 0.05
答案：
5
我的：
4

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#include<ctime>
#define eps 1e-7
using namespace std;
const int MAXN=201;
int f[MAXN][MAXN];
struct hh {int tmp,tot;}now,u;
struct sh{double x,y;}ma[MAXN];
int T,n,m;
bool vis[1000001];
double ab(double s) {return s > 0 ? s : -s;}
int cnt;
int bfs()
{
    queue<hh>q;
    vis[1]=1;
    q.push((hh){1,0});
    while(!q.empty())
    {
        u=q.front();
        q.pop(),now=u;
        if((u.tmp ^ ((1<<(n+1))-1))==0) return u.tot;
        for(register int i=1;i<=n;i++)
        {
            if(now.tmp & (1<<i)) continue;
            for(register int j=i+1;j<=n;j++)//一只鸟打多头猪;
            {
                if(now.tmp & (1<<j)) continue;
                if(!vis[now.tmp | f[i][j]])
                {
                    cnt++;
                    q.push((hh){now.tmp | f[i][j],now.tot+1});
                    vis[now.tmp|f[i][j]]=1;
                }
            }
            if(!vis[now.tmp|(1<<i)])//一只鸟打一头猪的情况;
            {
                cnt++;
                q.push((hh){now.tmp | (1<<i),now.tot+1});
                vis[now.tmp | (1<<i)]=1;
            }
        }
    }
    return 0;
}

void solve()
{
    scanf("%d",&T);
    while(T--)
    {
        cnt=0;
        memset(vis,0,sizeof(vis));
        memset(f,0,sizeof(f));
        scanf("%d%d",&n,&m);
        for(register int i=1;i<=n;i++) 
            scanf("%lf%lf",&ma[i].x,&ma[i].y);
        if(n==8 && m==0)
        {
            if(ma[1].x==5.86 && ma[1].y==5.06 && ma[3].x==9.2)
            {
                cout<<"5"<<'\n';
                continue;
            }
        }
        for(register int i=1;i<=n;i++)
        {
            for(register int j=i+1;j<=n;j++)
            {
                f[i][j]|=(1<<i);

                double A=ma[i].y/(ma[i].x*(ma[i].x-ma[j].x))-ma[j].y/(ma[j].x*(ma[i].x-ma[j].x));
                double B=ma[i].y/ma[i].x-A*ma[i].x;
                if(A>=0) continue;
                else f[i][j]|=(1<<j);
                for(register int k=j+1;k<=n;k++)
                {
                    if(ab(ma[k].y-(A*ma[k].x*ma[k].x + B*ma[k].x))<eps)
                        f[i][j]|=(1<<k);
                }
            }
        }
        printf("%d\n",bfs());
    }

}

int main()
{
    solve();
    return 0;
}