bzoj3680: 吊打XXX
来源:互联网 发布:战舰世界详细数据 编辑:程序博客网 时间:2024/05/22 01:35
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3680
思路:puts("nan nan");得AC
爬山就够了,膜你退火也可以。
膜你退火就是在膜你一个退火的过程,他和爬山的区别就在于,它多了一个温度参数。
我们可以发现,越到后面,我们就越接近。
所以我们应该把修改的范围越改越小,接受较劣解的可能性也应该调小。
于是我们引入一个温度变量T,膜你退火的过程,温度逐渐下降。
下面给出模拟退火的流程:
设定初始较高的温度T
while 温度>设定的最低值
随机得到一个新解(温度越高,新解与旧解的差异越大)
更新答案
根据新解与旧解之间的优劣关系,以一定概率接受新解(越优越有可能)
以一定方式降温
endwhile
为了保证答案的质量,我们可以在最优解附近再进行随机,并更新答案
#include<cmath>#include<cstdio>#include<cstdlib>#include<cstring>#include<iostream>#include<algorithm>const int maxn=10010;using namespace std;int n;double mins=1e100;struct poi{double x,y,g;}p[maxn],ans;double dis(const poi &a,const poi &b){return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}double check(poi a){double res=0.0;for (int i=1;i<=n;i++) res+=dis(a,p[i])*p[i].g;if (res<mins) ans=a,mins=res;return res;}double getrd(){return rand()%10000/10000.0;}void SA(double T){poi now=ans;while (T>0.001){poi newp;newp.x=now.x+T*(getrd()*2-1.0);newp.y=now.y+T*(getrd()*2-1.0);double delta=check(now)-check(newp);if (delta>0||exp(delta/T)>getrd()) now=newp;T*=0.991;}for (int i=1;i<=1000;i++){poi newp;newp.x=ans.x+T*(getrd()*2-1.0);newp.y=ans.y+T*(getrd()*2-1.0);check(newp);}}int main(){srand(8226312452);scanf("%d",&n);for (int i=1;i<=n;i++){scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].g);ans.x+=p[i].x,ans.y+=p[i].y;}ans.x/=(double)n,ans.y/=(double)n;SA(100000.0);printf("%.3lf %.3lf\n",ans.x,ans.y);return 0;}
0 0
- 【BZOJ3680】吊打XXX
- [BZOJ3680]吊打XXX
- bzoj3680: 吊打XXX
- BZOJ3680 吊打XXX
- bzoj3680 吊打XXX
- [BZOJ3680]吊打XXX && 模拟退火
- 【bzoj3680】【吊打XXX】【模拟退火】
- bzoj3680 吊打XXX 爬山算法
- 模拟退火法(吊打XXX)Bzoj3680
- 【BZOJ3680】吊打XXX 广义费马点 模拟退火
- BZOJ3680 吊打XXX(模拟退火算法)
- 【模拟退火\爬山算法】[HYSBZ/BZOJ3680]吊打XXX(吊打GTY)
- [bzoj3680] 吊打XXX:模拟退火 or 模拟力学情景
- 【BZOJ3680】吊打XXX 计算几何 广义费马点+模拟退火(爬山算法)
- bzoj-3680 吊打XXX
- 3680: 吊打XXX
- 【bzoj 3680】吊打XXX
- BZOJ 3680: 吊打XXX
- 浏览器存储
- python引入PIL做验证码,发现字体不支持的解决办法
- LeetCode97—Interleaving String
- Codeforces Round #345 (Div. 2)——A. Joysticks(模拟+特判)
- 搭建linux服务器之Debian安装Vim7.4
- bzoj3680: 吊打XXX
- sougou
- Android高阶UI之加入购物车选择颜色、尺码实现
- 番茄工作法学习笔记
- LA 2191(p247)----Potentiometers
- 总结&备忘:ShareSDK中微信分享报-6错误
- Codeforces Round #345 (Div. 2)
- java多线程断点下载
- 数列排序(选择排序法和冒泡排序法)