【洛谷1337】[JSOI2004] 吊打XXX（模拟退火经典题）

点此看题面

大致题意： 一个平面上有 n n n个点，每个点有1个权值，现在要选择平面上的一个点，使这 n n n个点的权值乘上到达选定点的距离之和最小。

模拟退火

我们可以用模拟退火来做这道题。

L i n k Link Link

模拟退火详见博客模拟退火，随机化下的贪心

先将 ( 0 , 0 ) (0,0) (0,0)设定为答案，随后不断选取一个新的坐标，比较选择该点时的代价与当前答案的代价。若小于当前答案的代价，则更新答案，否则，将有一定概率更新答案（更新坐标的幅度随时间的增大而减小，更新答案的概率随时间的增大以及两个代价之差的增大而减小）。

只要多模拟退火几遍，~~或者保持一颗虔诚的心，~~就能过了。

代码

#include<bits/stdc++.h>
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)<(y)?(x):(y))
#define LL long long
#define swap(x,y) (x^=y,y^=x,x^=y)
#define delta 0.99
#define N 1000
using namespace std;
int n,x[N+5],y[N+5],w[N+5];
inline char tc()
{static char ff[100000],*A=ff,*B=ff;return A==B&&(B=(A=ff)+fread(ff,1,100000,stdin),A==B)?EOF:*A++;
}
inline void read(int &x)
{x=0;int f=1;char ch;while(!isdigit(ch=tc())) f=ch^'-'?1:-1;while(x=(x<<3)+(x<<1)+ch-'0',isdigit(ch=tc()));x*=f;
}
inline void write(int x)
{if(x<0) putchar('-'),x=-x;if(x>9) write(x/10);putchar(x%10+'0');
}
inline double dis(double nx,double ny)//计算出n个点的权值乘上到选定点的距离之和
{register int i;double res=0.0;for(i=1;i<=n;++i) res+=(double)sqrt((x[i]-nx)*(x[i]-nx)+(y[i]-ny)*(y[i]-ny))*w[i];return res;
}
inline void SA(double &X,double &Y)//Simulated Annealing，模拟退火
{double tt=3000,res=dis(X,Y);//tt表示变化量，res表示当前代价while(tt>0.000000000000001){double nx=X+(rand()*2-RAND_MAX)*tt,ny=Y+(rand()*2-RAND_MAX)*tt,new_res=dis(nx,ny);//计算出新的坐标以及新的代价if(new_res<res||exp((res-new_res)/tt)*RAND_MAX>rand()) res=new_res,X=nx,Y=ny;//如果新的代价小于当前代价，或在一定的几率下，更新当前状态tt*=delta;//将变化量减小，是一个模拟物理学上的退火操作的过程}
}
int main()
{srand(time(NULL)),srand(rand()),srand(rand());register int i;for(read(n),i=1;i<=n;++i) read(x[i]),read(y[i]),read(w[i]);double ans_x=0.0,ans_y=0.0;for(i=1;i<=10;++i) SA(ans_x,ans_y);//模拟退火10次return printf("%.3lf %.3lf",ans_x,ans_y),0;
}

【洛谷1337】[JSOI2004] 吊打XXX（模拟退火经典题）相关推荐

洛谷 P1337 [JSOI2004]平衡点 / 吊打XXX
具体题目见洛谷 P1337 [JSOI2004]平衡点 / 吊打XXX 方法一:模拟退火思路:当绳子平衡时,系统的能量最小,则此时物体总的重力势能要最小,也就是物体重量一定的情况下绳长最长,即桌子上 ...
python刷题 NOI题库 python题解洛谷、牛客网、AcWing 刷题等
NOI题库 python题解-2022.01.07整理(1.1-1.3) NOI题库 python题解-2022.01.07整理(1.1-1.3)_dllglvzhenfeng的博客-CSDN博客 N ...
P1337 [JSOI2004]平衡点 / 吊打XXX 模拟退火
平衡点 / 吊打XXX solution 记录一下模拟退火的几个题目 code /*SiberianSquirrel*//*CuteKiloFish*/ #include <bits/stdc+ ...
洛谷、牛客网、AcWing 刷题（python版）
牛客网python专项练习整理(一) https://blog.csdn.net/weixin_41913008/article/details/87203468 牛客网剑指offer--python ...
洛谷P1337 [JSOI2004]平衡点 / 吊打XXX（模拟退火）
传送门先坑着,联赛活着回来的话我就写(意思就是我绝对不会写了) 1 //minamoto 2 #include<cstdio> 3 #include<cmath> 4 #in ...
BZOJ 3680 吊打XXX 模拟退火
首先这题应该改名叫吊打出题人题目大意:给定n个质点,求重心这n个质点的重心满足Σ(重心到点i的距离)*g[i]最小模拟退火的裸题尼玛交了两篇死活过不去各种改参数最后发现是我的INF不够大 ...
洛谷 P5713、P5719、P5730、P5737题题解（Java语言描述）
P5713 题目要求 P5713题目链接 AC代码(Java语言描述) import java.util.Scanner;public class Main {public static void m ...
大数高精加减乘除（洛谷P1601、P2142、P1303、P1480题题解，Java语言描述）
写在前面算法归算法,有现成的API能用为什么不用? (纯算法党别喷我,以后我再写高精算法的文) P1601 → A+B Problem 题目要求 P1601题目链接 AC代码(Java语言描述) i ...
洛谷P5733、P5734、P5735、P5739题题解（Java语言描述）
P5733题目要求 P5733题目链接 P5733分析 toUpperCase()就行了-- P5733AC代码(Java语言描述) import java.util.Scanner;public c ...

【洛谷1337】[JSOI2004] 吊打XXX（模拟退火经典题）

模拟退火

L i n k Link Link

代码

【洛谷1337】[JSOI2004] 吊打XXX（模拟退火经典题）相关推荐

最新文章

热门文章