JSOI2016 炸弹攻击Attack

炸弹攻击
【故事背景】
JYY 最近迷上了一款塔防游戏,在游戏里 JYY 除了建设建筑,还可以使用
炸弹对屏幕上的敌人进行范围杀伤。
【问题描述】
游戏地图可以简单认为是一个 2 维平面。
JYY 建造了 N 个建筑,每个建筑都是一个圆,其中第 i 个建筑的圆心位于
(x ” , y ” )且半径为r ” 。
地图上一共有 M 个敌人,一个敌人可以近似看成一个平面上的点,其中第 i
个敌人位于(p ” , q ” )。
JYY 可以使用一枚可以设置半径的炸弹,可以设置一个不超过 R的范围,然
后选择平面上的一个点引爆,范围内的所有敌人全部消灭。当然,由于炸弹威力
巨大,如果爆炸范围接触到 JYY 的建筑,那么 JYY 的建筑也会受到损伤。
(注:
如果炸弹的爆炸范围仅接触到了 JYY 建筑的边界,则不会对 JYY 的建筑造成损
伤;如果敌人出现在了爆炸范围的边界,则该敌人被消灭)
JYY 可以自由控制炸弹的爆炸地点和爆炸半径。作为一个保守的玩家,他希
望在保证自己建筑毫发无损的情况下,消灭尽量多的敌人。
【输入格式】
从文件 attack.in 中读入数据。
第一行包含三个非负整数,分别为 N,M,R。
接下来 N 行,每行 3 个整数,其中第 i 的为 x ” , y ” , r ” ,表示第 i 个建筑的位置
和半径。数据保证所有建筑不相交(但是有可能边界接触)。
接下来 M 行,每行 2 个整数,其中第 i 行为 p ” , q ” ,表示第 i 个敌人的位置。
【输出格式】
输出到文件 attack.out 中。
输出一行一个整数,表示 JYY 最多可以消灭的敌人数量。
【输入样例 1】
1 5 3
0 0 1
3 3
-3 3
3 -3
3 0
0 3
第 6 页
共 7 页江苏省省队第一轮选拔赛
第一试 炸弹攻击
【输出样例 1】
3
【输入样例 2】
4 10 100
0 0 3
10 0 3
10 10 3
0 10 3
0 4
0 5
0 6
5 3
5 -3
5 5
6 7
3 6
10 4
8 4
【输出样例 2】
5
【样例说明】
第一个样例中,最佳攻击选择应将炸弹在 3,3 引爆并将半径设置为 3。
【数据规模】
对于 20%的数据满足 M = 2。
对于另外 20%的数据,满足 N = 0;
对于另外 20%的数据,满足 M ≤ 50;
对 于 100% 的 数 据 , 满 足 0 ≤ N ≤ 10, 0 < M ≤ 10 D , 1 ≤ R, r ” ≤ 20000 ,
p ” , q ” , x ” , y ” ≤ 20000。

---

计算几何太弱了QwQ 只会写模拟退火> <
对着数据调到了100…不对着数据大概可以50-80吧？

#include <bits/stdc++.h>#define N 20#define M 1010#define INF INT_MAX#define eps 1e-7using namespace std;int n,m,R,ans(0);double Dis(INF);struct circle{    int x,y,r;}c[N];struct point{    double x,y;    point(){}    point(int d){        x=y=d;    }}e[M],cur,p,tmp,l(INF),r(-INF);inline double Rand(){    return (double)rand()/RAND_MAX;}inline double calc(point x){    double r=R,dis(0),t; int cnt=0;    for (int i=1;i<=n;i++)        r=min((double)sqrt((c[i].x-x.x)*(c[i].x-x.x)+(c[i].y-x.y)*(c[i].y-x.y))-c[i].r,r);    if (r<0) return 0;    for (int i=1;i<=m;i++)        if ((t=(double)sqrt((e[i].x-x.x)*(e[i].x-x.x)+(e[i].y-x.y)*(e[i].y-x.y))-r)<eps)            cnt++;        else            dis+=t;    if (cnt>ans || cnt==ans && dis<Dis)        ans=cnt,p=x,Dis=dis;    return dis;}void SA(double T){    cur.x=l.x+rand()%(int)(r.x-l.x+1);    cur.y=l.y+rand()%(int)(r.y-l.y+1);    double de;    for (calc(cur);T>0.01;T*=0.9998){        tmp.x=cur.x+T*(Rand()*2-1);        tmp.y=cur.y+T*(Rand()*2-1);        if (tmp.x<l.x || tmp.y<l.y || tmp.x>r.x || tmp.y>r.y)            continue;        de=calc(cur)-calc(tmp);        if (de>0 || Rand()<exp(de/T))            cur=tmp;    }    T=1000;    for (int i=1;i<=100000;i++){        tmp.x=p.x+T*(Rand()*2-1);        tmp.y=p.y+T*(Rand()*2-1);        calc(tmp);    }}int main(){    freopen("attack.in","r",stdin);    freopen("attack.out","w",stdout);    srand(20010823);    scanf("%d%d%d",&n,&m,&R);    for (int i=1;i<=n;i++)        scanf("%d%d%d",&c[i].x,&c[i].y,&c[i].r);    for (int i=1;i<=m;i++){        scanf("%lf%lf",&e[i].x,&e[i].y);        l.x=min(e[i].x-R,l.x),l.y=min(l.y,e[i].y-R);        r.x=max(e[i].x+R,r.x),r.y=max(r.y,e[i].y+R);    }    SA(max(r.x-l.x,r.y-r.x));    printf("%d\n",ans);    return 0;}