[luogu2831][noip d2t3]愤怒的小鸟_状压dp

愤怒的小鸟 noip-d2t3 luogu-2831

    题目大意:给你n个点,问最少需要多少条经过原点的抛物线将其覆盖。

    注释:1<=点数<=18,1<=数据组数<=30。且规定抛物线是开口向下的。

      想法:其实一开始的想法是很偏的,就是设dp[i][j][k]表示在状态k下建立$i_{th}$和$j_{th}$的抛物线的最少条数,然后向后转移。这显然是错误的,错误原因在于... 我日,没个转移。然后看了一下lijinnn的题解...啊?切了。

      是这样的,我们通过记录每条抛物线所能覆盖的点,将其记录在数组str中,不分先后顺序。然后,我们考虑状态

        dp[s]表示达到s状态的最少条数。

      转移:

        dp[s]=min(dp[ s ] , dp[ s ^ ( s & str[ i ] ) ] + 1);

    最后,附上丑陋的代码... ...

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef double db;
db ax[30],ay[30];
int n,m;
db ka,kb;
int str[1000100],cnt,dp[300100];
bool visit[1003000];
void dispose(int x,int y)
{
db tmp1=ax[x]*ax[y]*(ax[x]-ax[y]);//------------------------
db tmp2=ay[x]*ax[y]-ay[y]*ax[x];// |
if(tmp1<0)// |
{// |
tmp1=-tmp1,tmp2=-tmp2;// |
}// |
if(fabs(tmp1)<1e-6)// |
{// |
     return;//我们在计算过这两个点和原点的抛物线解析式 |
}// |
tmp1=tmp2/tmp1;// |
if(tmp1>0)// |
{// |
return;// |
}// |
ka=tmp1;// |
kb=(ay[x]-ka*ax[x]*ax[x])/ax[x];//--------------------------
int s=0;
for(int i=1;i<=n;i++)//枚举所有的点,计算该点是否在当前枚举的抛物线之内
{
db tmp=ax[i]*ax[i]*ka+kb*ax[i];
if(fabs(tmp-ay[i])<1e-6)
{
s+=(1<<(i-1));
}
}
if(!visit[s])
{
visit[s]=1;
str[++cnt]=s;//统计出一条抛物线能够杀死的pig的状态
}
return;
}
void original()
{
memset(visit,0,sizeof visit);
memset(dp,0x3f,sizeof dp);
cnt=0;
}
int main()
{
int cases;
scanf("%d",&cases);
while(cases--)
{
original();
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&ax[i],&ay[i]);
}
for(int i=1;i<=n;i++)
{
str[++cnt]=(1<<(i-1));
for(int j=1;j<=i-1;j++)
{
dispose(i,j);
}
}
dp[0]=0;
for(int s=0;s<(1<<n);s++)
{
for(int i=1;i<=cnt;i++)
{
if(s&str[i])
{
dp[s]=min(dp[s],dp[s^(s&str[i])]+1);//转移方程
}
}
}
printf("%d\n",dp[(1<<n)-1]);
}
return 0;
}

    小结:状态的选取决定着动态规划的走势----某乎上的dalao说的

上一篇:webpack打包nodejs项目(前端代码)


下一篇:Android 用 res 中文件名获取资源 id 的方法