Poj1379 poj 2539（模擬退火

題目：在矩形中有若干個點，求乙個點使得所有點到該點的最小距離最大。

思路：這個是2023年顧研**上的例題，可以比較簡單地用模擬退火演算法求解。所謂模擬退火就是先隨機出若干個點，然後以某一特定步長嘗試周圍的解，而後逐漸縮小步長，知道步長小於特定值，跳出。這個演算法雖然簡單易行，但是其正確性並不是非常有保障，（不合理的隨機方法可能會降低到達正確解的概率），而且幾個主要引數的選取也主要看經驗，畢竟對於乙個隨機化演算法做出像**中一樣的嚴格評估還是比較困難的。。。。所以說即使實現完全正確，能否通過有時也是玄學的領域。。。。。尤其對於不太適合隨機化的題。

ps：poj不能用time，一呼叫就re，查了乙個多小時錯才發現。。。。以後不用time了。。。

/*
* @author: cwind
* * 蒟蒻只能做幾個水題。。
*///
#pragma comment(linker, "/stack:102400000,102400000")
#include #include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
using
namespace
std;
#define ios std::ios::sync_with_stdio (false);std::cin.tie(0)
#define pb push_back
#define pb pop_back
#define bk back()
#define fs first
#define se second
#define sq(x) (x)*(x)
#define eps (1e-3)
#define iinf (1<<29)
#define linf (1ll<<59)
#define finf (1e100)
#define inf 1000000000
const
double pi=acos(-1.0
);typedef 
long
long
ll;typedef unsigned 
long
long
ull;
typedef pair
pii;
typedef pair
p;inline 
double dis(double x1,double y1,double x2,double
y2)int
t;int
x,y,m;
const
int l=30
;double zx[1005],zy[1005
];int
main()
for(int i=0;i<30;i++)
}while(delta>eps)
if(md>d[i])}}
delta*=0.8
; }
double ans=0
; 
int ansp=0
; 
for(int i=0;i<30;i++)
}printf(
"the safest point is (%.1f, %.1f).\n
",pos[ansp][0],pos[ansp][1
]); }
return0;
}

view code

poj2539

題目：求使得一點到其他各點距離之和最小的點。

思路：和上題一樣的寫法，很順利的1a了。。

/*
* @author: cwind
* * 蒟蒻只能做幾個水題。。
*///
#pragma comment(linker, "/stack:102400000,102400000")
#include #include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
using
namespace
std;
#define ios std::ios::sync_with_stdio (false);std::cin.tie(0)
#define pb push_back
#define pb pop_back
#define bk back()
#define fs first
#define se second
#define sq(x) (x)*(x)
#define eps (1e-3)
#define iinf (1<<29)
#define linf (1ll<<59)
#define finf (1e100)
#define inf 1000000000
const
double pi=acos(-1.0
);typedef 
long
long
ll;typedef unsigned 
long
long
ull;
typedef pair
pii;
typedef pair
p;inline 
double dis(double x1,double y1,double x2,double
y2)const
int maxp=50
;const
int maxx=10000
;int
n;double c[105][2
];double p[maxp][2
],d[maxp];
double delta=5000
;int l=30
;double a=0.8
;int
main()
for(int i=0;i)
d[i]=sdis;
}while(delta>eps)
if(sdis}}
delta*=a;
}double ans=1e100;
for(int i=0;i)
printf(
"%.0f
",ans);
return0;
}

view code

Poj1379 poj 2539（模擬退火

poj 2420 poj1379 模擬退火

poj1379 模擬退火

poj 1379 模擬退火

Poj1379 poj 2539（模擬退火

poj 2420 poj1379 模擬退火

poj1379 模擬退火

poj 1379 模擬退火

相關推薦