DP - 状压DP - NOIP2016 - 愤怒的小鸟

DP - 状压DP - NOIP2016 - 愤怒的小鸟

题意：

$T组测试用例，每组包括n个点的坐标。$

$现需要用过原点的抛物线y=a{x}^2+bx(a<0)，来覆盖这n个点，问至少需要多少条抛物线。$

输入样例：
2
2 0
1.00 3.00
3.00 3.00
5 2
1.00 5.00
2.00 8.00
3.00 9.00
4.00 8.00
5.00 5.00
输出样例：
1
1

数据范围：

$$\begin{gathered} 1<=n<=18，0<x_i,y_i<10。 \\ 时/空限制：2s/128MB \end{gathered}$$

---

分析：

$对于方程y=a{x}^2+bx，任意两点可以确定参数a和b，n个点最多有n^2条抛物线。$

$问题转化为从n^2条抛物线中最少选择几条，能够覆盖n个点，这n^2条抛物线可以先暴力预处理出来。$

$$\begin{gathered} 接着可以用一个n位二进制数表示哪些点已经被覆盖了，对于没有被覆盖的点x，我们再枚举所有能够覆盖x的抛物线， \\ 若所有点均被覆盖了(n位全1)，就更新一下使用的抛物线条数。 \end{gathered}$$

$$\begin{gathered} 解参数a,b： \\ 对于(x_1,y_1)和(x_2,y_2)所确定的抛物线，有方程组\{\begin{array}{l}{y}_1=a{x}_1^2+b{x}_1\\ y_2=a{x}_2^2+b{x}_2\end{array}，即\{\begin{array}{l}\frac{y_1}{x_1}=a{x}_1+b\\ \frac{y_2}{x_2}=a{x}_2+b\end{array}, \\ 两式作差，得\{\begin{array}{l}a=\frac{\frac{y_1}{x_1}-\frac{y_2}{x_2}}{x_1-x_2}\\ \\ b=\frac{y_1}{x_1}-a{x}_1\end{array},(x_1\ne x_2且x_1\ne 0)。 \end{gathered}$$

具体落实：

$$\begin{gathered} f[state]:表示覆盖状态为state需要的抛物线条数。 \\ path[i][j]:表示点i和点j所确定的抛物线能够覆盖的二进制状态。 \end{gathered}$$

$①、预处理path数组，尤其注意两个点的横坐标不能相同，抛物线开口要朝下(a<0)。$

$$\begin{gathered} ②、DP求f数组，对于每一个状态i，先任选一个该状态还未覆盖到的点x， \\ 接着枚举所有覆盖x点的抛物线所能覆盖的状态path[x][j],j\in [0,n-1]， \\ i|path[x][j]即表示加上这条新的抛物线后的状态，f[i|path[x][j]]=min(f[i|path[x][j]],f[i]+1)。 \\ 即更新新状态的最小值。 \end{gathered}$$

注意：$$\begin{gathered} 函数abs(intx)，是求整数绝对值的。 \\ 而函数fabs(doublex)，是求小数绝对值的。 \end{gathered}$$

代码：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>

#define P pair<double,double>
#define x first
#define y second

using namespace std;

const int N=20,M=1<<18;
double eps=1e-8;

int n,m,T;
int path[N][N];
int f[M];
P v[N];

int cmp(double a,double b)
{
    if(fabs(a-b)<eps) return 0;
    if(a<b) return -1;
    return 1;
}

int main()
{
    cin>>T;
    while(T--)
    {
        cin>>n>>m;
        for(int i=0;i<n;i++) cin>>v[i].x>>v[i].y;
        
        memset(path,0,sizeof path);
        for(int i=0;i<n;i++)
        {
            path[i][i]=1<<i;
            for(int j=0;j<n;j++)
            {
                double x1=v[i].x,y1=v[i].y;
                double x2=v[j].x,y2=v[j].y;
                if(!cmp(x1,x2)) continue;
                double a=(y1/x1-y2/x2)/(x1-x2);
                double b=y1/x1-a*x1;
                if(cmp(a,0)>=0) continue;
                
                int State=0;
                for(int k=0;k<n;k++)
                {
                    double x=v[k].x,y=v[k].y;
                    if(!cmp(a*x*x+b*x,y))
                        State+=1<<k;
                }
                
                path[i][j]=State;
            }
        }

        memset(f,0x3f,sizeof f);
        f[0]=0;
        for(int i=0;i+1<1<<n;i++)
        {
            int x=0;
            for(int j=0;j<n;j++)
                if(!(i>>j&1))
                {
                    x=j;
                    break;
                }
                
            for(int j=0;j<n;j++)
                f[i|path[x][j]]=min(f[i|path[x][j]],f[i]+1);
        }
        
        cout<<f[(1<<n)-1]<<endl;
    }
    
    return 0;
}