洛谷 P5544 [JSOI2016]炸弹攻击1【模拟退火】

滑稽大佬

于 2020-10-25 20:45:55 发布

阅读量144

点赞数 1

分类专栏：模拟

本文链接：https://blog.csdn.net/qq_35786326/article/details/109278534

版权

模拟专栏收录该内容

23 篇文章 0 订阅

订阅专栏

233

题目：
题意：
分析：
代码：

题目：

传送门

题意：

我们的攻击可以看作是一个半径不超过 $R$ 的圆
地图上有许多圆形的建筑和一些敌人，敌人只是一个坐标点
问在不碰到建筑的情况下最多能杀多少敌人

分析：

我们随机攻击的圆心坐标，然后注意到 $n$ 和 $m$ 都很小，所以在计算贡献时可以直接根据题意算
比较特殊的是我们要求的是最大值，所以大小关系和平时的模拟退火相反

代码：

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
#include<cstdlib>
#include<ctime>
#define down 0.9976
#define LL long long
using namespace std;
inline LL read() {
    LL d=0,f=1;char s=getchar();
    while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}
    while(s>='0'&&s<='9'){d=d*10+s-'0';s=getchar();}
    return d*f;
}
double x[15],y[15],r[15],p[1005],q[1005];
double n=read(),m=read(),R=read(),ans=0;/*
double get(double xx,double yy)
{
	double rr=R;
	for(double i=1;i<=n;i++) rr=min(rr,sqrt((x[i]-xx)*(x[i]-xx)+(y[i]-yy)*(y[i]-yy))-r[i]);
	double s=0;
	for(double i=1;i<=m;i++) if(sqrt((p[i]-xx)*(p[i]-xx)+(q[i]-yy)*(q[i]-yy))<=rr) s++;
	return s;
}*/

double cnt(double x1,double y1,double x2,double y2){
    return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}
double get(double xx,double yy){
    double tot=0;
    double k=R;
    for(int i=1;i<=n;i++) k=min(k,cnt(xx,yy,x[i],y[i])-r[i]);
    for(int i=1;i<=m;i++) if(cnt(xx,yy,p[i],q[i])<=k) tot++;
    return tot;
}
double ansx,ansy;
void sa()
{
	//srand(rand());
	double t=20060.404,X=ansx,Y=ansy;
	while(t>1e-18)
	{
		double xx=X+(rand()*2-RAND_MAX)*t,yy=Y+(rand()*2-RAND_MAX)*t;
		double now=get(xx,yy);
		double dif=now-ans;
		if(dif>0) ansx=X=xx,ansy=Y=yy,ans=now;
		else if(exp(-dif/t)*RAND_MAX<rand()) X=xx,Y=yy;
		t*=down;
	}
	return;
}
int main()
{
	srand(114514);
	for(int i=1;i<=n;i++) x[i]=read(),y[i]=read(),r[i]=read();
	for(int i=1;i<=m;i++) ansx+=p[i]=read(),ansy+=q[i]=read();
	ansx/=m;ansy/=m;ans=0;
    while((double)clock()/CLOCKS_PER_SEC<=0.8)  sa();
	printf("%.0lf",round(ans));
	return 0;
}