POJ 1379 Run Away(基础模拟退火)

转载请注明出处,谢谢http://blog.csdn.net/acm_cxlove/article/details/7854526       by---cxlove 

题目:找出一点,距离所有所有点的最短距离最大

http://poj.org/problem?id=1379 

随机化步长,贪心调整即可。

其实用模拟退火需要证明的~~~~%#@&%……@%¥%@证明完毕

没啥好说的,时限和精度要求都不高,随便搞吧

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<string>
#include<vector>
#include<algorithm>
#include<map>
#include<set>
#include<ctime>
#define maxn 200005
#define eps 1e-8
#define inf 2000000000
#define LL long long
#define zero(a) fabs(a)<eps
#define MOD 19901014
#define N 1000005
#define pi acos(-1.0)
using namespace std;
int t,n;
double X,Y;
double best[50];
struct Point{
	double x,y;
	bool check(){
		if(x>-eps&&x<eps+X&&y>-eps&&y<eps+Y)
			return true;
		return false;
	}
}p[1005],tp[50];
double dist(Point p1,Point p2){
	return sqrt((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double slove(Point p0){
	double ans=inf;
	for(int i=0;i<n;i++)
		ans=min(ans,dist(p[i],p0));
	return ans;
}
int main(){
	scanf("%d",&t);
	srand(time(NULL));
	while(t--){
		scanf("%lf%lf%d",&X,&Y,&n);
		for(int i=0;i<n;i++)
			scanf("%lf%lf",&p[i].x,&p[i].y);
		for(int i=0;i<15;i++){
			tp[i].x=(rand()%1000+1)/1000.0*X;
			tp[i].y=(rand()%1000+1)/1000.0*Y;
			best[i]=slove(tp[i]);
		}
		double step=max(X,Y)/sqrt(1.0*n);
		while(step>1e-3){
			for(int i=0;i<15;i++){
				Point cur,pre=tp[i];
				for(int j=0;j<35;j++){
					double angle=(rand()%1000+1)/1000.0*2*pi;
					cur.x=pre.x+cos(angle)*step;
					cur.y=pre.y+sin(angle)*step;
					if(!cur.check()) continue;
					double tmp=slove(cur);
					if(tmp>best[i]){
						tp[i]=cur;
						best[i]=tmp;
					}
				}
			}
			step*=0.85;
		}
		int idx=0;
		double ans=0;
		for(int i=0;i<15;i++){
			if(best[i]>ans){
				ans=best[i];
				idx=i;
			}
		}
		printf("The safest point is (%.1f, %.1f).\n",tp[idx].x,tp[idx].y);
	}
	return 0;
}


  • 2
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,这个问题属于技术问题。以下是一个简单的Python模拟登录POJ提交代码并抓取评测结果的代码示例: ```python import requests # 登录POJ,获取cookie def login(username, password): s = requests.Session() login_url = "http://poj.org/login" login_data = { "user_id1": username, "password1": password, "B1": "login", "url": "/" } s.post(login_url, data=login_data) return s # 提交代码 def submit_code(s, problem_id, language, source_code): submit_url = "http://poj.org/submit" submit_data = { "problem_id": problem_id, "language": language, "source": source_code } s.post(submit_url, data=submit_data) # 获取评测结果 def get_result(s, run_id): status_url = "http://poj.org/status" params = { "user_id": "", "result": "", "language": "", "top": run_id } r = s.get(status_url, params=params) table_start = r.text.find("<table cellpadding=0 cellspacing=0 border=0 width=100%>") table_end = r.text.find("</table>", table_start) table_html = r.text[table_start:table_end + 8] return table_html # 使用示例 username = "your_username" password = "your_password" problem_id = "1000" language = "G++" source_code = """ #include <iostream> using namespace std; int main() { int a, b; cin >> a >> b; cout << a + b << endl; return 0; } """ s = login(username, password) submit_code(s, problem_id, language, source_code) table_html = get_result(s, "12345678") # 替换成实际提交的run id print(table_html) ``` 其中,`login`函数模拟登录POJ并返回一个`Session`对象,`submit_code`函数提交代码,`get_result`函数获取评测结果。你可以根据实际需要修改代码中的`username`、`password`、`problem_id`、`language`和`source_code`等参数,并替换`get_result`函数中的`run_id`。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值