noip2016 愤怒的小鸟 luogu 2831

最新推荐文章于 2022-03-20 09:58:33 发布
最新推荐文章于 2022-03-20 09:58:33 发布
阅读量158 收藏 1
点赞数
分类专栏： 动态规划
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/elijahqi/article/details/79690835
版权

http://www.elijahqi.win/archives/535
题目描述

Kiana最近沉迷于一款神奇的游戏无法自拔。

简单来说，这款游戏是在一个平面上进行的。

有一架弹弓位于(0,0)处，每次Kiana可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如y=ax2+bxy=ax^2+bxy=ax​2​​+bx的曲线，其中a,b是Kiana指定的参数，且必须满足a<0。

当小鸟落回地面(即x轴）时，它就会瞬间消失。

在游戏的某个关卡里，平面的第一象限中有n只绿色的小猪，其中第i只小猪所在的坐标为(xi,yi)。

如果某只小鸟的飞行轨迹经过了(xi,yi),那么第i只小猪就会被消灭掉，同时小鸟将会沿着原先的轨迹继续飞行；

如果一只小鸟的飞行轨迹没有经过(xi,yi)，那么这只小鸟飞行的全过程就不会对第i只小猪产生任何影响。

例如，若两只小猪分别位于(1,3)和(3,3)，Kiana可以选择发射一只飞行轨迹为y=−x2+4xy=-x^2+4xy=−x​2​​+4x的小鸟，这样两只小猪就会被这只小鸟一起消灭。

而这个游戏的目的，就是通过发射小鸟消灭所有的小猪。

这款神奇游戏的每个关卡对Kiana来说都很难，所以Kiana还输入了一些神秘的指令，使得自己能更轻松地完成这个游戏。这些指令将在【输入格式】中详述。

假设这款游戏一共有T个关卡，现在Kiana想知道，对于每一个关卡，至少需要发射多少只小鸟才能消灭所有的小猪。由于她不会算，所以希望由你告诉她。

输入输出格式

输入格式：

第一行包含一个正整数T，表示游戏的关卡总数。

下面依次输入这T个关卡的信息。每个关卡第一行包含两个非负整数n,m，分别表示该关卡中的小猪数量和Kiana输入的神秘指令类型。接下来的n行中，第i行包含两个正实数(xi,yi)，表示第i只小猪坐标为(xi,yi)。数据保证同一个关卡中不存在两只坐标完全相同的小猪。

如果m=0，表示Kiana输入了一个没有任何作用的指令。

如果m=1，则这个关卡将会满足：至多用⌈n3+1⌉\left \lceil \frac{n}{3} + 1 \right \rceil⌈​3​​n​​+1⌉只小鸟即可消灭所有小猪。

如果m=2,则这个关卡将会满足：一定存在一种最优解，其中有一只小鸟消灭了至少⌊n3⌋\left \lfloor \frac{n}{3} \right \rfloor⌊​3​​n​​⌋只小猪。

保证1<=n<=18，0<=m<=2，0< xi,yi<10，输入中的实数均保留到小数点后两位。

上文中，符号⌈x⌉\left \lceil x \right \rceil⌈x⌉和⌊x⌋\left \lfloor x \right \rfloor⌊x⌋分别表示对c向上取整和向下取整

输出格式：

对每个关卡依次输出一行答案。

输出的每一行包含一个正整数，表示相应的关卡中，消灭所有小猪最少需要的小鸟数量

输入输出样例

输入样例#1：

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
1
输入样例#2：

3
2 0
1.41 2.00
1.73 3.00
3 0
1.11 1.41
2.34 1.79
2.98 1.49
5 0
2.72 2.72
2.72 3.14
3.14 2.72
3.14 3.14
5.00 5.00
输出样例#2：

2
2
3
输入样例#3：

1
10 0
7.16 6.28
2.02 0.38
8.33 7.78
7.68 2.09
7.46 7.86
5.77 7.44
8.24 6.72
4.42 5.11
5.42 7.79
8.15 4.99
输出样例#3：

6
洛谷评测48ms 本蒟蒻感觉应该挺快的了吧。（洛谷看到还有0ms 的神犇qwq)

这题是状态压缩

放一个80分的程序 应该是有一定的问题的，按理说洛谷应该可以跑85才对，有一个点wa了 蒟蒻我看不出来ORZ

思路：

因为一共18只鸟所以我们把18只鸟用二进制表示最后十进制储存

1表示编号为这一位的鸟还存活0表示已经挂了 因为有小数 注意判断下精度

开结构体a[i][j]表示编号i j两只鸟 exist表示是否一条抛物线可以打下两只鸟，a,b分别表示二次函数里的a和b

b数组的下标存的是状态 若为1表示这种状态存在

正式开始记忆化搜索：for (;s1>=1;s1=s&(s1-1)) 这个式子十分巧妙的把当前状态所有存在的鸟被打后的状态一一枚举出来不重不漏，好好推导一下

因为s1比s一定要小 所以直接十进制减法就可以代替状态里的鸟的存活状态的加减运算

提示：状态压缩时，假如只能打一只鸟b数组也要认为是1 就是认为打一只鸟也可以

鸟的标号从0开始比较好做，要么就在状压的时候统一好

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 20
#define N1 270000                   
int b[N1],data[N1],n,m;
double x[N],y[N];
struct node{
    int exist;double a,b;
}a[N][N];
inline bool judge(double a,double b){
    if (std::abs(a-b)<=1e-10) return true;else return false;
}
inline void work(int i,int j){
    if (judge(x[i],x[j])) {
        a[i][j].exist=0;a[i][j].a=a[i][j].b=0;return ;
    }
    double tmpa=(y[i]*x[j]-y[j]*x[i])/(x[i]*x[j]*(x[i]-x[j]));
    double tmpb=y[i]/x[i]-(y[i]*x[j]-y[j]*x[i])/(x[j]*(x[i]-x[j]));
    a[i][j].a=tmpa;a[i][j].b=tmpb;
    if (tmpa>0){
        a[i][j].exist=0;return ;
    }
    a[i][j].exist=1;
}
//for (int s1=s;s1>=1;s1=s&(s1-1)) 
void f(int s){
    if (b[s]) {data[s]=1;return;}
    int s1=s;s1=s&(s1-1);
    for (;s1>=1;s1=s&(s1-1)){
        if (data[s1]==50) f(s1);
        if (data[s-s1]==50) f(s-s1);
        data[s]=std::min(data[s],data[s1]+data[s-s1]);
    }
}
int main(){
//  freopen("ex1.in","r",stdin);
//  freopen("ex1.out","w",stdout);
    int T;
    scanf("%d",&T);
    while (T--){
        scanf("%d%d",&n,&m);
        for (int i=0;i<n;++i) scanf("%lf%lf",&x[i],&y[i]);
        memset(b,false,sizeof(b));
        for (int i=0;i<n;++i)
            for (int j=0;j<n;++j) work(i,j); 
    //  for (int i=0;i<n;++i)
    //      for (int j=i+1;j<n;++j) printf("%f %f\n",a[i][j].a,a[i][j].b);
        for (int i=0;i<n;++i) b[1<<i]=1;
        for (int i=0;i<n;++i){  
            for (int j=i+1;j<n;++j){
                int tmp=0;
                if (a[i][j].exist){
                    tmp+=(1<<i)+(1<<j);
                    for (int z=j+1;z<n;++z){
                        if (i!=z&&j!=z)
                            if (a[z][i].exist&&a[i][z].a&&a[i][j].a&&judge(a[i][z].a,a[i][j].a)&&judge(a[i][j].b,a[i][z].b)) {
                                tmp+=(1<<z);//  if (tmp==25) printf("%d %d %d\n",i,j,z);
                            }
                    }

                    b[tmp]=true;
                }
            }
        }   
        for (int i=1;i<=(1<<n)+3;++i) data[i]=50;
        f((1<<n)-1);
        printf("%d\n",data[(1<<n)-1]);
    //  for (int i=1;i<=1050;++i) if (b[i]) printf("%d\n",i);
        //printf("asffd\n");
    }
    return 0;
}

这个是满分的程序：

b1数组表示i,j两只鸟之间连成的抛物线打完鸟后最终平面内还剩余和不剩余的鸟 这种状态转成十进制储存在 b1数组中

这种做法因为我们肯定每次钦定都要打第一只鸟，然后枚举

s&((1<


#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define N 20
#define N1 270000                   
int b1[N][N],b[N1],data[N1],n,m;
double x[N],y[N];
struct node{
    int exist;double a,b;
}a[N][N];
inline bool judge(double a,double b){
    if (std::abs(a-b)<=1e-10) return true;else return false;
}
inline void work(int i,int j){
    if (judge(x[i],x[j])) {
        a[i][j].exist=0;a[i][j].a=a[i][j].b=0;return ;
    }
    double tmpa=(y[i]*x[j]-y[j]*x[i])/(x[i]*x[j]*(x[i]-x[j]));
    double tmpb=y[i]/x[i]-(y[i]*x[j]-y[j]*x[i])/(x[j]*(x[i]-x[j]));
    a[i][j].a=tmpa;a[i][j].b=tmpb;
    if (tmpa>0){
        a[i][j].exist=0;return ;
    }
    a[i][j].exist=1;
}
//for (int s1=s;s1>=1;s1=s&(s1-1)) 
/*void f(int s){
    if (b[s]) {data[s]=1;return;}
    int s1=s;s1=s&(s1-1);
    for (;s1>=1;s1=s&(s1-1)){
        if (data[s1]==50) f(s1);
        if (data[s-s1]==50) f(s-s1);
        data[s]=std::min(data[s],data[s1]+data[s-s1]);
    }
}*/
void f(int s){
    if (s==0) {data[s]=0;return;}
    int rest=0;for (int i=0;i<n;++i) if(s&(1<<i))++rest;
    if (rest==1){data[s]=1;return;} 
    if (rest==2) if(b[s]) {data[s]=1;return;}else {data[s]=2;return;}
    int i=0;while(!(s&(1<<i))) ++i;
    for (int j=0;j<n;++j)
        if (b1[i][j]) {
            if (data[s&((1<<n)-1-b1[i][j])]==50) f(s&((1<<n)-1-b1[i][j]));
            data[s]=std::min(data[s],data[s&((1<<n)-1-b1[i][j])]+1);
        }
}
int main(){
    freopen("ex1.in","r",stdin);
//  freopen("ex1.out","w",stdout);
    int T;
    scanf("%d",&T);
    while (T--){
        scanf("%d%d",&n,&m);
        for (int i=0;i<n;++i) scanf("%lf%lf",&x[i],&y[i]);
        memset(b,false,sizeof(b));memset(b1,0,sizeof(b1));
        for (int i=0;i<n;++i)
            for (int j=0;j<n;++j) work(i,j); 
    //  for (int i=0;i<n;++i)
    //      for (int j=i+1;j<n;++j) printf("%f %f\n",a[i][j].a,a[i][j].b);
        for (int i=0;i<n;++i) b[1<<i]=1,b1[i][i]=1<<i;
        for (int i=0;i<n;++i){  
            for (int j=i+1;j<n;++j){
                int tmp=0;
                if (a[i][j].exist){
                    tmp+=(1<<i)+(1<<j);
                    for (int z=j+1;z<n;++z){
                        if (i!=z&&j!=z)
                            if (a[z][i].exist&&a[i][z].a&&a[i][j].a&&judge(a[i][z].a,a[i][j].a)&&judge(a[i][j].b,a[i][z].b)) {
                                tmp+=(1<<z);//  if (tmp==25) printf("%d %d %d\n",i,j,z);
                            }
                    }
                    b1[i][j]=tmp;
                    b[tmp]=true;
                }
            }
        }   
        for (int i=1;i<=(1<<n);++i) data[i]=50;
        f((1<<n)-1);
        printf("%d\n",data[(1<<n)-1]);
    //  for (int i=1;i<=1050;++i) if (b[i]) printf("%d\n",i);
        //printf("asffd\n");
    }
    return 0;
}

递推：不好想，我脑子不太好，先贴一个leoly (ORZ)的

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 20
#define eps 1e-10
using namespace std;
struct aa{double x, y;}a[N];
int T, n, temp, bin[N], x, w[N][N], f[262150];
double x1, y1, x2, y2, A, B;
inline double abs1(double a){return max(a, -a);}
int main(){
    scanf("%d", &T);
    bin[0]=1; for(int i=1; i<=18; i++)bin[i]=bin[i-1]<<1;
    for(int owo=1; owo<=T; owo++){
        scanf("%d%d", &n, &temp);
        for(int i=1; i<=n; i++)scanf("%lf%lf", &a[i].x, &a[i].y);
        for(int i=1; i<=n-1; i++)
            for(int j=i+1; j<=n; j++){
                x1=a[i].x; y1=a[i].y; x2=a[j].x; y2=a[j].y;
                if(abs1(x1-x2)<eps){w[i][j]=bin[n]-1; continue;}
                A=(x1*y2-x2*y1)/(x1*x2*(x2-x1)); B=(y1-A*x1*x1)/x1;
                if(A>-eps){w[i][j]=bin[n]-1; continue;}
                w[i][j]=bin[n]-1;
                for(int k=1; k<=n; k++)
                    if(abs1(A*a[k].x*a[k].x+B*a[k].x-a[k].y)<eps)w[i][j]-=bin[k-1];
            }
        f[0]=0;
        for(int i=1; i<=bin[n]-1; i++){
            for(int j=n; j; j--)if(i&bin[j-1]){x=j; break;}
            f[i]=f[i-bin[x-1]]+1;
            for(int j=1; j<=x-1; j++)f[i]=min(f[i], f[i&w[j][x]]+1);
        }
        printf("%d\n", f[bin[n]-1]);
    }
    return 0;
}

2017.11.1 update 用真题当作模拟赛 然而怎么也写不满 现在还是95分 先弃坑待填

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
#define N 20
#define eps 1e-6
#define inf 0x3f3f3f3f
struct node{
    double x,y;
}bird[N];int dp[1<<N],n,bin[N],b1[N][N],T,m;
void f(int s){
    if (s==0){dp[s]=0;return;}if (dp[s]!=inf) return;
    int num=0;
    for (int i=0;i<n;++i) if (bin[i]&s)++num;
    if (num==1){dp[s]=1;return;}
    for (int i=0;i<n;++i){
        if (s&bin[i]){
            for (int j=0;j<n;++j){
                int s1=b1[i][j];if (!s1) continue;int s2=s^s1;f(s2);
                dp[s]=min(dp[s],1+dp[s2]);
            }break;
        }
    }
}
int main(){
//  freopen("angrybirds.in","r",stdin);
    scanf("%d",&T);for (int i=0;i<N;++i) bin[i]=1<<i;
    while (T--){
        scanf("%d%d",&n,&m);
        for (int i=0;i<n;++i) scanf("%lf%lf",&bird[i].x,&bird[i].y);memset(dp,0x3f,sizeof(dp));memset(b1,0,sizeof(b1));
        for (int i=0;i<n;++i) b1[i][i]=1<<i;
        for (int i=0;i<n;++i){
            for (int j=i+1;j<n;++j){
                int x1=i,x2=j;
                if (abs(bird[x1].x-bird[x2].x)<=eps) continue;
                double a=(bird[x1].y*bird[x2].x-bird[x2].y*bird[x1].x)/(bird[x1].x*bird[x2].x*(bird[x1].x-bird[x2].x));
                double b=bird[x1].y/bird[x1].x-(bird[x1].y*bird[x2].x-bird[x2].y*bird[x1].x)/(bird[x2].x*(bird[x1].x-bird[x2].x));
                if (a>=0) continue;int s=bin[x1]|bin[x2];
                for (int k=j+1;k<n;++k){
                    if (k==x1||k==x2) continue;
                    if (abs(bird[k].y-a*bird[k].x*bird[k].x-b*bird[k].x)<=eps) s|=bin[k];
                }dp[s]=1;b1[i][j]=s;
            }
        }f(bin[n]-1);printf("%d\n",dp[bin[n]-1]);
    }
    return 0;
}
elijahqi
关注
专栏目录
【洛谷题解/NOIP2016提高组】P2831 愤怒的小鸟
lovely_cyqの小窝
08-22 410
NOIP2016提高组愤怒的小鸟题解。
Luogu P2831 [NOIP2016 提高组] 愤怒的小鸟
fturceke的博客
11-13 701
Luogu P2831 [NOIP2016 提高组] 愤怒的小鸟
524. 愤怒的小鸟
GitHub：https://github.com/myp020301
03-20 3864
Kiana 最近沉迷于一款神奇的游戏无法自拔。 简单来说，这款游戏是在一个平面上进行的。 有一架弹弓位于 (0,0) 处，每次 Kiana 可以用它向第一象限发射一只红色的小鸟， 小鸟们的飞行轨迹均为形如 y=ax2+bx 的曲线，其中 a,b 是 Kiana 指定的参数，且必须满足 a<0。 当小鸟落回地面（即 x 轴）时，它就会瞬间消失。 在游戏的某个关卡里，平面的第一象限中有 n 只绿色的小猪，其中第 i 只小猪所在的坐标为 (xi,yi)。 如果某只小鸟的飞行轨迹经过了 (xi, yi)，那么
NOIP2016 愤怒的小鸟——解题报告
风中的微尘的博客
08-14 492
题目描述: Kiana最近沉迷于一款神奇的游戏无法自拔。 简单来说，这款游戏是在一个平面上进行的。 有一架弹弓位于(0,0)处，每次Kiana可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如y=ax2+bx的曲线，其中a,b是Kiana指定的参数，且必须满足a&amp;amp;lt;0。 当小鸟落回地面（即x轴）时，它就会瞬间消失。 在游戏的某个关卡里，平面的第一象限中有n只绿色的小猪，其...
愤怒的小鸟（NOIP2016提高组Day2T3）
11-26 256
【题目描述】 Kiana最近沉迷于一款神奇的游戏无法自拔。 简单来说，这款游戏是在一个平面上进行的。 有一架弹弓位于（0,0)处，每次Kiana可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如y=ax^2+bx的曲线，其中a,b是Kiana指定的参数，且必须满足a<0。 当小鸟落回地面（即x轴）时，它就会瞬间消失。 在游戏的某个关卡里，平...
[LUOGU]2831 愤怒的小鸟 [NOIP]2016 D2T3
Le Petit Prince
09-19 678
题目描述 Kiana最近沉迷于一款神奇的游戏无法自拔。 简单来说，这款游戏是在一个平面上进行的。 有一架弹弓位于(0,0)处，每次Kiana可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如y=ax2+bxy=ax^2+bxy=ax​2​​+bx的曲线，其中a,b是Kiana指定的参数，且必须满足a 当小鸟落回地面(即x轴）时，它就会瞬间消失。 在游戏的某个关卡里，平面的第一
【NOIP2016】愤怒的小鸟 搜索
weixin_30653097的博客
10-06 326
题目描述 Kiana最近沉迷于一款神奇的游戏无法自拔。 简单来说，这款游戏是在一个平面上进行的。 有一架弹弓位于(0,0)(0,0)处，每次Kiana可以用它向第一象限发射一只红色的小鸟，小鸟们的飞行轨迹均为形如y=ax^2+bxy=ax2+bx的曲线，其中a,ba,b是Kiana指定的参数，且必须满足a < 0a<0，a,ba,b都是实数。 当小鸟落回...
【题解】LuoGu2831/noip2016：愤怒的小鸟
ModestCoder_的博客
09-24 294
原题传送门 【题解】 我先想到了dfs，估计了一下时间复杂度，估不出来，觉得悬，就“瞟了一眼”题解，发现第一篇就是dfs，然后信心满满的码好，莫名wa掉，找了很长时间错也没找出来。 然后我又想起一年前某位老师讲状压dp的时候讲过这道题，就改用了状压dp dp[i]表示猪的状态为i时，用掉的鸟的最小值。 对于i，比如i=1001,那么表示第1、4只猪死掉了 预处理：num[i][j]表示过i,j两只...
Luogu 2831 [NOIP2016] 愤怒的小鸟
dashu497731727的博客
08-25 142
第一眼看成爆搜的状压dp，膜Chester大神犇。 考虑到三个不在同一直线上的点可以确定一条抛物线，而固定点$(0, 0)$和不在同一直线上这两个条件是题目中给定的，所以我们只要枚举两个点然后暴力算抛物线，然后chk一遍观察一下多少点在这一条抛物线上就行了。 想到状压之后状态和方程就显然了。 注意判解出来的抛物线$a \leqslant 0$的情况， eps开小一点。 因为决策集...
【NOIP2016】洛谷2831 愤怒的小鸟
YHaX的博客
11-28 2288
状压dp
NOIP2016普及组与提高组试题：NOIP2016复赛普及组试题&NOIP2016复赛提高组试题
11-20
《NOIP2016普及组与提高组试题解析》 全国青少年信息学奥林匹克联赛（National Olympiad in Informatics in Provinces，简称NOIP）是中国计算机学会主办的一项旨在选拔优秀中学生进入更高层次信息学竞赛的活动。...
解题报告——noip 2016提高组
11-14
### NOIP 2016 提高组解题报告 #### DAY1 T1 玩具谜题 **题目来源**：洛谷1563 **吐槽**：本题中，作者对NOIP出题人将“魔法师”这一英文单词拼写错误表示了一种幽默的调侃。 **思路分析**： 此题可通过模拟来...
NOIP2016普及组初赛复赛真题及解答与测试数据
02-05
总的来说，这套资源为准备参加NOIP2016普及组比赛的初中生提供了一个全面的学习平台，不仅有实际的编程题目，还有相应的测试数据和参考答案，有助于参赛者在实践中提升编程技巧和解题策略。对于信息学教育工作者而言...
基于springboot个人公务员考试管理系统源码数据库文档.zip
最新发布
11-18
基于springboot个人公务员考试管理系统源码数据库文档.zip
bimdata_api_client-4.2.1-py3-none-any.whl
11-18
bimdata_api_client-4.2.1-py3-none-any.whl
numpy-1.20.2-cp39-cp39-linux_armv7l.whl
11-18
numpy-1.20.2-cp39-cp39-linux_armv7l.whl
matplotlib-3.3.2-cp39-cp39-linux_armv7l.whl
11-18
matplotlib-3.3.2-cp39-cp39-linux_armv7l.whl
bimdata_api_client-4.0.0-py3-none-any.whl
11-18
bimdata_api_client-4.0.0-py3-none-any.whl
ta_lib-0.5.1-cp312-cp312-win32.whl
11-18
ta_lib-0.5.1-cp312-cp312-win32.whl
NOIP2016初赛C++试题与解析
"这是第二十二届全国青少年信息学奥林匹克联赛初赛普及组的C++语言试题，包含单项选择题，涉及计算机基础知识，如软件识别、颜色编码、通信技术、硬件生产商、键盘输入状态、二进制运算、进制转换、计算机体系结构...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值