upc备战省赛组队训练赛第六场 D:Cat and Mice(二分+状压dp)

6 篇文章 0 订阅
1 篇文章 0 订阅

题意:初始你在(0,0),有n只老鼠,给出每只老鼠的位置和第几秒消失,问你猫的速度最小是多少可以保证在每只老鼠消失之前吃掉所有老鼠?

题解:首先可以想到二分速度,然后考虑如何check?一开始在想复杂度比较低的贪心,但感觉没有什么太好的决策,发现n是15,很小,便开始考虑暴力状压dp,dp[i][j]表示此时状态是i最后到达的是第j只老鼠的最小时间,当这个时间小于等于当前老鼠的消失时间再去更新dp的状态,最后看一下最终状态是否被更新即可.(Note! memset对于double类型 置最大值用 0x7f 用其他的数值可能奇奇怪怪)

#include<bits/stdc++.h>
using namespace std;
#define Sheryang main
#define TYPE int
const int maxn=2e6+7;
typedef long long ll;
const int mod=2333333;
const double eps=1e-4;
void Smax(TYPE &a,TYPE b){if(a<b) a=b;}
void Smin(TYPE &a,TYPE b){if(a>b) a=b;}
///#define getchar()(p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1 << 21, stdin), p1 == p2) ? EOF : *p1++)
///char buf[(1 << 21) + 1], *p1 = buf, *p2 = buf;
#define IO cin.tie(0),ios::sync_with_stdio(false);
#define pi acos(-1)
#define PII pair<ll,ll>
ll read(){ll c = getchar(),Nig = 1,x = 0;while(!isdigit(c) && c!='-')c = getchar();if(c == '-')Nig = -1,c = getchar();while(isdigit(c))x = ((x<<1) + (x<<3)) + (c^'0'),c = getchar();return Nig*x;}
#define read read()
/** keep hungry and keep calm! **/
  

struct node{
	double x,y,t;
}p[20];
double dp[1<<16][20],vv,dis[20][20];
int num[1<<16];
int n;

void init(){
	for(int i=0;i<n;i++){
		for(int j=i+1;j<n;j++){
			dis[i][j]=sqrt((p[i].x-p[j].x)*(p[i].x-p[j].x)+(p[i].y-p[j].y)*(p[i].y-p[j].y));
			dis[j][i]=dis[i][j];
		}
	}
	for(int i=1;i<(1<<n);i++){
		for(int j=0;j<n;j++){
			if(i&(1<<j)){
				num[i]++;
			}
		}
	}
}

int check(double v){
	memset(dp,100,sizeof dp);
	
	for(int i=1;i<(1<<n);i++){
		for(int j=0;j<n;j++){
			if(i&(1<<j)){
				if(i==(1<<j)){
					double tmp=sqrt(p[j].x*p[j].x+p[j].y*p[j].y)/v;
					if(tmp<=p[j].t) dp[i][j]=min(dp[i][j],tmp);
				}
				for(int k=0;k<n;k++){
					if(k!=j && i&(1<<k)){
						double tmp=dp[i^(1<<j)][k]+dis[k][j]/(v*pow(vv,num[i]-1));
						if(tmp<=p[j].t) dp[i][j]=min(dp[i][j],tmp);
					}
				}
			}
		}
	}
	double res=1e11;
	for(int i=0;i<n;i++){
		res=min(res,dp[(1<<n)-1][i]);
	}
	return res<1e10;
}

int Sheryang()
{
	n=read;
	
	for(int i=0;i<n;i++){
		scanf("%lf %lf %lf",&p[i].x,&p[i].y,&p[i].t);
	}
	scanf("%lf",&vv);
	
	init();
	double l=0,r=1e10;
    while(r-l>=eps)
    {
        double mid=(l+r)/2.0;
        check(mid)?r=mid:l=mid;
    }

	printf("%.6f\n",l);
    return 0;
}

 

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值