Description
魔法炮来到了帝都,除了吃特色菜之外,还准备去尝一尝著名的北京烤鸭。帝都一共有n(1<=1<=100)个烤鸭店,可以看成是二维平面内的点。不过由于魔法炮在吃烤鸭之前没有带钱,所以吃完烤鸭之后只能留下刷盘子。刷完盘子之后,魔法炮除了不用付饭费之外,老板还会奖励他d(1<=d<=10000)元钱。魔法炮是一个特么喜欢吃烤鸭的孩子,所以在去过一家烤鸭店之后,魔法炮还准备去其他的烤鸭店。但是由于帝都路费较贵,每单位长度需要花费1元钱,所以魔法炮可能去不了所有其他的烤鸭店。在到达下一家烤鸭店之前,魔法炮会花掉手里所有钱,以便于下一次接着吃霸王餐。另外,魔法炮对于自己刷过盘子的烤鸭店有着特殊的感情,所以他要求在某一家烤鸭店吃完烤鸭后,可以到达全部已经吃过去过的烤鸭店。那么问题来了,魔法炮想知道自己最多能去多少家烤鸭店,以及这些烤鸭店都是哪些。你能帮帮他吗?
题目大意:给定平面内的n个点,选出一个点集S,使得S里的所有点两两之间欧几里得距离不超过d,问|S|的最大值以及S里的点都有哪些。若答案有多种,输出任意一个。
第一行两个整数n和d,分别表示烤鸭店数和老板给魔法炮的路费。
数据范围:1<=n<=100 1<=d<=10000 -10000<=x,y<=10000
Solution
可以发现要求一个最大团,随机化一下可以屮过去
一开始看成了曼哈顿距离,老年选手打出gg
Code
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <vector>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define fill(x,t) memset(x,t,sizeof(x))
const int N=205;
std:: vector <int> prt,rec;
struct pos {int x,y,id;} d[N];
int n,m;
bool vis[N];
int read() {
int x=0,v=1; char ch=getchar();
for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
return x*v;
}
int sqr(int x) {
return x*x;
}
int get_dis(pos a,pos b) {
return sqr(a.x-b.x)+sqr(a.y-b.y);
}
int solve() {
fill(vis,0); int ret=0;
rec.clear();
rep(i,1,n) {
if (vis[i]) continue;
ret++; rec.push_back(d[i].id);
rep(j,i+1,n) {
int dis=get_dis(d[i],d[j]);
if (dis>m) vis[j]=1;
}
}
return ret;
}
int main(void) {
srand(19260817);
n=read(),m=read(); m=m*m; int ans=0;
rep(i,1,n) {
d[i].x=read(); d[i].y=read(); d[i].id=i;
}
rep(i,1,1000) {
std:: random_shuffle(d+1,d+n+1);
int ret=solve();
if (ret>ans) {
ans=ret;
prt=rec;
}
}
printf("%d\n", ans);
for (int i=0;i<prt.size();i++) printf("%d ", prt[i]);
return 0;
}