【算法学习】模拟退火

合集 - 【算法学习】(35)

1.【算法学习】排序2024-08-272.【算法学习】Manacher 马拉车2024-08-293.【算法学习】KMP 算法2024-09-014.LCA 最近公共祖先（树链和倍增）这次真有树链了！！！2024-08-305.线段覆盖问题2024-08-296.【算法学习】学换根dp有感2024-09-207.二分图最大匹配2024-09-178.【算法学习】01BFS2024-09-179.洛谷 P1892 [BOI2003] 团伙 种类并查集 扩展域并查集2024-09-1110.【算法学习】高斯消元法2024-09-1011.贝叶斯公式2024-09-0712.背包2024-09-05

13.【算法学习】模拟退火2024-09-05

14.【算法学习】基环树2024-10-2915.【算法学习】树链部分2024-10-2916.【算法学习】莫队2024-10-2917.【算法学习】分块九讲2024-10-2818.平衡树2024-10-2319.圆方树2024-10-2220.【算法学习】点分治2024-10-1121.公式2024-10-0922.【算法学习】笛卡尔树2024-10-0423.【算法学习】悬线法2024-10-0424.欧几里得算法与 EX2024-09-2325.【算法学习】逆元与求解2024-09-2326.【算法学习】费马定理2024-09-2327.三分2024-09-2328.裴蜀定理2024-09-2329.【算法学习】欧拉函数φ2024-09-2330.【算法学习】二维转一维问题2024-11-0331.【算法学习】扫描线2024-11-0132.【算法学习】矩阵乘法2024-10-3133.【算法学习】同余最短路2024-11-0434.【算法学习】组合数学2024-11-0535.【算法学习】反悔贪心2024-11-06

收起

P1337 吊打XXX

SA 无疑就这几个步骤。

1. 设初温 t;

2. 降低系数 down<1

3. 设定阈值

4. 更劣随机概率选择（温度越低，概率越低）

AI 评价：

#include<bits/stdc++.h>
using namespace std;
#define down 0.996 //降低系数 
int n;
struct node{
	int x,y,w;
}point[10006];
 
double ansx,ansy,answ;
 
double energy(double x,double y){
	double r=0,dx,dy;
	
	for(int i=1;i<=n;i++){//计算能量总和，总和越小越稳定 
		dx=x-point[i].x;
		dy=y-point[i].y;
		r+=sqrt(dx*dx+dy*dy)*point[i].w;
	}
	return r;
}
 
void sa(){
	double t=3000;//设置初温，是控制接受较差解的概率 
	while(t>1e-15){//设置温度阈值,当温度低于阈值时，停止降温过程
		double ex=ansx+(rand()*2-RAND_MAX)*t;//根据当前温度随机生成新的解 
		double ey=ansy+(rand()*2-RAND_MAX)*t; 
		double ew=energy(ex,ey);//计算新值 
		double de=ew-answ;//计算新值和旧值的差 
		if(de<0){//越小越稳定，更优 
			ansx=ex;
			ansy=ey;
			answ=ew;
		}
		else if(exp(-de/t)*RAND_MAX>rand()){//概率选择更劣的值，温度越小概率越小 
			ansx=ex;
			ansy=ey;
		}
		t*=down;//降低温度 
	}
}
 
int main(){
	ios::sync_with_stdio(false);
	cin>>n;
	
	for(int i=1;i<=n;i++){
		cin>>point[i].x>>point[i].y>>point[i].w;
		ansx+=point[i].x;//初始化 
		ansy+=point[i].y;
	}
	
	ansx/=n;//取平均值更合理 
	ansy/=n;
	answ=energy(ansx,ansy);//初始值 
	
	while ((double)clock()/CLOCKS_PER_SEC<0.975) sa();//卡时 
	
	printf("%.3lf %.3lf",ansx,ansy);
	
    return 0;
}