UOJ 265（状压dp）（NOIP2016）

最新推荐文章于 2018-09-09 17:59:54 发布

嘉伟森的猫

最新推荐文章于 2018-09-09 17:59:54 发布

阅读量266

点赞数

CC 4.0 BY-SA版权

分类专栏： OI-dp 文章标签：状压dp

本文链接：https://blog.youkuaiyun.com/KGV093/article/details/78133785

OI-dp 专栏收录该内容

83 篇文章

订阅专栏

本文提供了一种解决NOIP2016第二天第三题的方法，通过预处理每两只猪所在抛物线能打到的猪的集合，并采用状态压缩动态规划求解最优解。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

传送门
NOIP 2016 D2T3
预处理每两只猪所在抛物线能打的猪的集合，然后状压dp即可。
P.S.注意有些猪由于在一条覆盖大量目标的抛物线上，所以可能重复打，于是不能加那句被去掉的”if”。天啊加上一个错误的判断居然只丢了5分。。。
又一年了，想到去年这时候连什么是dp都不知道。。。今年要好好努力了(ง •̀_•́)ง

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define eps 1e-10
const int INF=0x3f3f3f3f;
int n,waste,f[1200000],st[300],tot;
struct MPro {
    double x,y;
}p[20];
inline void init() {
    memset(f,INF,sizeof(f));
    memset(st,0,sizeof(st));
    tot=f[0]=0;
    for (int i=1;i<=n;++i)
        for (int j=i;j<=n;++j) {
            ++tot;
            if (i^j) {
                double t=p[i].x/p[j].x,t2=p[j].x*p[i].x-p[i].x*p[i].x,t3=p[j].y*t-p[i].y;
                double a=t3/t2;
                if (a>-eps) continue;
                double b=-(a*p[i].x*p[i].x-p[i].y)/p[i].x;
                for (int k=1;k<=n;++k)
                    if (a*p[k].x*p[k].x+b*p[k].x<p[k].y+eps&&a*p[k].x*p[k].x+b*p[k].x>p[k].y-eps)
                        st[tot]+=1<<k-1;
            }
            else st[tot]=1<<i-1;
        }
}
int main() {
    int kase;
    scanf("%d",&kase);
    while (kase--)  {
        scanf("%d%d",&n,&waste);
        for (int i=1;i<=n;++i) scanf("%lf%lf",&p[i].x,&p[i].y);
        init();
        for (int i=0;i<(1<<n);++i)
            for (int j=1;j<=tot;++j)
                /*if (!(i&st[j]))*/f[i|st[j]]=min(f[i|st[j]],f[i]+1);
        printf("%d\n",f[(1<<n)-1]);
    }
    return 0;
}