noip2016 day2T3题解

又是一道游戏题，与2014的day1T3都是鸟。。。。

看一下题，哇还真是愤怒的小鸟。。。

好 吐槽完了。。。

先把每一个鸟看成一个点，设第i个鸟的坐标为(xi,yi)

然后从原点发射一条线，这个线是开口向下的抛物线

显然我们要把这条线求出来，但只有两个点，而二次函数要三个点

怎么办？ 仔细想一下，最优解是什么，一条线过所有点。。。

所以我们每次假设这条抛物线可以打两个点 A(xa,ya) B(xb,yb)
由于过原点 所以y=ax^2+bx
a=(yb-ya×xb/x1)/(xb^2-xa×xb);
b=(ya-a×xa^2)/x1;
然后注意开口向上的线要减掉
xa=xb的时候不用算直接减掉

然后就可以求出来所有的线了
线可能有重复的减掉(具体实现见代码)

再怎么办，看一下n的范围 n<=18

想到了什么 暴搜？状态压缩？

对 就是状压

我们用1表示这个点已经打了，0表示还没有打

然后我们把线也变成这种类似的形式
1表示能打，0表示不能打
计算某一条线是否能过i的一个点
这里要注意double的精度问题

dp(S)表示到当前解所用的线的条数
然后 不难发现从S用某条线的到s dp(s)=min(dp(S)+1,dp(S)) s=S|l
考虑达一开始S=0 这就不是求0的sssp吗

果断用SPFA
开一个队列 一开始只有0
然后枚举每一条线，如果从S到s的方法所求到的s比原来的s小，就将s压到队列里面
当队列空的的时候，答案就算出来了

其实还有更多的细节。。。
在代码里面

ps: aeb表示的是a 乘以 10的b次方

pps:其他神犇都是用dp做的。。。。。
好像就我一个用的SPFA。。。。。。

代码

#include<iostream>
#include<cstdio>
#include<queue>
using namespace std;
//
const int MAX=20;
const int MAXN=1<<20;
const double MIN=1e-8;
const int INF=MAXN;
//
int T;
int n,m;
int l[MAX*MAX+10];
double x[MAX+10],y[MAX+10];
int vis[MAX+10]; 
int dp[MAXN];
int cur;
//
queue<int> Q;
void initialize(void);
void read(void);
void solve(void);
void getl(void);
bool cmp(double,double);
void SPFA(void);
//
int main(){
    cin>>T;
    while(T--){
        initialize();
        read();
        solve();
    }
}
//
void initialize(void){
    cur=0;
    for(int i=0;i<=n;i++) vis[i]=0;
    for(int i=0;i<=n;i++) l[i]=0;
}
void getl(void){
    int ll;
    for(int i=1;i<=n;i++){
        for(int j=i+1;j<=n;j++){
            if(x[i]==x[j]) continue;
            double a,b;
            double x1=x[i],x2=x[j];double y1=y[i],y2=y[j];
            a=(y2-y1*x2/x1)/(x2*x2-x1*x2);
            b=(y1-a*x1*x1)/x1;            
            if(a>=0) continue;
            ll=0;    
            for(int t=1;t<=n;t++){
                if(cmp(x[t]*x[t]*a+x[t]*b,y[t])){
                    vis[t]=1;    
                    ll+=(1<<(t-1));
                }
            }
            int o=1;
            for(int i=0;i<cur;i++){
                if(ll==l[i]){
                    o=0;
                    break;
                }
            }
            if(o) l[cur]=ll;
            else continue;
            cur++;
        }
    }
    for(int i=1;i<=n;i++){
        if(!vis[i]){
            l[cur]=(1<<(i-1));
            cur++; 
        }
    }
}
void read(void){
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        cin>>x[i]>>y[i];
    }
}
void SPFA(void){
    dp[0]=0;
    Q.push(0);
    while(!Q.empty()){
        int S=Q.front();
        Q.pop();
        for(int i=0;i<cur;i++)
        if(dp[l[i]|S]>dp[S]+1){
            dp[l[i]|S]=dp[S]+1;
            Q.push(l[i]|S);
        }
    }
}
bool cmp(double q,double p){
    return q>p? q-p<MIN:p-q<MIN; 
}
void solve(void){
    for(int i=0;i<=(1<<(n));i++) dp[i]=INF;
    getl();
    SPFA();
    cout<<dp[(1<<(n))-1]<<endl;
}