MyException - 我的异常网
当前位置:我的异常网» 编程 » 【BZOJ3680】吊打XXX 狭义费马点 模拟退火

【BZOJ3680】吊打XXX 狭义费马点 模拟退火

www.MyException.Cn  网友分享于:2015-02-11  浏览:0次
【BZOJ3680】吊打XXX 广义费马点 模拟退火

#include <stdio.h>
int main()
{
	puts("转载请注明出处[vmurder]谢谢");
	puts("网址:blog.csdn.net/vmurder/article/details/43526909");
}


!!!其实我脸一点也不黑!我天天洗脸的!


题解:

我的姿势是先模拟退火,然后少少爬下山来取优。


参数什么的看代码就好。

那个种子的生成方式是[生日^名字首字母的hash]


代码:

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 10100
#define inf 23333333333333333ll
#define eps 1e-3
#define down 0.993
using namespace std;
double ln,rn,lm,rm;
struct Point // 点
{
	double x,y,z; // 坐标、该点作为结束点的答案。
	Point(double _x=0.0,double _y=0.0):x(_x),y(_y){}
	bool in()
	{
		if(x<ln||x>rn)return 0;
		if(y<lm||y>rm)return 0;
		return 1;
	}
}P[N],now,ans;
int p;
inline double Rand(){return ((double)(rand()%1000+1.0))/1000.0;}
inline void init()
{
	double tempx=0.0,tempy=0.0;
	for(int i=1;i<=p;i++)
	{
		tempx+=P[i].x,tempy+=P[i].y;
		ln=min(P[i].x,ln),rn=max(P[i].x,rn);
		lm=min(P[i].y,lm),rm=max(P[i].y,rm);
	}
	now=Point(tempx/p,tempy/p);
	ans.z=inf;
}
inline double calc(const Point &a,const Point &b)
{return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
inline double Calc(Point &a)
{
	a.z=0.0;
	for(int i=1;i<=p;i++)a.z+=calc(P[i],a)*P[i].z;
	if(a.z<ans.z)ans=a;
	return a.z;
}
Point newpos(double T,const Point &a)
{
	double alp=2.0*acos(-1.0)*Rand();
	return Point(a.x+T*cos(alp)*Rand(),a.y+T*sin(alp)*Rand());
}
void Std_Anne(double T=sqrt((rn-ln)*(rn-ln)+(rm-lm)*(rm-lm))*0.3)
{
	for(;T>eps;T*=down)
	{
		Point next=newpos(T,now);
		if(next.in())
		{
			double dis=Calc(now)-Calc(next);
			if(exp(dis/T)>=Rand())now=next;
		}
	}
}
void ioi()
{
	for(int i=1;i<=1000;i++)
	{
		Point next=newpos(0.001,ans);
		if(next.in())
		{
			double dis=Calc(next)-Calc(now);
			if(dis>0)now=next;
		}
	}
}
int main()
{
	freopen("test.in","r",stdin);
	srand(20134858);
	scanf("%d",&p);
	for(int i=1;i<=p;i++)scanf("%lf%lf%lf",&P[i].x,&P[i].y,&P[i].z);
	init();
	Std_Anne(1000000.0);
	ioi();
	printf("%.3lf %.3lf\n",ans.x,ans.y);
	return 0;
}


文章评论

不懂技术不要对懂技术的人说这很容易实现
不懂技术不要对懂技术的人说这很容易实现
Web伟德国际app者需具备的8个好习惯
Web伟德国际app者需具备的8个好习惯
科技史上最臭名昭著的13大罪犯
科技史上最臭名昭著的13大罪犯
初级 vs 高级伟德国际app者 哪个性价比更高?
初级 vs 高级伟德国际app者 哪个性价比更高?
团队中“技术大拿”并非越多越好
团队中“技术大拿”并非越多越好
当下全球最炙手可热的八位少年创业者
当下全球最炙手可热的八位少年创业者
如何成为一名黑客
如何成为一名黑客
我是如何打败拖延症的
我是如何打败拖延症的
Google伦敦新总部 犹如星级庄园
Google伦敦新总部 犹如星级庄园
那些性感的让人尖叫的程序员
那些性感的让人尖叫的程序员
如何区分一个程序员是“老手“还是“新手“?
如何区分一个程序员是“老手“还是“新手“?
为什么程序员都是夜猫子
为什么程序员都是夜猫子
程序员最害怕的5件事 你中招了吗?
程序员最害怕的5件事 你中招了吗?
程序员眼里IE浏览器是什么样的
程序员眼里IE浏览器是什么样的
2013年中国软件伟德国际app者薪资调查报告
2013年中国软件伟德国际app者薪资调查报告
程序员和编码员之间的区别
程序员和编码员之间的区别
聊聊HTTPS和SSL/TLS协议
聊聊HTTPS和SSL/TLS协议
旅行,写作,编程
旅行,写作,编程
中美印日四国程序员比较
中美印日四国程序员比较
Java 与 .NET 的平台发展之争
Java 与 .NET 的平台发展之争
那些争议最大的编程观点
那些争议最大的编程观点
看13位CEO、创始人和高管如何提高工作效率
看13位CEO、创始人和高管如何提高工作效率
Web伟德国际app人员为什么越来越懒了?
Web伟德国际app人员为什么越来越懒了?
程序员应该关注的一些事儿
程序员应该关注的一些事儿
“懒”出效率是程序员的美德
“懒”出效率是程序员的美德
60个伟德国际app者不容错过的免费资源库
60个伟德国际app者不容错过的免费资源库
鲜为人知的编程真相
鲜为人知的编程真相
什么才是优秀的用户界面设计
什么才是优秀的用户界面设计
每天工作4小时的程序员
每天工作4小时的程序员
“肮脏的”IT工作排行榜
“肮脏的”IT工作排行榜
十大编程算法助程序员走上高手之路
十大编程算法助程序员走上高手之路
我跳槽是因为他们的显示器更大
我跳槽是因为他们的显示器更大
2013年美国伟德国际app者薪资调查报告
2013年美国伟德国际app者薪资调查报告
10个调试和排错的小建议
10个调试和排错的小建议
程序员都该阅读的书
程序员都该阅读的书
做程序猿的老婆应该注意的一些事情
做程序猿的老婆应该注意的一些事情
Java程序员必看电影
Java程序员必看电影
漫画:程序员的工作
漫画:程序员的工作
一个程序员的时间管理
一个程序员的时间管理
程序员必看的十大电影
程序员必看的十大电影
 程序员的样子
程序员的样子
编程语言是女人
编程语言是女人
代码女神横空出世
代码女神横空出世
程序员的一天:一寸光阴一寸金
程序员的一天:一寸光阴一寸金
我的丈夫是个程序员
我的丈夫是个程序员
为啥Android手机总会越用越慢?
为啥Android手机总会越用越慢?
要嫁就嫁程序猿—钱多话少死的早
要嫁就嫁程序猿—钱多话少死的早
程序猿的崛起——Growth Hacker
程序猿的崛起——Growth Hacker
程序员的鄙视链
程序员的鄙视链
软件伟德国际app程序错误异常ExceptionCopyright © 2009-2015 MyException 版权所有